壓電陶瓷測量原理

1、壓電陶瓷及其測量原理近年來，壓電陶瓷的研究發(fā)展迅速，取得一系列重大成果，應(yīng)用范圍不斷擴大，已深入到國民經(jīng)濟和尖端技術(shù)的各個方面中，成為不可或缺的現(xiàn)代化工業(yè)材料之一。由于壓電材料這就使得壓電陶瓷材料的的各向異性，每一項性能參數(shù)在不同的方向所表現(xiàn)出的數(shù)值不同,性能參數(shù)比一般各向同性的介質(zhì)材料多得多。同時，壓電陶瓷的眾多的性能參數(shù)也是它廣泛 應(yīng)用的重要基礎(chǔ)。（一）壓電陶瓷的主要性能及參數(shù)（1）壓電效應(yīng)與壓電陶瓷在沒有對稱中心的晶體上施加壓力、張力或切向力時，則發(fā)生與應(yīng)力成比例的介質(zhì)極化,同時在晶體兩端將出現(xiàn)正負(fù)電荷，這一現(xiàn)象稱為正壓電效應(yīng)；反之，在晶體上施加電場時,則將產(chǎn)生與電場強度成比例的變形或機

2、械應(yīng)力,這一現(xiàn)象稱為逆壓電效應(yīng)。這兩種正、逆壓即晶電效應(yīng)統(tǒng)稱為壓電效應(yīng)。 晶體是否出現(xiàn)壓電效應(yīng)由構(gòu)成晶體的原子和離子的排列方式,體的對稱性所決定。在聲波測井儀器中， 發(fā)射探頭利用的是正壓電效應(yīng)，接收探頭利用的是 逆壓電效應(yīng)。（2）壓電陶瓷的主要參數(shù)1 、介質(zhì)損耗介質(zhì)損耗是包括壓電陶瓷在內(nèi)的任何電介質(zhì)的重要品質(zhì)指標(biāo)之一。在交變電場下，電介質(zhì)所積蓄的電荷有兩種分量：一種是有功部分（同相），由電導(dǎo)過程所引起；另一種為無功部分tan（異相），由介質(zhì)弛豫過程所引起。介質(zhì)損耗是異相分量與同相分量的比值，如圖1C為同相分量，Ir為異相分量，心與總電流I的夾角為 ，其正切值為 -Cr 其中3為交變電場的角頻率

3、，R為損耗電阻，C為介質(zhì)電容。1 CCR2、機械品質(zhì)因數(shù)機械品質(zhì)因數(shù)是描述壓電陶瓷在機械振動時，材料內(nèi)部能量消耗程度的一個參數(shù)，它也是衡量壓電陶瓷材料性能的一個重要參數(shù)。機械品質(zhì)因數(shù)越大， 能量的損耗越小。產(chǎn)生能量損耗的原因在于材料的內(nèi)部摩擦。機械品質(zhì)因數(shù)Qm的定義為:Q諧振時振子儲存的機械能CQm諧振時振子每周所 損失的機械能 2兀機械品質(zhì)因數(shù)可根據(jù)等效電路計算而得1s L 1式中Ri為等效電阻（Q）s R1C 1R1s為串聯(lián)諧振角頻率（Hz），Ci為振子諧振時的等效電容（F）, Li為振子諧振時的等效電感。Qm與其它參數(shù)之間的關(guān)系將在后續(xù)詳細(xì)推導(dǎo)。不同的壓電器件對壓電陶瓷材料的Qm值的要求

4、不同，在大多數(shù)的場合下（包括聲波測井的壓電陶瓷探頭），壓電陶瓷器件要求壓電陶瓷的 Qm值要高。3、壓電常數(shù)壓電陶瓷具有壓電性， 即在其外部施加應(yīng)力時能產(chǎn)生額外的電荷。其產(chǎn)生的電荷與施加D （單位面積的電荷）和應(yīng)的應(yīng)力成比例，對于壓力和張力來說，其符號是相反的，電位移力的關(guān)系表達式為：D Q drA式中Q為產(chǎn)生的電荷（C），A為電極的面積（n2），d為壓電應(yīng)變常數(shù)（C/N）。在逆 壓電效應(yīng)中，施加電場E時將成比例地產(chǎn)生應(yīng)變 S，所產(chǎn)生的應(yīng)變 S是膨脹還是收縮， 取 決于樣品的極化方向。S=dED S兩式中的壓電應(yīng)變常數(shù) d在數(shù)值上是相同的，即 d -E另一個常用的壓電常數(shù)是壓電電壓常數(shù)g，它表示

5、應(yīng)力與所產(chǎn)生的電場的關(guān)系，或應(yīng)變 與所引起的電位移的關(guān)系。常數(shù) g與d之間有如下關(guān)系： g -式中 為介電系數(shù)。在聲波測井儀器中， 壓電換能器希望具有較高的壓電應(yīng)變常數(shù)和壓 電電壓常數(shù)，以便能發(fā)射較大能量的聲波并且具有較高的接受靈敏度。4、機電耦合系數(shù)當(dāng)用機械能加壓或者充電的方法把能量加到壓電材料上時，由于壓電效應(yīng)和逆壓電效。這種轉(zhuǎn)換的強弱用機電耦合應(yīng)，機械能（或電能）中的一部分要轉(zhuǎn)換成電能（或機械能） 系數(shù)k來表示，它是一個量綱為一的量。機電耦合系數(shù)是綜合反映壓電材料性能的參數(shù),它表示壓電材料的機械能和電能的耦合效應(yīng)。機電耦合系數(shù)的定義為:k2=電能轉(zhuǎn)變?yōu)闄C械能輸入電能或者k2=機械能轉(zhuǎn)變?yōu)?/p>

6、電能輸入機械能機電耦合系數(shù)不但與材料參數(shù)有關(guān)，還與具體壓電材料的工作方式有關(guān)。對于壓電陶瓷 來說，它的大小還與極化程度相關(guān)。它只是反映機、電兩類能量通過壓電效應(yīng)耦合的強弱,當(dāng)制并不代表兩類能量之間的轉(zhuǎn)換效率。壓電材料的耦合系數(shù)在不同的場合有不同的要求, 作換能器時，希望機電耦合系數(shù)越大越好。（二）壓電換能器的等效電路壓電換能器的等效電路表示法，是利用電學(xué)網(wǎng)絡(luò)術(shù)語表示壓電陶瓷的機械振動特性,把某些力學(xué)量模擬為電學(xué)量的方法。把壓電換能器用等效電路來表示，有很多優(yōu)點：其一,可以把力學(xué)上復(fù)雜的振動問題有效地進行簡化;其二，為了得到換能器的各個參數(shù),從而定量地分析或篩選換能器；其三，實際應(yīng)用的需要，因為

7、在實際的應(yīng)用當(dāng)中，壓電換能器也是 接入到具體的電子線路中的，得到壓電換能器的等效電路能夠更好地對其外圍電路進行匹配 設(shè)計。由此可見，得到壓電換能器的等效電路是十分必要的。2.3壓電換能器的諧振特性將壓電換能器按照圖2-2所示線路連接。當(dāng)改變信號頻率時，可以發(fā)現(xiàn)，通過壓電陶瓷換能器的電流也隨著發(fā)生變化，其變化規(guī)律如圖2-3 （a）所示。從圖2-3（a）可以看出,當(dāng)信號為某一頻率 fm時，通過壓電陶瓷換能器的電流出現(xiàn)最大值Imax ;而當(dāng)信號變到另頻率fn時，傳輸電流出現(xiàn)最小值Imin。由流經(jīng)它的電流隨頻率的變化可以看出，壓電陶瓷2-3(b)所示。換能器的阻抗是隨頻率的變化而變化的，其變化規(guī)律同電

8、流相反，如圖圖2-2壓電陶瓷換能器諧振特性接線示意圖圖2-3 壓電陶瓷換能器電流、阻抗同頻率的關(guān)系曲線電流-頻率關(guān)系曲線(b) 阻抗-頻率關(guān)系曲線從圖中可以看出，當(dāng)信號頻率為fm時，通過壓電陶瓷換能器的電流最大，即其等效阻抗最小，導(dǎo)納最大；當(dāng)信號頻率為fn時，通過壓電陶瓷換能器的電流最小，即其等效阻抗最大，導(dǎo)納最小。因此把 fm稱為最大導(dǎo)納頻率或最小阻抗頻率；而把fn稱為最小導(dǎo)納頻率或最大阻抗頻率。而當(dāng)信號頻率繼續(xù)增大時，還會出現(xiàn)一系列的電流的極大值和極小值，如 圖2-4所示。圖2-4 壓電陶瓷換能器電流隨頻率變化示意圖（多諧振模式）2.2.4壓電換能器的等效電路根據(jù)交流電路相關(guān)知識，對于圖2

9、-5所示好的LC電路來說，其阻抗Z也隨著頻率的變化而變化。在圖 2-2所示的線路中，用LC電路代替壓電陶瓷換能器，可以發(fā)現(xiàn)，在壓電陶瓷換能器的諧振頻率處，只要選擇合適的Li、Ci、R和Co，通過LC電路的電流和LC電路的阻抗的絕對值隨頻率的變化曲線，分別同圖2-1中的（b）和（c）的關(guān)系曲線非能量損耗可用一個并聯(lián)電阻Rj來等效。所以其最終等效電路圖如圖2-6所示。常相似。也就是說，在串聯(lián)諧振頻率附近， 壓電陶瓷換能器的阻抗特性和諧振特性同LC電路的阻抗特性和頻率特性非常相似。因此，利用機電類比的方法，可以用一個LC電路來表示壓電陶瓷換能器的參數(shù)和特性，這個LC電路即為壓電陶瓷換能器的等效電路。

10、L,pnrnC.4Co圖2-5 LC 電路對壓電陶瓷換能器來說，在任何串聯(lián)諧振頻率附近，其電行為可以用圖2-3所示的LC電路來表示。在壓電陶瓷換能器的串聯(lián)諧振頻率附近，如果值存在一種振動模式，即沒有其它寄生響應(yīng),則在串聯(lián)諧振頻率附近很窄的頻率范圍內(nèi)，可以認(rèn)為壓電陶瓷換能器的等效參數(shù) R、Ci、Ri和Co與頻率無關(guān)。在實際中通過選擇合適的尺寸進行加工處理，是可以將所需要的振動模式同其他模式充分隔離開來的。另外，考慮到在實際中，在通電之后，壓電陶瓷換能器必然會存在能量的損耗，這一2111_1CIIIIRo1.1Cl圖2-6壓電陶瓷換能器等效電路圖圖中串聯(lián)支路中的Li稱為壓電陶瓷換能器的動態(tài)電感，C

11、i稱為動態(tài)電容，Ri稱為動態(tài)電阻。這三個參數(shù)用來表征壓電陶瓷換能器在工作（加電源激勵產(chǎn)生振動） 的情況下，振動部分所受到的力阻抗和介質(zhì)對振動的反作用的強弱。并聯(lián)電容 Co又稱靜態(tài)電容，表征壓電陶瓷換能器在未加激勵的情況下等效為一個純電容,它的值的大小與換能器的形狀有關(guān)?？紤]到換能器的靜態(tài)電阻并不為零，則實際中的導(dǎo)納圓不可能與縱軸相切，而是向橫軸的正聯(lián)電阻 Ro又稱靜態(tài)電阻，表征換能器的電損耗的大小。225壓電換能器的導(dǎo)納特性根據(jù)已得到的壓電換能器的等效電路圖，來進一步分析其導(dǎo)納特性。 為了簡化推導(dǎo),假定壓電陶瓷換能器沒有電損耗，即ROnO,此時其等效電路即為一個LC電路，如圖2-5所示。則Y

12、Yo Yi(2-1)式中：丫為換能器的總的導(dǎo)納值,YojBo j 0為并聯(lián)支路的導(dǎo)納值,丫 GijBi為串聯(lián)支路的導(dǎo)納值。先對串聯(lián)支路進行分析。Yo jBo j 0GiRi2i 2Ri( Li -)若令Ri( LiRiRi2,Bij(j *得到:(LiCi(Li 丄)2Ri (Jx Ri2LiCii)2(2-2)iLi ；）。由式（2-2）Ci可得：R1X G1Bi（X2R1 )2G12B1G 2 G1G1巨（1 ）2 加上京）1 2(G1 藥)Bi1 2肩(2-3)若以電導(dǎo)為橫坐標(biāo)，電納為縱坐標(biāo)，則式一 12-3 ）表示一個以（,0 ）為圓心，2R12R1為半徑的圓，也即是我們所說的導(dǎo)納圓。

13、如圖2-7中虛線所示圖2-7導(dǎo)納圓圖對于串聯(lián)支路進行分析，根據(jù)串聯(lián)諧振頻率的定義，令Bi=0，則由式（2-3 ）可得到Gi=0或Gi 1/Ri 。由于實際的壓電陶瓷換能器的動態(tài)電阻Ro不可能為零，根據(jù)式（2-2）中G1的表達式可以知道，只有 G11/R1滿足串聯(lián)諧振的條件。即:L120，所以可以得到串聯(lián)支路的諧振頻率（又稱機械共振頻率）TLCi (2-4)接著考慮加入靜態(tài)電容后的情況。由式（2-1 ）可知，考慮靜態(tài)電容后換能器的導(dǎo)納相當(dāng)于在串聯(lián)支路的電納（虛部）加上丫0。鑒于一般情況下，壓電陶瓷換能器的機械品質(zhì)因數(shù)都較大，也即在串聯(lián)諧振頻率s附近，丫0 j C0的值隨頻率的變化很小，可以近似認(rèn)

14、而橫坐2-7中點劃線所示。若再為是一個常數(shù)。因此，只需將串聯(lián)支路所得到的導(dǎo)納圓的縱坐標(biāo)向上平移一個常數(shù), 標(biāo)保持不變即可得到加入靜態(tài)電容后換能器的導(dǎo)納關(guān)系圖，如圖，如圖2-7中實線圓所示向平移一定的量（平移距離的大小取決于靜態(tài)電阻的阻值）對導(dǎo)納圓圖進行簡要的分析可知：當(dāng)f fs即s時，電納值大于零，當(dāng)fs即s時，電納值小于零。所以，隨著頻率的增加，導(dǎo)納圓是沿順時針方向變化的。另外,在串聯(lián)諧振頻率的附近， 還存在著兩個頻率點使得換能器總的電納為零,此時電源信號經(jīng)過換能器之后只有幅值的改變，而沒有相位的變化，也即電壓和電流信號同相位。這兩個頻率中，值較小的那個頻率fr稱為諧振頻率，較大 fa的稱為

15、反諧振頻率。另外還存在使得換能器的導(dǎo)納值取得最大的頻率 fm，導(dǎo)納值最小的頻率 fn。連接原點和串聯(lián)諧振頻率點,與導(dǎo)納圓的交點處的頻率 fp稱為并聯(lián)諧振頻率。另外，需要特別指出的是， 上述討論是在一個振動模態(tài)諧振頻率s附近較小的頻率變化范圍內(nèi)進行的，并且只有在導(dǎo)納圓的直徑遠(yuǎn)大于這個頻率范圍內(nèi)Co的變化時才是正確的，否則換能器的導(dǎo)納曲線將變得十分復(fù)雜，具有蔓葉曲線的特征。根據(jù)以上導(dǎo)納圓圖的推導(dǎo)過程， 下面介紹一下壓電陶瓷換能器等效電路中各個參數(shù)和導(dǎo)納圓圖的關(guān)系，并給出各自的計算公式。在換能器的導(dǎo)納圓圖中作平行于縱軸的直徑，交導(dǎo)納圓于兩點，分別記作f1、f2。在f1點處，串聯(lián)支路的動態(tài)電導(dǎo)和電納值

16、相等，即Gi Bi。由式（2-2 ）可得:2Ri (iLiR _十)1C11(iLi -C)1C1(2-5)R1( 1L1 十)21C1iCi在f2可得：2RRi(2L1結(jié)合式(2-5)Li再由式（2-4）2L112C1士）"（兒亡）2和式（2-6），可得:(2-7)(2-6)可得:點處，串聯(lián)支路的動態(tài)電導(dǎo)和電納值相等，但符號相反，即GiBi。由式（2-2）(2-8)機械品質(zhì)因數(shù):QmsGRsL1(2-9)結(jié)合式（2-7 ）和（2-8 ）可得:Qm1 f2式（2-5）和式（2-6 ）消去Ri得到:LiCi 則j f f所以：Qm 亠丄二 （2-10）f2 f1動態(tài)電阻的值可以通過導(dǎo)納

17、圓的直徑求得:(2-11)Ro 1/D靜態(tài)電容 Co的值也可由導(dǎo)納圓偏離橫軸的距離來確定:Co匹（2-12）s式中y0為圓心的縱坐標(biāo)。靜態(tài)電阻Ro的值可由導(dǎo)納圓偏離縱軸的距離（或圓心的橫坐標(biāo)）來確定:RoXo12R(2-13)式中Xo為圓心的橫坐標(biāo)。至此，我們已得到壓電陶瓷換能器等效電路中所有參數(shù)的計算公式。2.3測量原理也即在上一節(jié)中，得到的壓電陶瓷換能器等效電路參數(shù)的計算公式都是基于導(dǎo)納圓的,是基于各個頻率下的電導(dǎo)和電納值的，因此我們需要得到每個頻率點的導(dǎo)納值。為此采用圖 2-8所示的測量原理圖進行測量。lintVnNr圖2-8壓電陶瓷換能器測量原理示意圖* *圖2-8中，AC為頻率可控的

18、交流信號源，R表示源內(nèi)阻，Rm稱為精密電阻，Ua為加在壓電陶瓷換能器山的電壓信號Ub為經(jīng)過換能器之后的電壓信號。根據(jù)前面章節(jié)所介U B相對于U A會有一個UA UAmej，Ub UBmej()(2-14)其中：UAm，UBm分別表示兩路信號的幅值,為信號的角頻率，為信號的初始紹的壓電陶瓷的導(dǎo)納特性可以知道，在經(jīng)過換能器之后的電壓信號 幅度和相位的變化。不失一般性，在這里設(shè)定:相位，為兩路信號的相位差。按照習(xí)慣表達，先求壓電陶瓷換能器的阻抗,再取倒數(shù)得到導(dǎo)納。_ U Ua UbI Ub, /Rm將式（2-14 ）代入得U Ame U BmeUBmej(Rm 護 e jRmUBm1)=Rm (護

19、COSU Bm1)sin jX對應(yīng)得到：RRm (嚴(yán) cosU Bm1），Rm 詈 sinU Bm(2-15)再由導(dǎo)納和阻抗的關(guān)系可得Y 丄 一1一 G jBZ R jXcRc即: G R才，B R2 X 2R2 X2(2-16)由以上推導(dǎo)可以看出，換能器的導(dǎo)納和阻抗值僅與加在其兩端的電壓信號的幅值比和相位差有關(guān)，因此只需要得到兩路信號的幅值和相位信息即可得到換能器等效電路的各個參數(shù)。而實際中，只需要對兩路信號進行采樣，再通過對采樣所得數(shù)據(jù)進行處理便可得到幅值和相位信息。2.4正弦信號的測量方法根據(jù)上一節(jié)介紹的測量原理可知，要得到壓電換能器在測試頻率下的電導(dǎo)和電納值,需要測得其兩端正弦信號的幅

20、值比和相位差。但是實際中，硬件電路實現(xiàn)的僅是對兩路信號的A/D轉(zhuǎn)換采集，也即是得到的是兩路正弦信號的一系列的離散的點。在這一節(jié)中，將介紹從這些采集到的離散的點計算其幅值和相位的方法。2.4.1數(shù)字相關(guān)法隨著微處理器和大規(guī)模集成電路的迅速發(fā)展，在測試系統(tǒng)中，越來越多的傳統(tǒng)的測量方法被數(shù)字化測量方法所取代。近年來，由于相關(guān)函數(shù)法具有提高測試精度，減少或簡化硬件設(shè)計，能夠充分利用測試系統(tǒng)中的數(shù)據(jù)采集系統(tǒng)和微型計算機,提高測試系統(tǒng)的可靠性和可維護性的諸多優(yōu)點，使得相關(guān)技術(shù)原理在相位差的測量及數(shù)字信號處理中得到了廣泛應(yīng)用,并展現(xiàn)出良好的應(yīng)用前景1、相關(guān)函數(shù)法原理相關(guān)函數(shù)法利用兩同頻正弦信號的延時為零時的

21、互相關(guān)函數(shù)值與其相位差的余弦值成正比的原理獲得相位差。設(shè)兩路被測信號為:x(t) Asin(2 ft) Nx(t)，y(t)Bsin(2 ft ) Ny(t)(2-17)其中：A、B分別表示兩路信號的幅值，f表示信號的頻率，Nx(t)、 Ny(t)分別表示兩路信號的干擾噪聲信號，表示兩路信號的相位差。顯然，信號x(t)和y(t)是相關(guān)的,則兩路信號的互相關(guān)函數(shù)為:1 TRxy( ) T 0 x(t)y(t1 TRxy( ) T 0Asin(2)dtft) Nx(t) Bsin(2f(t)(2-18)Ny(t )dt1式中T為信號的周期，即T -=0時，有1 TRxy(0) T 0Asin(2

22、ft) Nx(t) Bsin(2 ft ) Ny(t)dt由于噪聲信號之間不相關(guān)，噪聲和信號之間也不相關(guān)，將上式進一步展開得:Rxy所以,1 T-0 ABsin(2 ft)sin(2 ft )dt1 T-0 ABsin(2 ft)sin(2 ft)cos1 T-0 ABsin(2 ft)sin(2 ft) cos dt可以得到相位差的計算公式：cos(2 ft) sin dt1AB cos22Rxy(0)arccos()AB(2-19)而信號幅值的大小可由信號的自相關(guān)函數(shù)求得1 -Rx( ) T TX(t)X(t2r Asin(2 ft)Asin(2 f(t )dt (2-20)T =0時，有

23、Rx(0) T -T 2tA2sin2(2 ft)dtA2所以可得信號幅值的計算公式:A J2FU0)B j2Ry(0)(2-21)將上式代入式（2-19），可得相位差計算公式的另一種表達式:Rxy）arccos(寸 Rx(0)Ry (0)(2-22)而在實際中，是沒有完整精確的信號的表達式的，有的是對信號的模數(shù)轉(zhuǎn)換所得到的離 散的數(shù)據(jù)，離散序列的自相關(guān)和互相關(guān)的計算公式如下:Rxy(O)Rx(0)Ry(O)1 n 1-x(i)y(i)n i 01 n 11 x(i)2n i 01 n 1 2 y(i) n i 0(2-23)式中：n表示采樣個數(shù)，i表示第i個采樣點，x（i）、y（i）分別表示

24、兩路信號的第i個點的轉(zhuǎn)換得到數(shù)值。由式（2-23 ）分別求出兩路信號的自相關(guān)和互相關(guān)函數(shù)值之后，再由式（2-21 ）和式（2-22）即可得到兩路信號各自的幅值和它們之間的相位差。但是，需要指出的是，由數(shù)字相關(guān)法求得的相位差，并不能區(qū)分是超前還是滯后， 這就需要采用其他方法來確定相位差符號的正負(fù)號。根據(jù)前面測量原理中的介紹，由式（2-15 ）和式（2-16 ）可知，壓電陶瓷換能器的電導(dǎo)值僅取決于兩路信號相位差的余弦值，而電納的值是在電導(dǎo)值取得最大的時候發(fā)生變號。此，可以先求得電導(dǎo)的值，再通過循環(huán)找其最大值， 并從使電導(dǎo)取得最大值時的相位差開始,2-9表示把相位差變號，得到新的相位差序列，再由新的

25、相位差序列求電納的值即可。圖其中（a）表示的是導(dǎo)納圓圖,的為采用數(shù)字相關(guān)法對一號壓電換能器測量數(shù)據(jù)的處理結(jié)果,（b）表示的是電導(dǎo)和電納值隨測試頻率的變化曲線。22046屯呈mS5 2 1 壬 lb22 5232362i(b)圖2-9數(shù)字相關(guān)法處理結(jié)果2、相關(guān)函數(shù)法的特點及誤差分析相關(guān)函數(shù)法測量信號的幅值通過上面對相關(guān)函數(shù)法測量原理的理論推導(dǎo)過程可以看出,和相位差與信號的頻率無關(guān)。也即是說相關(guān)函數(shù)法不受頻率的影響，可以用來測量未知頻率因此，它只能用的信號的相位差。同時，相關(guān)函數(shù)法測量原理的推導(dǎo)都是基于正弦函數(shù)的, 于測量正弦或余弦信號，并不能測量一般的周期信號。由于噪聲干擾信號和原信號并不相關(guān)，

26、所以相關(guān)函數(shù)法能夠有效的抑制噪聲干擾。但是,如果在系統(tǒng)中存在相關(guān)性較強的干擾信號，并且信噪比又比較低的情況下，相關(guān)函數(shù)法測量誤差就會比較大。由相關(guān)函數(shù)法離散序列的最終計算公式可以看出,其計算結(jié)果與采樣的點數(shù)有關(guān)，也即是說測量誤差的大小與采樣點數(shù)是相關(guān)的,采樣點數(shù)越大，計算結(jié)果越接近真實值，測量誤差也就越小。綜合以上對相關(guān)函數(shù)法的特點的分析，可知相關(guān)函數(shù)法對于采樣轉(zhuǎn)換信號中的直流偏移和噪聲等干擾具有很強的抑制能力，它的誤差主要是因為采用有限長度的樣本代替了高斯白噪聲和均勻分布的 A/D量化誤差，使得被檢正弦信號與噪聲信號并非完全不相關(guān)。所以,相關(guān)函數(shù)法的測量誤差與 A/D轉(zhuǎn)換的位數(shù)、信號的信噪比

27、和采集點數(shù)有關(guān)。242快速離散傅里葉變換法現(xiàn)代信號分析采用數(shù)字化方式實現(xiàn)，其核心是離散傅立葉變換，它完成了從時域到頻域的轉(zhuǎn)換，不僅可以實現(xiàn)線性譜分析，而且還是均方譜分析的關(guān)鍵。離散傅立葉變換(DFT)實現(xiàn)了信號首次在頻域表示的離散化，使得頻域也能夠用計算機進行處理，但由于用于實際時計算量太大而使應(yīng)用受到限制。直到1965年由Cooly和Tukey建立了一種快速傅立葉變換一一FFT時，DFT的應(yīng)用才成為現(xiàn)實1、FFT獲取正弦波幅值和相位的原理設(shè)采集正弦信號得到的離散序列為x(n) , n=1,2 ,KN則該序列的離散傅里葉變換為:N 1X(k) DFTx( n)x(n)en 0.2 k j nN

28、 12 kx(n)cos(nn 0N2 kjsin(n) ReX(k) ImX(k) Nk 0,1,2,N-1(2-24)則其初始相位為:arctan冊),k其中：fs是信號的采樣頻率,N是采樣長度。在對時域離散序列進行傅立葉變換之后，可以得到其離散的幅度譜和相位譜，在幅度譜 和相位譜中找到對應(yīng)時域波形的頻率的譜線就可以得到時域的正弦波形的幅值和相位信息。圖2-10所示的是采用快速離散傅里葉變換法對采集到的數(shù)據(jù)處理的結(jié)果?！?b044:J- ' * 1 LB JLL LBC匕J亠3,|lu1V!'1H H H HH. HHH H H H H 4 II - -. - - . -

29、A. - . J. 1+. iJJ102463ME瘞m1030 -4p22.52323.5隕丕山He24(a)(b)圖2-10快速離散傅里葉變換法處理結(jié)果2、FFT的特點及誤差分析通過傅里葉變換可以只提取基波參數(shù)，因此諧波的存在并不影響基波成分,所以諧波的存在對應(yīng)用這種方法測量相位差幾乎沒有影響;對于噪聲干擾，只有當(dāng)高斯白噪聲接近基波的頻率分量時才會影響到基波的相位，所以應(yīng)用FFT法測量相位差也能有效地抑制高斯白FFT對其進行譜分析時，必須截噪聲干擾。但是，實際上信號是連續(xù)的無限長的序列，用短形成有限長序列，再進行周期延拓，這樣就不可避免的造成信號頻譜的泄漏，由此便產(chǎn)生 了相位差測量誤差。誤差

30、現(xiàn)象主要是：混疊現(xiàn)象、柵欄效應(yīng)和截斷效應(yīng)。要想減小相位差測 量誤差，就必須提高譜分辨率。實際中可通過提高采樣頻率或者增加采樣數(shù)據(jù)長度來提高譜 分辨率，進而達到減小相位差測量誤差的目的。243正弦曲線參數(shù)擬合法設(shè)被測的正弦信號為:f(t)A0sin(2 ft ) D(2-25)其中：f表示信號頻率， Ao表示被測信號幅值,表示被測信號的初始相位角，D表故只需對測得的數(shù)據(jù)進行三參數(shù)的示被測信號的直流分量。 由于被測信號的頻率為已知的, 正弦曲線擬合，即可得到被測信號的幅值和相位信息。為此，進一步將上式展開可得:(2-26)f Asin(2 ft) Bcos(2 ft) D其中：A Ad cos ,

31、B A0sin從而將被測信號的幅值和初始相角轉(zhuǎn)化為對參數(shù)A、B的求取。其基本思想就是尋找合適的 A、B和D的值，使得其測量殘差的平方和取得最小。設(shè)1每個頻率下測量的時間序列為ti ,(i1,2,3 n) , n為測量數(shù)據(jù)的個數(shù)，fs為采樣頻fs率，fi為每個點的測量值。則測量殘差的表達式為:nf (t) Asi n(2 fti) Bcos(2 fti) D2(2-27)nsin(2 fti)Asin(2 ftBcos(2 fti)i 0ncos(2 fti)Asin(2 fti) Bcos(2 fti)i 0nAsin(2 fti) Bcos(2 ftj Di 0DDf(i)對于式(2-29

32、),構(gòu)造如下三個矩陣:sin(2 fti)f(i)i 0ncos(2 fti)f(i) (2-29)sin 2 fticos2 fti1sin2 ft2.sin2 ftncos2 ft2.cos2ftn1 1 1,X,Ff1(ti) f1(t2)f1(tn)i 1要使得上式取得最小值，可對其參數(shù)求偏導(dǎo)，并令其為零。即:nf (i) Asin(2 fti) Bcos(2 fti)D20A A i 1n一 一 f(i) Asin(2 fti) Bcos(2 fti)D20(2-28)B B i 1n f(i) Asin(2 fti) Bcos(2 fti)D20進一步化簡得到:i 1式(2-29

33、)可寫成如下矩陣形式:T X tF上式中X的解為:X ( T )1( TF)(2-30)則被測信號幅值的計算公式為A0Va2 B2(2-31)初始相角的計算公式為:Barcta n()A 0(2-32)A 0Aarctan(B)至此，得到了被測信號的幅值和相位信息。采用同樣的方法對第二路信號的采樣數(shù)據(jù)進行處理，即可得到第二路信號的幅值和相位信息,從而求出兩路信號的幅值比和相位差，進步便可得到每個頻率下?lián)Q能器的電導(dǎo)和電納值。圖2-11為參數(shù)擬合法的處理結(jié)果。0.6054i-/ :V 0.4 d護0.30.2密:0.1卜00.100204io<6 oe坦 3/rnS圖2-11正弦曲線參數(shù)擬合

34、法處理結(jié)果2.5導(dǎo)納圓的帶約束最小二乘曲線擬合通過以上章節(jié)的介紹， 我們已經(jīng)得到了各個測試頻率下壓電換能器的電導(dǎo)和電納值, 制出了導(dǎo)納圓圖，但這還是不夠的。由壓電換能器等效電路的各個參數(shù)的計算公式可以看出,我們還需要得到導(dǎo)納圓的圓心和半徑的值。為此，就需要對所得到的離散點進行圓曲線擬合。擬合圓的方法有很多種，常用的有平均值法、加權(quán)平均法和最小二乘法。平均值法的思想是分別計算各個離散點的橫、縱坐標(biāo)的平均值，作為圓心的橫、縱坐標(biāo),將圓心到各個離散點的距離的平均值作為半徑。這種方法計算簡單，適用于離散點分布較均 勻的情況，但對于分布不均的情況， 所計算的圓心位置會偏向離散點分布較密集的一側(cè),徑的計算值也會偏小，誤差較大。加權(quán)平均法是對平均值法的改進，它在計算圓心坐標(biāo)時加 入一個與兩相鄰點間弧長相關(guān)的系數(shù)，降低了離散點分布不均的影響，減小了誤差。但是,，所以其誤本文采取一種帶由于兩相鄰點間的弧長是無法精確得到的（實際