METHOD FOR MEASURING EQUIVALENT CIRCUIT PARAMETERS AND RESONANT FREQUENCY OF PIEZOELECTRIC RESONATOR

Abstract

A method for measuring equivalent circuit parameters and resonant frequency of a piezoelectric resonator, by which the phase-frequency curve of the piezoelectric resonator is measured, and the resonant frequency and the anti-resonant frequency are obtained. Then, the slopes of the phase-frequency curve at the resonant frequency and the anti-resonant frequency are respectively measured. The resonant angular frequency and the anti-resonant angular frequency are also calculated. Finally, the equivalent circuit parameters of the piezoelectric resonator are obtained by solving a system of nonlinear equations.

Description

CROSS REFERENCE TO RELATED APPLICATIONS

This application is based upon and claims priority to Chinese Patent Application No. 201811393865.3, filed on Nov. 21, 2018 the entire contents of which are incorporated herein by reference.

TECHNICAL FIELD

The present invention belongs to the field of electronic technology, and more specifically relates to a method for measuring equivalent circuit parameters and resonant frequency of a piezoelectric resonator.

BACKGROUND

All electrostatically driven piezoelectric resonators may be equivalent to BVD (Butterworth-Van Dyke) models (according to V. E. Bottom, Introduce to crystal unite design. York: Van Nostrand Reinhold Company, pp:82, 1987.). A variety of piezoelectric materials, including quartz crystal, lithium tantalate, and piezoelectric ceramic lead zirconate titanate (PZT), can be equivalent to the BVD model. The piezoelectric resonators of some kinds, such as bulk acoustic wave (BAW) resonators, surface acoustic wave (SAW) resonators, quartz crystal microbalance (QCM), etc., can also be equivalent to the BVD model.

There are many methods which can be used to measure the parameters of the BVD equivalent model. For example, CI-miter method, impedance meter, measurement of quartz crystal device parameters by zero phase technology, and measurement of quartz crystal device parameters based on resonant frequency, anti-resonant frequency, load resonant frequency, and load anti-resonant frequency (according to Chinese patent No. CN201610556735.1, entitled Method For Measuring Equivalent Circuit Parameters Of Quartz Crystal Resonator, inventors Yang Wang, Hai Yun Shen, Dong Liu), etc.

Instruments for measuring equivalent circuit parameters currently available on the market include impedance meters, vector network analyzers, and the like, such as KH1800/KH1820 manufactured by a Hong Kong Kolinker Industrial Equipments Co. Ltd.

At present, measurement methods based on the IEC60444 standard are mainly used, where the motional resistance are replaced with the resonant resistance. It is reported from the simulations that a more than 20% error exists between the motional resistance and the resonant resistance when the resonator has a low Q value (according to Dong Liu, Research on Temperature Compensation Technology of Crystal Oscillator [D]. University of Electronic Science and Technology, 2017.).

The method proposed in the Chinese patent No. CN201610556735.1 (Yang Wang, Hai Yun Shen, Dong Liu, Method For Measuring equivalent Circuit Parameters Of Quartz Crystal Resonator, CN201610556735.1), does not have the above error theoretically. However, it has difficulty in practical operations. Therefore, as the frequency drifts with the temperature, a small temperature change would lead to great error in the calculation results. Liu Dong et al. proposed a method for calculating the equivalent parameters based on resonant frequency, anti-resonant frequency, load resonant frequency, load anti-resonant frequency, derivative of phase-frequency curve at the resonant frequency, and derivative of phase-frequency curve at the anti-resonant frequency. (D. Liu, X. H. Huang, J. G. Hu, Y. L. Tang, and Y. Wang, “Measurement of quartz crystal unit parameters based on the precise derivation of zero phase frequency” Electronics Letters, vol. 53, p. 2, Feb 2017.). According to the method, the phase-frequency curve of the resonator with and without load capacitance are needed. Therefore, this method requires complicated operations, and it is difficult been integrated in the current scientific instruments as software modules. (e.g. vector network analyzer and impedance meter, etc.).

All the methods described above fail to take the measurement errors of the resonant frequency, the anti-resonant frequency, the load resonant frequency, and the load anti-resonant frequency into account. For example, the measurement accuracy of the resonant frequency measured by the IEC60444-based method is susceptible to the variations of wire length and load. FIG. 1 demonstrates the error in related to the variations of wire length and load. Such an error will be passed to the parameter measurement, thereby affecting the accuracy of the parameter measurement.

SUMMARY

It is therefore an objective of the present invention to overcome the drawbacks of the prior art and provide a method for measuring equivalent circuit parameters and resonant frequency of a piezoelectric resonator without the need of impedance matching or a peripheral circuit such as π-type network or load capacitance, so as to measure the resonant frequency, anti-resonant frequency and equivalent circuit parameters.

To this end, the method for measuring equivalent circuit parameters and resonant frequency of the piezoelectric resonator provided by the present invention includes the following steps:

(1) measuring a phase-frequency curve of the piezoelectric resonator and finding the zero phase point on the phase-frequency curve, namely, obtaining a resonant frequency ƒ_r and an anti-resonant frequency ƒ_a;

(2) measuring a slope of the phase-frequency curve at the resonant frequency ƒ_r, recording the slope as Δ₁; and measuring a slope of the phase-frequency curve at the anti-resonant frequency ƒ_a, recording the slope as Δ₂;

(3) calculating a resonant angular frequency with a formula ω_r=2πƒ_r, and calculating an anti-resonant angular frequency with a formula ω_a=2πƒ_a;

(4) substituting ω_r, ω_a, Δ₁, and Δ₂ into the following system of nonlinear equations and solving the nonlinear equations to obtain the equivalent electrical parameters of the piezoelectric resonator: static capacitance C₀, motional capacitance C₁, motional inductance L₁, and motional resistance R₁;

$\{\begin{array}{c}\frac{{\omega }_a^2+{\omega }_r^2}{2}=\frac{1}{L_1C_1}+\frac{1}{2L_1C_0}-\frac{R_1^2}{2L_1^2}\\ {\omega }_r^2{\omega }_a^2=\frac{1}{L_1^2C_1^2}\left(1+\frac{C_1}{C_0}\right)\\ \frac{\partial \left(\mathrm{Pashe}\left(f\right)+{\mathrm{Phase}}_{\mathrm{offset}}\right)}{\partial f}{|}_{f=f_r}={\Delta }_1\\ \frac{\partial \left(\mathrm{Pashe}\left(f\right)+{\mathrm{Phase}}_{\mathrm{offset}}\right)}{\partial f}{|}_{f=f_a}={\Delta }_2\end{array}\hspace{1em}$

where Phase_offset is a phase offset value, Pashe(ƒ) is a phase-frequency function, and ƒ is a frequency point on the phase-frequency curve.

The objective of the present invention is realized according to the following processes.

According to the method for measuring the equivalent circuit parameters and the resonant frequency of the piezoelectric resonator provided by the present invention, the resonant frequency and the anti-resonant frequency are obtained by measuring the phase-frequency curve of the piezoelectric resonator; subsequently, the slopes of the phase-frequency curve at the resonant frequency and anti-resonant frequency are respectively measured, and then the resonant angular frequency and the anti-resonant angular frequency are calculated, and last the equivalent circuit parameters of the piezoelectric resonator are obtained by solving the system of the nonlinear equations.

At the same time, the method for measuring the equivalent circuit parameters and the resonant frequency of the piezoelectric resonator according to the present invention has the following advantages:

(1) the present invention neither need additional operations such as impedance matching, instrument zero phase point adjustment, nor peripheral circuit such as the π-type network and the load capacitance; in contrary, only the phase-frequency curve is required to calculate the four equivalent circuit parameters of the piezoelectric resonator;

(2) the accuracy of the resonant frequency measurement is susceptible to the errors induced by the variations of wire length and load; the present invention avoids the influence of the errors, and can directly calculate the zero phase point drift caused by the variations of wire length and load, so that the frequency measurement is more accurate and the frequency measurement error will not be passed to the measurement of parameters;

(3) the calculated parameters highly fit the measured phase-frequency curve, so the parameter measurement accuracy is high;

(4) the R1 measured according to the present invention is not affected by the variation of the Q value of the resonator, so the measurement of the motional resistance R1 is more accurate;

(5) while the present invention does not adopt an approximate model, thus theoretically the motional capacitance C1 is not affected by the M value, and thereby the measurement of the motional capacitance C1 is not affected by the M value of the resonator; and

(6) the static capacitance C0 is a capacitance value of a resonant frequency segment, while in the IEC-based method, the measured static capacitance C0 is higher than the capacitance value of the resonant frequency segment. The static capacitance C0 measured with the method according to the present invention is exactly the capacitance value of the resonant frequency segment, thus the capacitance value is more accurate.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a diagram showing a relation between a measurement error of a resonant frequency and variations of wire length and load according to an IEC-based method;

FIG. 2 is a diagram showing a BVD equivalent model of a piezoelectric crystal;

FIG. 3 is a comparison block diagram of three measurement solutions;

FIG. 4 is a schematic diagram of an ADS software (Advanced Design System) simulation;

FIG. 5 is a diagram showing a simulation result of the ADS software; and

FIG. 6 is a diagram showing an actually measured phase-frequency curve of a QCM (Quartz Crystal Microbalance) subjected to a water load and the phase-frequency curve derived from the calculated parameters.

DETAILED DESCRIPTION OF THE EMBODIMENTS

The specific embodiments of the present invention will be described hereinafter with reference to the drawings, so that those skilled in the art can better understand the present invention. It should be particularly noted that in the following description, when a detailed description of the well-known functions and designs may downplay the main content of the present invention, these descriptions will be omitted herein.

EMBODIMENT

FIG. 2 is a diagram showing a BVD equivalent model of a piezoelectric crystal.

In the present embodiment, the piezoelectric resonator may be a high-Q value crystal piezoelectric resonator, a quartz crystal microbalance piezoelectric resonator, or a microelectromechanical system piezoelectric resonator. As shown in FIG. 2, the piezoelectric crystal is connected to a vector network analyzer to measure S21 characteristics, and the measurement mode is switched to a phase measurement mode. Besides, the piezoelectric resonator may be made of quartz, lithium tantalate, lanthanum gallium silicate, piezoelectric ceramic lead zirconate titanate (PZT), and aluminum nitride (AlN).

The method for measuring equivalent circuit parameters and resonant frequency of a piezoelectric resonator of the present invention includes the following steps:

S1: the phase-frequency curve of the piezoelectric resonator is actually measured through the vector network analyzer and the zero phase point on the phase-frequency curve is found out, i.e., the resonant frequency ƒ_r and the anti-resonant frequency ƒ_a are obtained;

S2: a slope of the phase-frequency curve at the resonant frequency ƒ_r is measured and recorded as Δ₁; and a slope of the phase-frequency curve at the anti-resonant frequency ƒ_ais measured and recorded as Δ₂;

S3: a resonant angular frequency is calculated with a formula ω_r=2πƒ_r, and an anti-resonant angular frequency is calculated with a formula ω_a=2πƒ_a; and

S4: ω_r, ω_a, Δ₁, and Δ₂ are substituted into the following system of nonlinear equations and the nonlinear equations are solved to obtain the equivalent circuit parameters of the piezoelectric resonator: static capacitance C₀, motional capacitance C₁, motional inductance L₁, and motional resistance R₁;

$\{\begin{array}{c}\frac{{\omega }_a^2+{\omega }_r^2}{2}=\frac{1}{L_1C_1}+\frac{1}{2L_1C_0}-\frac{R_1^2}{2L_1^2}\\ {\omega }_r^2{\omega }_a^2=\frac{1}{L_1^2C_1^2}\left(1+\frac{C_1}{C_0}\right)\\ \frac{\partial \left(\mathrm{Pashe}\left(f\right)+{\mathrm{Phase}}_{\mathrm{offset}}\right)}{\partial f}{|}_{f=f_r}={\Delta }_1\\ \frac{\partial \left(\mathrm{Pashe}\left(f\right)+{\mathrm{Phase}}_{\mathrm{offset}}\right)}{\partial f}{|}_{f=f_a}={\Delta }_2\end{array}\hspace{1em}$

where Phase_offset is a phase offset value, Pashe(ƒ) is a phase-frequency function value, and ƒ is a frequency point on the phase-frequency curve.

In the present embodiment, assuming there are N frequency points on the phase-frequency curve, then ƒ may be represented as ƒ₁, ƒ₂, . . . , ƒf_N;

then the phase-frequency function value Pashe(ƒ_i) on the i_th frequency point meets the following equation:

$\mathrm{Pashe}\left(f_i\right)=\frac{180}{\pi }·\mathrm{atan}\left(\frac{\mathrm{imag}\left(G\left(f_i\right)\right)}{\mathrm{real}\left(G\left(f_i\right)\right)}\right)$ $\mathrm{where},G\left(f_i\right)=\frac{1}{Z\left(f_i\right)+R_0},$

Z(ƒ_i) is an impedance of the piezoelectric crystal, imag() represents an imaginary part, real()represents a real part, and R₀ is an internal resistance of the vector network analyzer;

then a phase difference of each frequency point is:

Pashe_i=Pashe(ƒ_i)+Phase_offset−Pashe_measure(ƒ_i)

where Pashe_measure(ƒ_i) represents the actual phase value of the i_th frequency point ƒ_i;

then a root mean square error of the phase differences is calculated as follows:

$\mathrm{RMSE}=\sqrt{\frac{1}{N}\sum _1^N\Delta {\mathrm{Pashe}}_i}$

The RMSE is minimized by continuously changing the Phase_offset, at that time, the calculated equivalent circuit parameters, resonant frequency and anti-resonant frequency are final results.

FIG. 3 is a comparison block diagram of three measurement solutions.

In the present embodiment, the IEC-based measurement method is shown in FIG. 3(a). The measurement method proposed by Dong Liu in 2017 is shown in FIG. 3(b), and the measurement method of the present invention is shown in FIG. 3(c).

The IEC-based method requires a pi-type network to facilitate impedance matching and has a relatively greater error.

The measurement method proposed by Dong Liu in 2017 requires load capacitance and series resistance, and the measurement error is related to the resistance. This method also needs to switch the state between with and without load-capacitance, which is difficult been integrated in the current scientific instruments as software modules.

The technical solution of the present invention requires no external circuit and no impedance matching, and is directly connected in series to the vector network analyzer. This solution is easy to operate, and has no theoretical error and small measurement error. Simply by measuring four points of the phase-frequency curve, so the resonant frequency, the anti-resonant frequency, and four equivalent circuit parameters can be calculated. The phase-frequency curve inversed from the parameters fits with the actually measured curve at more than 800 points, and the root mean square error of the 800 points is below 0.1107.

Embodiment 1

The simulation experiment was carried out with the ADS (Advanced Design System) software. The schematic diagram is shown in FIG. 4. The phase-frequency curve of the simulation is shown in FIG. 5. The results obtained are shown in Table 1.

	TABLE 1
		R1(Ω)	L1(mH)	C1(fF)	C0(pF)
	Set value	14.0000	75.0000	3.3600	3.0000
	Calculated value	14.0050	75.0000	3.3600	3.0000

Embodiment 2

The solution of the embodiment 2 is similar to embodiment 1. With varying set values of the motional resistance R₁, the calculated motional resistance R₁, which is obtain from the method proposed by the present invention, is shown in Table 2 respectively. The resonant resistance Rr, is also shown in Table 2. The resonant resistance Rr can be obtain from substitution methods similar to IEC.

TABLE 2
	Rr calculated	R1 calculated
	by the	by the
Set value	IEC-based	proposed
of R1	method	method
14	14.0094	14.0050
50	55.4998	50.0019
70	90.4585	70.0037

As shown in Table 2, the motional resistance calculated according to the method proposed by the present invention is closer to the set value. The resonant resistance is close to the set value, when the set R1 is small. However, the resonant resistance Rr is different from the motional resistance R1, when the set R1 is bigger.

Embodiment 3

The parameters of a QCM wafer loaded by a drop of water are measured with the vector network analyzer N9913A. The phase-frequency curve measured with the N9913A is shown as a solid line in FIG. 6. According to our nonlinear equations and a MATLAB numerical solution thereof, the obtained equivalent parameters are shown in Table 3.

TABLE 3
parameter	Rr (Ω)	R1(Ω)	L1(mH)	C1(fF)	C0(pF)
QCM water	380	321.8300	12.8864	20.0005	4.2462
load

The phase-frequency curve inverted from the parameters measured by the method of the present invention is shown as a dotted line in FIG. 6, and it can be seen that the inversion curve and the measured curve are almost coincident.

The resonant resistance Rr measured by the substitution method, which is similar to the IEC-based method, is also given in Table 3. It can be concluded that the R₁ measured by the method of the present invention is smaller than the Rr which is consistent with the simulation results. Also, it can be concluded from Table 3 that when QCM is loaded with liquid, a large gap exists between Rr and R₁.

Embodiment 4

The phase offsets caused by different loads are measured. The phase offsets of the same QCM wafer with no load, water load, or oil load are measured by the method of the present invention and shown in Table 4.

	TABLE 4
	QCM load
	situation	No-load	Water load	Oil load
	phase offset	−15.0441	−15.6539	−17.0441
	value (deg.)

It can be concluded that the present invention can directly measure the zero phase point drifts caused by the variations of load. Therefore, the present invention does not have the measurement errors caused by the variations of load as shown in FIG. 1.

Embodiment 5

The phase offsets caused by different wire lengths are measured. Two quartz crystals with brackets of different lengths are numbered as 1# and 2#, respectively. The phase offsets measured by the method of the present invention are shown in Table 5.

	TABLE 5
	Crystal number	Crystal 1#	Crystal 2#
	phase offset value	−45.0005	−35.1400
	(deg.)

It can be concluded that the method of the present invention can directly measure the phase offset caused by the variations of load or the wire length without introducing the measurement errors caused by the variations of load and the wire length as shown in FIG. 1.

The specific embodiments of the present invention are described above to facilitate the understanding of those skilled in the art. It should be noted, however, that the present invention is not limited to the scope of the specific embodiments. As for those of ordinary skill in the art, as long as the various variations fall within the spirit and scope of the present invention as defined and determined by the appended claims, these variations are obvious. All inventions derived from the idea of the present invention are covered by the present invention.

Claims

1. A method for measuring equivalent circuit parameters and a resonant frequency of a piezoelectric resonator comprising: measuring a phase-frequency curve of the piezoelectric resonator and finding a zero phase point on the phase-frequency curve, namely, obtaining the resonant frequency ƒr and an anti-resonant frequency ƒa;measuring a slope of the phase-frequency curve at the resonant frequency fr, recording the slope as Δ1; and measuring a slope of the phase-frequency curve at the anti-resonant frequency ƒa, recording the slope as Δ2;calculating a resonant angular frequency with a formula ωr=2πƒr, and calculating an anti-resonant angular frequency with a formula ωa=2πƒa; andsubstituting ω, ωa, Δ1, and Δ2 into the following system of nonlinear equations and solving the nonlinear equations to obtain the equivalent circuit parameters of the piezoelectric resonator including a static capacitance C0, a motional capacitance C1, a motional inductance L1, and a motional resistance R1;

2. The method for measuring the equivalent circuit parameters and the resonant frequency of the piezoelectric resonator according to claim 1, wherein the Pashe(ƒ) meets the following formula:

3. The method for measuring the equivalent circuit parameters and the resonant frequency of the piezoelectric resonator according to claim 1, wherein the piezoelectric resonator comprises a high-Q value crystal piezoelectric resonator, a quartz crystal microbalance piezoelectric resonator, and a microelectromechanical system piezoelectric resonator.

4. The method for measuring the equivalent circuit parameters and the resonant frequency of the piezoelectric resonator according to claim 1, wherein the piezoelectric resonator is made of quartz, lithium tantalate, lanthanum gallium silicate, piezoelectric ceramic lead zirconate titanate (PZT), or aluminum nitride (AlN).