基于M orlet小波函數的IDT建模與仿真

文理為，曾志超，牟慶豐，王代強

（貴州大學 大數據與信息工程學院，貴州 貴陽 550025）

基于M orlet小波函數的IDT建模與仿真

文理為，曾志超，牟慶豐，王代強

（貴州大學 大數據與信息工程學院，貴州 貴陽 550025）

為了得到更精確的IDT結構模型，本文提出了基于Morlet小波函數的IDT數學模型，通過Matlab對其頻響特性進行仿真分析，并與改進型δ函數模型IDT結構進行比較。結果表明，在相同的參數情況下，校正后的Morlet小波式IDT結構主瓣保持得很好，具有很好的通帶性和主瓣波形，相對于改進型δ函數模型IDT結構，校正后的Morlet小波式IDT結構旁瓣電平得到大大的抑制，聲電轉換效率更高。

Morlet小波函數；IDT；Matlab；δ函數模型；頻率響應

隨著聲表面波技術的快速發展，聲表面波器件在電子信息等領域的應用越來越廣泛，與此同時，為了滿足現代通信技術高頻化、微型化等特點，對聲表面波器件的要求也越來越高：更高的響應頻率、更高的帶寬、更低的插入損耗以及更小的尺寸等［1-2］。叉指換能器（IDT）作為聲表面波器件最重要的部分之一，要改善聲表面波器件的特性，提高其性能，就必須從改進IDT結構入手［3-4］。首先建立IDT數學模型，再對其進行模擬仿真，并與傳統模型仿真結果進行對比，改進模型參數以滿足IDT結構優化的目的［5］。Morlet小波函數是一種單頻復正弦調制高斯波，具有較好的時頻特性對稱性，在頻域具有很好的線性相位特性，且其時頻域局部性非常好，是最常用的復值小波，是分析信號的強有力工具［6］。文中是在IDT脈沖響應模型和改進型δ函數模型的基礎上選擇Morlet小波函數對IDT進行數學建模，再利用傳統數學模型對同一IDT結構進行仿真對比，結果表明：校正后的Morlet小波函數的IDT模型比改進型的δ函數IDT模型更精確，達到預期結果。

1　M orlet小波函數型IDT的建立

1.1M orlet小波函數

小波函數有正交性和非正交性之分，Morlet小波函數是一種非正交波，不存在尺度函數，其數學表達式為［7］：

其中，ω0為小波的中心頻率。對（1）式進行傅里葉變換，得其頻率響應：

由于小波函數是處理信號的工具，SAW是變化緩慢的信號，又由于小波函數的時頻窗是靈活伸縮的，因此用此特性就能靈活的設計IDT結構［8］。

1.2改進型IDT的δ函數模型

貴州大學王代強教授提出了一種改進型δ函數模型［9］，如圖1（b）所示，在IDT結構中用邊長為aik的若干正方形來表示一根指條，并將每個正方形看作一個獨立的聲波激發源，假設有k個正方形，每個聲波激發源用基本的δ函數單元表示，再將整根叉指指條激發的SAW進行歸一化處理。

為了表達得更明確，選取一個規則型等間距的IDT結構進行分析。在IDT結構的中心取一坐標原點x=0，如圖1（a）所示。同時，由于IDT結構可看作橫向濾波器，則對其場源分布進行傅里葉變換可求出整個IDT結構的頻率響應為［10］：

上式中，ai表示第i根叉指電極的寬度，di表示第i根叉指電極的長度，xi表示第i根叉指電極的位置，vs表示SAW的傳播速度，Ii表示第i根叉指電極的聲波幅值。

圖1　改進型IDT結構δ函數模型

1.3M orlet小波式IDT結構模型

文中是在IDT的脈沖響應模型和改進型δ函數模型的基礎上提出的基于Morlet小波函數的IDT模型。構造Morlet小波函數IDT的基本思想是IDT可看作橫向濾波器，因此，首先將想求解的的頻率響應函數H（f）進行變換，得到相應的脈沖響應函數h（t）［11］。根據抽樣定理，以1/4f0對脈沖響應函數h（t）均勻抽樣，然后用小波函數φ（t）截斷，在小波函數窗口內得到有限長度φ（tn）·h（tn）。為了對IDT最終頻響進行精度調整，要求計算的頻率響應比常規方法計算的精度更高，因此采用時域上對脈沖響應直接奇偶分解與小波窗函數截斷相結合的方法研究IDT結構［12-13］。根據網絡脈沖響應函數理論，因為IDT的脈沖響應與它的結構之間有著密切的關系，很自然地想到用脈沖響應函數去描述分析IDT結構。在濾波器理論的基礎上，脈沖響應和頻率響應之間互為傅里葉變換對，可由需要的頻響函數得到相應的脈沖響應函數［14］。

當Morlet小波函數的伸縮因子c=1，平移因子d=0時，（4）式中的脈沖響應函數h（t）對應的IDT結構為Morlet母小波式IDT結構模型。（4）式中的脈沖響應表達式是頻率響應函數校正后的結果，若不進行頻率響應函數的校正，那么h（t）結果為：

我們知道，在改進型的δ函數模型中，為了計算更加簡單方便，選取IDT結構的中心為一坐標原點x=0。同時，對于一個電極間隔b均勻的IDT結構，不管IDT結構是否被加權，它的頻率響應都是對稱的，因此為了計算簡單方便，也取Morlet母小波式IDT中心為坐標原點，則整個Morlet小波式IDT的脈沖響應為：

采用1/4f0對式（6）抽樣，其抽樣序列為：

由上式可知，整個Morlet小波式IDT結構的脈沖響應函數是奇偶分量與校正函數之和，當中心頻率為f0時，只剩下偶分量和頻率響應的調整函數。根據脈沖響應函數與頻率響應函數之間的關系，當IDT結構的聲孔徑長度按照Morlet母小波函數設計時，小波式IDT結構脈沖響應函數偶分量所對應的頻率響應為：

式中，Bi為第i根叉指電極的校正函數幅度，可以取正負值，fi為第i根叉指電極校正函數的調整點頻率。將（9）與（10）兩式相加，即：

2　IDT的M orlet小波函數模型仿真

2.1IDT的M orlet小波函數模型結構參數

為了能夠設計出IDT結構，求解出小波式IDT的叉指電極位置以及相鄰叉指間隔很關鍵，同時還需確定IDT結構的叉指電極數和中心頻率f0。而叉指換能器的指條位置與基片材料的選取以及時間有關，其關系如下：

在IDT結構中，基片材料確定后，SAW的傳播速度vs也就確定了。因此中心頻率f0可由下式確定：

對于均勻的IDT結構，通常取ai=bi。我們知道叉指電極數n也是決定IDT結構性能好壞的參數之一，叉指電極數不但影響著聲波的激發強弱，也影響著帶寬，在IDT結構中要根據實際情況來確定其個數。

2.2仿真分析

在IDT結構中，設每一根叉指電極有k個基本δ信號源，而且每根叉指電極產生的信號源強度與其聲孔徑Wi（t）成正比，由于SAW的激發強度僅與聲孔徑有關，因此在本文中，將其聲孔徑按照Morlet小波函數設計時，聲孔徑是可變化的，這樣就能夠得出不同的小波式IDT結構［15］。最后結合改進型δ函數模型以及IDT結構的脈沖響應模型得出其小波式IDT結構的頻率響應。

根據改進型δ函數對規則IDT結構的建模，采用Matlab數學工具，其中取叉指寬度a=7.5μm、叉指間隔b=7.5μm、叉指對數N=45，根據規則型IDT的頻率響應表達式（3），得其結構的頻率響應，如圖2所示。同時保持參數值不變，根據（6）式，用Matlab數學工具，仿真得到Morlet母小波式IDT的頻率響應圖形，如圖3所示。

圖2　規則型IDT結構頻響仿真圖

對比圖2和圖3，可見Morlet小波式IDT在與規則的δ函數型IDT具有同樣的叉指電極時，它們的帶寬大體一樣。但是Morlet小波式IDT結構的通帶性及主瓣波形較好，相對于規則型IDT結構，Morlet小波式IDT結構的第一旁瓣電平抑制從38.5 dB提高到41.9 dB，第二旁瓣電平抑制從42.1 dB提高到43.7 dB，第三旁瓣電平抑制從43.6 dB提高到46.9 dB。

圖3　Morlet小波式IDT頻率響應仿真圖

若對 Morlet小波式 IDT結構頻響進行校正，根據式（11），所取參數值與上面的相同，用Matlab作為數學工具將頻率響應仿真出來，可以得到校正后的Morlet母小波式IDT結構的頻率響應如圖4所示。

圖4　Morlet小波式IDT頻率響應校正后的仿真圖

對比圖3、圖4，經過校正后的Morlet小波式IDT結構與規則型IDT結構在叉指電極對數N相同、指寬以及間距都保持一樣時，帶寬能夠保持大體相同，同時，其主瓣可以保存的更好，衰減更小，有較好的通帶性，而且旁瓣電平得到了更好的抑制，第一旁瓣電平抑制從38.5 dB提高到57.1 dB，第二旁瓣電平抑制從42.1 dB提高到62.7 dB，第三旁瓣電平抑制從43.6 dB提高到66.9 dB。由于在IDT結構中其頻率響應的諧振峰高度與叉指電極對數N成正比，寬度與叉指電極對數N成反比，可見校正后的Morlet小波式IDT結構具有更高的聲電轉換效率，頻率響應特性更好。在生產中，可以根據實際情況，根據帶寬、中心頻率、實際濾波的需要有選擇的挑選所需要的IDT結構。

3　結　論

文中是在IDT的脈沖響應模型和改進型的δ函數模型的基礎上提出了Morlet小波函數的IDT數學模型，采用Matlab對IDT的Morlet小波函數模型以及改進型δ函數模型的頻響特性進行仿真分析、比較。結果表明，改進型δ函數型IDT結構和校正后的Morlet小波式IDT結構在叉指電極對數N相同、指寬以及間距都保持一樣時，可以取得相同的帶寬，同時，校正后的Morlet小波式IDT結構主瓣保持得更好，衰減更小，通帶性也很好，相對于改進型δ函數型IDT結構，校正后的Morlet小波式IDT結構旁瓣電平得到大大的抑制，其第一旁瓣電平抑制提高了18.6 dB，第二旁瓣電平抑制提高了20.6 dB，第三旁瓣電平抑制提高了23.3 dB，這一結果說明Morlet小波式IDT結構具有比改進型δ函數型IDT結構更高的精確度。

［1］徐希嬪，王代強，劉橋.階梯型IDT結構的仿真研究［J］.現代電子技術，2009（6）:176-177.

［2］王代強，姚祖明，陳雨青，等.圓形IDT結構聲表面波濾波器的設計與實現［J］.電子元件與材料，2013，32（4）:50-53.

［3］王代強，陳雨青，徐希嬪，等.變跡加權三角形IDT結構的SAW濾波器設計與實現［J］.西南大學學報：自然科學版，2011，33（7）:31-35.

［4］王代強，陳雨青，姚祖明，等.硅基AlN橢圓形IDT-SAW濾波器的設計［J］.傳感技術學報，2013，26（12）:1691-1694.

［5］魏培永，朱長純，劉君華.用聲表面波器件實現小波變換［J］.壓電與聲光，2001，23（3）:173-175.

［6］羅光坤.Morlet小波變換理論與應用研究及軟件實現［D］.南京：南京航空航天大學，2007.

［7］何嶺松，李巍華.用Morlet小波進行包絡檢波分析［J］.振動工程學報，2002，15（1）:119-121.

［8］LI Yuan-yuan，LU Wen-ke，GONG Lin-jia.The research of wavelet transformation technology by using surface acoustic wave devices［J］.IEEE Conference Publications，2011，3:926-929.

［9］王代強，陳雨青，徐希嬪，等.SAW器件叉指換能器δ函數模型的改進［J］.電子元件與材料，2010，29（9）:62-65.

［10］Hashimoto K，Yamaguchi M.Derivation of coupling-of-modes parametersforSAWdeviceanalysisbymeansof boundaryelment method［J］.IEEE Ultrason.Symp.，1991，1: 21-26.

［11］張曉慶.高頻聲表面波濾波器的納米壓印技術制備研究［D］.武漢：華中科技大學，2011.

［12］趙義禎.聲表面波濾波器的設計及其衍射補償的研究［D］.杭州：浙江工業大學，2007.

［13］董輝平.聲表面波傳感器模擬與仿真［D］.太原：中北大學，2008.

［14］馮亞菲.聲表面波器件的制備及其性能研究［D］.哈爾濱：黑龍江大學，2008.

［15］董溯前.聲表面波式小波變換器件的研究［D］.上海：東華大學，2008.

Modeling and simulation of IDT based on Morletwavelet

WEN Li-wei，ZENG Zhi-chao，MOU Qing-feng，WANG Dai-qiang
（College of Big Data and Information Engineering，GuiZhou University，Guiyang 550025，China）

In order to obtain a more accurate IDT（Inter Digital Transducer）structure，a new IDT mathematicalmodel is presented，based on Morlet wavelet function.Compared with the modifiedδfunction model IDT of its frequency response characteristics by Matlab.The results shows，in the case of the same parameters，the adjusted Morletwavelet IDT'smain lobe kept well，with good passband and the main lobe of the waveform.Compared with the modifiedδfunction model IDT，the sidelobe levelwasgreatly suppressed，and ithasa higheracoustic-electric conversion efficiency.

Morletwavelet function；IDT；Matlab；δfunctionmodel；frequency response

TN65

A

1674-6236（2016）19-0007-04

2016-03-14稿件編號：201603158

國家自然科學基金項目（11564005）

文理為（1992—），男，湖南醴陵人，碩士研究生。研究方向：聲表面波器件，電路與系統。