簡(jiǎn)介聯(lián)合目標(biāo)雷達(dá)特征估計(jì)與目標(biāo)檢測(cè)的寬帶認(rèn)知雷達(dá)波形優(yōu)化是一篇面向雷達(dá)信號(hào)處理與認(rèn)知雷達(dá)領(lǐng)域研究人員、工程師的研究論文。論文針對(duì)寬帶認(rèn)知雷達(dá)中目標(biāo)雷達(dá)特征TRS隨目標(biāo)姿態(tài)與時(shí)間變化的問(wèn)題提出在確保TRS估計(jì)精度的前提下通過(guò)最大化接收回波平均信號(hào)與雜波噪聲比來(lái)設(shè)計(jì)波形并建立凸優(yōu)化模型求解同時(shí)設(shè)計(jì)快速分層策略滿(mǎn)足實(shí)時(shí)處理需求基于數(shù)值仿真驗(yàn)證了檢測(cè)性能的提升。資源為1個(gè)PDF文件壓縮包大小1.55MB包含論文全文從問(wèn)題建模、理論推導(dǎo)到算法實(shí)現(xiàn)與仿真分析均有完整呈現(xiàn)適合作為WCR波形優(yōu)化研究的參考文獻(xiàn)。目前已吸引130人瀏覽學(xué)習(xí)適合具備雷達(dá)基礎(chǔ)知識(shí)的讀者深入研讀有助于把握認(rèn)知雷達(dá)自適應(yīng)波形設(shè)計(jì)的核心方法。1. 寬帶認(rèn)知雷達(dá)為什么要做檢測(cè)和估計(jì)的聯(lián)合波形優(yōu)化先厘清目標(biāo)到底在哪寬帶雷達(dá)一上來(lái)很多朋友還在沿用窄帶那套“目標(biāo)就是一個(gè)點(diǎn)、先檢測(cè)后測(cè)量、兩步走”的思路結(jié)果一到實(shí)測(cè)就翻車(chē)——目標(biāo)回波跨了幾十個(gè)距離單元檢測(cè)器明明把目標(biāo)框住了后面的特征提取卻在旁瓣和雜波里把參數(shù)估得亂七八糟。這個(gè)問(wèn)題的根子不在接收機(jī)而在發(fā)射波形根本沒(méi)有為“既要發(fā)現(xiàn)目標(biāo)、又要看清目標(biāo)”這件事做過(guò)優(yōu)化。聯(lián)合目標(biāo)雷達(dá)特征估計(jì)和目標(biāo)檢測(cè)的寬帶認(rèn)知雷達(dá)波形優(yōu)化就是要在發(fā)射端把這兩件事當(dāng)成一個(gè)整體來(lái)設(shè)計(jì)讓波形本身既帶檢測(cè)友好性又帶估計(jì)友好性。我最早接觸這個(gè)方向是在某高校實(shí)驗(yàn)室做某雷達(dá)系統(tǒng)的預(yù)研項(xiàng)目當(dāng)時(shí)導(dǎo)師給的一句話(huà)我一直留著“波形優(yōu)化不是挑一個(gè)窗函數(shù)是把你對(duì)環(huán)境和目標(biāo)的認(rèn)知全部壓進(jìn)發(fā)射信號(hào)里?！闭J(rèn)知雷達(dá)的閉環(huán)長(zhǎng)這樣上一幀接收數(shù)據(jù)先做感知估計(jì)出目標(biāo)散射中心分布、雜波協(xié)方差等先驗(yàn)然后用這些先驗(yàn)設(shè)計(jì)下一幀的發(fā)射波形接收端再針對(duì)這個(gè)發(fā)射波形做匹配或自適應(yīng)濾波。檢測(cè)和估計(jì)的聯(lián)合優(yōu)化就落在“設(shè)計(jì)下一幀波形”這個(gè)環(huán)節(jié)上。這篇文章不聊科普直接講清楚建模、準(zhǔn)則構(gòu)造、求解流程和那些讓我反復(fù)折騰的坑新手能順著跑通仿真熟手能對(duì)照著調(diào)參數(shù)。2. 寬帶場(chǎng)景下的信號(hào)模型與認(rèn)知閉環(huán)目標(biāo)不再是一個(gè)點(diǎn)波形優(yōu)化才有意義2.1 多散射中心模型下的回波表達(dá)為什么窄帶模型在寬帶下會(huì)失效寬帶雷達(dá)的帶寬大距離分辨單元小一個(gè)真實(shí)的飛機(jī)、艦船、車(chē)輛目標(biāo)在距離維上會(huì)占據(jù)多個(gè)分辨單元每個(gè)單元里有一個(gè)或多個(gè)強(qiáng)散射中心。這時(shí)候目標(biāo)不能再建模成單個(gè)時(shí)延副本而應(yīng)該用一個(gè)目標(biāo)沖激響應(yīng)序列來(lái)表達(dá)。設(shè)發(fā)射波形為s(t)目標(biāo)沖激響應(yīng)為h(t)則理想回波是兩者卷積的結(jié)果再加上噪聲和雜波r(t) s(t) * h(t) c(t) n(t)離散化后寫(xiě)成向量形式會(huì)更直觀r S h c n其中S是由發(fā)射波形s構(gòu)造的卷積矩陣h是目標(biāo)散射中心的復(fù)幅度向量維度等于目標(biāo)占據(jù)的距離單元數(shù)。窄帶模型下h只有一個(gè)非零元素寬帶模型下h可能有幾十個(gè)非零元素。這一下就把問(wèn)題復(fù)雜度抬上來(lái)了波形優(yōu)化不再是選一個(gè)波形讓匹配濾波輸出峰值最高而是要保證整個(gè)h向量都能被感知到并且估計(jì)誤差足夠小。我一般會(huì)直接把h的先驗(yàn)協(xié)方差R_h當(dāng)作已知量放進(jìn)設(shè)計(jì)里這樣做出來(lái)的波形對(duì)目標(biāo)結(jié)構(gòu)有明確的針對(duì)性。散射中心模型還有個(gè)好處它天然支持檢測(cè)和估計(jì)的聯(lián)合——檢測(cè)關(guān)注的是S h在接收端的能量匯聚程度估計(jì)關(guān)注的是h各分量的可辨識(shí)度這兩個(gè)訴求在同一組波形參數(shù)里被統(tǒng)一起來(lái)。2.2 認(rèn)知閉環(huán)里真正能優(yōu)化的自由度發(fā)射波形、接收濾波器和調(diào)度參數(shù)認(rèn)知雷達(dá)與傳統(tǒng)自適應(yīng)雷達(dá)最大的區(qū)別在于發(fā)射端可編程。傳統(tǒng)雷達(dá)的發(fā)射波形是出廠定死的接收端再花哨也只能被動(dòng)適配。認(rèn)知雷達(dá)的發(fā)射機(jī)支持在幀間快速切換不同的相位編碼、頻率編碼或調(diào)頻斜率所以?xún)?yōu)化變量首先落在波形序列上。恒模約束是必須考慮的因?yàn)榇蠖鄶?shù)發(fā)射機(jī)功放工作在飽和區(qū)幅度調(diào)制會(huì)帶來(lái)嚴(yán)重的非線(xiàn)性失真。除了波形本身的碼片相位接收濾波器也是聯(lián)合優(yōu)化的一半變量。常見(jiàn)做法是波形設(shè)計(jì)后配套一個(gè) MVDR 類(lèi)的空時(shí)濾波器兩者交替迭代求解整體性能遠(yuǎn)比單獨(dú)設(shè)計(jì)一端要好。調(diào)度層面的自由度也容易被忽略波形是在相干處理間隔CPI級(jí)切換還是脈沖級(jí)切換決定了感知信息和波形設(shè)計(jì)的銜接顆粒度。我習(xí)慣按 CPI 為單位做閉環(huán)因?yàn)閱蚊}沖的信噪比太低估計(jì)出的目標(biāo)沖激響應(yīng)方差太大直接用這個(gè)結(jié)果去優(yōu)化波形很容易被噪聲牽著走。一個(gè) CPI 內(nèi)的多脈沖回波做相干積累后目標(biāo)沖激響應(yīng)的估計(jì)可信度高一個(gè)數(shù)量級(jí)這時(shí)候再算波形優(yōu)化才有意義。設(shè)計(jì)維度整理如下自由度典型取值范圍對(duì)優(yōu)化問(wèn)題的影響碼片相位連續(xù)相位或有限相位集決定頻譜形狀和距離旁瓣碼片數(shù) N642048N 越大距離分辨率越高卷積矩陣越大恒模約束必須滿(mǎn)足非線(xiàn)性約束投影處理的關(guān)鍵頻譜模板頻帶內(nèi)平坦、帶外抑制約束波形頻譜不能越界接收濾波器權(quán)值與波形聯(lián)合迭代把干擾抑制和信號(hào)匹配分離2.3 感知先驗(yàn)的獲取目標(biāo)沖激響應(yīng)估計(jì)與雜波協(xié)方差矩陣的實(shí)用性處理波形優(yōu)化的輸入是感知層給出的先驗(yàn)。目標(biāo)沖激響應(yīng)h的估計(jì)可以用上一幀的高分辨距離像做幅相提取也可以用多幀回波做參數(shù)化估計(jì)。關(guān)鍵是這個(gè)估計(jì)一定是有誤差的而且誤差會(huì)直接傳導(dǎo)到波形優(yōu)化結(jié)果里。雜波協(xié)方差矩陣R_c的估計(jì)則需要用不含目標(biāo)的參考距離單元做樣本協(xié)方差估計(jì)樣本數(shù)要大于雜波自由度數(shù)的兩倍否則矩陣病態(tài)。這個(gè)數(shù)字是我在仿真里對(duì)比出來(lái)的樣本數(shù)不夠時(shí)優(yōu)化出來(lái)的波形會(huì)把功率集中到雜波樣本的隨機(jī)凹陷處看起來(lái)目標(biāo)函數(shù)的理論值很好看實(shí)測(cè)一塌糊涂。處理這些先驗(yàn)工程上有一個(gè)穩(wěn)健化手段——對(duì)角加載。我一般會(huì)在估計(jì)出的協(xié)方差矩陣上加一個(gè)正則項(xiàng)R R_hat epsilon * Iepsilon 取矩陣最大特征值的 0.01 到 0.1 倍。這個(gè)操作會(huì)犧牲一點(diǎn)理論最優(yōu)性但換來(lái)了對(duì)估計(jì)誤差的魯棒性。認(rèn)知雷達(dá)的閉環(huán)最容易死在哪就死在感知層給出的先驗(yàn)里帶著感知層的錯(cuò)誤優(yōu)化層再把這錯(cuò)誤放大成波形里的定向缺陷。所以現(xiàn)在我做波形優(yōu)化前第一件事不是調(diào)算法而是先把先驗(yàn)的置信度量化出來(lái)誤差大的維度直接在優(yōu)化目標(biāo)里降權(quán)。3. 聯(lián)合目標(biāo)函數(shù)的構(gòu)造讓檢測(cè)準(zhǔn)則和估計(jì)準(zhǔn)則在一個(gè)數(shù)學(xué)問(wèn)題里共存3.1 檢測(cè)性能的量化擴(kuò)展目標(biāo)下的輸出 SINR 表達(dá)式檢測(cè)性能的度量不能直接從檢測(cè)概率P_d出發(fā)因?yàn)镻_d是高度非線(xiàn)性的門(mén)限函數(shù)直接放進(jìn)優(yōu)化問(wèn)題里求梯度很痛苦。工程上普遍采用輸出信干噪比作為代理指標(biāo)。對(duì)擴(kuò)展目標(biāo)接收機(jī)將多個(gè)距離單元的相干積累結(jié)果合并輸出 SINR 寫(xiě)成:SINR_out (w^H R_s w) / (w^H R_in w)其中R_s是目標(biāo)信號(hào)分量經(jīng)過(guò)波形調(diào)制后的協(xié)方差矩陣R_in是噪聲加雜波的協(xié)方差矩陣w是接收濾波器權(quán)向量。R_s里包含發(fā)射波形通過(guò)卷積矩陣對(duì)目標(biāo)沖激響應(yīng)的調(diào)制效應(yīng)所以改變波形會(huì)直接改變SINR_out。為什么用 SINR 而不是直接算P_d因?yàn)閷?duì)非起伏目標(biāo)檢測(cè)概率是 SINR 的單調(diào)函數(shù)最大化 SINR 等價(jià)于最大化檢測(cè)概率而且 SINR 對(duì)波形變量的梯度計(jì)算方便。對(duì) Swerling 起伏目標(biāo)也有近似關(guān)系只是要把目標(biāo)起伏方差計(jì)入R_s的構(gòu)造。我在做聯(lián)合優(yōu)化時(shí)檢測(cè)項(xiàng)一律用 SINR 的歸一化形式避免它和估計(jì)項(xiàng)在數(shù)值量級(jí)上打架。3.2 估計(jì)性能的量化CRB 和互信息怎么選估計(jì)性能的度量有兩個(gè)主流選擇——克拉美羅界CRB和互信息MI。CRB 給的是目標(biāo)參數(shù)估計(jì)方差的下界物理含義清晰直接對(duì)應(yīng)測(cè)距、測(cè)速、散射中心幅度的估計(jì)精度。構(gòu)造 CRB 需要計(jì)算 Fisher 信息矩陣?yán)锩姘l(fā)射波形對(duì)目標(biāo)參數(shù)的導(dǎo)數(shù)項(xiàng)所以 CRB 自然成為波形的函數(shù)。互信息衡量的是接收回波與目標(biāo)沖激響應(yīng)之間的統(tǒng)計(jì)依賴(lài)程度最大化互信息等價(jià)于讓接收回波攜帶盡可能多的目標(biāo)信息。兩者不是完全等價(jià)的CRB 是參數(shù)估計(jì)視角MI 是信息論視角。當(dāng)目標(biāo)模型是確定參數(shù)加高斯噪聲時(shí)最大化 MI 往往也同時(shí)逼近最優(yōu) CRB但數(shù)值上 CRB 更容易算尤其是散射中心個(gè)數(shù)超過(guò)幾十個(gè)時(shí)MI 的矩陣行列式計(jì)算量會(huì)漲得很快。我個(gè)人的選型習(xí)慣是目標(biāo)模型明確、主要關(guān)心測(cè)距和測(cè)速精度時(shí)用 CRB目標(biāo)結(jié)構(gòu)復(fù)雜、說(shuō)不清多少個(gè)散射中心時(shí)用 MI。聯(lián)合優(yōu)化里估計(jì)項(xiàng)的歸一化基準(zhǔn)一般取“參考波形下的 CRB 值”這樣估計(jì)項(xiàng)是無(wú)量綱的相對(duì)量和檢測(cè)項(xiàng)的相對(duì) SINR 可以放在同一個(gè)量級(jí)上比較。CRB 計(jì)算還有一個(gè)細(xì)節(jié)當(dāng)目標(biāo)沖激響應(yīng)的某些分量幅度接近零時(shí)Fisher 信息矩陣容易接近奇異此時(shí)要做對(duì)角加載或直接刪掉這些弱散射中心對(duì)應(yīng)的維度。這個(gè)處理在仿真里不做的話(huà)優(yōu)化器會(huì)瘋狂地把能量投向噪聲級(jí)的弱散射中心整個(gè)波形被帶偏。3.3 聯(lián)合準(zhǔn)則的歸一化設(shè)計(jì)加權(quán)求和與分層約束兩條路線(xiàn)把檢測(cè)和估計(jì)放進(jìn)同一個(gè)目標(biāo)函數(shù)最簡(jiǎn)單的形式是加權(quán)求和max (1 - λ) * SINR_norm λ * MI_norm 約束 恒模約束 |s_n| 1 頻譜模板約束λ 在 0 到 1 之間體現(xiàn)系統(tǒng)對(duì)檢測(cè)和估計(jì)的偏好。λ 0 就是純檢測(cè)優(yōu)化λ 1 就是純估計(jì)優(yōu)化。這個(gè)形式簡(jiǎn)單直接但有一個(gè)大坑如果SINR_norm和MI_norm的數(shù)值范圍差一個(gè)數(shù)量級(jí)優(yōu)化結(jié)果會(huì)被量級(jí)大的一方主導(dǎo)λ 的白調(diào)節(jié)。所以我每次構(gòu)造目標(biāo)函數(shù)前都會(huì)先固定一個(gè)初始波形通常用線(xiàn)性調(diào)頻算出兩者的初始值再做歸一化。另一條路線(xiàn)是分層約束先保證檢測(cè)性能不低于門(mén)限再把估計(jì)誤差最小化。形式是:min CRB_norm s.t. SINR_out SINR_threshold |s_n| 1 頻譜模板約束這種約束型設(shè)計(jì)更貼近實(shí)際雷達(dá)系統(tǒng)的設(shè)計(jì)流程——檢測(cè)是底線(xiàn)估計(jì)是優(yōu)化目標(biāo)。缺點(diǎn)是可行域的搜索比加權(quán)型復(fù)雜需要結(jié)合懲罰函數(shù)法處理。我的經(jīng)驗(yàn)是系統(tǒng)指標(biāo)里對(duì)P_d有硬性要求的場(chǎng)景用分層約束追求整體最優(yōu)平衡的場(chǎng)景用加權(quán)求和。λ 的取值可以跑一組 Pareto 曲線(xiàn)來(lái)確定后面第 6 章會(huì)給出具體做法。4. 波形優(yōu)化求解與仿真實(shí)現(xiàn)交替迭代框架下的完整流程4.1 交替迭代求解框架與核心偽代碼聯(lián)合優(yōu)化問(wèn)題里發(fā)射波形和接收濾波器是耦合的最優(yōu)接收濾波器依賴(lài)于發(fā)射波形最優(yōu)發(fā)射波形也依賴(lài)于接收濾波器。把兩者同時(shí)求解是非凸的直接丟給通用優(yōu)化器又慢又不穩(wěn)定。我用的框架是交替迭代——固定波形解接收濾波器再固定濾波器解波形來(lái)回迭代。這個(gè)方法的工程依據(jù)是給定波形后最優(yōu)接收濾波器有閉式解即 MVDR 類(lèi)加權(quán)給定濾波器后波形優(yōu)化變成有恒模約束的二次優(yōu)化可以用投影梯度法處理。% 聯(lián)合波形與接收濾波器的交替迭代優(yōu)化 % 輸入: h_hat 目標(biāo)沖激響應(yīng)估計(jì), R_in 雜波加噪聲協(xié)方差, s0 初始波形 % 輸出: s_opt 優(yōu)化波形, w_opt 接收濾波器權(quán)向量 s s0; % 初始波形, 建議用LFM N_iter 80; % 最大迭代次數(shù) tol 1e-4; % 目標(biāo)函數(shù)相對(duì)變化容差 mu 0.02; % 波形投影梯度步長(zhǎng) J_old -inf; for iter 1:N_iter % 1) 固定當(dāng)前波形, 構(gòu)造卷積矩陣 S convmtx(s, length(h_hat)); R_s S * (h_hat * h_hat) * S; % 2) 求解最優(yōu)接收濾波器 (MVDR類(lèi), 歸一化無(wú)偏) R_x R_s R_in; w (R_x \ S * h_hat) / (h_hat * S * (R_x \ S * h_hat)); % 3) 固定接收濾波器, 計(jì)算聯(lián)合目標(biāo)函數(shù)對(duì)波形的梯度 grad compute_gradient(s, w, h_hat, R_in, lambda); % 4) 梯度上升并做約束投影 s_tmp s mu * real(grad); s project_constant_modulus(s_tmp); % 恒模投影 s project_spectrum_mask(s, mask); % 頻譜模板投影 % 5) 收斂判定 J_new joint_objective(s, w, h_hat, R_in, lambda); if abs(J_new - J_old) / abs(J_old) tol break; end J_old J_new; end s_opt s; w_opt w;邏輯很簡(jiǎn)單卷積矩陣S把波形和目標(biāo)沖激響應(yīng)的關(guān)系線(xiàn)性化R_s是目標(biāo)經(jīng)過(guò)波形調(diào)制后的信號(hào)協(xié)方差R_x是信號(hào)加干擾的總協(xié)方差。接收濾波器閉式解本質(zhì)上是維納濾波的歸一化形式它同時(shí)完成信號(hào)匹配和干擾抑制。compute_gradient這一步要按目標(biāo)函數(shù)的具體形式展開(kāi)如果是加權(quán)型目標(biāo)函數(shù)梯度就是 SINR 項(xiàng)的梯度加上 MI 項(xiàng)的梯度先用數(shù)值差分驗(yàn)證一階解析梯度的正確性再做解析化。我在這里吃過(guò)虧解析梯度一個(gè)符號(hào)寫(xiě)錯(cuò)迭代曲線(xiàn)前幾十次看著在上升后期突然發(fā)散定位了好久才發(fā)現(xiàn)是梯度項(xiàng)的轉(zhuǎn)置寫(xiě)漏了。4.2 恒模與頻譜約束的投影實(shí)現(xiàn)參數(shù)怎么設(shè)才不崩恒模約束是所有發(fā)射機(jī)都繞不開(kāi)的物理限制。投影操作本身很簡(jiǎn)單把每個(gè)碼片的幅度歸一化為同一值相位保持不變。但要注意投影和梯度步長(zhǎng)之間的配合。步長(zhǎng)mu太大會(huì)導(dǎo)致投影后相位亂跳目標(biāo)函數(shù)劇烈振蕩太小則收斂極慢。我一般把mu設(shè)置為 0.01 到 0.05 之間配合歸一化的梯度方向使用。頻譜模板投影比恒模投影復(fù)雜做法是先把波形做 FFT 到頻域把超過(guò)模板約束的譜線(xiàn)拉回邊界內(nèi)再 IFFT 回時(shí)域。但直接這樣做會(huì)破壞恒模性質(zhì)所以我用交替投影先在頻域限幅再投影回恒模多迭代幾輪。實(shí)際參數(shù)參考下表參數(shù)推薦值調(diào)整說(shuō)明交替迭代次數(shù)50100超過(guò) 100 次不收斂基本是梯度或初值問(wèn)題梯度步長(zhǎng) mu0.010.05目標(biāo)函數(shù)超調(diào)時(shí)減半收斂容差 tol1e-4相對(duì)目標(biāo)函數(shù)變化量恒模投影方式逐一碼片取模歸一注意保持總能量不變頻譜模板裕量模板限值下留 0.5 dB防止限幅后的離散化溢出對(duì)角加載系數(shù)0.010.1 倍最大特征值先驗(yàn)誤差大時(shí)取上限恒定模投影在迭代后期有個(gè)容易忽略的問(wèn)題如果波形受到頻譜投影的反復(fù)拉扯相位會(huì)出現(xiàn)高頻抖動(dòng)導(dǎo)致距離旁瓣升高。解決的辦法是頻譜投影的頻點(diǎn)間隔不要取得太細(xì)一般讓頻點(diǎn)數(shù)和碼片數(shù)一致即可。更激進(jìn)的做法是改變約束層級(jí)先保頻譜后保恒模交替迭代時(shí)給頻譜投影更高的優(yōu)先級(jí)因?yàn)轭l譜外泄對(duì)系統(tǒng)的影響通常比恒模偏差更致命。4.3 仿真驗(yàn)證指標(biāo)檢測(cè)概率、估計(jì)誤差和波形特性三張表要一起看仿真不能只看目標(biāo)函數(shù)曲線(xiàn)指標(biāo)要分開(kāi)驗(yàn)證。檢測(cè)側(cè)要跑蒙特卡洛計(jì)算檢測(cè)概率門(mén)限用 CA-CFAR 或者理想門(mén)限都可以但要在同一門(mén)限設(shè)置下對(duì)比不同波形才有意義。估計(jì)側(cè)計(jì)算 CRB 和實(shí)測(cè) RMSE 的比值這個(gè)比值越接近 1 說(shuō)明估計(jì)器效率越高如果遠(yuǎn)大于 1 說(shuō)明波形和估計(jì)器不匹配。波形側(cè)看自相關(guān)旁瓣電平PSLR和積分旁瓣電平ISLR因?yàn)榧词鼓繕?biāo)函數(shù)優(yōu)化得很好波形自身的旁瓣特性差接收端的弱小目標(biāo)會(huì)被主目標(biāo)旁瓣淹沒(méi)。我在仿真里固定一個(gè)指標(biāo)清單每次實(shí)驗(yàn)跑完就生成三張表指標(biāo)類(lèi)別具體指標(biāo)對(duì)比基準(zhǔn)檢測(cè)性能輸出 SINR、P_d、P_faLFM 基線(xiàn)波形估計(jì)性能CRB、RMSE、RMSE/CRB 比純檢測(cè)優(yōu)化波形波形質(zhì)量PSLR、ISLR、頻譜帶外抑制原始波形穩(wěn)健性對(duì)角線(xiàn)加載后的性能保持率無(wú)加載波形仿真場(chǎng)景設(shè)置也要注意目標(biāo)散射中心的位置和幅度要隨機(jī)化多次取平均不能只測(cè)一組固定目標(biāo)。信噪比從低到高掃 5 到 10 個(gè)點(diǎn)畫(huà)出檢測(cè)曲線(xiàn)和估計(jì)誤差曲線(xiàn)。只測(cè)一個(gè)高信噪比點(diǎn)得出來(lái)的結(jié)論不具代表性——低信噪比下聯(lián)合優(yōu)化的優(yōu)勢(shì)最大因?yàn)榇藭r(shí)傳統(tǒng)兩步走方案丟失的信息最多而這恰恰是整個(gè)方向最有價(jià)值的應(yīng)用場(chǎng)景。5. 波形優(yōu)化的常見(jiàn)坑位先驗(yàn)失配、量綱失衡與迭代翻車(chē)的對(duì)癥處理5.1 場(chǎng)景一感知先驗(yàn)失配優(yōu)化波形比固定 LFM 還差現(xiàn)象仿真里感知模塊給的目標(biāo)沖激響應(yīng)估計(jì)和真實(shí)值有偏差偏差大概 20% 以上優(yōu)化出來(lái)的波形性能反而低于不優(yōu)化的 LFM 基線(xiàn)。原因波形優(yōu)化把能量定向投向了估計(jì)錯(cuò)誤的散射中心位置真實(shí)目標(biāo)在那幾個(gè)距離單元上是空的能量全部浪費(fèi)。解決先驗(yàn)不確定度大的維度在目標(biāo)函數(shù)里直接降權(quán)或者退到穩(wěn)健優(yōu)化即把目標(biāo)沖激響應(yīng)用一個(gè)協(xié)方差矩陣而不是單點(diǎn)向量來(lái)描述優(yōu)化目標(biāo)變成在協(xié)方差范圍內(nèi)平均性能最優(yōu)。我后來(lái)養(yǎng)成的習(xí)慣是感知模塊輸出的每一個(gè)先驗(yàn)量都附帶一個(gè)置信度標(biāo)簽置信度低于閾值的分量在波形優(yōu)化里直接凍結(jié)不參與能量分配。5.2 場(chǎng)景二檢測(cè)項(xiàng)和估計(jì)項(xiàng)的量綱失衡聯(lián)合優(yōu)化變成單目標(biāo)優(yōu)化現(xiàn)象λ 取了 0.5看起來(lái)是五五開(kāi)實(shí)際優(yōu)化結(jié)果跟純檢測(cè)優(yōu)化完全一樣。查了目標(biāo)函數(shù)的數(shù)值才發(fā)現(xiàn)SINR 的數(shù)值在 20 dB 量級(jí)MI 的數(shù)值是 0.1 量級(jí)加權(quán)后檢測(cè)項(xiàng)直接壓垮了估計(jì)項(xiàng)。這不是 λ 的問(wèn)題是歸一化沒(méi)做好。解決每一輪迭代前都用參考波形計(jì)算 SINR 和 MI 的初始值然后各自除以初始值做歸一化再進(jìn)入加權(quán)。歸一化之后 λ 才有真實(shí)含義。我仿真里還遇到過(guò)反向問(wèn)題歸一化后估計(jì)項(xiàng)偶爾會(huì)劇烈波動(dòng)導(dǎo)致 λ 0.9 時(shí)波形在純檢測(cè)和純估計(jì)之間反復(fù)橫跳原因是對(duì)角加載系數(shù)太小CRB 矩陣有時(shí)接近奇異數(shù)值上撐出尖峰把加載系數(shù)從 0.01 提高到 0.1 后曲線(xiàn)平滑了。5.3 場(chǎng)景三恒模投影后頻譜外泄超標(biāo)發(fā)射機(jī)帶外輻射打不過(guò)去現(xiàn)象優(yōu)化后的波形頻譜旁瓣比模板要求高了 3 dB整個(gè)設(shè)計(jì)被射頻前端否決。原因恒模投影每次只做幅度歸一完全不考慮頻譜形狀而梯度方向并沒(méi)有天然地抑制帶外分量。換句話(huà)說(shuō)恒模投影這個(gè)操作本身就產(chǎn)生新的頻譜擴(kuò)展。解決投影順序調(diào)整——先做頻域限幅投影再做恒模投影循環(huán)往復(fù) 20 到 30 次讓兩個(gè)約束交替滿(mǎn)足。同時(shí)頻譜模板不要設(shè)得太緊給投影迭代留一點(diǎn)收斂緩沖域工程上我會(huì)把模板邊界放寬 0.5 dB。這個(gè)方法不保證兩個(gè)約束同時(shí)嚴(yán)格滿(mǎn)足但能讓誤差落在一個(gè)可接受的范圍內(nèi)。真正要求嚴(yán)格滿(mǎn)足時(shí)可以把頻譜模板轉(zhuǎn)成懲罰項(xiàng)放進(jìn)目標(biāo)函數(shù)用乘法子權(quán)重調(diào)節(jié)。5.4 場(chǎng)景四波形更新速率過(guò)快目標(biāo)檢測(cè)框沿距離維跳變現(xiàn)象幀與幀之間波形切換后同樣一個(gè)目標(biāo)在距離像上的位置偏移了幾個(gè)距離單元顯示端的目標(biāo)軌跡出現(xiàn)抖動(dòng)。原因波形變化導(dǎo)致匹配濾波輸出的峰值位置發(fā)生了輕微偏移而目標(biāo)沖激響應(yīng)的強(qiáng)散射中心相對(duì)時(shí)延沒(méi)變波形本身的群延遲特性變了。解決把波形更新的決策機(jī)制加滯后——只有當(dāng)新波形相對(duì)當(dāng)前波形的目標(biāo)函數(shù)增益超過(guò) 1 dB 時(shí)才切換否則保持上一幀波形。我一般還會(huì)約束連續(xù)幀之間波形相位變化的平滑性用一階低通濾波對(duì)碼片相位做時(shí)間平滑。這個(gè)處理雖然讓性能收斂慢了一點(diǎn)但系統(tǒng)穩(wěn)定性提升非常明顯尤其對(duì)跟蹤場(chǎng)景來(lái)說(shuō)軌跡平滑比單幀信噪比重要得多。5.5 場(chǎng)景五迭代從隨機(jī)相位初值出發(fā)每次結(jié)果都不一樣現(xiàn)象同一組參數(shù)跑十次得到十個(gè)不同的最優(yōu)波形目標(biāo)函數(shù)值還都差不多。原因問(wèn)題是高度非凸的隨機(jī)初值會(huì)收斂到不同的局部最優(yōu)點(diǎn)。仿真報(bào)告里如果只報(bào)一次結(jié)果別人復(fù)現(xiàn)時(shí)對(duì)不上數(shù)據(jù)可信度直接崩盤(pán)。解決統(tǒng)一用 LFM 波形做初值因?yàn)?LFM 相位是確定性的任何人在任何機(jī)器上跑都能得到同樣的初始狀態(tài)。想要跳出局部最優(yōu)就用多重啟動(dòng)策略——LFM 初值跑一輪隨機(jī)相位跑十輪取目標(biāo)函數(shù)最高的結(jié)果。這樣既保證可復(fù)現(xiàn)性又兼顧全局搜索。我自己的實(shí)踐里L(fēng)FM 初值雖然不一定是最優(yōu)的但它的物理意義明確優(yōu)化軌跡穩(wěn)定最適合做性能對(duì)比實(shí)驗(yàn)的統(tǒng)一基線(xiàn)。6. 怎么驗(yàn)證這套波形優(yōu)化真的有用一張圖和三組對(duì)照實(shí)驗(yàn)的復(fù)現(xiàn)習(xí)慣無(wú)論目標(biāo)函數(shù)寫(xiě)得多漂亮最終要落在“聯(lián)合優(yōu)化比傳統(tǒng)方法好”這個(gè)結(jié)論上。我的驗(yàn)證流程固定四步。第一步固定一個(gè)仿真場(chǎng)景目標(biāo)設(shè)置 5 個(gè)散射中心幅度按 1:0.8:0.5:0.3:0.2 遞減排列信噪比從 -5 dB 掃到 15 dB雜波協(xié)方差用 AR(4) 模型生成。第二步準(zhǔn)備四組對(duì)照波形LFM 基線(xiàn)、純檢測(cè)優(yōu)化波形、純估計(jì)優(yōu)化波形、聯(lián)合優(yōu)化波形保證四組波形帶寬和發(fā)射能量完全一致。第三步每組波形跑 1000 次蒙特卡洛統(tǒng)計(jì)輸出 SINR、檢測(cè)概率、測(cè)距 RMSE 和 CRB 比值。第四步出一張 Pareto 散點(diǎn)圖橫軸是歸一化 CRB 的倒數(shù)越大越好縱軸是歸一化 SINR越大越好把四個(gè)波形的點(diǎn)標(biāo)上去。這張圖是最有說(shuō)服力的驗(yàn)證——如果聯(lián)合優(yōu)化波形的點(diǎn)位于純檢測(cè)波形右上方的區(qū)域說(shuō)明它在不犧牲檢測(cè)性能的前提下確實(shí)換來(lái)了更好的估計(jì)精度如果聯(lián)合優(yōu)化的點(diǎn)落在純檢測(cè)和純估計(jì)兩點(diǎn)的連線(xiàn)路徑上說(shuō)明它是在兩目標(biāo)之間做權(quán)衡而不是真正同時(shí)收益。我做過(guò)的仿真里聯(lián)合優(yōu)化波形通常在中低信噪比區(qū)表現(xiàn)出明顯優(yōu)勢(shì)高信噪比下和純檢測(cè)優(yōu)化的差距收窄這符合直覺(jué)——信噪比足夠高時(shí)檢測(cè)已經(jīng)毫無(wú)壓力估計(jì)器喂給它的數(shù)據(jù)質(zhì)量也不差波形優(yōu)化的邊際收益自然變小。所以匯報(bào)結(jié)論時(shí)一定要注明信噪比區(qū)間避免夸大。另一個(gè)值得養(yǎng)成的習(xí)慣是參數(shù)敏感性分析。每次收斂后把感知模塊的目標(biāo)沖激響應(yīng)估計(jì)手動(dòng)加噪聲噪聲功率從 0% 加到 50%觀察優(yōu)化波形的 SINR 和 CRB 劣化曲線(xiàn)。如果 10% 的估計(jì)誤差就導(dǎo)致性能崩潰說(shuō)明這個(gè)波形對(duì)先驗(yàn)過(guò)度擬合工程上沒(méi)有實(shí)用價(jià)值。我會(huì)把 30% 誤差下的性能保持率作為一個(gè)硬性指標(biāo)低于 80% 就回爐調(diào)對(duì)角加載系數(shù)。最后提醒一句跑完所有仿真后把初始波形、中間迭代波形、最終波形和每輪目標(biāo)函數(shù)值全部落盤(pán)。這個(gè)習(xí)慣救過(guò)我很多次——某次仿真結(jié)果異常我翻記錄才發(fā)現(xiàn)是初始波形被上一組實(shí)驗(yàn)的變量污染了。整理好復(fù)現(xiàn)數(shù)據(jù)再下結(jié)論希望幫到你。本文還有配套的精品資源點(diǎn)擊獲取