压电陶瓷测量原理

1、压电陶瓷及其测量原理近年来，压电陶瓷的研究发展迅速，取得一系列重大成果，应用范围不断扩大，已深入到国民经济和尖端技术的各个方面中，成为不可或缺的现代化工业材料之一。由于压电材料的各向异性，每一项性能参数在不同的方向所表现出的数值不同，这就使得压电陶瓷材料的性能参数比一般各向同性的介质材料多得多。同时，压电陶瓷的众多的性能参数也是它广泛应用的重要基础。（一）压电陶瓷的主要性能及参数（1）压电效应与压电陶瓷在没有对称中心的晶体上施加压力、张力或切向力时，则发生与应力成比例的介质极化，同时在晶体两端将出现正负电荷，这一现象称为正压电效应；反之，在晶体上施加电场时，则将产生与电场强度成比例的变形或机械

2、应力，这一现象称为逆压电效应。这两种正、逆压电效应统称为压电效应。晶体是否出现压电效应由构成晶体的原子和离子的排列方式，即晶体的对称性所决定。在声波测井仪器中，发射探头利用的是正压电效应，接收探头利用的是逆压电效应。（2）压电陶瓷的主要参数1、介质损耗介质损耗是包括压电陶瓷在内的任何电介质的重要品质指标之一。在交变电场下，电介质所积蓄的电荷有两种分量：一种是有功部分（同相），由电导过程所引起；另一种为无功部分（异相），由介质弛豫过程所引起。介质损耗是异相分量与同相分量的比值，如图1所示，IC为同相分量，*为异相分量，IC与总电流I的夹角为8,其正切值为tan8-Lr-1ICRC其中s为父变电场

3、的角频率，R为损耗电阻，C为介质电容。图1交流电路中电压-电流矢量图（有损耗时）2、机械品质因数机械品质因数是描述压电陶瓷在机械振动时，材料内部能量消耗程度的一个参数，它也是衡量压电陶瓷材料性能的一个重要参数。机械品质因数越大，能量的损耗越小。产生能量损耗的原因在于材料的内部摩擦。机械品质因数0的定义为：mQ谐振时振子储存的机械能°nm=谐振时振子每周所损失的机械能X2机械品质因数可根据等效电路计算而得Q_m3L3Rc=s111式中R1为等效电阻（Q）,°s为串联谐振角频率（Hz）,Ci为振子谐振时的等效电容（F）,L1为振子谐振时的等效电感。Q与其它参数之间的关系将在后续

4、详细推导。1m不同的压电器件对压电陶瓷材料的Qm值的要求不同，在大多数的场合下（包括声波m测井的压电陶瓷探头），压电陶瓷器件要求压电陶瓷的Q值要高。m3、压电常数压电陶瓷具有压电性，即在其外部施加应力时能产生额外的电荷。其产生的电荷与施加的应力成比例，对于压力和张力来说，其符号是相反的，电位移D（单位面积的电荷）和应力b的关系表达式为：D=drA式中Q为产生的电荷（C）,A为电极的面积（m2）,d为压电应变常数（C/N）。在逆压电效应中，施加电场E时将成比例地产生应变S,所产生的应变S是膨胀还是收缩，取决于样品的极化方向。S=dEDS两式中的压电应变常数d在数值上是相同的，即d=bE另一个常用

5、的压电常数是压电电压常数g,它表示应力与所产生的电场的关系，或应变d与所引起的电位移的关系。常数g与d之间有如下关系：g=-£式中£为介电系数。在声波测井仪器中，压电换能器希望具有较高的压电应变常数和压电电压常数，以便能发射较大能量的声波并且具有较高的接受灵敏度。4、机电耦合系数当用机械能加压或者充电的方法把能量加到压电材料上时，由于压电效应和逆压电效应，机械能（或电能）中的一部分要转换成电能（或机械能）。这种转换的强弱用机电耦合系数k来表示，它是一个量纲为一的量。机电耦合系数是综合反映压电材料性能的参数，它表示压电材料的机械能和电能的耦合效应。机电耦合系数的定义为：殳2=

6、电能转变为机械能或者上2二机械能转变为电能输入电能输入机械能机电耦合系数不但与材料参数有关，还与具体压电材料的工作方式有关。对于压电陶瓷来说，它的大小还与极化程度相关。它只是反映机、电两类能量通过压电效应耦合的强弱，并不代表两类能量之间的转换效率。压电材料的耦合系数在不同的场合有不同的要求，当制作换能器时，希望机电耦合系数越大越好。（二）压电换能器的等效电路压电换能器的等效电路表示法，是利用电学网络术语表示压电陶瓷的机械振动特性，即把某些力学量模拟为电学量的方法。把压电换能器用等效电路来表示，有很多优点：其一，可以把力学上复杂的振动问题有效地进行简化；其二，为了得到换能器的各个参数，从而定量地

7、分析或筛选换能器；其三，实际应用的需要，因为在实际的应用当中，压电换能器也是接入到具体的电子线路中的，得到压电换能器的等效电路能够更好地对其外围电路进行匹配设计。由此可见，得到压电换能器的等效电路是十分必要的。2.3压电换能器的谐振特性将压电换能器按照图2-2所示线路连接。当改变信号频率时，可以发现，通过压电陶瓷换能器的电流也随着发生变化，其变化规律如图2-3（a）所示。从图2-3（a）可以看出,当信号为某一频率fm时，通过压电陶瓷换能器的电流出现最大值1max；而当信号变到另频率/时，传输电流出现最小值1。由流经它的电流随频率的变化可以看出，压电陶瓷nmin换能器的阻抗是随频率的变化而变化的

8、，其变化规律同电流相反，如图2-3(b)所示。图2-2压电陶瓷换能器谐振特性接线示意图fnf®tb)图2-3压电陶瓷换能器电流、阻抗同频率的关系曲线(a)电流-频率关系曲线(b)阻抗-频率关系曲线从图中可以看出，当信号频率为f时，通过压电陶瓷换能器的电流最大，即其等效阻m抗最小，导纳最大；当信号频率为f时，通过压电陶瓷换能器的电流最小，即其等效阻抗n最大，导纳最小。因此把/称为最大导纳频率或最小阻抗频率；而把f称为最小导纳频率mn或最大阻抗频率。而当信号频率继续增大时，还会出现一系列的电流的极大值和极小值，如图2-4所示。图2-4压电陶瓷换能器电流随频率变化示意图（多谐振模式）2.2

9、.4压电换能器的等效电路根据交流电路相关知识，对于图2-5所示好的LC电路来说，其阻抗Z也随着频率的变化而变化。在图2-2所示的线路中，用LC电路代替压电陶瓷换能器，可以发现，在压电陶瓷换能器的谐振频率处，只要选择合适的L、C1、R和C0，通过LC电路的电流和LC电路的阻抗的绝对值随频率的变化曲线，分别同图2-1中的（b）和（c）的关系曲线非常相似。也就是说，在串联谐振频率附近，压电陶瓷换能器的阻抗特性和谐振特性同LC电路的阻抗特性和频率特性非常相似。因此，利用机电类比的方法，可以用一个LC电路来表示压电陶瓷换能器的参数和特性，这个LC电路即为压电陶瓷换能器的等效电路。LClRl1=16II图

10、2-5LC电路对压电陶瓷换能器来说，在任何串联谐振频率附近，其电行为可以用图2-3所示的LC电路来表示。在压电陶瓷换能器的串联谐振频率附近，如果值存在一种振动模式，即没有其它寄生响应，则在串联谐振频率附近很窄的频率范围内，可以认为压电陶瓷换能器的等效参数R、C1、R和C0与频率无关。在实际中通过选择合适的尺寸进行加工处理，是可以将所需要的振动模式同其他模式充分隔离开来的。另外，考虑到在实际中，在通电之后，压电陶瓷换能器必然会存在能量的损耗，这一能量损耗可用一个并联电阻R0来等效。所以其最终等效电路图如图2-6所示。LiClRiIIQG.1|图2-6压电陶瓷换能器等效电路图图中串联支路中的L1称

11、为压电陶瓷换能器的动态电感，C1称为动态电容，R1称为动态电阻。这三个参数用来表征压电陶瓷换能器在工作（加电源激励产生振动）的情况下，振动部分所受到的力阻抗和介质对振动的反作用的强弱。并联电容Co又称静态电容，表征压电陶瓷换能器在未加激励的情况下等效为一个纯电容，它的值的大小与换能器的形状有关。并联电阻Ro又称静态电阻，表征换能器的电损耗的大小。2.2.5压电换能器的导纳特性根据已得到的压电换能器的等效电路图，来进一步分析其导纳特性。为了简化推导，先假定压电陶瓷换能器没有电损耗，即Ro=0，此时其等效电路即为一个LC电路，如图2-5所示。则Y二Yo+Y1（2-1）式中：Y为换能器的总的导纳值，

12、Yo二j叮巴为并联支路的导纳值，Y1=G1=jB为串联支路的导纳值。先对串联支路进行分析。_1_Z"11R+j®L+11j6C1Rj©L-厶)11®C=产得到:R2+LF)11®C1G1=R11R2+(®L)2113C1©L-丄)B_1XB1R2+(®L)2113C1(2-2)1若令R2+(L1C）2二fxR12二（七云C）。由式（2-2）可得:11x=RGiB21GxR2)2=1(xR2)=G2+1R111G-1iR1所以GB2+G21二011R1边同时加上（2R）2,可得111(G)2+B2=()2(2_3)i

13、2Ri2R(2-3)11若以电导为横坐标，电纳为纵坐标，则式2-3）表示一个以112r,0)为圆心，2r11为半径的圆，也即是我们所说的导纳圆。如图2-7中虚线所示图2-7导纳圆图对于串联支路进行分析，根据串联谐振频率的定义，令B1=0,则由式（2-3）可得到Gi=0或=1R。由于实际的压电陶瓷换能器的动态电阻R0不可能为零,根据式（2-2）中G1的表达式可以知道，只有=1R满足串联谐振的条件。即：£示C二0，所以11可以得到串联支路的谐振频率（又称机械共振频率）：°s二7（2-4）接着考虑加入静态电容后的情况。由式（2-1）可知，考虑静态电容后换能器的导纳相当于在串联支路

14、的电纳（虚部）加上Y0。鉴于一般情况下，压电陶瓷换能器的机械品质因数都较大，也即在串联谐振频率®s附近，Y0二jC0的值随频率的变化很小，可以近似认s00为是一个常数。因此，只需将串联支路所得到的导纳圆的纵坐标向上平移一个常数，而横坐标保持不变即可得到加入静态电容后换能器的导纳关系图，如图2-7中点划线所示。若再考虑到换能器的静态电阻并不为零，则实际中的导纳圆不可能与纵轴相切，而是向横轴的正向平移一定的量（平移距离的大小取决于静态电阻的阻值），如图2-7中实线圆所示对导纳圆图进行简要的分析可知：当f<f即<叫时，电纳值大于零，当f>f即sss°>&#

15、176;s时，电纳值小于零。所以，随着频率的增加，导纳圆是沿顺时针方向变化的。另外,在串联谐振频率的附近，还存在着两个频率点使得换能器总的电纳为零，此时电源信号经过换能器之后只有幅值的改变，而没有相位的变化，也即电压和电流信号同相位。这两个频率中，值较小的那个频率f称为谐振频率，较大f的称为反谐振频率。另外还存在使得换ra能器的导纳值取得最大的频率f，导纳值最小的频率f。连接原点和串联谐振频率点，mn与导纳圆的交点处的频率fp称为并联谐振频率。另外，需要特别指出的是，上述讨论是在一个振动模态谐振频率®s附近较小的频率变化范围内进行的，并且只有在导纳圆的直径远大于这个频率范围内

16、6;Co的变化时才是正确的，否则换能器的导纳曲线将变得十分复杂，具有蔓叶曲线的特征。根据以上导纳圆图的推导过程，下面介绍一下压电陶瓷换能器等效电路中各个参数和导纳圆图的关系，并给出各自的计算公式。在换能器的导纳圆图中作平行于纵轴的直径，交导纳圆于两点，分别记作。在f点处，串联支路的动态电导和电纳值相等，即G二B。由式（2-2）可得:R2(°L_-4)111°C11-(°L-丄)11°CR2+(°L-)2111°C11(2-5)在f2点处，串联支路的动态电导和电纳值相等，但符号相反，即G1=-B1。由式（2-2）可得:R11R2+(&#

17、176;L)121°C21°L-R2+(°12°CL-°C)2(2-6)结合式（2-5）和式（2-6），可得：(2-7)R1°-°21再由式（2-4）可得：(2-8)C=-1-1w2Ls1机械品质因数：QmCRs11wL=s1R1(2-9)结合式（2-7）和（2-8）可得：Qm=w:smww21式（2-5）和式（2-6）消去R得到：而"12二LC211s所以：Q=%：mf2-f2-10)动态电阻的值可以通过导纳圆的直径求得：R二1.D02-11)静态电容C0的值也可由导纳圆偏离横轴的距离来确定：CoW(212)s式

18、中y0为圆心的纵坐标。静态电阻R0的值可由导纳圆偏离纵轴的距离（或圆心的横坐标）来确定:R=-0x02R12-13)式中x0为圆心的横坐标。至此，我们已得到压电陶瓷换能器等效电路中所有参数的计算公式。2.3测量原理在上一节中，得到的压电陶瓷换能器等效电路参数的计算公式都是基于导纳圆的，也即是基于各个频率下的电导和电纳值的，因此我们需要得到每个频率点的导纳值。为此采用图2-8所示的测量原理图进行测量。压峽瓷-:iI山Ur如Ml图2-8压电陶瓷换能器测量原理示意图图2-8中，AC为频率可控的交流信号源，R表示源内阻，Rm称为精密电阻，UA为加在压电陶瓷换能器山的电压信号UB为经过换能器之后的电压信

19、号。根据前面章节所介绍的压电陶瓷的导纳特性可以知道，在经过换能器之后的电压信号UB相对于UA会有一个幅度和相位的变化。不失一般性，在这里设定：Ua二匕e网申，U二Uej(叫+0)(2-14)AmBBm其中：U，U分别表示两路信号的幅值，®为信号的角频率，0为信号的初始AmBm相位，0为两路信号的相位差。按照习惯表达，先求压电陶瓷换能器的阻抗，再取倒数得到导纳。7UU-UZmBBrm将式(2-14)代入得R(Ame-j01)mUBmUUejra+0ZAmBmUej(哪+e)/Bmrm二R(AmCOSsin0R+jX=mUUBmBm2-15)对应得到：RRm(UA"cos01)

20、，X=RmUA-Sin9BmBm再由导纳和阻抗的关系可得YIR+jXG+jB1G二RBx即：r+x，R+X(2-16)2222由以上推导可以看出，换能器的导纳和阻抗值仅与加在其两端的电压信号的幅值比和相位差有关，因此只需要得到两路信号的幅值和相位信息即可得到换能器等效电路的各个参数。而实际中，只需要对两路信号进行采样，再通过对采样所得数据进行处理便可得到幅值和相位信息。2.4正弦信号的测量方法根据上一节介绍的测量原理可知，要得到压电换能器在测试频率下的电导和电纳值，就需要测得其两端正弦信号的幅值比和相位差。但是实际中，硬件电路实现的仅是对两路信号的A/D转换采集，也即是得到的是两路正弦信号的一

21、系列的离散的点。在这一节中，将介绍从这些采集到的离散的点计算其幅值和相位的方法。2.4.1数字相关法随着微处理器和大规模集成电路的迅速发展，在测试系统中，越来越多的传统的测量方法被数字化测量方法所取代。近年来，由于相关函数法具有提高测试精度，减少或简化硬件设计，能够充分利用测试系统中的数据采集系统和微型计算机，提高测试系统的可靠性和可维护性的诸多优点，使得相关技术原理在相位差的测量及数字信号处理中得到了广泛应用，并展现出良好的应用前景1、相关函数法原理相关函数法利用两同频正弦信号的延时为零时的互相关函数值与其相位差的余弦值成正比的原理获得相位差。设两路被测信号为：x(t)=Asin(2兀)+N

22、(t)，y(t)=Bsin(2劝+*)+N(t)(217)xy其中：A、B分别表示两路信号的幅值，f表示信号的频率，N(t)、N(t)分别表示xy两路信号的干扰噪声信号，0表示两路信号的相位差。显然，信号x(t)和y(t)是相关的，则两路信号的互相关函数为：R(T)=丄Tx(t)y(t+t)dt218)xyT0R(t)=丄JTAsin(2矽)+N(t)-Bsin(2砒(t+t)+申)+N(t+t)dtxyT0xy式中T为信号的周期，即T二4=0时，有R(0)=JTAsin(2f)+N(t)-Bsin(2f+申)+N(t)dtxyT0xy由于噪声信号之间不相关，噪声和信号之间也不相关，将上式进一

23、步展开得R=丄JTABsin(2妙)sin(2妙+申)dtxyT0=JTABsin(2矽)sin(2矽)cos申+cos(2矽)sin申dtT011=JTABsin(2矽)sin(2矽)cos申dt=ABcos申T02所以，可以得到相位差的计算公式：2R(0)(2-19)申=arccos(xy)AB而信号幅值的大小可由信号的自相关函数求得R(T)=xx(t)x(t+T)dt=TJ_2tAsin(2劝)Asin(2吋(t+t)dt(2-20)当T=0时，有1.zA2R(0)=-J2A2sin2(2；ft)dt=一xT_z2_2所以可得信号幅值的计算公式：A=x(2-21)B=QR(0)y将上式代

24、入式(2-19)，可得相位差计算公式的另一种表达式R(0)(2-22)申=arccos(JR(0)R(0)xy而在实际中，是没有完整精确的信号的表达式的，有的是对信号的模数转换所得到的离散的数据，离散序列的自相关和互相关的计算公式如下：R(0)=1艺x(i)y(i)xyni=0R(0)=1艺x(i)2(2-23)xni=0R(0)=-艺y(i)2yni=0式中：n表示采样个数，i表示第i个采样点，x(i)、y(i)分别表示两路信号的第i个点的转换得到数值。由式(2-23)分别求出两路信号的自相关和互相关函数值之后，再由式(2-21)和式(2-22)即可得到两路信号各自的幅值和它们之间的相位差。

25、但是，需要指出的是，由数字相关法求得的相位差，并不能区分是超前还是滞后，这就需要采用其他方法来确定相位差符号的正负号。根据前面测量原理中的介绍，由式(2-15)和式(2-16)可知，压电陶瓷换能器的电导值仅取决于两路信号相位差的余弦值，而电纳的值是在电导值取得最大的时候发生变号。由此，可以先求得电导的值，再通过循环找其最大值，并从使电导取得最大值时的相位差开始，把相位差变号，得到新的相位差序列，再由新的相位差序列求电纳的值即可。图2-9表示的为采用数字相关法对一号压电换能器测量数据的处理结果，其中(a)表示的是导纳圆图，(b)表示的是电导和电纳值随测试频率的变化曲线。(a)(b)图2-9数字相

26、关法处理结果2、相关函数法的特点及误差分析通过上面对相关函数法测量原理的理论推导过程可以看出，相关函数法测量信号的幅值和相位差与信号的频率无关。也即是说相关函数法不受频率的影响，可以用来测量未知频率的信号的相位差。同时，相关函数法测量原理的推导都是基于正弦函数的，因此，它只能用于测量正弦或余弦信号，并不能测量一般的周期信号。由于噪声干扰信号和原信号并不相关，所以相关函数法能够有效的抑制噪声干扰。但是，如果在系统中存在相关性较强的干扰信号，并且信噪比又比较低的情况下，相关函数法测量误差就会比较大。由相关函数法离散序列的最终计算公式可以看出，其计算结果与采样的点数有关，也即是说测量误差的大小与采样

27、点数是相关的，采样点数越大，计算结果越接近真实值，测量误差也就越小。综合以上对相关函数法的特点的分析，可知相关函数法对于采样转换信号中的直流偏移和噪声等干扰具有很强的抑制能力，它的误差主要是因为采用有限长度的样本代替了高斯白噪声和均匀分布的A/D量化误差，使得被检正弦信号与噪声信号并非完全不相关。所以，相关函数法的测量误差与A/D转换的位数、信号的信噪比和采集点数有关。2.4.2快速离散傅里叶变换法现代信号分析采用数字化方式实现，其核心是离散傅立叶变换，它完成了从时域到频域的转换，不仅可以实现线性谱分析，而且还是均方谱分析的关键。离散傅立叶变换(DFT)实现了信号首次在频域表示的离散化，使得频

28、域也能够用计算机进行处理，但由于用于实际时计算量太大而使应用受到限制。直到1965年由Cooly和Tukey建立了一种快速傅立叶变换FFT时，DFT的应用才成为现实1、FFT获取正弦波幅值和相位的原理设采集正弦信号得到的离散序列为x(n),n=1,2,駅。则该序列的离散傅里叶变换为:N-1.2kX(k)二DFTx(n)二£x(n)e-nnn=0=£1x(n)cos(2kn)-jsin(2n)=ReX(k)+ImX(k)NNn=0k=0,1,2,N-1(2-24)则其初始相位为:9=arctan(ReX(k),k=f+1ImX(k)fs其中：fs是信号的采样频率，N是采样长度

29、。在对时域离散序列进行傅立叶变换之后，可以得到其离散的幅度谱和相位谱，在幅度谱和相位谱中找到对应时域波形的频率的谱线就可以得到时域的正弦波形的幅值和相位信息图2-10所示的是采用快速离散傅里叶变换法对采集到的数据处理的结果。图2-10快速离散傅里叶变换法处理结果2、FFT的特点及误差分析通过傅里叶变换可以只提取基波参数，因此谐波的存在并不影响基波成分，所以谐波的存在对应用这种方法测量相位差几乎没有影响；对于噪声干扰，只有当高斯白噪声接近基波的频率分量时才会影响到基波的相位，所以应用FFT法测量相位差也能有效地抑制高斯白噪声干扰。但是，实际上信号是连续的无限长的序列，用FFT对其进行谱分析时，必

30、须截短形成有限长序列，再进行周期延拓，这样就不可避免的造成信号频谱的泄漏，由此便产生了相位差测量误差。误差现象主要是：混叠现象、栅栏效应和截断效应。要想减小相位差测量误差，就必须提高谱分辨率。实际中可通过提高采样频率或者增加采样数据长度来提高谱分辨率，进而达到减小相位差测量误差的目的。2.4.3正弦曲线参数拟合法设被测的正弦信号为：f(t)二Asin(2硏+Q)+D(2-25)0其中：f表示信号频率，A°表示被测信号幅值，©表示被测信号的初始相位角，D表示被测信号的直流分量。由于被测信号的频率为已知的，故只需对测得的数据进行三参数的正弦曲线拟合，即可得到被测信号的幅值和相位

31、信息。为此，进一步将上式展开可得：f=Asin(2硏)+Bcos(2妙)+D(2-26)苴中.A=AcosQ,B=AsinQ00从而将被测信号的幅值和初始相角转化为对参数A、B的求取。其基本思想就是寻找合适的A、B和D的值，使得其测量残差的平方和取得最小。设1每个频率下测量的时间序列为ti二，（i二1,2,3.）,n为测量数据的个数，f为采样频ifss率，f为每个点的测量值。则测量残差的表达式为：iE=2f（t）-Asin（2胡）-Bcos（2劝）-D2（227）iii=1要使得上式取得最小值，可对其参数求偏导，并令其为零。即真=斗工f(i)-Asin(2矽)-Bcos(2矽)-0cAcAI1

32、1Ji=1<=舟工f(i)-Asin(2Kft)-Bcos(2)-DUL0cBcBIiii=1H=CD目f(i)-Asin(2叫-Bcos(2叫-D2=0IDDIi=1(228)进一步化简得到：工sin(2妙)Asin(2胡)+Bcos(2胡)+D=工sin(2胡)f(i)iiiii=0i=0工cos(2妙)Asin(2胡)+Bcos(2胡)+D=工cos(2胡)f(i)iiii<i=0工Asin(2矽)+Bcos(2妙)+D=工f(i)iii=0i=0对于式（229），构造如下三个矩阵：i=0rf1(t)sin2兀ft.sin2矽A"12nf1(t)cos2胡.cos2

33、妙n,X=B,F=2111D_f1(t)_sin2矽cos2妙n式（229）可写成如下矩阵形式：屮t屮X=屮tF上式中X的解为：X=（屮T屮）-1（屮TF）（230）则被测信号幅值的计算公式为:IA=；'A2+B2（231）0(229)初始相角的计算公式为：Barctan()A>0(2-32)ABarctan()+兀A<0A至此，得到了被测信号的幅值和相位信息。采用同样的方法对第二路信号的采样数据进行处理，即可得到第二路信号的幅值和相位信息，从而求出两路信号的幅值比和相位差，进一步便可得到每个频率下换能器的电导和电纳值。图2-11为参数拟合法的处理结果。-Q.1a11+H-

34、T.一u+1+.-丄:*r-c02040$032343.6审宰出HrSE_盖图2-11正弦曲线参数拟合法处理结果2.5导纳圆的带约束最小二乘曲线拟合通过以上章节的介绍，我们已经得到了各个测试频率下压电换能器的电导和电纳值，绘制出了导纳圆图，但这还是不够的。由压电换能器等效电路的各个参数的计算公式可以看出我们还需要得到导纳圆的圆心和半径的值。为此，就需要对所得到的离散点进行圆曲线拟合拟合圆的方法有很多种，常用的有平均值法、加权平均法和最小二乘法。平均值法的思想是分别计算各个离散点的横、纵坐标的平均值，作为圆心的横、纵坐标将圆心到各个离散点的距离的平均值作为半径。这种方法计算简单，适用于离散点分布较均匀的情况，但对于分布不均的情况，所计算的圆心位置会偏向离散点分布较密集的一侧，半径的计算值也会偏小，误差较大。加权平均法是对平均值法的改进，它在计算圆心坐标时加入一个与两相邻点间弧长相关的系数，降低了离散点分布不均的影响，减小了误差。但是由于两相邻点间的弧长是无法精确得到的（实际中采用两点间的距离来代替的），所以其误差仍较大。相比之下，最小二乘法具