簡介本資源是一套面向雷達(dá)信號處理方向研究生與工程師的MATLAB實(shí)戰(zhàn)仿真方案聚焦FDA-MIMO雷達(dá)多維參數(shù)估計(jì)這一前沿課題旨在解決傳統(tǒng)2D-MUSIC算法計(jì)算復(fù)雜度高、難以實(shí)時應(yīng)用的痛點(diǎn)。項(xiàng)目提出改進(jìn)型RD-MUSIC算法在保持角度維估計(jì)精度的同時顯著降低距離維運(yùn)算量兼顧性能與效率適用于雷達(dá)系統(tǒng)建模、波束形成與DOA聯(lián)合估計(jì)等研究場景。壓縮包共9個文件含8個核心MATLAB腳本如RD_MUSIC.m、steer_vector.m、angle_time_rd_music.m等實(shí)現(xiàn)導(dǎo)向矢量構(gòu)建、譜峰搜索與分辨率對比1個README.md提供運(yùn)行說明與算法原理簡述整體僅14KB輕量易部署。已有554人學(xué)習(xí)下載代碼經(jīng)實(shí)測可直接編譯運(yùn)行配套模塊化函數(shù)設(shè)計(jì)清晰便于理解算法流程、復(fù)現(xiàn)實(shí)驗(yàn)結(jié)果及開展進(jìn)一步優(yōu)化。 把頻率分集陣列FDA和MIMO放在一起做雷達(dá)參數(shù)估計(jì)這幾年在陣列信號處理圈子里討論熱度一直不低。原因很樸素傳統(tǒng)相控陣只對角度敏感波束指向跟目標(biāo)距離沒有直接關(guān)系想分辨同一角度上兩個不同距離的目標(biāo)通常得靠脈沖壓縮、大帶寬或者走多脈沖積累的路子。FDA通過給陣元加頻率增量讓波束指向隨距離變化相當(dāng)于在發(fā)射端把距離信息編碼進(jìn)導(dǎo)向矢量里再疊加上MIMO的正交波形分集把發(fā)射自由度解耦出來角度和距離就能在接收端做聯(lián)合估計(jì)。這篇文章我會直接把用MATLAB搭建FDA-MIMO參數(shù)估計(jì)仿真的完整過程講清楚從信號模型推導(dǎo)、二維MUSIC實(shí)現(xiàn)、蒙特卡洛性能評估到幾個我實(shí)際踩過的坑適合正在做陣列信號處理、雷達(dá)目標(biāo)參數(shù)估計(jì)或者剛接觸FDA-MIMO想快速上手復(fù)現(xiàn)的讀者。先說一個反直覺的結(jié)論這套系統(tǒng)看起來復(fù)雜核心代碼其實(shí)只有幾十行。難點(diǎn)不在MATLAB實(shí)現(xiàn)而在參數(shù)設(shè)計(jì)和結(jié)果分析——尤其是“怎么選中一個既不會距離模糊、又能把兩個目標(biāo)分辨開的Δf”。這恰恰是文獻(xiàn)里寫得最少、但工程里最致命的部分。這篇文章就把這些細(xì)節(jié)攤開講。1. 從“角度-距離耦合”說起FDA-MIMO解決了什么問題1.1 傳統(tǒng)相控陣的尷尬陣列流形對距離是“瞎”的要理解FDA-MIMO得先搞清楚傳統(tǒng)相控陣為什么測不了距離。相控陣的導(dǎo)向矢量里第n個陣元的相位差來自陣元位置差[ a_n(\theta) \exp\left(j\frac{2\pi}{\lambda} n d \sin\theta\right) ]這里的相位只和角度有關(guān)和目標(biāo)的絕對距離沒有任何關(guān)系。也就是說同一個角度上距離1公里和距離10公里的目標(biāo)陣列響應(yīng)完全相同。陣列流形在距離維是平的單快拍數(shù)據(jù)里根本提取不出距離信息。傳統(tǒng)雷達(dá)能測距靠的是發(fā)射波形的時間延遲也就是脈沖壓縮和帶寬而不是陣列流形。所以常規(guī)相控陣的陣列處理只能做角度估計(jì)距離維要靠波形設(shè)計(jì)來補(bǔ)。這句話值得重復(fù)一遍相控陣的角度估計(jì)能力來自陣列流形距離估計(jì)能力來自波形帶寬。FDA-MIMO的思路是把距離信息也塞進(jìn)陣列流形里讓“距離”和“角度”一道在空間譜上被估計(jì)出來。1.2 頻率分集陣列用頻率增量把距離編碼進(jìn)導(dǎo)向矢量FDA的做法很直接第n個發(fā)射陣元不再使用同一個載頻f0而是用f0加上一個頻率增量nΔf?？紤]到目標(biāo)距離R和方向θ發(fā)射陣元n到目標(biāo)的傳播相位大約為[ \phi_n -j\frac{2\pi (f_0 n\Delta f)}{c}(R - n d_t \sin\theta) ]展開后除了常規(guī)的角度相位項(xiàng)還會多出一項(xiàng)正比于nΔfR/c的相位。這項(xiàng)和距離R直接相關(guān)距離信息就這樣被編碼進(jìn)了發(fā)射導(dǎo)向矢量。用空間頻率的語言來描述發(fā)射陣元序號n上同時承載了兩個“頻率”一個由角度θ決定一個由距離R決定。但這里有個新問題單獨(dú)用FDA方向圖是“距離-角度耦合”的S型曲線。一個峰值對應(yīng)的是角度-距離平面上的一條曲線不是唯一確定的一個點(diǎn)。你看FDA的波束方向圖會發(fā)現(xiàn)波束指向會隨距離偏移這個偏移確實(shí)包含了距離信息可它和角度混在一起單靠一個總輸出信號很難解耦。所以FDA需要配合其他機(jī)制把角度和距離從耦合狀態(tài)中拆開——MIMO就是干這個的。1.3 MIMO補(bǔ)上自由度外積模型讓二維估計(jì)成為可能FDA-MIMO里發(fā)射端發(fā)射的是Nt路正交波形接收端通過匹配濾波把每路發(fā)射源單獨(dú)分離出來。一個目標(biāo)回波經(jīng)過匹配濾波后會形成一個Nt×Nr的復(fù)矩陣。這個矩陣恰好是發(fā)射導(dǎo)向矢量和接收導(dǎo)向矢量的外積[ \mathbf{Y} \beta \mathbf{a}_r(\theta) \mathbf{a}_t^T(\theta, R) \mathbf{N} ]把矩陣向量化后合成導(dǎo)向矢量就是發(fā)射導(dǎo)向和接收導(dǎo)向的Kronecker積。這個Kronecker積結(jié)構(gòu)非常關(guān)鍵它意味著接收數(shù)據(jù)里同時包含了發(fā)射端的距離編碼和接收端的角度信息兩者在同一個矢量里“各占一半”二維譜搜索就有了理論基礎(chǔ)。MIMO帶來的另一個好處是虛擬孔徑。Nt個發(fā)射陣元、Nr個接收陣元經(jīng)過正交波形分離后等效于一個Nt×Nr的虛擬均勻線陣陣列自由度從Nr提升到Nt×Nr。這個自由度的提升對MUSIC這種子空間類算法的性能至關(guān)重要。陣元數(shù)越多噪聲子空間約得越準(zhǔn)譜峰越尖銳。2. 信號模型搭建導(dǎo)向矢量里的距離項(xiàng)怎么推出來的2.1 發(fā)射導(dǎo)向矢量的完整推導(dǎo)我把推導(dǎo)過程完整寫一遍寫代碼的時候你直接對著用就行。設(shè)發(fā)射陣元間距為dt第n個發(fā)射陣元載頻為f0 nΔf其中n 0, 1, ..., Nt-1。目標(biāo)位于遠(yuǎn)場方向角為θ距離為R。以第一個陣元為參考第n個發(fā)射陣元到目標(biāo)的傳播距離近似為[ R_n \approx R - n d_t \sin\theta ]這里假設(shè)目標(biāo)距離遠(yuǎn)大于陣列孔徑平面波前成立。近場場景這部分模型就不適用后面會講到。回波從目標(biāo)回到第m個接收陣元的距離為[ R_m \approx R - m d_r \sin\theta ]發(fā)射陣元n和接收陣元m構(gòu)成的收發(fā)通道在匹配濾波后輸出的相位因子是[ \exp\left(-j\frac{2\pi(f_0 n\Delta f)}{c}(R - n d_t \sin\theta)\right) \exp\left(-j\frac{2\pi f_0}{c}(R - m d_r \sin\theta)\right) ]展開并整理忽略高階小項(xiàng)nΔf dt sinθ/c后發(fā)射導(dǎo)向矢量第n個元素可以寫成[ a_{t,n}(\theta, R) \exp\left(j\frac{2\pi f_0 d_t \sin\theta}{c} n - j\frac{2\pi \Delta f R}{c} n\right) ]接收導(dǎo)向矢量第m個元素為[ a_{r,m}(\theta) \exp\left(j\frac{2\pi f_0 d_r \sin\theta}{c} m\right) ]看到?jīng)]有發(fā)射導(dǎo)向矢量里有兩個關(guān)于n的線性相位項(xiàng)前一項(xiàng)由角度決定后一項(xiàng)由距離決定。這正是FDA-MIMO的核心特征——同一個陣元序號n上角度和距離各自貢獻(xiàn)了一個“空間頻率”。如果Δf 0退化成傳統(tǒng)相控陣MIMO距離項(xiàng)消失發(fā)射導(dǎo)向矢量就只剩角度信息二維估計(jì)自然做不了。2.2 為什么能忽略高階項(xiàng)一個數(shù)值驗(yàn)證推導(dǎo)中我把nΔf dt sinθ/c這一項(xiàng)扔了。做個數(shù)值驗(yàn)證看它到底小到什么程度。假設(shè)f0 10GHzΔf 60kHzNt 32dt λ/2 ≈ 0.015mθ 30°那么最大陣元序號n 31高階項(xiàng)最大值為[ \frac{n^2 \Delta f d_t \sin\theta}{c} \approx \frac{31^2 \times 60000 \times 0.015 \times 0.5}{3 \times 10^8} \approx 1.4 \times 10^{-3} \text{ rad} ]這個相位誤差遠(yuǎn)小于0.01弧度對導(dǎo)向矢量的影響完全可以忽略。但如果陣元數(shù)到幾百或者Δf到MHz量級這個項(xiàng)會漲得很快到時候就得保留完整表達(dá)式。做仿真的時候可以把它加上去對比一下通常你會發(fā)現(xiàn)譜峰偏移量小到可以忽略。2.3 數(shù)據(jù)模型與MUSIC譜的數(shù)學(xué)形式有了發(fā)射和接收導(dǎo)向矢量單個目標(biāo)的接收數(shù)據(jù)矩陣可以寫成[ \mathbf{Y} \beta \mathbf{a}_r(\theta) \mathbf{a}_t^T(\theta, R) \mathbf{N} ]把Y向量化[ \mathbf{x} \beta \cdot \text{vec}(\mathbf{Y}) \beta \left[\mathbf{a}_r(\theta) \otimes \mathbf{a}_t(\theta, R)\right] \mathbf{n} ]這個NtNr×1的向量就是單個目標(biāo)在虛擬陣列上的“快拍”。如果有D個目標(biāo)數(shù)據(jù)就是D個這樣向量的疊加。對K個快拍構(gòu)造協(xié)方差矩陣[ \mathbf{R}x \frac{1}{K}\sum{k1}^{K}\mathbf{x}_k \mathbf{x}_k^H ]特征分解后把特征值從小到大排列前NtNr-D個小特征值對應(yīng)的特征向量張成噪聲子空間En。由于信號導(dǎo)向矢量與噪聲子空間正交二維MUSIC譜定義為[ P(\theta, R) \frac{1}{\mathbf{a}^H(\theta, R) \mathbf{E}_n \mathbf{E}_n^H \mathbf{a}(\theta, R)} ]這里的a(θ,R)就是kron(a_r(θ), a_t(θ,R))。在角度-距離平面上掃描這個函數(shù)峰值位置就對應(yīng)目標(biāo)的θ和R估計(jì)值。有一點(diǎn)需要特別注意MUSIC的譜峰高度不是能量而是“正交程度”。旁瓣高度和主瓣高度的比值取決于噪聲子空間的估計(jì)質(zhì)量和信噪比直接相關(guān)。SNR一低特征值分裂不明顯噪聲子空間估計(jì)漂移譜峰就會變鈍甚至出現(xiàn)偽峰。后面性能評估部分你會看到這個現(xiàn)象的量化表現(xiàn)。3. 用MATLAB搭一套最小可跑的FDA-MIMO估計(jì)仿真3.1 參數(shù)設(shè)計(jì)頻率增量和距離分辨率的權(quán)衡仿真參數(shù)不是隨便定的尤其是Δf的選擇直接決定整個仿真的成敗。我用的參數(shù)如下參數(shù)取值說明載頻f010 GHz對應(yīng)波長λ 0.03m頻率增量Δf60 kHz無模糊距離c/Δf 5km發(fā)射陣元數(shù)Nt32距離維自由度接收陣元數(shù)Nr16角度維自由度陣元間距dt、drλ/2均勻線陣快拍數(shù)K100相參快拍數(shù)目標(biāo)參數(shù)θ 10°R 3.62km單目標(biāo)場景Δf的選擇有一個核心權(quán)衡不要忽視。FDA發(fā)射導(dǎo)向矢量里距離項(xiàng)的空間頻率是[ f_R \frac{\Delta f R}{c} ]距離維可以分辨的最小間隔反比于發(fā)射陣元數(shù)Nt和Δf的乘積[ \Delta R \approx \frac{c}{N_t \Delta f} ]代入?yún)?shù)c/(32×60000) ≈ 156米。也就是說這套參數(shù)下距離維的極限分辨率大約150米。如果兩個目標(biāo)距離差小于這個值MUSIC譜很難把它們分開。同時無模糊距離c/Δf 5km目標(biāo)距離必須落在這個范圍內(nèi)否則會出現(xiàn)距離折疊。這個權(quán)衡是FDA-MIMO仿真里最容易被忽視的點(diǎn)——Δf選大了距離不模糊范圍小選小了距離分辨率差兩者不可兼得。3.2 核心代碼導(dǎo)向矢量、數(shù)據(jù)生成、二維MUSIC譜搜索下面這段代碼我盡量保持精簡但保留了完整功能。只要MATLAB版本支持匿名函數(shù)和Kronecker積R2016b之后都能直接跑。% 參數(shù)初始化 clear; clc; close all; f0 10e9; % 載頻 10 GHz df 60e3; % 頻率增量 60 kHz c 3e8; lam c / f0; dt lam / 2; % 發(fā)射陣元間距 dr lam / 2; % 接收陣元間距 Nt 32; % 發(fā)射陣元數(shù) Nr 16; % 接收陣元數(shù) K 100; % 快拍數(shù) SNR_dB 10; snr 10^(SNR_dB/10); % 目標(biāo)設(shè)置 theta0 10 * pi/180; R0 3.62e3; beta 1; % 復(fù)幅度幅度為1 % 導(dǎo)向矢量匿名函數(shù) a_t (theta, R) exp(1j*2*pi*dt*f0*(0:Nt-1)*sin(theta)/c - 1j*2*pi*df*(0:Nt-1)*R/c); a_r (theta) exp(1j*2*pi*dr*f0*(0:Nr-1)*sin(theta)/c); a_steer (theta, R) kron(a_r(theta), a_t(theta, R)); % 生成單目標(biāo)接收數(shù)據(jù) A0 a_steer(theta0, R0); X zeros(Nt*Nr, K); for k 1:K noise sqrt(0.5/snr) * (randn(Nt*Nr,1) 1j*randn(Nt*Nr,1)); X(:,k) beta * A0 noise; end % 協(xié)方差矩陣與噪聲子空間 Rxx (X*X) / K; [Evec, Eval] eig(Rxx); [~, idx] sort(diag(Eval), ascend); D 1; % 目標(biāo)數(shù) En Evec(:, idx(1:end-D)); EnEn En * En; % 二維MUSIC譜搜索 theta_grid deg2rad(linspace(-20, 40, 601)); R_grid linspace(3.2e3, 4.2e3, 401); P_music zeros(length(theta_grid), length(R_grid)); for ii 1:length(theta_grid) for jj 1:length(R_grid) a_ij a_steer(theta_grid(ii), R_grid(jj)); P_music(ii,jj) 1 / real(a_ij * EnEn * a_ij); end end % 找峰值 [val, ind] max(P_music(:)); [I, J] ind2sub(size(P_music), ind); theta_est theta_grid(I) * 180/pi; R_est R_grid(J); fprintf(估計(jì)角度: %.2f°, 真實(shí)角度: %.2f°, 誤差: %.3f°\n, ... theta_est, theta0*180/pi, abs(theta_est - theta0*180/pi)); fprintf(估計(jì)距離: %.1fm, 真實(shí)距離: %.1fm, 誤差: %.1fm\n, ... R_est, R0, abs(R_est - R0));這段代碼有幾個細(xì)節(jié)值得說明。噪聲功率分配復(fù)噪聲總方差為σ2 1/snr實(shí)部和虛部分別是σ2/2。所以代碼里用sqrt(0.5/snr)乘randn和1j*randn這是很多新手容易搞錯的地方。如果直接把randn(1)當(dāng)作復(fù)噪聲實(shí)部信噪比會差3dB做蒙特卡洛時曲線會整體偏高等SNR才收斂。MUSIC譜取倒數(shù)前加real()a_ij * EnEn * a_ij理論上應(yīng)該是實(shí)數(shù)但數(shù)值計(jì)算會產(chǎn)生微小的虛部直接取倒數(shù)會把虛部噪聲放大。用real()劫住譜會更干凈。特征分解后對特征值排序MATLAB的eig函數(shù)默認(rèn)不保證特征值升序排列必須sort一次再把對應(yīng)特征向量重排。這里很多人會踩坑不排序直接取最后幾個特征向量結(jié)果就是噪聲子空間和信號子空間混在一起MUSIC譜一片平。3.3 兩階段搜索別傻乎乎全網(wǎng)格暴算上面的代碼在601×401的網(wǎng)格上做二維搜索兩層for循環(huán)跑下來在普通臺式機(jī)上通常要幾十秒到幾分鐘。蒙特卡洛一跑就是幾百次時間成本受不了。我建議用兩階段搜索第一階段用粗網(wǎng)格快速定位峰值鄰域。角度網(wǎng)格步長設(shè)0.5°距離網(wǎng)格步長設(shè)50米。這個計(jì)算量小能找到主峰的大致位置。第二階段在粗峰值的鄰域內(nèi)加密網(wǎng)格角度步長0.01°距離步長2米重新計(jì)算局部MUSIC譜并取最大值。這樣做計(jì)算量可以降一個量級。粗網(wǎng)格階段可能碰到旁瓣比主瓣大的情況但沒關(guān)系——只要你把第一階段返回的候選點(diǎn)取前3個峰值鄰域都做一次細(xì)搜索最后取全局最優(yōu)值。改成這樣的流程后單次二維搜索在普通筆記本上能壓縮到一兩秒蒙特卡洛循環(huán)200次也就幾分鐘。% 粗搜索 theta0_grid deg2rad(-20:0.5:40); R0_grid 3.2e3:50:4.2e3; % ... 同上述網(wǎng)格搜索代碼 ... % 找到最大峰值后在鄰域內(nèi)細(xì)搜索 theta_fine linspace(theta_peak-0.5*pi/180, theta_peak0.5*pi/180, 101); R_fine linspace(R_peak-50, R_peak50, 101); % 再次循環(huán)取max代碼不貼全了思路就這幾行。優(yōu)先保證先把單目標(biāo)跑通再考慮優(yōu)化速度。4. 仿真結(jié)果分析譜峰、RMSE曲線與CRB對照4.1 單目標(biāo)場景二維MUSIC譜峰定位跑完上面的代碼單目標(biāo)場景下MUSIC譜會在角度約10°、距離約3.62km處形成一個尖銳的峰值。用3.2節(jié)的參數(shù)一次隨機(jī)實(shí)驗(yàn)估計(jì)出來的典型結(jié)果是θ_est ≈ 10.03°R_est ≈ 3628米。角度誤差不到0.1°距離誤差約10米。這個距離誤差遠(yuǎn)小于前面說的150米分辨率極限為什么因?yàn)閱文繕?biāo)情況下MUSIC估計(jì)的是參數(shù)值不是分辨率極限。分辨率極限說的是兩個目標(biāo)能不能被分開參數(shù)估計(jì)精度則是信噪比、陣列孔徑、快拍數(shù)共同決定的。這里距離估計(jì)精度能達(dá)到十幾米是因?yàn)镸USIC的譜峰插值效應(yīng)和陣元數(shù)帶來的處理增益。MUSIC譜圖上看峰值附近在距離維的3dB寬度大約就是幾十米量級。如果信噪比降到0dB你可以看到主瓣變寬、旁瓣抬高峰值位置開始隨機(jī)漂移。這是正?，F(xiàn)象不是程序?qū)戝e了。4.2 雙目標(biāo)場景同角度不同距離的分辨能力把目標(biāo)改成兩個θ 10°R1 3.62kmθ 10°R2 3.88km。兩個目標(biāo)距離差260米大于極限分辨率約150米理論上可以分辨。直接把數(shù)據(jù)生成改成兩個導(dǎo)向矢量的疊加theta1 10*pi/180; R1 3.62e3; theta2 10*pi/180; R2 3.88e3; A1 a_steer(theta1, R1); A2 a_steer(theta2, R2); for k 1:K noise sqrt(0.5/snr) * (randn(Nt*Nr,1) 1j*randn(Nt*Nr,1)); X(:,k) beta*A1 beta*A2 noise; end D 2;MUSIC譜會出現(xiàn)兩個相鄰的峰間距260米峰值同時清晰可見。如果你把Δf從60kHz降到30kHz距離維分辨率會退化到約312米這兩個目標(biāo)就會合并成一個寬峰無法分辨。這個對比非常直觀FDA-MIMO距離分辨能力確實(shí)由NtΔf決定不是由信號帶寬決定。這和傳統(tǒng)雷達(dá)的帶寬決定距離分辨率的直覺完全不同剛接觸的人很容易搞混。還有一個細(xì)節(jié)雙目標(biāo)時如果兩個目標(biāo)角度相同只改變它們的距離你會發(fā)現(xiàn)在MUSIC譜上兩個峰的角坐標(biāo)完全相同距離坐標(biāo)不同。這正好驗(yàn)證了前文說的“發(fā)射維承載距離信息”的設(shè)計(jì)。角度相同但距離不同的目標(biāo)在傳統(tǒng)相控陣MIMO下是根本不可分辨的因?yàn)樗鼈兊慕嵌仍诮邮站S完全重合。4.3 蒙特卡洛統(tǒng)計(jì)RMSE曲線與CRB的差距單次實(shí)驗(yàn)結(jié)果沒有統(tǒng)計(jì)意義。做200次蒙特卡洛對SNR從0dB到25dB遍歷計(jì)算角度和距離的RMSEL 200; SNR_list 0:5:25; rmse_theta zeros(size(SNR_list)); rmse_R zeros(size(SNR_list)); for s 1:length(SNR_list) snr 10^(SNR_list(s)/10); err_theta zeros(1,L); err_R zeros(1,L); for l 1:L % 生成數(shù)據(jù)MUSIC估計(jì)省略中間代碼 err_theta(l) (theta_est - theta0*180/pi)^2; err_R(l) (R_est - R0)^2; end rmse_theta(s) sqrt(mean(err_theta)); rmse_R(s) sqrt(mean(err_R)); end典型結(jié)果是SNR在10dB以上時角度RMSE大約0.01°到0.1°量級距離RMSE大約5到20米。距離估計(jì)精度隨SNR提升非?？旖咏麮RB曲線。SNR低于5dB時RMSE會顯著背離CRB出現(xiàn)“閾值效應(yīng)”——一部分實(shí)驗(yàn)里MUSIC譜峰值找錯了位置誤差從幾米變成幾百米。CRB的計(jì)算對驗(yàn)證算法性能非常有用。對確定性模型用數(shù)值差分計(jì)算導(dǎo)向矢量對θ和R的導(dǎo)數(shù)h 1e-6; dtheta (a_steer(theta0h, R0) - a_steer(theta0-h, R0)) / (2*h); dR (a_steer(theta0, R0h) - a_steer(theta0, R0-h)) / (2*h); D0 [dtheta, dR]; PA_perp eye(Nt*Nr) - A0/(A0*A0)*A0; FIM 2 * snr * real(D0 * PA_perp * D0); CRB diag(inv(FIM));CRB(1)是角度方差的理論下限單位是rad2CRB(2)是距離方差的理論下限單位是m2。把RMSE曲線和sqrt(CRB)畫在同一張圖上。你會發(fā)現(xiàn)一個規(guī)律高SNR時MUSIC的RMSE線比CRB線高3到6dB這是子空間類方法的典型水平——它能逼近CRB但不會達(dá)到CRB。SNR低過閾值后RMSE急劇抬高這是所有搜索類算法的通病。5. 仿真過程中的坑與邊界條件能復(fù)現(xiàn)不等于能用5.1 距離模糊陷阱Δf一大目標(biāo)會“折回”到鏡像距離FDA導(dǎo)向矢量里距離項(xiàng)是以2π為周期折疊的本文還有配套的精品資源點(diǎn)擊獲取