用改進(jìn)的Chebyshev函數(shù)提高干涉圖切趾效果*

江峰，武高衛(wèi)，盛文，豐驍，馮至

(空軍預(yù)警學(xué)院， 湖北 武漢 430019)

用改進(jìn)的Chebyshev函數(shù)提高干涉圖切趾效果*

江峰，武高衛(wèi)，盛文，豐驍，馮至

(空軍預(yù)警學(xué)院， 湖北 武漢 430019)

針對Chebyshev切趾函數(shù)在低旁瓣衰減時功率泄漏率高以及旁瓣衰減的等波動性問題，提出了一種改進(jìn)的Chebyshev切趾函數(shù)。首先通過仿真分析指出Chebyshev切趾函數(shù)在光譜復(fù)原中存在的缺陷；然后選取Blackman-Harris函數(shù)對其進(jìn)行改進(jìn)，并與改進(jìn)前進(jìn)行了對比分析；最后，用該函數(shù)對單色光干涉圖進(jìn)行切趾處理并仿真驗(yàn)證，仿真結(jié)果表明，改進(jìn)后Chebyshev切趾函數(shù)能夠較好地克服功率泄露率高及旁瓣衰減等波動性的問題，其應(yīng)用能力也更加靈活。

Chebyshev；切趾函數(shù)；Blackman-Harris；干涉圖；光譜復(fù)原；仿真分析

0 引言

在傅里葉變換光譜儀光譜復(fù)原的過程中，存在傅里葉變換需要無限的積分區(qū)間而實(shí)際光譜儀只能提供有限光程差的矛盾，這個矛盾會造成復(fù)原光譜存在旁瓣效應(yīng)進(jìn)而影響鄰近的、微弱信號的測定，此時就需要將光譜儀的干涉圖進(jìn)行切趾處理，以消除或者減弱虛假的旁瓣，增強(qiáng)有用信號的被檢測率[1-2]。切趾的普遍做法是在信號處理中用切趾函數(shù)乘以干涉圖，迄今學(xué)界已經(jīng)提出了眾多的切趾函數(shù)，包括三角函數(shù)、Hamming函數(shù)、Blackman函數(shù)、KaiserBessel函數(shù)以及各種函數(shù)族[3]。切趾函數(shù)的選取一般要考慮以下原則[4]：

(1) 切趾函數(shù)的形式盡量簡單，切趾計算過程方便，計算量小；

(2) 切趾處理后復(fù)原光譜主瓣盡量窄，以保證系統(tǒng)分辨率；

(3) 切趾處理后復(fù)原光譜旁瓣盡量低，旁瓣衰減率盡量高。

實(shí)際應(yīng)用中，用形式簡單的切趾函數(shù)處理干涉圖容易丟失有用信號，在抑制旁瓣的同時，系統(tǒng)分辨率也會降低，因此，要綜合考慮切趾函數(shù)性質(zhì)以及系統(tǒng)要求，選取合適的切趾函數(shù)。在給定旁瓣高度下，Chebyshev切趾函數(shù)的主瓣寬度最小，滿足切趾函數(shù)的最大震幅比準(zhǔn)則[5]，同時其旁瓣衰減量可調(diào)節(jié)范圍寬，是一種性能優(yōu)異的切趾函數(shù)；但是當(dāng)干涉圖信號功率較弱，需要小的旁瓣衰減情形時，其主瓣功率泄露嚴(yán)重，在頻域上表現(xiàn)為相位誤差增大，并且由于旁瓣衰減的等波動性，旁瓣衰減速率為零，不利于臨近強(qiáng)干擾信號的抑制。文獻(xiàn)[6]對基于傳統(tǒng)Chebyshev函數(shù)設(shè)計的濾波器相位特性較差的缺點(diǎn)進(jìn)行了改進(jìn)，犧牲幅度特性提升其相位特性；文獻(xiàn)[7]針對其旁瓣衰減的等波動性，提出離散傅立葉變換后對旁瓣進(jìn)行相位旋轉(zhuǎn)，利用疊加校正的方法使旁瓣干擾降到最低，然而，這種方法在增加計算量、引入誤差的同時，計算精度得不到保證。因此，本文提出一種改進(jìn)的Chebyshev切趾函數(shù)，在保留Chebyshev切趾函數(shù)原有優(yōu)點(diǎn)的基礎(chǔ)上，通過引入Blackman-harris函數(shù)改善其旁瓣衰減的等波動性、低旁瓣衰減時主瓣功率泄露嚴(yán)重的問題，同時，新切趾函數(shù)第一旁瓣衰減的可調(diào)節(jié)性也更加靈活。

1 Chebyshev切趾函數(shù)及其分析

Chebyshev切趾函數(shù)是由Chebyshev多項(xiàng)式在單位圓上作等間隔抽樣，再作DFT反變換得到的[8]。Chebyshev多項(xiàng)式定義為

(1)

式(1)可以通過迭代產(chǎn)生：Tk+1(x)=2xTk(x)-Tk-1(x),T0(x)=1,T1(x)=x，其中x∈Z,k∈R且k≥1，令函數(shù)長度N=2M+1，則Chebyshev切趾函數(shù)在時域的表達(dá)式為

wC(n,γ)=

-M≤n≤M,

(2)

(3)

式中：γ為用分?jǐn)?shù)表示的旁瓣與主瓣幅度的比值[9-10]，則第一旁瓣相對于主瓣衰減值m=20lgγ(dB)，在實(shí)際應(yīng)用中，僅需調(diào)整γ的值就可以獲得不同的主瓣寬度、旁瓣衰減等指標(biāo)，Chebyshev切趾函數(shù)的性質(zhì)如表1所示，與其他典型切趾函數(shù)的幅頻響應(yīng)對比如圖1所示。

表1 Chebyshev切趾函數(shù)在不同旁瓣衰減時的性質(zhì)(N=128)

圖1 Chebyshev切趾函數(shù)與兩種典型切趾函數(shù)在相同旁瓣衰減下的對比(N=128)Fig.1 Comparison of Chebyshev apodization function with two typical apodization functions at the same side-lobe attenuation(N=128)

圖1中，a)是Chebyshev切趾函數(shù)與Hamming函數(shù)(形式簡單、切趾性能不可調(diào))的對比，b)是Chebyshev切趾函數(shù)與Gausswin函數(shù)(形式復(fù)雜、切趾性能可調(diào))的對比。通過圖1可以發(fā)現(xiàn)，在相同第一旁瓣衰減情況下，Chebyshev切趾函數(shù)具有更窄的主瓣寬度，這使得干涉圖經(jīng)過Chebyshev切趾后光譜分辨率更高；同時，Chebyshev切趾函數(shù)第一旁瓣衰減可調(diào)節(jié)范圍寬，具有更強(qiáng)的工程實(shí)用價值。然而通過表1可以發(fā)現(xiàn)，當(dāng)旁瓣衰減較小時，Chebyshev切趾函數(shù)主瓣功率泄露嚴(yán)重，并且其旁瓣衰減率為零，不利于臨近強(qiáng)干擾信號的抑制。為了更好地利用Chebyshev切趾函數(shù)的優(yōu)點(diǎn)，對Chebyshev切趾函數(shù)的改進(jìn)顯得十分有必要。

2 改進(jìn)的Chebyshev切趾函數(shù)及其分析

通過第1節(jié)對Chebyshev切趾函數(shù)的分析，在眾多切趾函數(shù)中選取Blackman-Harris函數(shù)用于對Chebyshev切趾函數(shù)的改進(jìn)。Blackman-Harris切趾函數(shù)是一種具有良好旁瓣性能的四項(xiàng)系數(shù)三階余弦函數(shù)，與其他切趾函數(shù)相比，它具有以下優(yōu)點(diǎn)：第1，作為一種余弦函數(shù)，其切趾運(yùn)算簡單易行[11]；第2，在同步誤差較小時，旁瓣抑制能力強(qiáng)，通過在Matlab中調(diào)用Wvtool函數(shù)對其幅頻響應(yīng)進(jìn)行分析，發(fā)現(xiàn)其旁瓣電平相對于主瓣可降低-92 dB[12]，能較好的適應(yīng)Chebyshev切趾函數(shù)旁瓣衰減調(diào)節(jié)范圍寬的特點(diǎn)；第3，Blackman-Harris函數(shù)有一定的旁瓣衰減率，對Chebyshev切趾函數(shù)旁瓣衰減率的提高具有改善作用。Blackman-Harris函數(shù)表達(dá)式可由余弦窗的一般表達(dá)式推導(dǎo)得出：

(4)

式(4)為余弦窗的一般表達(dá)式，其中0≤n≤N-1，k為項(xiàng)數(shù)。當(dāng)k=0時，式(4)為矩形切趾函數(shù)，當(dāng)k=1時，為漢寧切趾函數(shù)，當(dāng)k=2時，為Blackman切趾函數(shù)，k=3時，為Blackman-Harris函數(shù)。長度為N的Blackman-Harris函數(shù)其時域表達(dá)式為[12-13]

(5)

式中：a0=0.358 75,a1=0.488 29,a2=0.141 28,a3=0.011 68。為此，在原有可調(diào)參數(shù)γ的基礎(chǔ)上，增加了可調(diào)參數(shù)L(0≤L≤N)，對Chebyshev切趾函數(shù)進(jìn)行改進(jìn)：

wC_B_H(n,γ,N,L)=

(6)

式(6)為改進(jìn)后的Chebyshev切趾函數(shù)時域表達(dá)式，其波形如圖2所示(N=128,L=80)，改進(jìn)的Chebyshev切趾函數(shù)同時擁有Chebyshev與Blackman-Harris函數(shù)的優(yōu)點(diǎn)，通過Matlab計算，其性質(zhì)如表2所示。

圖2 改進(jìn)的Chebyshev切趾函數(shù)Fig.2 Improved Chebyshev apodization function

表2 改進(jìn)的Chebyshev切趾函數(shù)在不同旁瓣衰減時的性質(zhì)

Table 2 Property of improved Chebyshev apodization function at different side-lobe attenuation

m/dBL實(shí)際第一旁瓣衰減/dB功率泄露系數(shù)(%)歸一化主瓣寬度旁瓣衰減率/(dB/otc)-2510-48．20．010．02930-6-2560-18．61．640．02340-12-25120-21．71．060．01563-10-5010-51．300．02930-6-5060-23．70．440．02340-12-50120-41．80．020．01953-10-802-80．800．02930-6-80126-77．300．02539-10

通過表1和表2的對比可以發(fā)現(xiàn)在實(shí)際旁瓣衰減較小時，改進(jìn)后切趾函數(shù)功率泄露系數(shù)得到了明顯的改善，同時通過調(diào)節(jié)參數(shù)L，其仍可以保留Chebyshev切趾函數(shù)歸一化主瓣寬度較窄、第一旁瓣衰減調(diào)節(jié)范圍寬等優(yōu)點(diǎn)，并且旁瓣衰減率得到了增加。改進(jìn)的切趾函數(shù)與原Chebyshev切趾函數(shù)的幅頻響應(yīng)對比如圖3所示。

圖3中，a)為改進(jìn)的切趾函數(shù)與Chebyshev切趾函數(shù)在第一旁瓣衰減43 dB時的比較，其中L=127，m=-46 dB,改進(jìn)的切趾函數(shù)功率泄露系數(shù)為0.07%，原Chebyshev切趾函數(shù)為0.23%。b)為改進(jìn)的切趾函數(shù)與原Chebyshev函數(shù)在第一旁瓣衰減25 dB時的比較，其中L=127，m=-27 dB，改進(jìn)的切趾函數(shù)功率泄露系數(shù)為5.89%，原Chebyshev切趾函數(shù)為16.01%。c)為改進(jìn)的切趾函數(shù)與原Chebyshev函數(shù)在第一旁瓣衰減25 dB時的比較，其中L=101，m=-38 dB, 改進(jìn)的切趾函數(shù)功率泄露系數(shù)為0.4%，原Chebyshev切趾函數(shù)為16.01%，改進(jìn)后函數(shù)旁瓣衰減率達(dá)到了12 dB/otc。

圖3 改進(jìn)的切趾函數(shù)與原Chebyshev切趾函數(shù)在相同旁瓣衰減下的比較(N=128)Fig.3 Comparison between the improved apodization function and the original Chebyshev apodization function under the same side-lobe attenuation (N=128)

從圖3a)，b)，c)的結(jié)果可以看出，通過對參數(shù)γ及參數(shù)L的調(diào)節(jié)，新的切趾函數(shù)在功率泄露系數(shù)、旁瓣衰減率方面得到了很大的改善，從圖3b),c)可以看出，不同的γ及L值可以得到相同的第一旁瓣衰減，但功率泄露系數(shù)、旁瓣衰減率和主瓣寬度會有相應(yīng)的增減，這使得改進(jìn)的Chebyshev切趾函數(shù)可以適應(yīng)不同的切趾要求，應(yīng)用能力更加靈活。在一般情況下，其主瓣寬度較原Chebyshev函數(shù)有所增寬，但是當(dāng)取L=N時(圖3d)，改進(jìn)的切趾函數(shù)恢復(fù)原Chebyshev函數(shù)的效果，主瓣寬度最窄。

3 仿真與討論

為進(jìn)一步驗(yàn)證改進(jìn)型切趾函數(shù)的性能，選取余弦波模擬單色光的干涉圖，首先在Matlab中生成余弦信號，對此信號采樣(采樣頻率100 Hz，采樣點(diǎn)數(shù)512點(diǎn))，然后進(jìn)行切趾、光譜還原處理。則系統(tǒng)理想干涉圖為

IR(x)=100 cos(2π·4x).

(7)

限制于干涉儀只能提供有限的光程差，實(shí)際截斷干涉圖為

(8)

式(8)經(jīng)過反傅里葉變換，即可得到矩形函數(shù)切趾后的光譜圖[14]：

(9)

式中：σ為光譜波數(shù)；BR(σ)表示波數(shù)σ處的光譜強(qiáng)度。

然后用改進(jìn)的切趾函數(shù)和4種典型切趾函數(shù)乘以式(8)，對實(shí)際干涉圖切趾后進(jìn)行式(9)的變換，得復(fù)原光譜如圖4所示。

圖4 幾種典型切趾函數(shù)對單色光干涉圖切趾還原后光譜圖效果對比Fig.4 Contrast effect of restored spectrum of monochromatic light from several typical apodization functions

圖4中，a)為模擬原始干涉圖及其未經(jīng)切趾處理的復(fù)原光譜，b)為矩形切趾及三角切趾后的復(fù)原光譜，實(shí)際應(yīng)用中，限制于光譜儀只能提供有限的光程差，矩形切趾后的干涉圖即為一般意義上的原始干涉圖[15]，從圖b)可以看出，經(jīng)矩形切趾后的復(fù)原光譜旁瓣衰減不足，經(jīng)三角切趾后的復(fù)原光譜旁瓣雖然得到了一定的抑制，但是主瓣寬度增加較多，光譜分辨率下降明顯；c)為布萊克曼與漢明切趾后的復(fù)原光譜，旁瓣抑制能力進(jìn)一步增強(qiáng)，但仍然存在主瓣寬度增加明顯的問題，切趾處理后光譜還原效果并不理想；d)中，2種切趾函數(shù)在相同的第一旁瓣衰減下(-25 dB)，改進(jìn)的Chebyshev切趾函數(shù)旁瓣抑制能力較強(qiáng)，旁瓣衰減率得到了提高的同時，主瓣寬度增加并不明顯，因此，改進(jìn)后的Chebyshev函數(shù)是一種性能優(yōu)異的切趾函數(shù)。

4 結(jié)束語

本文對傅里葉變換光譜儀干涉圖切趾函數(shù)進(jìn)行了研究，提出了改進(jìn)的Chebyshev切趾函數(shù)，利用Chebyshev切趾函數(shù)切趾處理后光譜圖主瓣寬度窄、旁瓣衰減可調(diào)節(jié)的優(yōu)點(diǎn)與Blackman-Harris切趾函數(shù)旁瓣衰減能力強(qiáng)，具有一定旁瓣衰減率的特點(diǎn)相結(jié)合，對Chebyshev切趾函數(shù)進(jìn)行了改進(jìn)。將改進(jìn)后的切趾函數(shù)進(jìn)行分析與仿真驗(yàn)證，結(jié)果表明：①改進(jìn)后的Chebyshev切趾函數(shù)較正了原切趾函數(shù)旁瓣衰減較低時，功率泄露嚴(yán)重的問題；②改進(jìn)后的Chebyshev切趾函數(shù)較正了原切趾函數(shù)旁瓣等波動性，旁瓣衰減率為0的問題；③改進(jìn)后的Chebyshev切趾函數(shù)相比于原切趾函數(shù)更加靈活，可以在相同第一旁瓣衰減量下得到不同的功率泄露效果、旁瓣衰減率和主瓣寬度。

[1] 肖青，柳欽火，李小文，等.高分辨率機(jī)載遙感數(shù)據(jù)的交叉輻射影響及其校正[J].遙感學(xué)報，2005(6): 3-11. XIAO Qing, LIU Qin-huo, LI Xiao-wen,et al. Analysis and Correction of Atmospheric Cross Radiation for High Geometric Resolution Airborne Remote Sensing Data[J].Journal of Remote Sensing, 2005(6): 3-11.

[2] 李蘇寧，朱日紅，李建欣，等.傅里葉干涉成像光譜技術(shù)中的重構(gòu)方法[J].應(yīng)用光學(xué),2009,30(2):268-272. LI Su-ning, ZHU Ri-hong, LI Jian-xin,et al. Method of Reconstruction on Fourier Transform Spectroscopy[J].Journal of Applied Optics, 2009,30(2):268-272.

[3] 景娟娟，相里斌，呂群波，等.干涉光譜數(shù)據(jù)處理技術(shù)研究進(jìn)展[J].光譜學(xué)與光譜分析,2011,31(4):865-870. JING Juan-juan, XIANG Li-bin, Lü Qun-bo,et al.Advance in Interferogram Data Processing Technique[J].Spectroscopy and Spectral Analysis,2011,31(4): 865-870.[4] 楊琨. 傅里葉變換紅外光譜儀若干核心技術(shù)研究及其應(yīng)用[D]. 武漢：武漢大學(xué)，2010: 20-21. YANG Kun. Research and Application on a Certain Number of Core Technologies of Fourier Transform Infrared Spectrometer[D]. Wuhan: Wuhan University, 2010: 20-21.

[5] 《數(shù)學(xué)辭?！肪庉嬑瘑T會. 數(shù)學(xué)辭海(第三卷)[M].北京:中國科學(xué)技術(shù)出版社，2002:373. Mathematic Dictionary Editorial Board. Mathematic Dictionary(Volume III)[M]. Beijing: China Science and Technology Press,2002:373.

[6] 王曉凱，李鋒.Chebyshev 過渡型濾波器設(shè)計和分析[J].信息與電子工程，2009,7(2):115-118. WANG Xiao-kai, LI Feng. Design and Analysis of Chebyshev Transitional Filters[J]. Information and Electronic Engineering, 2009,7(2):115-118.

[7] 孫仲民，黃俊，楊鍵維，等. 基于切比雪夫窗的電力系統(tǒng)諧波/間諧波高精度分析方法[J]. 電力系統(tǒng)自動化，2015, 39(7)：117-122. SUN Zhong-min, HUANG Jun, YANG Jian-wei,et al. A High Accuracy Analysis Method for Harmonics and Inter-Harmonics in Power Systems Based on Dolph-Chebyshev Windows[J]. Automation of Electric Power Systems, 2015, 39(7)：117-122.

[8] 謝俊好，楊曦.天線陣角分辨率準(zhǔn)則的比較研究[J].雷達(dá)科學(xué)與技術(shù)，2010,8(4):382-389. XIE Jun-hao, YANG Xi. Comparative Research on Angular Resolution Criteria for Antenna Array[J]. Rafar Science and Technology, 2010,8(4):382-389.

[9] OPPENHEIM A V，SCHAFER R W，BUCK J R. 離散時間信號處理[M]. 2版.劉樹棠，黃建國,譯. 西安：西安交通大學(xué)出版社,2001. OPPENHEIM A V,SCHAFER R W,BUCK J R.Discrete-Time Signal Processing[M].2nd ed.LIU Shu-tang,HUANG Jian-guo,translated.Xi’an Jiaotong Uni versity Press,2001.

[10] MITRA S K. 數(shù)字信號處理——基于計算機(jī)的方法[M].3版. 孫洪，譯. 北京: 電子工業(yè)出版社，2006:435-436. MITRA S K.Digital Signal Processing——A Computer-Based Aproach[M]. 3rd ed. SUN Hong，Translated，Beijing: Elecronic Industry Press,2006:435-436.

[11] 趙文春，馬偉明，胡安.電機(jī)測試中諧波分析的高精度FFT算法[J].中國電機(jī)工程學(xué)報，2001, 9(12) : 83-87. ZHAO Wen-chun, MA Wei-ming, HU An. FFT Algorithm with High Accuracy for Harmonic Analysis in the Electric Machine[J]. Proceedings of the CSEE, 2001, 9(12) : 83-87.

[12] 劉海升，付志紅，張淮清，等.基于Blackman-Harris相位差校正信號諧波分析方法[J].數(shù)據(jù)采集與處理，2011,26(4):468-472. LIU Hai-sheng, FU Zhi-hong, ZHANG Huai-qing,et al. Harmonic Analysis Based on Blackman-Harris Window and Phase Difference Correction[J]. Journal of Data Acquisition and Processing,2011, 26(4):468-472.

[13] ZHANG F, GENG Z, YUAN W. The Algorithm of Interpolation Windowed FFT for Harmonic Analysis of Electric Power System[J]. IEEE Transactions on Power System, 2001,16(2):160-164.

[14] 尚博. 高光譜干涉成像重構(gòu)技術(shù)研究[D]. 南京：南京理工大學(xué)，2008:23-24. SHANG Bo. The Study of Hyperspectral Imaging Interferometer Reconstruction[D]. Nanjing: Nanjing University of Technology and Engeneering, 2008:23-24.

[15] 韓剛. 干涉光譜復(fù)原算法研究[D]. 西安:中國科學(xué)院研究生院，2012:9-12. HAN Gang. Research on the Algorithm of Inteferogram Spectrum Reconstruction[D]. Xi’an: Chinese Academy of Sciences of Master of Engineering, 2012:9-12.

Improved Chebyshev Function to Raise the Effect of Interferogram Apodization

JIANG Feng, WU Gao-wei, SHENG Wen, FENG Xiao, FENG Zhi

(Air Force Early Warning Academy, Hubei Wuhan 430019, China)

In order to rectify the shortcomings of high power leakage rate at low side-lobe attenuation and identical fluctuation of side-lobe attenuation of Chebyshev apodization function, an improved Chebyshev apodization function is put forward. Firstly, the defects of Chebyshev apodization function in spectrum restoration is analyzed by simulation; and then, the Blackman-Harris function is selected to improve it and compared with the previous one; at last, the function is used for monochromatic interferogram apodization process and simulation. The simulation results show that the new Chebyshev apodization function can solve the problems of high power leakage rate at low side-lobe attenuation and identical fluctuation of side-lobe attenuation, and its application ability is more flexible.

Chebyshev; apodization function; Blackman-Harris; interferogram; spectrum reconstruction; simulation analysis

2016-04-05；

2016-07-10 作者簡介：江峰(1992-)，男，四川眉山人。碩士生，主要研究方向?yàn)榧t外光譜復(fù)原技術(shù)研究。

10.3969/j.issn.1009-086x.2017.02.025

TN911.1；TH744.1

A

1009-086X(2017)-02-0160-06

通信地址：430019 湖北省武漢市江岸區(qū)黃浦大街288號空軍預(yù)警學(xué)院研究生管理大隊21隊E-mail：969656256@qq.com