压电陶瓷测量基本知识

压电陶瓷及其测量原理
近年来，压电陶瓷的研究发展迅速，取得一系列重大成果，应用范围不断扩大，已深入到国民经济和尖端技术的各个方面中，成为不可或缺的现代化工业材料之一。

由于压电材料的各向异性，每一项性能参数在不同的方向所表现出的数值不同，这就使得压电陶瓷材料的性能参数比一般各向同性的介质材料多得多。

同时，压电陶瓷的众多的性能参数也是它广泛应用的重要基础。

（一）压电陶瓷的主要性能及参数
（1）压电效应与压电陶瓷
在没有对称中心的晶体上施加压力、张力或切向力时，则发生与应力成比例的介质极化，同时在晶体两端将出现正负电荷，这一现象称为正压电效应；反之，在晶体上施加电场时，则将产生与电场强度成比例的变形或机械应力，这一现象称为逆压电效应。

这两种正、逆压电效应统称为压电效应。

晶体是否出现压电效应由构成晶体的原子和离子的排列方式，即晶体的对称性所决定。

在声波测井仪器中，发射探头利用的是正压电效应，接收探头利用的是逆压电效应。

（2）压电陶瓷的主要参数
1、介质损耗
介质损耗是包括压电陶瓷在内的任何电介质的重要品质指标之一。

在交变电场下，电介质所积蓄的电荷有两种分量：一种是有功部分（同相），由电导过程所引起；另一种为无功部分（异相），由介质弛豫过程所引起。

介质损耗是异相分量与同相分量的比值，如图1 所示，C I为同相分量，R I为异相分量，C I与总电流I 的夹角为 ，其正切值为
CR I I C R ωδ1tan == 其中ω 为交变电场的角频率，R 为损耗电阻，C 为介质电容。

图 1 交流电路中电压-电流矢量图（有损耗时）
2、机械品质因数 机械品质因数是描述压电陶瓷在机械振动时，材料内部能量消耗程度的一个参数，它也是衡量压电陶瓷材料性能的一个重要参数。

机械品质因数越大，能量的损耗越小。

产生能量损耗的原因在于材料的内部摩擦。

机械品质因数m Q 的定义为：
π2的机械能
谐振时振子每周所损失能谐振时振子储存的机械⨯=m Q
机械品质因数可根据等效电路计算而得 11111R L C R Q s s m ωω==
式中1R 为等效电阻（Ω），s ω 为串联谐振角频率（Hz ），1C 为振子谐振时的等效电容（F ），1L 为振子谐振时的等效电感。

m Q 与其它参数之间的关系将在后续详细推导。

不同的压电器件对压电陶瓷材料的m Q 值的要求不同，在大多数的场合下（包括声波测井的压电陶瓷探头），压电陶瓷器件要求压电陶瓷的m Q 值要高。

3、压电常数
压电陶瓷具有压电性，即在其外部施加应力时能产生额外的电荷。

其产生的电荷与施加的应力成比例，对于压力和张力来说，其符号是相反的，电位移 D （单位面积的电荷）和
应力σ 的关系表达式为：dr A
Q D == 式中 Q 为产生的电荷（C ），A 为电极的面积（m ²），d 为压电应变常数（C/N ）。

在逆压电效应中，施加电场 E 时将成比例地产生应变 S ，所产生的应变 S 是膨胀还是收缩，取决于样品的极化方向。

S=dE
两式中的压电应变常数 d 在数值上是相同的，即E
S D
d ==σ 另一个常用的压电常数是压电电压常数 g ，它表示应力与所产生的电场的关系，或应变与所引起的电位移的关系。

常数 g 与 d 之间有如下关系： εd
g =
式中ε为介电系数。

在声波测井仪器中，压电换能器希望具有较高的压电应变常数和压电电压常数，以便能发射较大能量的声波并且具有较高的接受灵敏度。

4、机电耦合系数
当用机械能加压或者充电的方法把能量加到压电材料上时，由于压电效应和逆压电效应，机械能（或电能）中的一部分要转换成电能（或机械能）。

这种转换的强弱用机电耦合系数 k 来表示，它是一个量纲为一的量。

机电耦合系数是综合反映压电材料性能的参数，它表示压电材料的机械能和电能的耦合效应。

机电耦合系数的定义为：
k ²=输入电能电能转变为机械能或者k ²=输入机械能
机械能转变为电能 机电耦合系数不但与材料参数有关，还与具体压电材料的工作方式有关。

对于压电陶瓷来说，它的大小还与极化程度相关。

它只是反映机、电两类能量通过压电效应耦合的强弱，并不代表两类能量之间的转换效率。

压电材料的耦合系数在不同的场合有不同的要求，当制作换能器时，希望机电耦合系数越大越好。

（二）压电换能器的等效电路
压电换能器的等效电路表示法，是利用电学网络术语表示压电陶瓷的机械振动特性，即把某些力学量模拟为电学量的方法。

把压电换能器用等效电路来表示，有很多优点：其一，可以把力学上复杂的振动问题有效地进行简化；其二，为了得到换能器的各个参数，从而定量地分析或筛选换能器；其三，实际应用的需要，因为在实际的应用当中，压电换能器也是接入到具体的电子线路中的，得到压电换能器的等效电路能够更好地对其外围电路进行匹配设计。

由此可见，得到压电换能器的等效电路是十分必要的。

2.3 压电换能器的谐振特性
将压电换能器按照图 2-2 所示线路连接。

当改变信号频率时，可以发现，通过压电陶瓷换能器的电流也随着发生变化，其变化规律如图 2-3（a ）所示。

从图2-3（a ）可以看出，当信号为某一频率m f 时，通过压电陶瓷换能器的电流出现最大值m ax I ；而当信号变到另一频率n f 时，传输电流出现最小值m in I 。

由流经它的电流随频率的变化可以看出，压电陶瓷换能器的阻抗是随频率的变化而变化的，其变化规律同电流相反，如图 2-3（b ）所示。

图 2-2 压电陶瓷换能器谐振特性接线示意图
图 2-3 压电陶瓷换能器电流、阻抗同频率的关系曲线
(a)电流-频率关系曲线 (b)阻抗-频率关系曲线 从图中可以看出，当信号频率为m f 时，通过压电陶瓷换能器的电流最大，即其等效阻抗最小，导纳最大；当信号频率为n f 时，通过压电陶瓷换能器的电流最小，即其等效阻抗最大，导纳最小。

因此把m f 称为最大导纳频率或最小阻抗频率；而把n f 称为最小导纳频率或最大阻抗频率。

而当信号频率继续增大时，还会出现一系列的电流的极大值和极小值，如图 2-4 所示。

图 2-4 压电陶瓷换能器电流随频率变化示意图（多谐振模式）
2.2.4 压电换能器的等效电路
根据交流电路相关知识，对于图 2-5 所示好的 LC 电路来说，其阻抗 Z 也随着频率的变化而变化。

在图 2-2 所示的线路中，用 LC 电路代替压电陶瓷换能器，可以发现，在压电陶瓷换能器的谐振频率处，只要选择合适的1L 、1C 、1R 和 0C ，通过 LC 电路的电流
和LC 电路的阻抗的绝对值随频率的变化曲线，分别同图2-1中的（b）和（c）的关系曲线非常相似。

也就是说，在串联谐振频率附近，压电陶瓷换能器的阻抗特性和谐振特性同LC 电路的阻抗特性和频率特性非常相似。

因此，利用机电类比的方法，可以用一个LC 电路来表示压电陶瓷换能器的参数和特性，这个LC 电路即为压电陶瓷换能器的等效电路。

图2-5 LC 电路
对压电陶瓷换能器来说，在任何串联谐振频率附近，其电行为可以用图2-3所示的LC 电路来表示。

在压电陶瓷换能器的串联谐振频率附近，如果值存在一种振动模式，即没有其它寄生响应，则在串联谐振频率附近很窄的频率范围内，可以认为压电陶瓷换能器的等效参数1R、1C、1R和0C与频率无关。

在实际中通过选择合适的尺寸进行加工处理，是可以将所需要的振动模式同其他模式充分隔离开来的。

另外，考虑到在实际中，在通电之后，压电陶瓷换能器必然会存在能量的损耗，这一能量损耗可用一个并联电阻
R来等效。

所以其最终等效电路图如图2-6所示。

图2-6 压电陶瓷换能器等效电路图
图中串联支路中的1L称为压电陶瓷换能器的动态电感，1C称为动态电容，1R称为动态电阻。

这三个参数用来表征压电陶瓷换能器在工作（加电源激励产生振动）的情况下，振动
部分所受到的力阻抗和介质对振动的反作用的强弱。

并联电容 0C 又称静态电容，表征压电陶瓷换能器在未加激励的情况下等效为一个纯电容，它的值的大小与换能器的形状有关。

并联电阻 0R 又称静态电阻，表征换能器的电损耗的大小。

2.2.5 压电换能器的导纳特性
根据已得到的压电换能器的等效电路图，来进一步分析其导纳特性。

为了简化推导，先假定压电陶瓷换能器没有电损耗，即 0R =0，此时其等效电路即为一个 LC 电路，如图 2-5 所示。

则 10Y Y Y += (2-1)
式中：Y 为换能器的总的导纳值，000ωj jB Y ==为并联支路的导纳值，111jB G Y ==为串联支路的导纳值。

先对串联支路进行分析。

000ωj jB Y ==
)1()1(11111211111111C L R C L j R C j L j R Z ωωωωωω-+-
-=++==得到： 211111)1(2C L R R G ωω-+=,2111111)1()1(2C L R C L B ωωωω-+-
-= (2-2) 若令x C L R =-+2111)1(2ωω则)1(1
112C L R x ωω-=-。

由式（2-2）可得: 1
1G R x =,111111*********)()(R G G R x R G R x R G B +-=-=--= 所以，0111122=-+R G G B 两边同时加上21)21(R ，可得21
1211)21()21(2R B R G =+-(2-3) 若以电导为横坐标，电纳为纵坐标，则式（2-3）表示一个以（121R ,0）为圆心，1
21R
为半径的圆，也即是我们所说的导纳圆。

如图 2-7 中虚线所示
图 2-7 导纳圆图 对于串联支路进行分析，根据串联谐振频率的定义，令 1B =0，则由式（2-3）可得到 1G =0 或 111R G =。

由于实际的压电陶瓷换能器的动态电阻 0R 不可能为零，根据式
（2-2）中1G 的表达式可以知道，只有111R G =满足串联谐振的条件。

即：
0111=-C L ωω，所以可以得到串联支路的谐振频率（又称机械共振频率）：1
11C L s =ω (2-4) 接着考虑加入静态电容后的情况。

由式（2-1）可知，考虑静态电容后换能器的导纳相当于在串联支路的电纳（虚部）加上 0Y 。

鉴于一般情况下，压电陶瓷换能器的机械品质因数都较大，也即在串联谐振频率s ω 附近，00C j Y ω=的值随频率的变化很小，可以近似认为是一个常数。

因此，只需将串联支路所得到的导纳圆的纵坐标向上平移一个常数，而横坐标保持不变即可得到加入静态电容后换能器的导纳关系图，如图 2-7 中点划线所示。

若再考虑到换能器的静态电阻并不为零，则实际中的导纳圆不可能与纵轴相切，而是向横轴的正向平移一定的量（平移距离的大小取决于静态电阻的阻值），如图 2-7 中实线圆所示
对导纳圆图进行简要的分析可知：当s f f <即s ωω<时，电纳值大于零，当s f f >即s ωω>时，电纳值小于零。

所以，随着频率的增加，导纳圆是沿顺时针方向变化的。

另外，在串联谐振频率的附近，还存在着两个频率点使得换能器总的电纳为零，此时电源信号经过
换能器之后只有幅值的改变，而没有相位的变化，也即电压和电流信号同相位。

这两个频率中，值较小的那个频率r f 称为谐振频率，较大a f 的称为反谐振频率。

另外还存在使得换能器的导纳值取得最大的频率m f ，导纳值最小的频率n f 。

连接原点和串联谐振频率点，与导纳圆的交点处的频率p f 称为并联谐振频率。

另外，需要特别指出的是，上述讨论是在一个振动模态谐振频率s ω 附近较小的频率变化范围内进行的，并且只有在导纳圆的直径远大于这个频率范围内0C ω的变化时才是正确的，否则换能器的导纳曲线将变得十分复杂，具有蔓叶曲线的特征。

根据以上导纳圆图的推导过程，下面介绍一下压电陶瓷换能器等效电路中各个参数和导纳圆图的关系，并给出各自的计算公式。

在换能器的导纳圆图中作平行于纵轴的直径，交导纳圆于两点，分别记作1f 、2f 。

在1f 点处，串联支路的动态电导和电纳值相等，即11B G =。

由式（2-2）可得：
21111211111111111)1()1()1(2C L R C L C L R R ωωωωωω-+-
-=- (2-5) 在2f 点处，串联支路的动态电导和电纳值相等，但符号相反，即11B G -=。

由式（2-2）
可得:21212211212121211)1(1
)1(2C L R C L C L R R ωωωωωω-+-
=-+ (2-6) 结合式（2-5）和式（2-6），可得：
121
1ωω-=R L (2-7)
再由式（2-4）可得：
1121
L C s ω= (2-8)
机械品质因数：1
11111111C L R R L R C Q s s M ===ωω (2-9) 结合式（2-7）和（2-8）可得：121
2f f f Q s s
m -=-=ωωω 式（2-5）和式（2-6）消去1R 得到：1121211C L s ==ωωω 则21f f f s = 所以：1
221f f f f Q m -= （2-10） 动态电阻的值可以通过导纳圆的直径求得：
D R 10= （2-11）
静态电容 0C 的值也可由导纳圆偏离横轴的距离来确定：s y C ω0
0=
（2-12） 式中0y 为圆心的纵坐标。

静态电阻0R 的值可由导纳圆偏离纵轴的距离（或圆心的横坐标）来确定： 100211
R x R -= （2-13） 式中0x 为圆心的横坐标。

至此，我们已得到压电陶瓷换能器等效电路中所有参数的计算公式。

2.3 测量原理
在上一节中，得到的压电陶瓷换能器等效电路参数的计算公式都是基于导纳圆的，也即是基于各个频率下的电导和电纳值的，因此我们需要得到每个频率点的导纳值。

为此采用图 2-8 所示的测量原理图进行测量。

图 2-8 压电陶瓷换能器测量原理示意图
图 2-8 中，AC 为频率可控的交流信号源，R 表示源内阻，m R 称为精密电阻，A U 为加在压电陶瓷换能器山的电压信号B U 为经过换能器之后的电压信号。

根据前面章节所介绍的压电陶瓷的导纳特性可以知道，在经过换能器之后的电压信号B U 相对于A U 会有一个幅度和相位的变化。

不失一般性，在这里设定：
ωϕj Am A e U U =，)(θωϕ+=j Bm
B e U U （2-14） 其中： Am U ，Bm U 分别表示两路信号的幅值，ω 为信号的角频率，φ为信号的初始相位，θ为两路信号的相位差。

按照习惯表达，先求压电陶瓷换能器的阻抗，再取倒数得到导纳。

m
B B A R U U U I U Z -== 将式（2-14）代入得
m
j Bm j Bm j Am R e U e U e U Z )(θωϕθωωϕ++-==)1(--θj Bm Am m e U U R =jX R U U j U U R Bm
Am Bm Am m +=--]sin )1cos [(θθ 对应得到：)1cos (
-=θBm Am m U U R R ，θsin Bm Am m U U R X -= （2-15） 再由导纳和阻抗的关系可得
jB G jX
R Z Y +=+==11 即：22X R R G +=，22X R X B +-= （2-16） 由以上推导可以看出，换能器的导纳和阻抗值仅与加在其两端的电压信号的幅值比和相位差有关，因此只需要得到两路信号的幅值和相位信息即可得到换能器等效电路的各个参数。

而实际中，只需要对两路信号进行采样，再通过对采样所得数据进行处理便可得到幅值和相位信息。

2.4 正弦信号的测量方法
根据上一节介绍的测量原理可知，要得到压电换能器在测试频率下的电导和电纳值，就需要测得其两端正弦信号的幅值比和相位差。

但是实际中，硬件电路实现的仅是对两路信号的 A/D 转换采集，也即是得到的是两路正弦信号的一系列的离散的点。

在这一节中，将介绍从这些采集到的离散的点计算其幅值和相位的方法。

2.4.1 数字相关法
随着微处理器和大规模集成电路的迅速发展，在测试系统中，越来越多的传统的测量方法被数字化测量方法所取代。

近年来，由于相关函数法具有提高测试精度，减少或简化硬件设计，能够充分利用测试系统中的数据采集系统和微型计算机，提高测试系统的可靠性和可维护性的诸多优点，使得相关技术原理在相位差的测量及数字信号处理中得到了广泛应用，并展现出良好的应用前景
1、相关函数法原理
相关函数法利用两同频正弦信号的延时为零时的互相关函数值与其相位差的余弦值成正比的原理获得相位差。

设两路被测信号为：
)()2sin()(t N ft A t x x +=π，)()2sin()(t N ft B t y y ++=ϕπ （2-17）
其中：A 、B 分别表示两路信号的幅值，f 表示信号的频率，)(t N x 、 )(t N y 分别表示两路信号的干扰噪声信号，φ表示两路信号的相位差。

显然，信号 x(t)和y(t)是相关的，则两路信号的互相关函数为：
dt t N t f B t N ft A T R dt t y t x T R T y x xy T xy ⎰⎰++++⋅+=+=
00)]())(2sin([)]()2sin([1)()()(1)(τϕτππτττ （2-18） 式中 T 为信号的周期，即f
T 1=
当 τ =0时，有 dt t N ft B t N ft A T
R y T x xy )]()2sin([)]()2sin([1)0(0++⋅+=⎰ϕππ 由于噪声信号之间不相关，噪声和信号之间也不相关，将上式进一步展开得：
⎰+=
T xy dt ft ft AB T
R 0)2sin()2sin(1ϕππ ⎰+=T dt ft ft ft AB T 0
]sin )2cos(cos )2)[sin(2sin(1ϕπϕππ ⎰==T AB dt ft ft AB T 0cos 21cos )2sin()2sin(1ϕϕππ 所以，可以得到相位差的计算公式：
))
0(2arccos(AB R xy =ϕ (2-19)
而信号幅值的大小可由信号的自相关函数求得
dt t f A ft A T dt t x t x T R T T T T x ))(2sin()2sin(1)()(1)(2222
τππττ+=+=⎰⎰-- (2-20) 当 τ =0时，有
⎰-==222
222
)2(sin 1)0(T T x A dt ft A T R π 所以可得信号幅值的计算公式：
)0(2)
0(2y x R B R A == (2-21)
将上式代入式（2-19），可得相位差计算公式的另一种表达式：
))0()0()
0(arccos(y x xy R R R =ϕ (2-22)
而在实际中，是没有完整精确的信号的表达式的，有的是对信号的模数转换所得到的离散的数据，离散序列的自相关和互相关的计算公式如下：
∑∑∑-=-=-====1
2
1
2
1
)(1)0()(1)0()()(1)0(n i y n i x n i xy i y n R i x n R i y i x n R (2-23)
式中：n 表示采样个数，i 表示第 i 个采样点，x(i)、y(i)分别表示两路信号的第 i 个点的转换得到数值。

由式（2-23）分别求出两路信号的自相关和互相关函数值之后，再由式（2-21）和式（2-22）即可得到两路信号各自的幅值和它们之间的相位差。

但是，需要指出的是，由数字相关法求得的相位差，并不能区分是超前还是滞后，这就需要采用其他方法来确定相位差符号的正负号。

根据前面测量原理中的介绍，由式（2-15）和式（2-16）可知，压电陶瓷换能器的电导值仅取决于两路信号相位差的余弦值，而电纳的值是在电导值取得最大的时候发生变号。

由此，可以先求得电导的值，再通过循环找其最大值，并从使电导取得最大值时的相位差开始，把相位差变号，得到新的相位差序列，再由新的相位差序列求电纳的值即可。

图2-9 表示的为采用数字相关法对一号压电换能器测量数据的处理结果，其中（a ）表示的是导纳圆图，（b ）表示的是电导和电纳值随测试频率的变化曲线。

(a) (b)
图2-9 数字相关法处理结果
2、相关函数法的特点及误差分析
通过上面对相关函数法测量原理的理论推导过程可以看出，相关函数法测量信号的幅值和相位差与信号的频率无关。

也即是说相关函数法不受频率的影响，可以用来测量未知频率的信号的相位差。

同时，相关函数法测量原理的推导都是基于正弦函数的，因此，它只能用于测量正弦或余弦信号，并不能测量一般的周期信号。

由于噪声干扰信号和原信号并不相关，所以相关函数法能够有效的抑制噪声干扰。

但是，如果在系统中存在相关性较强的干扰信号，并且信噪比又比较低的情况下，相关函数法测量误差就会比较大。

由相关函数法离散序列的最终计算公式可以看出，其计算结果与采样的点数有关，也即是说测量误差的大小与采样点数是相关的，采样点数越大，计算结果越接近真实值，测量误差也就越小。

综合以上对相关函数法的特点的分析，可知相关函数法对于采样转换信号中的直流偏移和噪声等干扰具有很强的抑制能力，它的误差主要是因为采用有限长度的样本代替了高斯白噪声和均匀分布的A/D 量化误差，使得被检正弦信号与噪声信号并非完全不相关。

所以，相关函数法的测量误差与A/D 转换的位数、信号的信噪比和采集点数有关。

2.4.2 快速离散傅里叶变换法
现代信号分析采用数字化方式实现，其核心是离散傅立叶变换，它完成了从时域到频域的转换，不仅可以实现线性谱分析，而且还是均方谱分析的关键。

离散傅立叶变换（DFT ）实现了信号首次在频域表示的离散化，使得频域也能够用计算机进行处理，但由于用于实际时计算量太大而使应用受到限制。

直到1965 年由 Cooly 和 Tukey 建立了一种快速傅立叶变换——FFT 时，DFT 的应用才成为现实
1、FFT 获取正弦波幅值和相位的原理
设采集正弦信号得到的离散序列为 x(n)，n=1,2，…KN 。

则该序列的离散傅里叶变换为：
∑-=-==102)()]([)(N n n N k j e
n x n x DFT k X π
∑-=+=-=10)](Im[)](Re[)]2sin()2)[cos(
(N n k X k X n N
k j n N k n x ππ 1-N ,2,1,0⋯⋯=k （2-24）
则其初始相位为：
1),)](Im[)](Re[arctan(0+==s
f f k k X k X θ 其中：s f 是信号的采样频率，N 是采样长度。

在对时域离散序列进行傅立叶变换之后，可以得到其离散的幅度谱和相位谱，在幅度谱和相位谱中找到对应时域波形的频率的谱线就可以得到时域的正弦波形的幅值和相位信息。

图 2-10 所示的是采用快速离散傅里叶变换法对采集到的数据处理的结果。

（a ） (b)
图 2-10 快速离散傅里叶变换法处理结果
2、FFT 的特点及误差分析
通过傅里叶变换可以只提取基波参数，因此谐波的存在并不影响基波成分，所以谐波的存在对应用这种方法测量相位差几乎没有影响；对于噪声干扰，只有当高斯白噪声接近基波的频率分量时才会影响到基波的相位，所以应用 FFT 法测量相位差也能有效地抑制高斯白噪声干扰。

但是，实际上信号是连续的无限长的序列，用 FFT 对其进行谱分析时，必须截短形成有限长序列，再进行周期延拓，这样就不可避免的造成信号频谱的泄漏，由此便产生了相位差测量误差。

误差现象主要是：混叠现象、栅栏效应和截断效应。

要想减小相位差测量误差，就必须提高谱分辨率。

实际中可通过提高采样频率或者增加采样数据长度来提高谱分辨率，进而达到减小相位差测量误差的目的。

2.4.3 正弦曲线参数拟合法
设被测的正弦信号为：
D ft A t f ++=)2sin()(0ϕπ (2-25)
其中： f 表示信号频率，0A 表示被测信号幅值，φ表示被测信号的初始相位角，D 表示被测信号的直流分量。

由于被测信号的频率为已知的，故只需对测得的数据进行三参数的正弦曲线拟合，即可得到被测信号的幅值和相位信息。

为此，进一步将上式展开可得：
D ft B ft A f ++=)2cos()2sin(ππ (2-26)
其中：ϕϕsin ,cos 00A B A A ==
从而将被测信号的幅值和初始相角转化为对参数 A 、B 的求取。

其基本思想就是寻找合适的 A 、B 和 D 的值，使得其测量残差的平方和取得最小。

设每个频率下测量的时间序列为)n 3,2,1(,1⋯==i f t s
i ，n 为测量数据的个数，s f 为采样频率，i f 为每个点的测量值。

则测量残差的表达式为：
∑=---=n
i i i D ft B ft A t f 12])2cos()2sin()([ππε (2-27)
要使得上式取得最小值，可对其参数求偏导，并令其为零。

即： }}}
⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧⎩⎨⎧=---∂∂=∂∂⎩⎨⎧=---∂∂=∂∂⎩⎨⎧=---∂∂=∂∂∑∑∑===n i i i n i i i n i i i D ft B ft A i f D D D ft B ft A i f B B D ft B ft A i f A A 121
21
20])2cos()2sin()([0])2cos()2sin()([0])2cos()2sin()([ππεππεππε (2-28)
进一步化简得到：
⎪⎪⎪⎩⎪⎪⎪⎨⎧=++=++=++∑∑∑∑∑∑======n
i n i i i n
i n i i i i i n
i n i i i i i i f D ft B ft A i f ft D ft B ft A ft i f ft D ft B ft A ft 0
00
00
0)
(])2cos()2sin([)()2cos(])2cos()2sin()[2cos()
()2sin(])2cos()2sin()[2sin(ππππππππππ (2-29)
对于式（2-29），构造如下三个矩阵：
⎥⎥⎥⎥
⎦⎤⎢⎢⎢⎢⎣
⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=)(1...)(1)(1,,11112cos ...2cos 2cos 2sin ...2sin 2sin 212121n n n T t f t f t f F D B A X ft ft ft ft ft ft ππππππψ 式（2-29）可写成如下矩阵形式：
F X T T ψψψ=
上式中 X 的解为：
)()(1F X T T ψψψ-= (2-30)
则被测信号幅值的计算公式为: 220
B A A += (2-31)
初始相角的计算公式为：
⎪⎪⎩
⎪⎪⎨⎧<+≥=0)arctan(0)arctan(A A B A A B πϕ (2-32)
至此，得到了被测信号的幅值和相位信息。

采用同样的方法对第二路信号的采样数据进行处理，即可得到第二路信号的幅值和相位信息，从而求出两路信号的幅值比和相位差，进一步便可得到每个频率下换能器的电导和电纳值。

图 2-11为参数拟合法的处理结果。

图 2-11 正弦曲线参数拟合法处理结果
2.5 导纳圆的带约束最小二乘曲线拟合
通过以上章节的介绍，我们已经得到了各个测试频率下压电换能器的电导和电纳值，绘制出了导纳圆图，但这还是不够的。

由压电换能器等效电路的各个参数的计算公式可以看出，我们还需要得到导纳圆的圆心和半径的值。

为此，就需要对所得到的离散点进行圆曲线拟合。

拟合圆的方法有很多种，常用的有平均值法、加权平均法和最小二乘法。

平均值法的思想是分别计算各个离散点的横、纵坐标的平均值，作为圆心的横、纵坐标，将圆心到各个离散点的距离的平均值作为半径。

这种方法计算简单，适用于离散点分布较均匀的情况，但对于分布不均的情况，所计算的圆心位置会偏向离散点分布较密集的一侧，半径的计算值也会偏小，误差较大。

加权平均法是对平均值法的改进，它在计算圆心坐标时加入一个与两相邻点间弧长相关的系数，降低了离散点分布不均的影响，减小了误差。

但是，由于两相邻点间的弧长是无法精确得到的（实际中采用两点间的距离来代替的），所以其误差仍较大。

相比之下，最小二乘法具有较高的精度。

在查阅相关文献之后，本文采取一种带 约束的最小二乘法进行导纳圆的拟合
设计算出的电导和电纳值的离散点为：
N i y x i i ,...,2,1),(=
要得到的圆心坐标记为:),(000y x P 半径为 r 。

则构造一个目标函数：
{}N i y x y x r y y x x d i i ,...,2,1),(),(,)(2)(00=∈--+-= (2-33)
从而，问题转化为寻找合适的0x ，0y 和 r ，使得目标函数的模值取得最小。

由于上式中含有根号，不便于处理和计算，对其进行变形得到：
])()[(2122020r y y x x r
d --+-=
=)22(21222220000r y x yy x x y x r -++--+ (2-34) 式（2-34）和式(2-33)对圆的几何特征参数的估计是一致的。

令)(21,,,2120202021r y x r c b b B y x X r a -+=⎥⎦
⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡== 则式(2-34)可表示为：
c X B X aX X F
d T T T ++==)(
同时有：142221=-+ac b b （2-35）。