基于pri變換的單通道多分量sfm信號(hào)參數(shù)估計(jì)方法
【專利摘要】本發(fā)明請(qǐng)求保護(hù)一種基于PRI變換的單通道多分量SFM信號(hào)參數(shù)估計(jì)方法，屬于信號(hào)處理領(lǐng)域技術(shù)。本方法對(duì)接收端離散信號(hào)進(jìn)行FFT變換，提取頻譜中有效譜線，轉(zhuǎn)換為脈沖序列后，根據(jù)PRI變換估計(jì)得到接收端信號(hào)中信號(hào)分量個(gè)數(shù)與各分量信號(hào)的調(diào)制頻率；然后調(diào)整離散SFM基函數(shù)，再對(duì)接收信號(hào)與調(diào)整后離散SFM基函數(shù)的乘積進(jìn)行FFT變換；通過(guò)峰值搜索估計(jì)得到各分量信號(hào)的載波頻率和相應(yīng)調(diào)制系數(shù)，最后構(gòu)造相應(yīng)分量信號(hào)與接收端信號(hào)共軛相乘得到其幅度估計(jì)。本方法相對(duì)計(jì)算量較低，在低信噪比環(huán)境下能夠精確的估計(jì)出多分量SFM信號(hào)的信號(hào)分量個(gè)數(shù)、各分量信號(hào)的調(diào)制頻率、載波頻率、調(diào)制系數(shù)和幅度，因此具有良好的應(yīng)用前景。
【專利說(shuō)明】
基于PR I變換的單通道多分量SFM信號(hào)參數(shù)估計(jì)方法
技術(shù)領(lǐng)域
[0001]本發(fā)明涉及雷達(dá)通信信號(hào)處理，具體為一種基于脈沖重復(fù)間隔變換（Pulse repetition intervals，PRI)的單通道多分量SFM信號(hào)參數(shù)估計(jì)方法。
【背景技術(shù)】
[0002] 正弦調(diào)頻信號(hào)（Sinusoidal Frequency Modulation)是一種特殊的非平穩(wěn)時(shí)頻信 號(hào)，它的頻率隨時(shí)間呈正弦變化。因?yàn)樾盘?hào)的截獲率低，在雷達(dá)目標(biāo)檢測(cè)、近程探測(cè)、引信抗 干擾等領(lǐng)域得到廣泛應(yīng)用。所以在這些領(lǐng)域中單通道多分量正弦調(diào)頻信號(hào)的參數(shù)估計(jì)的方 法研究具有重要的意義。
[0003] 時(shí)頻分析方法用于非平穩(wěn)信號(hào)參數(shù)的估計(jì)，最常見(jiàn)的是把Wigner-Ville分布 (WVD)和陣列信號(hào)處理結(jié)合起來(lái)，然而在多分量信號(hào)的情況下這種方法的計(jì)算非常復(fù)雜，并 且對(duì)信號(hào)的采樣率要求很高，還存在交叉項(xiàng)干擾問(wèn)題，這些都導(dǎo)致此類方法準(zhǔn)確性、實(shí)用性 的降低，局限性大。黃浩等提出基于循環(huán)自相關(guān)的方法估計(jì)SFM信號(hào)參數(shù)，但是該方法只適 用于單分量SFM且抗噪性能較差。熊輝等人提出基于卡森準(zhǔn)則的SFM信號(hào)參數(shù)估計(jì)方法也只 適用于單分量SFM且抗噪性能較差。陳晶等提出基于離散正弦調(diào)頻變換的多分量正弦調(diào)頻 信號(hào)參數(shù)估計(jì)方法，適用于多分量基帶SFM信號(hào)，但是由于其相關(guān)參數(shù)正整數(shù)設(shè)置影響了參 數(shù)估計(jì)的精度，對(duì)抗噪性能也有一定的影響。朱航等提出基于改進(jìn)自適應(yīng)分解法的單通道 雷達(dá)引信混合信號(hào)分離，該方法將五參量的chirplet信號(hào)作為基函數(shù)，對(duì)混合信號(hào)進(jìn)行自 適應(yīng)分解，但是參量增加帶來(lái)巨大的計(jì)算量且適用于較高信噪比環(huán)境。
[0004] 因此本發(fā)明提出基于PRI變換的單通道多分量SFM信號(hào)參數(shù)估計(jì)方法。

【發(fā)明內(nèi)容】

[0005] 本發(fā)明所要解決的技術(shù)問(wèn)題，針對(duì)現(xiàn)有技術(shù)存在的單通道多分量SFM信號(hào)參數(shù)估 計(jì)中計(jì)算量大，低信噪比下估計(jì)性能差等缺陷，提出一種基于PRI變換的估計(jì)方法，解決這 一難題。該方法相對(duì)計(jì)算量較低，在低信噪比環(huán)境下能夠精確的估計(jì)出多分量SFM信號(hào)的 信號(hào)分量個(gè)數(shù)、各分量信號(hào)的調(diào)制頻率、載波頻率、調(diào)制系數(shù)和幅度。
[0006] 本發(fā)明解決上述技術(shù)問(wèn)題的技術(shù)方案是:一種基于PRI變換的估計(jì)方法，其步驟在 于，對(duì)接收端離散信號(hào)進(jìn)行FFT變換，提取頻譜中有效譜線，轉(zhuǎn)換為脈沖序列后，根據(jù)PRI變 換估計(jì)得到接收端信號(hào)中信號(hào)分量個(gè)數(shù)與各分量信號(hào)的調(diào)制頻率;然后調(diào)整離散SFM基函 數(shù)，再對(duì)接收信號(hào)與調(diào)整后離散SFM基函數(shù)的乘積進(jìn)行FFT變換;通過(guò)峰值搜索估計(jì)得到各 分量信號(hào)的載波頻率和相應(yīng)調(diào)制系數(shù)，最后構(gòu)造相應(yīng)分量信號(hào)與接收端信號(hào)相乘得到對(duì)應(yīng) 幅度估計(jì)。
【附圖說(shuō)明】
[0007] 圖1本發(fā)明多分量SFM信號(hào)參數(shù)估計(jì)方法流程圖；
[0008] 圖2本發(fā)明經(jīng)50次頻域累加平均后混合信號(hào)頻譜圖及譜線門限；
[0009] 圖3本發(fā)明消除干擾后經(jīng)PRI變換估計(jì)調(diào)制頻率圖；
[0010] 圖4本發(fā)明SFM信號(hào)經(jīng)調(diào)整SFM基匹配后再進(jìn)行FFT變換的三維圖；
[0011] 圖5本發(fā)明SFM信號(hào)經(jīng)調(diào)整SFM基匹配后再進(jìn)行FFT變換的三維圖；
[0012] 圖6本發(fā)明SFM信號(hào)經(jīng)調(diào)整SFM基匹配后再進(jìn)行FFT變換的三維圖；
[0013] 圖7本發(fā)明不同累加次數(shù)不同信噪比條件下信號(hào)分量調(diào)制頻率檢測(cè)估計(jì)性能圖；
[0014] 圖8本發(fā)明載波頻率估計(jì)的性能圖；
[0015] 圖9本發(fā)明調(diào)制系數(shù)估計(jì)的性能圖；
[0016] 圖10本發(fā)明幅度估計(jì)的性能圖；
【具體實(shí)施方式】
[0017] 以下結(jié)合附圖和具體實(shí)例，對(duì)本發(fā)明的實(shí)施作進(jìn)一步的描述。
[0018]
[0019]單通道信號(hào)的數(shù)學(xué)模型如式（1):
[0020] Y=AS+ff (1)
[0021] 其中，Y是I XN的向量，表示接收端觀測(cè)信號(hào)，N表示信號(hào)長(zhǎng)度;A是I XM的向量，表 示混合系數(shù)，M是源信號(hào)個(gè)數(shù);S是MX N的矩陣;W也是I X N向量，表示高斯白噪聲。
[0022]本文針對(duì)的是多分量SFM信號(hào)，對(duì)于高斯白噪聲環(huán)境中Q分量SFM信號(hào)混合模型如 式⑵：
[0023] (2)
[0024]
[0025] (3)
[0026] 其中y[n]=y[n · Ts]，Sq[n] = Sq[n · Ts]，w[n]=w[n · Ts]，TS為采樣周期;Aq、fCq、 fmq、mfj別對(duì)應(yīng)第q(q=l，2, . . .Q)個(gè)分量信號(hào)的幅度、載波頻率、調(diào)制頻率和調(diào)頻系數(shù)。
[0027] 脈沖重復(fù)間隔變換法是估計(jì)雷達(dá)脈沖序列重復(fù)間隔的有效方法，能抑制信號(hào)的脈 沖重復(fù)間隔及其整數(shù)倍同時(shí)存在的現(xiàn)象。令tP，p = 0，1，…，P-I為脈沖到達(dá)時(shí)間，其中P是采 樣脈沖數(shù)，如果只考慮到達(dá)時(shí)間唯一參數(shù)，采樣脈沖串就可以模型化為單位沖激函數(shù)的和， 即
[0028]
(4)
[0029]其中δ( ·)是Dirac函數(shù)。
[0030] g(t)的積分變換公式為：
[0031；
(5)：
[0032]其中τ>0，該算法稱為PRI變換。因?yàn)镮 D( τ) I給出了一種PRIS的譜圖，在代表真PRI 值的地方將出現(xiàn)峰值。將式(4)代入式(5 )，可得
[0033；
(6)
[0034] 因?yàn)镾FM信號(hào)頻譜上的特性，這里PRI變換運(yùn)用于SFM信號(hào)調(diào)制頻率的估計(jì)。具體理 論分析如下：
[0039] (8)
[0035] 對(duì)其分量SFM信號(hào)，由雅可bh屏開(kāi)忒可知，
[0036]
[0037]
[0038]則該分量信號(hào)的傅里葉變換為：
[0040]
[0041] (9)
[0042] 其頻譜相當(dāng)于多路間隔為各自調(diào)制頻率的譜線混合在一起。由此可提取信號(hào)頻譜 中的譜線，對(duì)其歸一化后，轉(zhuǎn)化為脈沖序列，即可應(yīng)用PRI變換估計(jì)其調(diào)制頻率。
[0043]但是對(duì)于高斯白噪聲環(huán)境下的多分量正弦調(diào)頻信號(hào)，其頻譜受到噪聲的影響，直 接提取相關(guān)信息進(jìn)行PRI轉(zhuǎn)換不能準(zhǔn)確的估計(jì)調(diào)制頻率。所以為了消除噪聲的影響，需要設(shè) 置一個(gè)合理門限，提取門限值以上的譜線，門限的設(shè)置既要保證譜線的數(shù)量足夠又要確定 沒(méi)有噪聲的干擾。因?yàn)樵肼暤碾S機(jī)性，接收端混合信號(hào)的頻譜不穩(wěn)定，這為門限的設(shè)置帶 來(lái)了巨大的障礙，所以這里采用頻域累加平均的方式降低噪聲隨機(jī)性在頻域的干擾。
[0044] 而對(duì)于門限選取，根據(jù)頻譜中特征兩端信號(hào)的譜線數(shù)目較少，所以可以選擇前段 后端多點(diǎn)求取平均值作為門限值，考慮到一些突發(fā)點(diǎn)的影響，在本文中門限值選取為
[0045]
(10)
[0046] 其中1.1 1.4,21為求取平均值的總點(diǎn)數(shù)，YY = abs(FFT(Y))。
[0047] 雖然混合信號(hào)中SFM對(duì)應(yīng)的PRI值較大，但是還是存在其他的干擾，所以需要設(shè)置 門限，這里門限設(shè)置為：
[0048] thred2=max( ID(x) I )/2 (11)
[0049] 如此即可應(yīng)用PRI變換獲得各分量信號(hào)調(diào)制頻率的估計(jì)。
[0050] 在獲得調(diào)制頻率的估計(jì)之后，接下來(lái)對(duì)多分量SFM信號(hào)的調(diào)制系數(shù)、載頻進(jìn)行估 計(jì)。對(duì)于1^點(diǎn)離散SFM信號(hào)，其DSFMT(Discrete Sinusoidal Frequence Modulation Transform，DSFMT)變換定義為
[0051 ]
(12)
[0052]其中0 < i < N-I，0〈k < L-I，0〈1 < L-I，且k和I為整數(shù)。
[0053]由DSFMT的定義可知，SFM信號(hào)的DSFMT變換在信號(hào)與離散正弦調(diào)頻基函數(shù)最為匹 配，即調(diào)制頻率與調(diào)制頻偏匹配時(shí)，信號(hào)與離散正弦調(diào)頻基函數(shù)的能量幅度取得最大值，此 時(shí)可以通過(guò)譜峰搜索得到調(diào)制頻率與調(diào)制頻偏的估計(jì)，即可轉(zhuǎn)化得到本文中調(diào)制頻率與調(diào) 制系數(shù)的估計(jì)。但是k，l的整數(shù)取值極大的限制了參數(shù)估計(jì)的精度。所以引入兩個(gè)因子a、i3， 其取值均為正數(shù)，有V =ak，V =m，aj的取值影響參數(shù)估計(jì)的精度和運(yùn)算量，可根據(jù)計(jì)算 量和參數(shù)估計(jì)精度的需求調(diào)整。
[0054] 結(jié)合本文的信號(hào)模型表達(dá)形式，將接收端離散信號(hào)與離散正弦調(diào)頻基函數(shù)相乘：
[0055]
[0056] 其中mf，fm分別為調(diào)制系數(shù)與調(diào)制頻率，其取值精度可通過(guò)兩個(gè)因子放縮。
[0057] 利用PRI獲得了各分量信號(hào)調(diào)制頻率的估計(jì)/胃，對(duì)應(yīng)調(diào)制頻率其估計(jì)誤差較小時(shí)
[0058]
(14)
[0059] 由此可以調(diào)整SFM基函數(shù)，只需對(duì)調(diào)制系數(shù)進(jìn)行搜索，進(jìn)一步降低計(jì)算量，DSFMy則 可表示為：
[0060]
(15)
[0061] 利用PRI變換獲得了各分量信號(hào)調(diào)制頻率的估計(jì)九7后，只需分量信號(hào)中調(diào)制系數(shù) mfq與基函數(shù)中調(diào)制系數(shù)Hlf最為匹配，對(duì)于某分量信號(hào)調(diào)制系數(shù)Hlfq與基函數(shù)中調(diào)制系數(shù)Hlf最 為匹配時(shí)，DSFMy可化簡(jiǎn)為：
[0062] (16) ci.=<j.q^e
[0063] 對(duì)式（16)進(jìn)行FFT變換，可知在其頻譜內(nèi)出現(xiàn)最大值而所對(duì)應(yīng)位置即為載波頻率 的估計(jì)I，且心=，，即
[0064]
(17)
[0065]所以將接收端離散信號(hào)與離散正弦調(diào)頻基函數(shù)乘積的結(jié)果進(jìn)行FFT變換，對(duì)FFT變 換的結(jié)果，搜索其最大峰值，即可獲得相應(yīng)分量的調(diào)制系數(shù)mfq和載波頻率&。
[0066] 在得到各分量信號(hào)的調(diào)制頻率/_、載波頻率&、調(diào)制系數(shù)"\的估計(jì)后，即可重構(gòu) 相應(yīng)分量
與混合信號(hào)共輒相乘，對(duì)各分量對(duì)應(yīng)幅度Aq進(jìn)行估
計(jì)：
[0067] (IB)
[0068] 利用仿真實(shí)驗(yàn)對(duì)本發(fā)明算法的理論推導(dǎo)進(jìn)行驗(yàn)證，實(shí)驗(yàn)參數(shù)設(shè)置:在_5dB高斯白 噪聲環(huán)境下，信號(hào)長(zhǎng)度為1024,米樣頻率為2048。三分量SFM信號(hào)，各分量信號(hào)幅度均為1，其 他參數(shù):SFMl載波頻率240Hz，調(diào)制系數(shù)5.8，調(diào)制頻率9. OHz ; SFM2載波頻率256Hz，調(diào)制系數(shù) 5.1，調(diào)制頻率19. OHz; SFM3載波頻率220Hz，調(diào)制系數(shù)6.9，調(diào)制頻率12. OHz。
[0069] 對(duì)接收端混合信號(hào)進(jìn)行50次頻域累加平均，其累加后頻譜如圖2。提取門限值以上 的譜線，轉(zhuǎn)化為脈沖序列后，進(jìn)行PRI變換，估計(jì)的調(diào)制頻率參數(shù)圖如圖3;將接收端混合信 號(hào)分解在調(diào)整后的SFM基函數(shù)后再進(jìn)行FFT變換，通過(guò)譜峰搜索獲得載頻、調(diào)制系數(shù)估計(jì)，三 分量SFM信號(hào)相應(yīng)三維圖如圖4、5、6。
[0070] 此時(shí)放縮因子為0.1，根據(jù)峰值搜索得到各分量信號(hào)的載頻和調(diào)制系數(shù)估計(jì)，在最大 值處可以得到分量SFMl信號(hào)的載頻為240、調(diào)制系數(shù)為5.8，分量SFM2信號(hào)的載頻為256，調(diào) 制系數(shù)為5.0,分量SFM3信號(hào)的載頻為220,調(diào)制系數(shù)為6.9。由此即可重構(gòu)相應(yīng)分量信號(hào)如：
分 另IJ與接收端信號(hào)共輒相乘，獲得相應(yīng)幅值A(chǔ)i、A2、A3估計(jì)。
[0071] 綜合實(shí)驗(yàn)數(shù)據(jù)，-5dB高斯白噪聲環(huán)境下估計(jì)的三分量SFM信號(hào)參數(shù)SFMl信號(hào)的幅 度為0.977，載頻為240、調(diào)制系數(shù)為5.8，調(diào)制頻率9.01;分量SFM2信號(hào)的幅度為0.949，載頻 為256，調(diào)制系數(shù)為5.0，調(diào)制頻率19.01;分量SFM3信號(hào)的幅度為1.07，載頻為220，調(diào)制系數(shù) 為6.9,調(diào)制頻率12.01。經(jīng)相似系數(shù)計(jì)算可得原信號(hào)和估計(jì)信號(hào)的相似系數(shù)分別Sp 14 = 0.9930、P25 = 0.9922、P36 = 0.99 6 2。
[0072] 當(dāng)不累加時(shí)由于噪聲隨機(jī)性的影響，對(duì)門限的設(shè)置帶來(lái)了極大的干擾，只能在較 高性噪比條件下實(shí)現(xiàn)SFM信號(hào)參數(shù)的估計(jì)。這里分別對(duì)累加次數(shù)為1=10、20、30、50，信噪比 從-15dB到OdB間隔IdB變化的情況做500次蒙特卡洛仿真，調(diào)制頻率檢測(cè)并估計(jì)正確率如圖 7。由圖7可知，累加次數(shù)的增加在一定程度上能夠改善估計(jì)的正確率;隨著信噪比的增加， 對(duì)調(diào)制頻率檢測(cè)、估計(jì)的正確率也越來(lái)越高，直至最佳。在J = 2 0時(shí)，信噪比為-6 d B的情況 下，可以檢測(cè)并準(zhǔn)確的估計(jì)各分量調(diào)制頻率。
[0073] 考慮調(diào)制頻率的估計(jì)誤差，各分量信號(hào)載頻、調(diào)制系數(shù)、幅度估計(jì)性能圖如圖8、9、 10。由圖8、9、10可知隨著信噪比的增加，對(duì)載頻、調(diào)制系數(shù)、幅度的估計(jì)性能越來(lái)越好，各分 量載頻估計(jì)均方誤差可以忽略不計(jì)，在信噪比為_(kāi)5dB時(shí)各分量調(diào)制頻率估計(jì)、幅度估計(jì)的 均方誤差都在1〇_ 3量級(jí)，誤差較小即可以有效的估計(jì)各分量信號(hào)的載頻、調(diào)制系數(shù)、幅度。
【主權(quán)項(xiàng)】
1. 一種基于PRI變換的單通道多分量SFM信號(hào)參數(shù)估計(jì)方法，其步驟在于，對(duì)接收端離 散信號(hào)進(jìn)行FFT變換，提取頻譜中有效譜線，轉(zhuǎn)換為脈沖序列后，根據(jù)PRI變換估計(jì)得到接收 端信號(hào)中信號(hào)分量個(gè)數(shù)與各分量信號(hào)的調(diào)制頻率;然后調(diào)整離散SFM基函數(shù)，再對(duì)接收信號(hào) 與調(diào)整后離散SFM基函數(shù)的乘積進(jìn)行FFT變換;通過(guò)峰值捜索估計(jì)得到各分量信號(hào)的載波頻 率和相應(yīng)調(diào)制系數(shù)，最后構(gòu)造相應(yīng)分量信號(hào)與接收端信號(hào)共輛相乘得到其幅度估計(jì)。2. 根據(jù)權(quán)利要求1所述的估計(jì)方法，其特征在于對(duì)接收信號(hào)進(jìn)行傅里葉變換，提取并轉(zhuǎn) 化處理得到含有調(diào)制頻率信息的脈沖序列;經(jīng)PRI變換得到接收信號(hào)中的分量個(gè)數(shù)、各個(gè)調(diào) 制頻率參數(shù);獲得分量信號(hào)的調(diào)制頻率參數(shù)，調(diào)整離散SFM基函數(shù)，再對(duì)接收信號(hào)與調(diào)整后 離散SFM基函數(shù)的乘積進(jìn)行FFT變換，通過(guò)最大峰值捜索獲得相應(yīng)載波頻率和調(diào)制系數(shù)的估 計(jì)，重構(gòu)相應(yīng)分量信號(hào)哪以2'了 4" +抑,戶(站么"；)]與接收端信號(hào)共輛相乘獲得相應(yīng)幅值估 計(jì)。
【文檔編號(hào)】G01R23/16GK105842534SQ201610151338
【公開(kāi)日】2016年8月10日
【申請(qǐng)日】2016年3月16日
【發(fā)明人】張?zhí)祢U, 廖暢, 葉飛, 張剛, 羅忠濤
【申請(qǐng)人】重慶郵電大學(xué)