CN115060355A - Harmonic oscillator quality factor measuring method based on linear frequency modulation pulse - Google Patents
Harmonic oscillator quality factor measuring method based on linear frequency modulation pulse Download PDFInfo
- Publication number
- CN115060355A CN115060355A CN202210380357.1A CN202210380357A CN115060355A CN 115060355 A CN115060355 A CN 115060355A CN 202210380357 A CN202210380357 A CN 202210380357A CN 115060355 A CN115060355 A CN 115060355A
- Authority
- CN
- China
- Prior art keywords
- frequency
- quality factor
- range
- value
- harmonic oscillator
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000001228 spectrum Methods 0.000 claims abstract description 18
- 238000005070 sampling Methods 0.000 claims abstract description 7
- 230000003595 spectral effect Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000005284 excitation Effects 0.000 claims description 4
- 238000010408 sweeping Methods 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 2
- 230000001131 transforming effect Effects 0.000 abstract 1
- 230000010355 oscillation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H13/00—Measuring resonant frequency
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention provides a harmonic oscillator quality factor measuring method based on a linear frequency modulation pulse signal(ii) a For response signalSampling, FFT transforming, modulus taking and normalizing to obtain frequency spectrum(ii) a Cut-select swept range spectrumCalculating to obtain the estimated value of the vibration frequency(ii) a Establishing a Lorentz curve model and a frequency range of the first Lorentz fittingQ value range 1 Q M,Performing a first estimation of the figure of merit; combining the first estimateIn frequency sweep universeRange of Q valueAnd carrying out second Lorentz fitting to obtain a final quality factor estimation value. The method not only improves the larger error of the half-power bandwidth method in the low signal-to-noise ratio environment, but also saves the operation time through the two-step fitting process, and particularly has obvious advantages for some harmonic oscillators with high Q values.
Description
Technical Field
The invention relates to the field of resonant device parameter measurement, in particular to a harmonic oscillator quality factor measuring method based on linear frequency modulation pulse.
Background
The quality factor is an important parameter of the vibration characteristic of the resonant device, and plays an important role in drive control and error mechanism analysis. The quality factor measuring method is mainly divided into a time domain method and a frequency domain method. The time domain method is based on the amplitude attenuation characteristic during free oscillation, and is calculated by utilizing the attenuation degree of the oscillation amplitude in specific time and the corresponding relation between the attenuation degree and the Q value, the peak value or the amplitude envelope of transient response needs to be detected, and the method is relatively complex. The frequency domain method is commonly used, and is generally calculated by using a resonance frequency and a pass band width based on an amplitude-frequency characteristic curve, and is called a half-power bandwidth method. However, when the signal-to-noise ratio is low, the reliability of the half-power bandwidth method is low, and the error is large.
Disclosure of Invention
In order to solve the problems, the invention discloses a harmonic oscillator quality factor measuring method based on linear frequency modulation pulses, which aims to solve the problem of large quality factor measuring error under the condition of low signal-to-noise ratio.
In order to achieve the above purpose, the invention adopts the following scheme:
a harmonic oscillator quality factor measuring method based on linear frequency modulation pulse comprises the following steps:
step 1: carrying out frequency sweeping excitation on the harmonic oscillator by using a chirp signal to obtain a response signal s (t);
step 2: sampling, FFT (fast Fourier transform), modulus taking and normalization are carried out on the response signal S (t) to obtain a frequency spectrum | S' (n) |; interception and selection sweep range spectrum | S 1 ' (n) | calculation obtains the estimated value of the vibration frequency
And step 3: and establishing a Lorentz curve model, sequentially carrying out Lorentz fitting twice aiming at different frequency ranges, and combining the resonant frequency value obtained in the previous step to obtain the quality factor.
The step 3 comprises the following steps:
step 3.1: establishing Lorentz curve model Is the estimated resonant frequency, Q is the quality factor, and the Q max is set to M. If Q ═ m, drawing L (f) non-calculation Q=m And normalizing to obtain L' (f) dichotomy Q=m ;
Step 3.2: first Lorentzian fitting frequency rangeThe Q value range is more than or equal to 1 and less than or equal to M. Subtracting the amplitude of each data point in the actual frequency range from the amplitude of the frequency point corresponding to the fitting curve point by point and accumulating to obtain a residual error (m);
step 3.3: index residual error sequence error (i)) The point with the smallest value, the sequence number is denoted index0, and the figure of merit of the first estimation is obtained
Step 3.4: the second Lorentz fitting f range is in the whole sweep frequency domain interval f 1 ≤f≤f 2 Selecting the range of Q value asRepeating the step 3.3 to obtain a residual error sequence error (i);
step 3.5: the point with the minimum value in the index residual error sequence error (i) is marked as index1 to obtain the final estimation value of the quality factor
The invention has the beneficial effects that: the method not only improves the larger error of the half-power bandwidth method in the low signal-to-noise ratio environment, but also saves the operation time through the two-step fitting process, and particularly has obvious advantages for some harmonic oscillators with high Q values.
Drawings
Fig. 1 is a flowchart of a harmonic oscillator quality factor measurement method based on chirp pulses provided in the present invention;
FIG. 2 is a comparison of simulation results of the half-power bandwidth method and the method of the present invention with respect to error.
Detailed Description
The present invention will be further illustrated with reference to the accompanying drawings and specific embodiments, which are to be understood as merely illustrative of the invention and not as limiting the scope of the invention. It should be noted that the terms "front," "back," "left," "right," "upper" and "lower" used in the following description refer to directions in the drawings, and the terms "inner" and "outer" refer to directions toward and away from, respectively, the geometric center of a particular component.
As shown in fig. 1-2, the method for measuring the quality factor of a harmonic oscillator based on chirp of this embodiment sets chirpParameters of the signal: starting frequency f 1 End frequency f 2 Sweep frequency time T, signal amplitude A, generating chirp signal x (T):
Carrying out frequency sweeping excitation on the harmonic oscillator by using a chirp signal to obtain a response signal s (t);
by the sampling rate f s Sampling the response signal S (T), then performing FFT (fast Fourier transform) and modulus to obtain a frequency spectrum | S (N) |, the length is N, the resolution ratio delta f of the frequency spectrum is 1/T, and normalizing the frequency spectrum to obtain | S' (N) |.
|S(n)|=|X(n)|·|G(n)| (2)
Wherein | x (n) | is the frequency spectrum of chirp in the frequency domain; | G (n) | is the mode of harmonic oscillator amplitude-frequency characteristics; n is 1, 2, …, N;
the sweep frequency range in interception | S' (n) | is f 1 ≤f≤f 2 Inner spectrum | S 1 ′(n)|,f 1 Corresponding spectral line number N 1 、f 2 Corresponding spectral line number N 2 I.e. | S 1 ′(n)|=|S′(n)|,N 1 ≤n≤N 2 :
To | S 1 ' (n) | is subjected to spectral peak indexing, and the serial number of a spectral line corresponding to a spectral peak is k 0 Peak amplitude of | S 1 ′(k 0 ) Obtaining an estimate of the resonant frequency
Establishing a Lorentz curve model L (f):
If Q ═ m, drawing L (f) & gtnon-calculation Q=m Is normalized to obtain L' (f) luminance Q=m ;
The first Lorentzian fitting frequency f is in the rangeThe Q value range is more than or equal to 1 and less than or equal to M.
The frequency spectrum sequence corresponding to the frequency range is that the amplitude of each data point in the actual frequency range and the frequency point corresponding to the fitting curve are subtracted point by point and accumulated to obtain a residual error (m):
the point with the smallest value in the index residual error sequence error (i), with the sequence number marked as index0, obtains the first estimated quality factor value
The second Lorentz fitting f range is in the whole sweep frequency domain interval f 1 ≤f≤f 2 Selecting the range of Q value asRepeating the step 3.3 to obtain a residual error sequence error (i):
the point with the minimum value in the index residual error sequence error (i) is marked as index1 to obtain the final estimation value of the quality factor
Because the half-power bandwidth method only utilizes three points of the response signal frequency spectrum to calculate, and Lorentz fitting utilizes all the points in the effective frequency range, the error can be reduced for the signal with low signal-to-noise ratio; compared with the Lorentz fitting only once, the two-step fitting can save the operation time, especially for some harmonic oscillators with high Q values.
The technical means disclosed in the invention scheme are not limited to the technical means disclosed in the above embodiments, but also include the technical scheme formed by any combination of the above technical features.
Claims (4)
1. A harmonic oscillator quality factor measuring method based on chirp is characterized by comprising the following steps:
step 1: carrying out frequency sweeping excitation on the harmonic oscillator by using a chirp signal to obtain a response signal s (t);
step 2: sampling, FFT (fast Fourier transform), modulus taking and normalization are carried out on the response signal S (t) to obtain a frequency spectrum | S' (n) |; interception and selection sweep range spectrum | S 1 ' (n) | calculation obtains the estimated value of the vibration frequency
And step 3: establishing a Lorentz curve model, and obtaining a frequency range of the first Lorentz fitting by using a half-power bandwidth methodQ is more than or equal to 1 and less than or equal to M, and the quality factor is estimated for the first time; combining the first estimateAt swept frequency universe f 1 ≤f≤f 2 ,And carrying out second Lorentz fitting to obtain a final quality factor estimation value.
2. The method of claim 1, wherein the harmonic oscillator quality factor is measured based on a chirp signal,
setting parameters of the chirp signals in the step 1: starting frequency f 1 End frequency f 2 Sweep frequency time T, signal amplitude A, generating chirp signal x (T):
3. The method for measuring the quality factor of the harmonic oscillator based on the chirp according to claim 1, wherein the step 2 is to use a sampling rate f s Sampling the response signal S (T), then performing FFT (fast Fourier transform) and modulus to obtain a frequency spectrum | S (N) |, the length is N, the resolution ratio delta f of the frequency spectrum is 1/T, and normalizing the frequency spectrum to obtain | S' (N) |.
|S(n)|=|X(n)|·|G(n)| (2)
Where | x (n) | is the frequency spectrum of chirp in the frequency domain; | G (n) | is the mode of harmonic oscillator amplitude-frequency characteristics; n is 1, 2, …, N;
the sweep frequency range in the interception and selection | S' (n) | is f 1 ≤f≤f 2 Inner frequency spectrum | S 1 ′(n)|,f 1 Corresponding spectral line number N 1 、f 2 Corresponding spectral line number N 2 I.e. | S 1 ′(n)|=|S′(n)|,N 1 ≤n≤N 2 :
4. The method for measuring the quality factor of the harmonic oscillator based on the chirp according to claim 1, wherein the step 3 comprises:
step 3.1: establishing Lorentz curve model Is the estimated resonant frequency, Q is the quality factor, and the maximum value of Q is set to M; if Q ═ m, drawing L (f) & gtnon-calculation Q=m And normalizing to obtain L' (f) dichotomy Q=m ;
Step 3.2: the first Lorentzian fitting frequency f is in the rangeThe Q value range is more than or equal to 1 and less than or equal to M. Subtracting the amplitude of each data point in the actual frequency range from the amplitude of the frequency point corresponding to the fitting curve point by point and accumulating to obtain a residual error (m);
step 3.4: the point with the smallest value in the index residual error sequence error (i), with the sequence number marked as index0, obtains the first estimated quality factor value
Step 3.5: the second Lorentz fitting f range is in the whole sweep frequency domain interval, and the Q value range is selected asRepeating the step 3.3 to obtain a residual error sequence error (i);
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210380357.1A CN115060355B (en) | 2022-04-12 | 2022-04-12 | Harmonic oscillator quality factor measurement method based on linear frequency modulation pulse |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210380357.1A CN115060355B (en) | 2022-04-12 | 2022-04-12 | Harmonic oscillator quality factor measurement method based on linear frequency modulation pulse |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115060355A true CN115060355A (en) | 2022-09-16 |
CN115060355B CN115060355B (en) | 2024-03-26 |
Family
ID=83196876
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210380357.1A Active CN115060355B (en) | 2022-04-12 | 2022-04-12 | Harmonic oscillator quality factor measurement method based on linear frequency modulation pulse |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115060355B (en) |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5257544A (en) * | 1992-01-22 | 1993-11-02 | The Board Of Trustees Of The Leland Stanford Junior University | Resonant frequency method for bearing ball inspection |
US20040098215A1 (en) * | 2002-11-19 | 2004-05-20 | Kevin Gee | Method for the rapid estimation of figures of merit for multiple devices based on nonlinear modeling |
US20060251264A1 (en) * | 2003-02-27 | 2006-11-09 | Toa Corporation | Dip filter frequency characteristic decision method |
US20100324852A1 (en) * | 2009-06-23 | 2010-12-23 | Exxonmobil Research And Engineering Company | Determining the resonance parameters for mechanical oscillators |
CN202562949U (en) * | 2012-04-27 | 2012-11-28 | 南京信息工程大学 | Resonant type micro-accelerometer based on static rigidity |
US20150059482A1 (en) * | 2013-08-28 | 2015-03-05 | Robert Bosch Gmbh | Mems component |
EP2894521A1 (en) * | 2014-01-13 | 2015-07-15 | Ecole Polytechnique Federale de Lausanne (EPFL) | Isotropic harmonic oscillator and associated time base without escapement or simplified escapement |
CN105258786A (en) * | 2015-11-03 | 2016-01-20 | 中国科学院半导体研究所 | Rapid measurement of resonant frequency and quality factor of high-frequency harmonic oscillator |
CN108507556A (en) * | 2018-03-19 | 2018-09-07 | 中国人民解放军国防科技大学 | Method and device for correcting and adjusting non-uniform quality factors of gyro harmonic oscillators with cylindrical shells |
CN109490853A (en) * | 2017-09-10 | 2019-03-19 | 北京遥感设备研究所 | Spectral line value determines method at a kind of chirp pulse signal centre frequency |
CN109490852A (en) * | 2017-09-10 | 2019-03-19 | 北京遥感设备研究所 | A kind of chirp pulse signal chirp rate polarity determination method |
CN110823530A (en) * | 2019-11-13 | 2020-02-21 | 南京大学 | Method for obtaining quality factor of micro-resonant cavity |
CN111272193A (en) * | 2020-02-17 | 2020-06-12 | 东南大学 | MEMS gyroscope online frequency difference identification method based on noise power spectrum estimation |
CN112559956A (en) * | 2020-12-23 | 2021-03-26 | 江苏睿科大器机器人有限公司 | Wavelet threshold adaptive shrinkage method, system, electronic device and storage medium |
CN113218547A (en) * | 2020-02-06 | 2021-08-06 | 天津大学 | Structure of high-Q perforated flexible micro-ring resonant cavity for double sensing |
-
2022
- 2022-04-12 CN CN202210380357.1A patent/CN115060355B/en active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5257544A (en) * | 1992-01-22 | 1993-11-02 | The Board Of Trustees Of The Leland Stanford Junior University | Resonant frequency method for bearing ball inspection |
US20040098215A1 (en) * | 2002-11-19 | 2004-05-20 | Kevin Gee | Method for the rapid estimation of figures of merit for multiple devices based on nonlinear modeling |
US20060251264A1 (en) * | 2003-02-27 | 2006-11-09 | Toa Corporation | Dip filter frequency characteristic decision method |
US20100324852A1 (en) * | 2009-06-23 | 2010-12-23 | Exxonmobil Research And Engineering Company | Determining the resonance parameters for mechanical oscillators |
CN202562949U (en) * | 2012-04-27 | 2012-11-28 | 南京信息工程大学 | Resonant type micro-accelerometer based on static rigidity |
US20150059482A1 (en) * | 2013-08-28 | 2015-03-05 | Robert Bosch Gmbh | Mems component |
EP2894521A1 (en) * | 2014-01-13 | 2015-07-15 | Ecole Polytechnique Federale de Lausanne (EPFL) | Isotropic harmonic oscillator and associated time base without escapement or simplified escapement |
CN105258786A (en) * | 2015-11-03 | 2016-01-20 | 中国科学院半导体研究所 | Rapid measurement of resonant frequency and quality factor of high-frequency harmonic oscillator |
CN109490853A (en) * | 2017-09-10 | 2019-03-19 | 北京遥感设备研究所 | Spectral line value determines method at a kind of chirp pulse signal centre frequency |
CN109490852A (en) * | 2017-09-10 | 2019-03-19 | 北京遥感设备研究所 | A kind of chirp pulse signal chirp rate polarity determination method |
CN108507556A (en) * | 2018-03-19 | 2018-09-07 | 中国人民解放军国防科技大学 | Method and device for correcting and adjusting non-uniform quality factors of gyro harmonic oscillators with cylindrical shells |
CN110823530A (en) * | 2019-11-13 | 2020-02-21 | 南京大学 | Method for obtaining quality factor of micro-resonant cavity |
CN113218547A (en) * | 2020-02-06 | 2021-08-06 | 天津大学 | Structure of high-Q perforated flexible micro-ring resonant cavity for double sensing |
CN111272193A (en) * | 2020-02-17 | 2020-06-12 | 东南大学 | MEMS gyroscope online frequency difference identification method based on noise power spectrum estimation |
CN112559956A (en) * | 2020-12-23 | 2021-03-26 | 江苏睿科大器机器人有限公司 | Wavelet threshold adaptive shrinkage method, system, electronic device and storage medium |
Non-Patent Citations (6)
Title |
---|
BOGDAN SAPIŃSKI: "Electrical harmonic oscillator with MR damper and energy harvester operating as TMD: Experimental study", MECHATRONICS * |
CHENG, HJ;: "Study on the Output of LD-Pumped Passively Q-switched Subnanosecond Microlasers", INTERNATIONAL SYMPOSIUM ON PHOTOELECTRONIC DETECTION AND IMAGING 2013: HIGH POWER LASERS AND APPLICATIONS * |
LIN, TY;: "Analyzing OAM mode purity in optical fibers with CNN-based deep learning", CHINESE OPTICS LETTERS * |
吕正: "一种高Q值硅微陀螺谐振频率快速锁定方法", 测控技术 * |
马雪琳: "基于多普勒测振的谐振子关键参数测量方法研究", 中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑) * |
高晓峰: "光谱分辨率增强方法品质因子研究", 光子学报 * |
Also Published As
Publication number | Publication date |
---|---|
CN115060355B (en) | 2024-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6512788B1 (en) | RF output spectrum measurement analyzer and method | |
US10998926B2 (en) | Method for compensating gain flatness of transceiver | |
US7161511B2 (en) | Linearization system and method | |
KR101294771B1 (en) | Filter equalization using magnitude measurement data | |
CN112462356B (en) | Method for improving radar ranging accuracy based on Gaussian interpolation | |
EP3586128B1 (en) | Probing a structure of concrete by means of electromagnetic waves | |
CN111693136B (en) | Acoustic surface wave resonator frequency estimation algorithm adopting echo signal autocorrelation phase spectrum | |
CN109682492B (en) | Frequency estimation method based on frequency domain Gaussian fitting | |
JP2007509328A (en) | Method for determining humidity and density of dielectric materials | |
CN114374446B (en) | Amplitude-phase characteristic measurement method based on pulse signals | |
CN115060355A (en) | Harmonic oscillator quality factor measuring method based on linear frequency modulation pulse | |
CN109655775A (en) | A kind of multiple dimensioned calibration method of arbitrary waveform generator amplitude frequency sweep and device | |
CN106483445B (en) | A kind of built-in measurement method and device of wideband circuit phase nonlinear distortion | |
Verspecht et al. | The vector component analyzer: A new way to characterize distortions of modulated signals in high-frequency active devices | |
CN113406386B (en) | Signal frequency accurate estimation method based on digital down-conversion | |
US8023558B2 (en) | Method and apparatus for measuring the input frequency response of a digital receiver | |
JPH05273335A (en) | Method and device for transmitting and receiving pulsed wave | |
CN112883787B (en) | Short sample low-frequency sinusoidal signal parameter estimation method based on spectrum matching | |
CN107872283B (en) | Frequency sweeping method and device | |
CN114966198A (en) | Sinusoidal signal-to-noise ratio measuring method | |
Wang et al. | Radar emitter signal intra-pulse feature extraction based on empirical mode decomposition | |
Keller et al. | A new method of emission measurement | |
CN110687362A (en) | Real-time phase noise testing method | |
CN100473175C (en) | Method for determining envelope lines of modulation signals | |
JP2009500593A (en) | Delay calculation apparatus and method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |