本發(fā)明屬于通信技術(shù)和數(shù)字信號(hào)處理領(lǐng)域，具體涉及一種基于循環(huán)譜的無線通信信號(hào)調(diào)制識(shí)別方法以及一種基于循環(huán)譜的無線通信信號(hào)調(diào)制識(shí)別與參數(shù)估計(jì)方法。

背景技術(shù)：

常用的無線通信信號(hào)包括2ASK(2Amplitude Shift Keying二進(jìn)制幅度鍵控)、2FSK(2Frequency Shift Keying二進(jìn)制頻移鍵控)、BPSK(Binary Phase Shift Keying二進(jìn)制相移鍵控)、QPSK(Quadrature Phase Shift Keying正交相移鍵控)、MSK(Minimum Shift Keying最小頻移鍵控)、16QAM(16Quadrature Amplitude Modulation16正交幅度調(diào)制)信號(hào)。在很多非協(xié)作通信中，譬如單點(diǎn)對多點(diǎn)網(wǎng)絡(luò)通信、無線電監(jiān)控與管理、通信系統(tǒng)故障檢測、認(rèn)知無線電等領(lǐng)域，甚至于軍事上的通信偵查電子對抗等相關(guān)領(lǐng)域，接收端在對信號(hào)進(jìn)行處理的過程中總要面臨著一項(xiàng)或多項(xiàng)未知信息，甚至根本不知道任何相關(guān)發(fā)射信息。在這種情況下，所接收信號(hào)的各種參數(shù)對于接收機(jī)都是未知的，接收機(jī)不能采用數(shù)據(jù)輔助算法來對信號(hào)的調(diào)制方式進(jìn)行識(shí)別和參數(shù)進(jìn)行估計(jì)。

目前國內(nèi)外已發(fā)表的有關(guān)調(diào)制識(shí)別與參數(shù)估計(jì)的文獻(xiàn)中，對信號(hào)調(diào)制識(shí)別的算法主要分類兩類，一類是基于最大似然決策理論的方法，將調(diào)制識(shí)別過程作為一個(gè)復(fù)合假設(shè)檢驗(yàn)問題，利用最大似然準(zhǔn)則進(jìn)行判決。該方法識(shí)別范圍小，計(jì)算復(fù)雜，不易實(shí)現(xiàn)。第二類是基于信號(hào)特征提取與分析的識(shí)別方法，這類方法包括三個(gè)步驟：信號(hào)預(yù)處理、特征提取和分類識(shí)別，根據(jù)所選用的特征不同，調(diào)制識(shí)別也分為不同的方向，主要包括瞬時(shí)特征、統(tǒng)計(jì)特征以及變換域特征等。

同時(shí)信號(hào)參數(shù)估計(jì)主要包括兩部分，分別為載頻估計(jì)和符號(hào)率估計(jì)。對于載頻的估計(jì)，在協(xié)作通信中常通過在信號(hào)中添加導(dǎo)頻或者訓(xùn)練序列來完成，而在非協(xié)作通信中，利用非線性變換來對載頻進(jìn)行估計(jì)是一種常見的方法。使用非線性變換方法使得這類載頻估計(jì)受信噪比影響較大，在低信噪比環(huán)境下性能較差。在頻譜監(jiān)測中，經(jīng)常通過對接收信號(hào)的頻譜進(jìn)行分析得到載頻信息，常見的有Welch周期圖法，這類方法可以對載頻值進(jìn)行粗略估計(jì)。在對信號(hào)的符號(hào)率估計(jì)上，常見的有基于包絡(luò)譜，小波變換，循環(huán)自相關(guān)和循環(huán)譜的符號(hào)速率估計(jì)方法?？偟膩碚f，對信號(hào)調(diào)制識(shí)別與參數(shù)估計(jì)的研究一般劃分為兩個(gè)獨(dú)立的研究方向，但隨著非協(xié)作通信的發(fā)展，調(diào)制識(shí)別與參數(shù)估計(jì)越來越發(fā)展為一個(gè)整體，需要一種聯(lián)合方法能夠?qū)ξ粗ㄐ判盘?hào)的調(diào)制識(shí)別與參數(shù)估計(jì)都能夠完成。

技術(shù)實(shí)現(xiàn)要素：

發(fā)明目的：針對上述現(xiàn)有技術(shù)存在的問題與不足，本發(fā)明提出了一種基于循環(huán)譜的無線通信信號(hào)調(diào)制識(shí)別及參數(shù)估計(jì)方法，能夠同時(shí)完成對未知信號(hào)的調(diào)制識(shí)別與參數(shù)估計(jì)的方法。

技術(shù)方案：為實(shí)現(xiàn)上述發(fā)明目的，本發(fā)明中基于循環(huán)譜的無線通信信號(hào)調(diào)制識(shí)別方法，其特征在于，包括以下步驟：

1)信號(hào)帶寬粗估計(jì)：接收端在接收未知信號(hào)后，對該未知信號(hào)的帶寬進(jìn)行粗略估計(jì)，得到帶寬估計(jì)值和中心頻率；

2)帶通濾波：根據(jù)得到的信號(hào)中心頻率和帶寬估計(jì)值，對其進(jìn)行帶通濾波；

3)計(jì)算信號(hào)循環(huán)譜：利用快速傅里葉變換累加算法計(jì)算信號(hào)的循環(huán)譜；

4)提取特征參數(shù)：根據(jù)計(jì)算得到的循環(huán)譜在其不同截面的譜峰值特點(diǎn)上提取特征參數(shù)，這些特征參數(shù)包括：截面峰值個(gè)數(shù)NR、截面譜峰下降度D、截面最大值M、截面最大值與次大值比R₁和截面最大值與截面最大值比值R₂；

5)調(diào)制方式識(shí)別：根據(jù)得到的信號(hào)特征參數(shù)，采用基于二叉樹分類方式對未知信號(hào)進(jìn)行調(diào)制識(shí)別，二叉樹分類方式以預(yù)先設(shè)定的閾值作為比較判決的依據(jù)。

其中，步驟3)中計(jì)算信號(hào)循環(huán)譜包括以下步驟：

將帶通濾波后的數(shù)據(jù)x(n)每隔L點(diǎn)進(jìn)行一次抽取，每次抽取的數(shù)據(jù)長度記為N′，總的抽取次數(shù)為P，對每次抽取的數(shù)據(jù)都要進(jìn)行加窗操作以減小頻譜泄露，再對加窗后的每塊數(shù)據(jù)進(jìn)行N'點(diǎn)的FFT變換，傅里葉變換的輸出為對應(yīng)頻率成分的復(fù)包絡(luò)；

將復(fù)信號(hào)向基帶進(jìn)行頻移，頻移距離為k/N′，并求出頻移后復(fù)包絡(luò)的共軛相關(guān)，將結(jié)果逐行放入大小為N′²的乘法矩陣中，通過不斷重復(fù)得到一個(gè)P×N′²大小的矩陣；

對得到的矩陣求P點(diǎn)FFT變換，對FFT后的結(jié)果進(jìn)行坐標(biāo)變換得到最終的信號(hào)循環(huán)譜。

其中，步驟5)中能識(shí)別的無線通信信號(hào)的調(diào)制方式包括：2ASK、2FSK、BPSK、QPSK、MSK或16QAM。

進(jìn)一步地，步驟5)中二叉樹分類方式具體為：比較截面峰值個(gè)數(shù)NR與第一預(yù)設(shè)閾值Thr1，當(dāng)NR＝Thr1時(shí)，則調(diào)制方式為2FSK；否則，比較截面譜峰下降度D與第二預(yù)設(shè)閾值Thr2，若D＞Thr2，則調(diào)制方式2ASK；否則，比較截面最大值M與第三預(yù)設(shè)閾值Thr3，若M＜Thr3，則調(diào)制方式為BPSK或MSK，需做進(jìn)一步識(shí)別判斷，若M≥Thr3，則調(diào)制方式為QPSK或16QAM，需做進(jìn)一步識(shí)別判斷；比較截面最大值與次大值比R₁與第四預(yù)設(shè)閾值Thr4，若R₁＜Thr4，則調(diào)制方式為MSK，否則為BPSK；比較截面最大值與截面最大值比R₂與第五預(yù)設(shè)閾值，若R₂＜Thr5，則調(diào)制方式為16QAM，否則為QPSK。

更進(jìn)一步地，所述第一預(yù)設(shè)閾值Thr1為4，第二預(yù)設(shè)閾值Thr2為0.3，第三預(yù)設(shè)閾值Thr3為0.9，第四預(yù)設(shè)閾值Thr4為2.5，第五預(yù)設(shè)閾值Thr5為0.6。

本發(fā)明還公開了一種基于循環(huán)譜的無線通信信號(hào)調(diào)制識(shí)別與參數(shù)估計(jì)方法，在利用上述基于循環(huán)譜的無線通信信號(hào)調(diào)制識(shí)別方法獲取到信號(hào)的調(diào)制方式后，進(jìn)一步獲取信號(hào)的載波頻率估計(jì)值和符號(hào)率估計(jì)值載波頻率估計(jì)值的獲取具體如下：

若調(diào)制方式為2ASK或BPSK，則載波頻率估計(jì)值若調(diào)制方式為2FSK，則載波頻率估計(jì)值和若調(diào)制方式為MSK，則載波頻率估計(jì)值若調(diào)制方式為QPSK，則載波頻率估計(jì)值若調(diào)制方式為16QAM，則載波頻率估計(jì)值其中，f_s為采樣頻率，N′為循環(huán)譜計(jì)算中的FFT點(diǎn)數(shù)，I_max1為截面上譜峰最大值所對應(yīng)的位置，I_{max_square}為二次方循環(huán)譜截面上譜線最大值所對應(yīng)的坐標(biāo)，I_{max_forth}為四次方循環(huán)譜截面上譜線最大值所對應(yīng)的坐標(biāo)。

符號(hào)率估計(jì)值的獲取具體如下：

找到原始信號(hào)循環(huán)譜截面的最大值f_c的坐標(biāo)，對循環(huán)譜截面進(jìn)行搜索，記錄除中心點(diǎn)外譜峰最大值所在的坐標(biāo)I_max3；

若調(diào)制方式為2ASK或BPSK，則符號(hào)率估計(jì)值若調(diào)制方式為2FSK，則符號(hào)率估計(jì)值若調(diào)制方式為MSK，則符號(hào)率估計(jì)值若調(diào)制方式為QPSK，則符號(hào)率估計(jì)值若調(diào)制方式為16QAM，則符號(hào)率估計(jì)值

有益效果：本發(fā)明中基于循環(huán)譜的無線通信信號(hào)調(diào)制識(shí)別方法，通過快速傅里葉變換獲取信號(hào)的循環(huán)譜，并利用循環(huán)譜提取特征參數(shù)，進(jìn)一步利用不同調(diào)制方式特征參數(shù)的區(qū)別，通過預(yù)設(shè)閾值進(jìn)行比較，最終確定信號(hào)的調(diào)制方式，能夠在無任何先驗(yàn)信息的情況下對通信信號(hào)調(diào)制方式進(jìn)行準(zhǔn)確識(shí)別；在此基礎(chǔ)上進(jìn)一步提出基于循環(huán)譜的無線通信信號(hào)調(diào)制識(shí)別與參數(shù)估計(jì)方法，無需獲取更多的參數(shù)，便可實(shí)現(xiàn)對載波頻率和符號(hào)率進(jìn)行有效估計(jì)，使得調(diào)制識(shí)別與參數(shù)估計(jì)越來越發(fā)展為一個(gè)整體，有益于非協(xié)作通信的發(fā)展。

附圖說明

圖1為本發(fā)明中基于循環(huán)譜的無線通信信號(hào)調(diào)制識(shí)別及參數(shù)估計(jì)方法的整體框架圖；

圖2為本發(fā)明中基于循環(huán)譜的無線通信信號(hào)調(diào)制識(shí)別及參數(shù)估計(jì)方法的流程圖；

圖3為循環(huán)譜計(jì)算的流程示意圖；

圖4為本發(fā)明中二叉樹分類方式的調(diào)制識(shí)別流程圖；

圖5為本發(fā)明中特征參數(shù)與信噪比變化關(guān)系示意圖，圖5(a)為截面峰值個(gè)數(shù)NR與信噪比變化關(guān)系示意圖，圖5(b)為截面譜峰下降度D與信噪比變化關(guān)系示意圖，圖5(c)為截面最大值M與信噪比變化關(guān)系示意圖，圖5(d)為截面最大值與次大值比R₁與信噪比變化關(guān)系示意圖，圖5(e)為截面與截面最大值比值R₂與信噪比變化關(guān)系示意圖；

圖6為本發(fā)明中的載波頻率估計(jì)和符號(hào)率估計(jì)的均方誤差性能圖；圖6(a)為載波頻率估計(jì)與信噪比的變化關(guān)系圖，圖6(b)為載波頻率估計(jì)與數(shù)據(jù)點(diǎn)數(shù)的變化關(guān)系圖，圖6(c)為載波頻率估計(jì)與采樣率的變化關(guān)系圖；圖6(d)為符號(hào)率估計(jì)與信噪比的變化關(guān)系圖，圖6(e)為符號(hào)率估計(jì)與數(shù)據(jù)點(diǎn)數(shù)的變化關(guān)系圖，圖6(f)為符號(hào)率估計(jì)與采樣率的變化關(guān)系圖。

具體實(shí)施方式

下面結(jié)合附圖對本發(fā)明的技術(shù)方案作進(jìn)一步解釋。

如圖1所示為基于循環(huán)譜的調(diào)制識(shí)別與參數(shù)估計(jì)整體結(jié)構(gòu)圖，首先對模擬信號(hào)的帶寬進(jìn)行粗略估計(jì)，利用估計(jì)的帶寬值和中心頻率對信號(hào)進(jìn)行模擬下變頻到中頻頻段，經(jīng)過帶通濾波和數(shù)模轉(zhuǎn)換后對數(shù)字中頻信號(hào)的循環(huán)譜進(jìn)行計(jì)算，計(jì)算的結(jié)果用于調(diào)制識(shí)別與參數(shù)估計(jì)。

相應(yīng)地，圖2中基于循環(huán)譜的無線通信信號(hào)調(diào)制識(shí)別及參數(shù)估計(jì)方法包括以下步驟：

1)信號(hào)帶寬粗估計(jì)：接收端在接收未知信號(hào)后，對該未知信號(hào)的帶寬進(jìn)行粗略估計(jì)，得到帶寬和中心頻率；

2)帶通濾波：根據(jù)得到的信號(hào)中心頻率和帶寬估計(jì)值，對其進(jìn)行帶通濾波；

3)計(jì)算信號(hào)循環(huán)譜：利用快速傅里葉變換累加算法計(jì)算信號(hào)的循環(huán)譜密度函數(shù)，根據(jù)公式獲得信號(hào)循環(huán)譜密度值其中為信號(hào)x(t)在處理時(shí)間T長度內(nèi)的短時(shí)傅里葉變換，相應(yīng)地，信號(hào)x(t)的循環(huán)譜密度函數(shù)為

4)提取特征參數(shù)：根據(jù)計(jì)算得到的循環(huán)譜在其不同截面的譜峰值特點(diǎn)上提取特征參數(shù)，這些特征參數(shù)具體包括：截面峰值個(gè)數(shù)NR、截面譜峰下降度D、截面最大值M、截面最大值與次大值比R₁和截面與截面最大值比值R₂；

5)調(diào)制方式識(shí)別：根據(jù)得到的信號(hào)特征參數(shù)，采用基于二叉樹分類器的分類方法對未知信號(hào)進(jìn)行調(diào)制識(shí)別，二叉樹分類采用基于閾值比較判決的方法，根據(jù)前期仿真首先設(shè)定相關(guān)閾值的大小作為比較判決的依據(jù)，本發(fā)明方法適用于2ASK、2FSK、BPSK、QPSK、MSK、16QAM這幾種調(diào)制方式的識(shí)別；

6)載波頻率估計(jì)：根據(jù)調(diào)制方式識(shí)別結(jié)果，當(dāng)確定未知信號(hào)的調(diào)制方式時(shí)，若調(diào)制方式為2ASK/2FSK/BPSK/MSK，則載波頻率估計(jì)值由原始信號(hào)循環(huán)譜的截面計(jì)算出；若調(diào)制方式為QPSK時(shí)，載波頻率估計(jì)值由二次方循環(huán)譜的截面計(jì)算出；若調(diào)制方式為16QAM時(shí)，載波頻率估計(jì)值由四次方循環(huán)譜的截面計(jì)算出；

7)符號(hào)率估計(jì)：找到原始信號(hào)循環(huán)譜截面的最大值f_c的坐標(biāo)，符號(hào)率估計(jì)值由原始信號(hào)循環(huán)譜的截面計(jì)算出。

上述步驟3)中計(jì)算信號(hào)循環(huán)譜，如圖3所示，可以分為三個(gè)階段，分別為：

第一階段：將帶通濾波后的數(shù)據(jù)x(n)每隔L點(diǎn)進(jìn)行一次抽取，每次抽取的數(shù)據(jù)長度記為N′，總的抽取次數(shù)為P，對每次抽取的數(shù)據(jù)都要進(jìn)行加窗操作以減小頻譜泄露，再對加窗后的每塊數(shù)據(jù)進(jìn)行N'點(diǎn)的FFT變換，傅里葉變換的輸出為對應(yīng)頻率成分的復(fù)包絡(luò)；

第二階段：將復(fù)信號(hào)向基帶進(jìn)行頻移，頻移距離為k/N′，并求出頻移后復(fù)包絡(luò)的共軛相關(guān)，將結(jié)果逐行放入大小為N′²的乘法矩陣中，通過不斷重復(fù)得到一個(gè)P×N′²大小的矩陣；

第三階段：對得到的矩陣求P點(diǎn)FFT變換，對FFT后的結(jié)果進(jìn)行坐標(biāo)變換得到最終的信號(hào)循環(huán)譜。

步驟4)中提取特征參數(shù)，中各參數(shù)具體定義如下：

截面峰值個(gè)數(shù)NR：大于預(yù)設(shè)閾值的極大值個(gè)數(shù)為NR；截面譜峰下降度的最大值D：其中Δf為循環(huán)譜的頻率分辨率；截面最大值M：循環(huán)譜截面的譜峰最大值；截面最大值與次大值比R₁：其中secmax()表示取次大值的大??；截面與截面最大值比值R₂：

利用如下步驟獲取上述各特征參數(shù)：對循環(huán)譜截面進(jìn)行搜索，確定譜峰極大值、譜峰差值的最大值以及譜峰最大值所在的坐標(biāo)，通過將譜峰極大值與預(yù)設(shè)閾值進(jìn)行比較，當(dāng)極大值大于閾值時(shí)，特征參數(shù)NR的值加1，否則NR的值保持不變，通過譜峰差值的最大值確定特征參數(shù)D的值，并且記錄下譜峰最大值所在的坐標(biāo)f_c；

對循環(huán)譜截面進(jìn)行搜索，確定譜峰最大值與其對應(yīng)的坐標(biāo)I_max1，譜峰次大值與其對應(yīng)的坐標(biāo)I_max2，通過譜峰最大值確定特征參數(shù)M的值，通過最大值與次大值的比值確定特征參數(shù)R₁的值，通過截面的最大值與截面最大值的比值確定特征參數(shù)R₂的值。

上述步驟5)中采用二叉樹分類的方式完成調(diào)制方式的識(shí)別，利用特征參數(shù)NR從待識(shí)別的六種信號(hào)中將2FSK信號(hào)識(shí)別出來，即將信號(hào)分為兩組：{2FSK}和{2ASK、BPSK、QPSK、MSK、16QAM}；特征參數(shù)D將{2ASK、BPSK、QPSK、MSK、16QAM}分為{2ASK}和{BPSK、QPSK、MSK、16QAM}；特征參數(shù)M將{BPSK、QPSK、MSK、16QAM}分為{BPSK、MSK}和{QPSK、16QAM}；特征參數(shù)R₁將{BPSK、MSK}分為{BPSK}和{MSK}；特征參數(shù)R₂將{QPSK、16QAM}分為{QPSK}和{16QAM}。

如圖4所示，當(dāng)NR＝Thr1時(shí)將2FSK信號(hào)從{2ASK、2FSK、BPSK、QPSK、MSK、16QAM}里識(shí)別出來，當(dāng)D＞Thr2時(shí)將2ASK信號(hào)從{2ASK、BPSK、QPSK、MSK、16QAM}里識(shí)別出來，當(dāng)M＜Thr3時(shí)將{BPSK、MSK}從{BPSK、QPSK、MSK、16QAM}里識(shí)別出來，反之當(dāng)M≥Thr3時(shí)將{QPSK、16QAM}從{BPSK、QPSK、MSK、16QAM}里識(shí)別出來，當(dāng)R₁＜Thr4時(shí)將MSK信號(hào)從{BPSK、MSK}里識(shí)別出來，反之當(dāng)R₁≥Thr4將BPSK信號(hào)從{BPSK、MSK}里識(shí)別出來，當(dāng)R₂＜Thr5時(shí)將16QAM信號(hào)從{QPSK、16QAM}里識(shí)別出來，反之當(dāng)R₂≥Thr5時(shí)將QPSK信號(hào)從{QPSK、16QAM}里識(shí)別出來，其中Thr1、Thr2、Thr3、Thr4、Thr5分別表示設(shè)定的不同門限值。

圖5為載波頻率為500Hz，符號(hào)率為250bps，采樣頻率f_s＝4000Hz，考慮到實(shí)際短距離無線通信環(huán)境，設(shè)置信噪比變化范圍為-10dB到20dB條件下五個(gè)特征參數(shù)的變化圖。圖5(a)中可以看出隨著信噪比的增大，不同調(diào)制類型信號(hào)的NR趨于一個(gè)穩(wěn)定值，對于2FSK信號(hào)當(dāng)SNR＞5dB時(shí)其NR趨近于4，而對于其他5種調(diào)制信號(hào)當(dāng)SNR＞5dB時(shí)NR趨近于2；圖5(b)看出2ASK信號(hào)有著較大的譜峰下降度，在SNR＞-5dB時(shí)其譜峰下降度的最大值明顯大于其他信號(hào)；圖5(c)看出對于BPSK和MSK信號(hào)在SNR＞0dB且不斷增加時(shí)，其譜峰的歸一化幅度值趨近于1，而對于QPSK和16QAM信號(hào)，其譜峰歸一化幅度值一直保持在一個(gè)較小的水平；圖5(d)看出隨著信噪比的增加BPSK信號(hào)和MSK信號(hào)的特征值R₁大小呈現(xiàn)除了較大的差異；圖5(e)看出QPSK信號(hào)的特征參數(shù)R₂隨著信噪比的增加越來越趨近0.7，而16QAM信號(hào)的特征參數(shù)R₂隨著信噪比的增加不斷趨近于0.5。因此，上述流程中的門限值可分別設(shè)置為Thr1＝4，Thr2＝0.3，Thr3＝0.9，Thr4＝2.5，Thr5＝0.6，從圖5各附圖中可以看出在SNR＝10dB時(shí)識(shí)別率能達(dá)到95％以上。

上述步驟7)中符號(hào)率估計(jì)具體包括：

找到原始信號(hào)循環(huán)譜截面的最大值f_c的坐標(biāo)，對循環(huán)譜截面進(jìn)行搜索，記錄除中心點(diǎn)外譜峰最大值所在的坐標(biāo)I_max3；

根據(jù)步驟5)中調(diào)制識(shí)別的結(jié)果，結(jié)合步驟4)得到的坐標(biāo)值，對于2ASK/BPSK調(diào)制方式，載波頻率的估計(jì)值符號(hào)率估計(jì)值對于2FSK調(diào)制方式，載波頻率估計(jì)值和符號(hào)率估計(jì)值對于MSK調(diào)制方式，載波頻率估計(jì)值符號(hào)率估計(jì)值對于QPSK調(diào)制方式，載波頻率估計(jì)值符號(hào)率估計(jì)值對于16QAM調(diào)制方式，載波頻率估計(jì)值符號(hào)率估計(jì)值其中，f_s為采樣頻率，N′為循環(huán)譜計(jì)算中的FFT點(diǎn)數(shù)，I_max1為截面上譜峰最大值所對應(yīng)的位置，I_max3表示本步驟中記錄的坐標(biāo)值，I_{max_square}為二次方循環(huán)譜截面上譜線最大值所對應(yīng)的坐標(biāo)，I_{max_forth}為四次方循環(huán)譜截面上譜線最大值所對應(yīng)的坐標(biāo)。

如圖6所示為從信噪比、數(shù)據(jù)點(diǎn)數(shù)、采樣率三個(gè)方面分析的載波頻率估計(jì)和符號(hào)率估計(jì)的均方誤差變化圖，其中均方誤差定義為和仿真參數(shù)與圖5一致，可以看到總的來看，在SNR＞10dB，數(shù)據(jù)點(diǎn)數(shù)大于1024，采樣率的偏差Δf_s/f_b＜0.5時(shí)，參數(shù)估計(jì)的均方誤差均能夠達(dá)到10^-4，性能較好。

以上所述，僅為本發(fā)明較佳的具體實(shí)施方式，但本發(fā)明的保護(hù)范圍不僅限于此。本發(fā)明的保護(hù)范圍以權(quán)利要求的保護(hù)范圍為準(zhǔn)。