本發(fā)明涉及雷達(dá)技術(shù)領(lǐng)域，特別涉及一種頻率捷變信號(hào)的時(shí)域調(diào)制域參數(shù)聯(lián)合測(cè)量方法，可用于雷達(dá)信號(hào)處理。

背景技術(shù)：

頻率捷變信號(hào),是指雷達(dá)信號(hào)相鄰脈沖的載頻在一定頻帶內(nèi)隨機(jī)快速改變的一種信號(hào)類(lèi)型。測(cè)量頻率捷變信號(hào)的參數(shù)是雷達(dá)信號(hào)處理工程領(lǐng)域的重點(diǎn)之一。傳統(tǒng)的測(cè)量方法僅針對(duì)時(shí)域或者調(diào)制域，未把時(shí)域和調(diào)制域結(jié)合在一起進(jìn)行參數(shù)測(cè)量，而且測(cè)量頻率捷變信號(hào)調(diào)制域參數(shù)需要進(jìn)行時(shí)頻分析，常用的時(shí)頻分析方法，如短時(shí)傅里葉變換STFT和Wigner-Vile分布，它們?cè)趯?duì)頻率捷變信號(hào)時(shí)頻分析時(shí)存在著諸多缺陷。STFT存在著窗長(zhǎng)和窗函數(shù)選擇困難的問(wèn)題，由于使用了傅里葉變換，運(yùn)算量較大，而且不能解決時(shí)間分辨率和頻率分辨率的矛盾。Wigner-Vile分布雖然不使用窗函數(shù)，不存在STFT的問(wèn)題，但是對(duì)多分量信號(hào)進(jìn)行時(shí)頻分布時(shí)存在著“交叉項(xiàng)”干擾，時(shí)頻分布將變得模糊，對(duì)頻率估計(jì)不準(zhǔn)，而且它采用了積分運(yùn)算，運(yùn)算量較大。

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

本發(fā)明的目的在于針對(duì)上述已有技術(shù)的不足，提出一種頻率捷變信號(hào)的時(shí)域調(diào)制域參數(shù)聯(lián)合測(cè)量方法，以減小運(yùn)算量和誤差，獲得準(zhǔn)確的捷變頻率調(diào)制域參數(shù)，提高測(cè)量精度。

實(shí)現(xiàn)本發(fā)明目的的技術(shù)思路是聯(lián)合時(shí)域、調(diào)制域測(cè)量，時(shí)域分析時(shí)采用統(tǒng)計(jì)脈沖信號(hào)包絡(luò)數(shù)值的方法估計(jì)判決門(mén)限，實(shí)現(xiàn)脈沖個(gè)數(shù)、脈沖寬度等時(shí)域參數(shù)的測(cè)量，在調(diào)制域分析時(shí)采用相位差分法來(lái)得到信號(hào)的瞬時(shí)頻率曲線，獲得準(zhǔn)確的捷變頻率等調(diào)制域參數(shù)，其技術(shù)步驟包括如下：

(1)測(cè)量頻率捷變信號(hào)的離散數(shù)值；

1a)用信號(hào)分析儀讀取得到待測(cè)量的頻率捷變信號(hào)實(shí)部I(t)和虛部Q(t)的數(shù)值大小；

1b)設(shè)定采樣率f_s，抽樣選取1秒時(shí)間長(zhǎng)度的待測(cè)量的頻率捷變信號(hào)實(shí)部和虛部，得到此頻率捷變信號(hào)實(shí)部離散取值I(n)和虛部離散取值Q(n)；

(2)通過(guò)離散取值得到頻率捷變信號(hào)的時(shí)域參數(shù)：

(2a)根據(jù)頻率捷變信號(hào)的實(shí)、虛部離散取值I(n)和Q(n)，計(jì)算頻率捷變信號(hào)包絡(luò)：

(2b)根據(jù)頻率捷變信號(hào)包絡(luò)y(n)估計(jì)出頻率捷變信號(hào)的近似頂值A(chǔ)'_top和近似底值A(chǔ)'_ba，得到近似幅度測(cè)量值為：A'＝A'_top-A'_ba；

(2c)用10％的近似幅度測(cè)量值10％A'作為上升沿和下降沿的近似檢測(cè)門(mén)限，檢測(cè)得到頻率捷變信號(hào)的脈沖個(gè)數(shù)n，用每個(gè)脈沖上升沿的近似檢測(cè)門(mén)限所在的點(diǎn)數(shù)和下降沿的近似檢測(cè)門(mén)限所在的點(diǎn)數(shù)分別除以采樣率，得到相對(duì)應(yīng)的每個(gè)脈沖的上升沿的近似判決時(shí)刻t′_升和下降沿的近似判決時(shí)刻t′_降，計(jì)算脈沖信號(hào)各個(gè)脈沖的脈沖寬度為：τ＝t′_升-t′_降；

(2d)在每個(gè)脈沖的脈沖寬度內(nèi)，運(yùn)用密度分布平均法，計(jì)算每個(gè)脈沖精確的頂值A(chǔ)_top和底值A(chǔ)_ba，并求得每個(gè)脈沖精確幅度值A(chǔ)為：A＝A_top-A_ba；

(2e)根據(jù)每個(gè)脈沖的最大值V_max和最小值V_min，計(jì)算每個(gè)脈沖的過(guò)沖S_over和下沖S_under：

S_over＝V_max-A_top，S_under＝A_ba-V_min，

(2f)根據(jù)頻率捷變信號(hào)的每個(gè)脈沖精確的頂值A(chǔ)_top、底值A(chǔ)_ba和精確幅度值A(chǔ)，計(jì)算頻率捷變信號(hào)的脈沖幅度上參考線M_上和脈沖幅度下參考線M_下：

M_上＝A_ba+90％×A，M_下＝A_ba+10％×A；

(2g)根據(jù)頻率捷變信號(hào)的脈沖幅度上參考線M_上、脈沖幅度下參考線M_下和脈沖寬度τ，在脈沖寬度內(nèi)分別查找上升沿中脈沖幅度上參考線對(duì)應(yīng)的時(shí)刻值t_r上、脈沖幅度下參考線對(duì)應(yīng)的時(shí)刻值t_r下和下降沿中脈沖幅度上參考線對(duì)應(yīng)的時(shí)刻值t_f上、脈沖幅度下參考線對(duì)應(yīng)的時(shí)刻值t_f下，以及相鄰下一次上升沿中脈沖幅度上參考線所對(duì)應(yīng)的時(shí)刻t'_r下，計(jì)算頻率捷變信號(hào)脈沖的上升時(shí)間t_r、下降時(shí)間t_f，脈沖周期T，關(guān)閉時(shí)間t_off和占空比d_t：

t_r＝t_r上-t_r下，t_f＝t_f下-t_f上，T＝t'_r下-t_r下，t_off＝T-τ，d_t＝τ/T；

(3)通過(guò)離散取值得到頻率捷變信號(hào)的調(diào)制域參數(shù)：

(3a)根據(jù)頻率捷變信號(hào)的實(shí)、虛部離散取值I(n)和Q(n)，計(jì)算頻率捷變信號(hào)的瞬時(shí)相位：

(3b)根據(jù)頻率捷變信號(hào)的瞬時(shí)相位采用相位差分法對(duì)頻率捷變信號(hào)進(jìn)行時(shí)頻分析，得到頻率捷變信號(hào)的瞬時(shí)頻率：其中f_s為信號(hào)的采樣率；

(3c)根據(jù)頻率捷變信號(hào)的瞬時(shí)頻率f_c(n)，計(jì)算頻率捷變信號(hào)瞬時(shí)頻率的平均值F_avg；

(3d)根據(jù)頻率捷變信號(hào)瞬時(shí)頻率的平均值F_avg和脈沖個(gè)數(shù)n，統(tǒng)計(jì)頻率捷變信號(hào)瞬時(shí)頻率在平均值附近的個(gè)數(shù)N，若滿足90％n≤N≤n，則判斷該頻率捷變信號(hào)的捷變方式為脈內(nèi)捷變，否則判斷為脈間捷變；

(3e)根據(jù)頻率捷變信號(hào)實(shí)、虛部的離散取值I(n)和Q(n)，分別對(duì)每個(gè)脈沖內(nèi)的點(diǎn)的實(shí)部、虛部做快速傅里葉變換得到對(duì)應(yīng)的實(shí)部I(ω)和虛部Q(ω)，并計(jì)算其瞬時(shí)功率：

(3f)根據(jù)頻率捷變信號(hào)的瞬時(shí)功率spec(ω)，計(jì)算脈內(nèi)捷變方式或脈間捷變方式的捷變個(gè)數(shù)C和捷變頻率freq(i)，i＝1,2,…,C。

本發(fā)明與現(xiàn)有技術(shù)相比具有如下優(yōu)點(diǎn)：

1、本發(fā)明在估計(jì)脈沖的近似幅度的過(guò)程中，采用密度分布平均法來(lái)估計(jì)頻率捷變信號(hào)的近似頂值A(chǔ)'_top和近似底值A(chǔ)'_ba，與現(xiàn)有的密度分布眾數(shù)法相比精度更高，當(dāng)頻率捷變信號(hào)的采樣點(diǎn)數(shù)越多時(shí)，脈沖包絡(luò)最大值y_max和脈沖包絡(luò)最小值y_min之間均分的區(qū)間個(gè)數(shù)越多，所求得的近似底值A(chǔ)'_top和近似頂值A(chǔ)'_ba越精確，可以滿足更好的精度要求。

2、本發(fā)明進(jìn)行時(shí)頻分析時(shí)采用相位差分法，相對(duì)于STFT法，運(yùn)算量大大減小，而且沒(méi)有窗函數(shù)，也不受測(cè)不準(zhǔn)準(zhǔn)則的限制；相對(duì)于Wigner-Vile分布法，沒(méi)有“交叉項(xiàng)”的干擾問(wèn)題，相對(duì)于其他方法，具有測(cè)量參數(shù)運(yùn)算量小，誤差小等優(yōu)點(diǎn)，適合于工程實(shí)踐。

下面結(jié)合附圖，對(duì)本發(fā)明做進(jìn)一步詳細(xì)描述。

附圖說(shuō)明

圖1是本發(fā)明的實(shí)現(xiàn)流程圖；

圖2是用現(xiàn)有STFT方法對(duì)頻率捷變信號(hào)進(jìn)行時(shí)頻分析的時(shí)頻曲線圖；

圖3是用現(xiàn)有Wigner-Vile方法對(duì)頻率捷變信號(hào)進(jìn)行時(shí)頻分析的時(shí)頻曲線圖；

圖4是用本發(fā)明方法進(jìn)行時(shí)頻分析的瞬時(shí)相位曲線及時(shí)頻曲線圖。

具體實(shí)施方式

在雷達(dá)接收中，雷達(dá)回波的時(shí)域、調(diào)制域參數(shù)搭載了很多有用信息，而雷達(dá)回波的參數(shù)是未知的，并不能通過(guò)雷達(dá)接收天線直接獲取。對(duì)于不同類(lèi)型的雷達(dá)信號(hào)，其參數(shù)測(cè)量方法也不盡相同，頻率捷變信號(hào)是雷達(dá)常用的一種發(fā)射信號(hào)，測(cè)量頻率捷變信號(hào)參數(shù)是雷達(dá)信號(hào)處理工程領(lǐng)域的重點(diǎn)之一。頻率捷變信號(hào)由于脈內(nèi)調(diào)制復(fù)雜的特性使參數(shù)測(cè)量尤為困難，運(yùn)算量巨大，經(jīng)常出現(xiàn)測(cè)量不準(zhǔn)確的問(wèn)題。為了能夠準(zhǔn)確快速得測(cè)量頻率捷變信號(hào)的時(shí)域、調(diào)制域參數(shù)，必須采取一種運(yùn)算量小并且精度更高的測(cè)量方法。本發(fā)明就是通過(guò)信號(hào)分析儀讀取頻率捷變信號(hào)的數(shù)值大小，進(jìn)而采用聯(lián)合時(shí)域、調(diào)制域的思想，用密度分布平均法和相位差分法更加精確詳細(xì)地測(cè)量頻率捷變信號(hào)的具體參數(shù)。

參照?qǐng)D1，本發(fā)明的實(shí)施步驟如下：

步驟1：測(cè)量頻率捷變信號(hào)的離散數(shù)值。

1a)雷達(dá)通過(guò)接收天線接收頻率捷變信號(hào)形式的雷達(dá)回波，并將接收到的頻率捷變信號(hào)的實(shí)部I(t)和虛部Q(t)以數(shù)據(jù)形式保存到.dat文件中；

1b)將.dat文件上傳到信號(hào)分析儀，讀取得到頻率捷變信號(hào)的實(shí)部I(t)和虛部Q(t)的數(shù)值大小；

1c)設(shè)定采樣率f_s，抽取1秒時(shí)間長(zhǎng)度的待測(cè)量的頻率捷變信號(hào)，對(duì)其進(jìn)行離散采樣，得到此1秒長(zhǎng)度的頻率捷變信號(hào)的實(shí)部離散取值I(n)和虛部離散取值Q(n)。

步驟2：計(jì)算頻率捷變信號(hào)的時(shí)域參數(shù)。

2a)根據(jù)頻率捷變信號(hào)的實(shí)、虛部離散取值I(n)和Q(n)，計(jì)算頻率捷變信號(hào)包絡(luò)：

2b)估計(jì)頻率捷變信號(hào)的近似頂值A(chǔ)'_top和近似底值A(chǔ)'_ba，以及近似幅度A'；

2b1)查找得到信號(hào)包絡(luò)y(n)的最大值y_max和最小值y_min，其中n＝1,2,…,N，N為頻率捷變信號(hào)抽樣所得的總點(diǎn)數(shù)；

2b2)將信號(hào)包絡(luò)y(n)從小到大分為100個(gè)區(qū)間，其中第k個(gè)區(qū)間E(k)取值范圍的最小值為E_min(k)，最大值為E_max(k)：

2b3)統(tǒng)計(jì)信號(hào)包絡(luò)y(n)落在每個(gè)區(qū)間E(k)內(nèi)的個(gè)數(shù)，分別記為計(jì)數(shù)值C(k),k＝1,2…100；

2b4)在計(jì)數(shù)值C(k)中查找前50個(gè)計(jì)數(shù)值的最大值所在的位置，記為k₁，查找后50個(gè)計(jì)數(shù)值的最大值所在的位置，記為k₂；

2b5)計(jì)算區(qū)間E(k₁)的平均值，即為頻率捷變信號(hào)的近似頂值A(chǔ)'_top，計(jì)算區(qū)間E(k₂)的平均值，即為頻率捷變信號(hào)的近似底值A(chǔ)'_ba；

2b6)根據(jù)頻率捷變信號(hào)的近似底值A(chǔ)'_ba和近似頂值A(chǔ)'_top，得到頻率捷變信號(hào)的近似幅度A'＝A'_top-A'_ba；

2c)檢測(cè)頻率捷變信號(hào)的脈沖個(gè)數(shù)m，并記錄每個(gè)脈沖的近似上升沿判決時(shí)刻t'_升、近似下降沿的判決時(shí)刻t'_降；

2c1)將頻率捷變信號(hào)的脈沖個(gè)數(shù)記為m，并初始化m為0，從頻率捷變信號(hào)的起始位置開(kāi)始查找第一個(gè)近似上升沿，將查找結(jié)果記為標(biāo)志位flag，若查找成功，則將標(biāo)志位flag置1，并記錄本次近似上升沿判決時(shí)刻t'_升，否則，將標(biāo)志位flag置0；

2c2)查找到近似上升沿后，繼續(xù)往后查找近似下降沿，若查找成功，則將標(biāo)志位flag置1，并記錄本次近似下降沿判決時(shí)刻t'_降，否則，將標(biāo)志位flag置0；

2c3)查找到近似下降沿后，繼續(xù)往后查找近似上升沿，若查找成功，則將標(biāo)志位flag置1，記錄本次近似上升沿判決時(shí)刻t'_升，并令頻率捷變信號(hào)的脈沖個(gè)數(shù)m加1，否則，將標(biāo)志位flag置0；

2c4)檢測(cè)標(biāo)志位flag的當(dāng)前值，若標(biāo)志位flag＝1，則返回步驟2c2)，若標(biāo)志位flag＝0，則結(jié)束查找；

2d)計(jì)算頻率捷變信號(hào)各個(gè)脈沖的寬度：τ＝t'_降-t'_升；

2e)在每個(gè)脈沖的脈沖寬度內(nèi)，運(yùn)用密度分布平均法，計(jì)算得到每個(gè)脈沖精確的頂值A(chǔ)_top和底值A(chǔ)_ba，并求得每個(gè)脈沖精確幅度值A(chǔ)為：A＝A_top-A_ba；

2f)檢測(cè)頻率捷變信號(hào)每個(gè)脈沖的最大值V_max和最小值V_min，計(jì)算每個(gè)脈沖的過(guò)沖S_over＝V_max-A_top和下沖S_under＝A_ba-V_min；

2g)計(jì)算出脈沖的上升時(shí)間t_r、下降時(shí)間t_f，脈沖周期T，關(guān)閉時(shí)間t_off及占空比d_t：

2g1)設(shè)定頻率捷變信號(hào)的脈沖幅度上參考線M_上和脈沖幅度下參考線M_下：

M_上＝A_ba+90％×A，

M_下＝A_ba+10％×A；

2g2)在脈沖寬度τ內(nèi)查找脈沖幅度上參考線M_上和脈沖幅度下參考線M_下對(duì)應(yīng)的時(shí)刻值，該對(duì)應(yīng)的時(shí)刻值包括：上升沿中脈沖幅度上參考線所對(duì)應(yīng)的時(shí)刻t_r上、上升沿中脈沖幅度下參考線所對(duì)應(yīng)的時(shí)刻t_r下、下降沿中脈沖幅度上參考線所對(duì)應(yīng)的時(shí)刻t_f上、下降沿中脈沖幅度下參考線所對(duì)應(yīng)的時(shí)刻t_f下，及相鄰下一次上升沿中脈沖幅度上參考線所對(duì)應(yīng)的時(shí)刻t'_r下；

2g3)根據(jù)2g2)查找的出的時(shí)刻值，計(jì)算出脈沖的如下參數(shù)：

上升時(shí)間t_r：t_r＝t_r上-t_r下，

下降時(shí)間t_f：t_f＝t_f下-t_f上，

脈沖周期T：T＝t'_r下-t_r下，

關(guān)閉時(shí)間t_off：t_off＝T-τ，

占空比d_t：d_t＝τ/T。

步驟3：計(jì)算頻率捷變信號(hào)的調(diào)制域參數(shù)。

3a)根據(jù)頻率捷變信號(hào)的實(shí)、虛部離散取值I(n)和Q(n)，計(jì)算頻率捷變信號(hào)每個(gè)點(diǎn)的瞬時(shí)相位：n＝1,2,…,N，其中N為頻率捷變信號(hào)抽樣所得的總點(diǎn)數(shù)；

3b)根據(jù)瞬時(shí)相位采用相位差分法對(duì)頻率捷變信號(hào)進(jìn)行時(shí)頻分析，提取頻率捷變信號(hào)的瞬時(shí)頻率：n＝1,2,…,N，其中f_s為信號(hào)的采樣率，N為頻率捷變信號(hào)抽樣所得的總點(diǎn)數(shù)；

3c)根據(jù)頻率捷變信號(hào)各點(diǎn)的瞬時(shí)頻率f_c(n)，求得頻率捷變信號(hào)瞬時(shí)頻率的平均值；

3d)統(tǒng)計(jì)頻率捷變信號(hào)的瞬時(shí)頻率在平均值附近的個(gè)數(shù)M，若滿足90％m≤M≤m，則判斷該頻率捷變信號(hào)的捷變方式為脈內(nèi)捷變頻，否則判斷為脈間捷變頻，m為頻率捷變信號(hào)的脈沖個(gè)數(shù)；

3e)根據(jù)頻率捷變信號(hào)實(shí)、虛部的離散取值I(n)和Q(n)，分別對(duì)每個(gè)脈沖內(nèi)的點(diǎn)的實(shí)部、虛部做快速傅里葉變換得到對(duì)應(yīng)的實(shí)部I(ω)和虛部Q(ω)，并計(jì)算其瞬時(shí)功率：

3f)根據(jù)頻率捷變信號(hào)的瞬時(shí)功率spec(ω)，計(jì)算脈內(nèi)捷變方式的捷變個(gè)數(shù)C和捷變頻率freq(i)，i＝1,2,…,C；

3f1)統(tǒng)計(jì)得到瞬時(shí)功率spec(ω)的值域spec(j)，j＝1,2,…,N，其中N為頻率捷變信號(hào)抽樣所得的總點(diǎn)數(shù)；

3f2)計(jì)算得到瞬時(shí)均值avg：avg＝(spec(1)+spec(2)+…+spec(j)+…+spec(N))/N；

3f3)依次對(duì)每個(gè)值域spec(j)的值進(jìn)行判斷，若滿足spec(j)>avg并且spec(j-2)<spec(j-1)，spec(j-1)<spec(j)，spec(j)>spec(j+1)，spec(j+1)>spec(j+2)，則將該點(diǎn)記為一個(gè)頻譜峰值pSpec(ω)；

3f4)統(tǒng)計(jì)查找出的所有的頻譜峰值pSpec(ω)，依次比較每個(gè)頻譜峰值pSpec(ω)的大小，確定出最大頻譜峰值pSpec_max(ω)以及最大頻譜峰值的個(gè)數(shù)，最大譜峰值個(gè)數(shù)即為捷變個(gè)數(shù)C；

3f5)將每個(gè)脈沖的脈沖寬度分為捷變個(gè)數(shù)段，分別對(duì)每一段計(jì)算其擬合頻率即為捷變頻率freq(i)，i＝1,2,…,C；

3g)根據(jù)頻率捷變信號(hào)的瞬時(shí)功率spec(ω)，計(jì)算脈間捷變方式的捷變個(gè)數(shù)C和捷變頻率freq(i)，i＝1,2,…,C；

3g1)分別對(duì)頻率捷變信號(hào)的每個(gè)脈沖計(jì)算其擬合頻率freq(k),k＝1,2,…,m，其中m為頻率捷變信號(hào)的脈沖個(gè)數(shù)；

3g2)查找得出擬合頻率freq(k)中的最大值freq_max及其所在脈沖的位置；

3g3)根據(jù)頻率捷變信號(hào)擬合頻率的最大值freq_max所在的脈沖位置，計(jì)算相鄰兩個(gè)最大值所在脈沖位置之差，該差值即為捷變個(gè)數(shù)C；

3g4)取前捷變個(gè)數(shù)C的擬合頻率，即為捷變頻率freq(i),i＝1,2,…,C。

本發(fā)明對(duì)頻率捷變信號(hào)的時(shí)域、調(diào)制域參數(shù)測(cè)量的效果可通過(guò)以下仿真實(shí)驗(yàn)進(jìn)一步驗(yàn)證。

實(shí)驗(yàn)1，用本發(fā)明方法對(duì)頻率捷變信號(hào)的時(shí)域參數(shù)測(cè)量進(jìn)行Matlab仿真，參數(shù)設(shè)置為：捷變方式為脈間捷變，脈沖寬度為10us，脈沖重復(fù)周期為60us，完整脈沖個(gè)數(shù)為6，采樣率為20MHz，帶寬20MHz。測(cè)得的時(shí)域參數(shù)如表1所示：

表1頻率捷變信號(hào)時(shí)域參數(shù)表

從表1可以看出，測(cè)量的脈沖個(gè)數(shù)為6，脈沖寬度為10us，重復(fù)周期為60us，上升時(shí)間為0.05us，下降時(shí)間為0.05us，該測(cè)量參數(shù)與實(shí)際設(shè)定的參數(shù)一致。

實(shí)驗(yàn)2，用本發(fā)明方法對(duì)頻率捷變信號(hào)的調(diào)制域參數(shù)測(cè)量進(jìn)行Matlab仿真，參數(shù)設(shè)置為：帶寬20Mhz，脈內(nèi)最大頻率、最小頻率為10.2e+008和9.8e+008。測(cè)量所得各脈內(nèi)調(diào)制參數(shù)如表2所示:

表2頻率捷變信號(hào)脈內(nèi)調(diào)制參數(shù)表

從表2可以看出，捷變方式為脈間捷變，捷變個(gè)數(shù)為3個(gè)，捷變頻率分別為0MHz，1.5MHz，3MHz，測(cè)量結(jié)果準(zhǔn)確，效果很好。

實(shí)驗(yàn)3，分別用現(xiàn)有STFT方法和Wigner-Vile方法對(duì)頻率捷變信號(hào)進(jìn)行時(shí)頻分析得到相應(yīng)時(shí)頻曲線，測(cè)量結(jié)果如圖2、圖3所示。其中圖2為STFT方法測(cè)得的頻率捷變信號(hào)的時(shí)頻分布圖，圖3為Wigner-Vile方法測(cè)得的頻率捷變信號(hào)的時(shí)頻分布圖。

從圖2中可以看出，用現(xiàn)有STFT方法得到的時(shí)頻分布圖的分辨率較差。

從圖3中可以看出，采用現(xiàn)有Wigner-Vile分布方法，頻率捷變信號(hào)的時(shí)頻分布出現(xiàn)了“虛假”，有用的時(shí)變譜圖變得模糊，因此STFT和Wigner-Vile分布都不適合分析頻率捷變信號(hào)。

實(shí)驗(yàn)4，用本發(fā)明方法，對(duì)頻率捷變信號(hào)進(jìn)行時(shí)頻分析得到瞬時(shí)相位曲線和時(shí)頻曲線，測(cè)量結(jié)果如圖4所示，其中圖4(a)為頻率捷變信號(hào)的瞬時(shí)相位曲線，圖4(b)為頻率捷變信號(hào)的時(shí)頻曲線，圖4(c)為頻率捷變信號(hào)的瞬時(shí)相位曲線在頻率為1.5MHz時(shí)的放大圖，圖4(d)為頻率捷變信號(hào)的瞬時(shí)相位曲線在頻率為3MHz時(shí)的放大圖。

從圖4仿真結(jié)果可看出，而本發(fā)明的頻率捷變信號(hào)參數(shù)的測(cè)量方法，具有算法簡(jiǎn)單，運(yùn)算量小等特點(diǎn)，不存在“交叉項(xiàng)”干擾，特別適合于工程實(shí)踐上大數(shù)據(jù)量的測(cè)量。

綜上所述，本發(fā)明的頻率捷變信號(hào)參數(shù)的測(cè)量方法能夠更為精確地對(duì)頻率捷變信號(hào)進(jìn)行時(shí)頻分析，更加快速的得到各個(gè)時(shí)域參數(shù)和調(diào)制域參數(shù)。