The disclosure herein belongs to the field of high-frequency signal test, and more particularly relates to a Fourier analysis method with a variable sampling frequency.
The most common Fourier analysis is discrete Fourier transform (DFT) and fast Fourier transform (FFT). The DFT is an important harmonic analysis tool, can perform mathematical transformation on a sampling sequence of complex signals to separate fundamental signals from each harmonic signal. Generally, in order to ensure a frequency resolution, a sampling sequence length N needs to be increased. When an N value is greater, the DFT needs N2-time complex multiplication operation, and the required time is too long. The requirements on hardware may be greatly increased to ensure good real-time performance.
The FFT algorithm adopts a butterfly operation mode, and can realize harmonic detection in a short time, but the frequency resolution is low, and synchronous sampling and integer period truncation are required. If accurate synchronous sampling can be guaranteed, the available measurement accuracy of the FFT on harmonic waves is very high. However, the measurement of inter-harmonics depends on the frequency resolution. Generally, window lengths and sampling frequencies set by commercial power analyzers in compliance with IEC standards can only meet the frequency resolution in a range of 1 to 10 Hz, but it is often insufficient for inter-harmonic measurement of non-integer frequencies. Under the condition of non-synchronous sampling, errors of amplitude value and frequency measurement may be greatly increased due to inherent spectrum leakage and fence effect of the Fourier method, but it is difficult to realize strict synchronous sampling in practical engineering application, so how to reduce the spectrum leakage and fence effect is a research focus of scholars worldwide.
Analysis results of the Fourier method are greatly influenced by the spectrum leakage and fence effect, and they are complementary to each other. Only when a measured frequency component just coincides with a frequency axis unit, accurate analysis results can be obtained. Generally, existing harmonic analysis instruments and analysis methods may realize the precondition by using dual limitations of synchronous sampling and fence effect in the low-frequency field, but in fact, a strong constraint is added to the sampling condition, and an effective action space of the sampling frequency is reduced, so that the analysis in the high-frequency field may be limited by the sampling frequency.
In order to make the measured frequency coincide with a frequency axis unit point as much as possible, the common harmonic analysis algorithm are required to meet two conditions of synchronous sampling and integer period truncation. For the strict synchronous sampling, the sampling frequency is integer times of all frequency components. Otherwise, once the sampling frequency forms non-synchronous sampling of a certain harmonic, all spectral line results of Fourier analysis may be influenced. In practical engineering, the signal contains harmonics and inter-harmonics, the kinds of frequency components are various, and the frequency of each component is unknown, so it is difficult to achieve strict synchronous sampling. At the same time, integer period truncation needs wave filtering by instruments at an earlier stage to determine the fundamental period so as to calculate the truncation time window length. Under the condition that PWM is used or noise exists, large amplitude vibration at a zero crossing point of a waveform may cause measurement inaccuracy of the fundamental period. If a filter is added and used, the amplitude value of a main component may be reduced to a certain extent
A necessary and sufficient condition that no spectrum leakage occurs is that the measured frequency coincides with the frequency unit, and synchronous sampling and integer period truncation are derived conditions for realizing the precondition. For high-speed motors with operating frequencies up to hundreds or even thousands of Hertz, a limit of the existing hardware sampling frequency may be exceeded if a least common multiple of each harmonic number is solved to achieve synchronous sampling.
Therefore, by aiming at this condition, the disclosure provides a novel Fourier decomposition method with a variable sampling frequency, so as to improve the precision of the Fourier decomposition method, and effectively reduce the spectrum leakage and fence effect on high-frequency signals, especially harmonic signals.
In order to solve the problems of analysis acquisition and analysis of high-frequency signals of a high-speed motor, the disclosure provides a Fourier analysis method with a variable sampling frequency so as to improve the precision of a Fourier decomposition method and effectively reduce the spectrum leakage and fence effect on high-frequency signals, especially harmonic signals.
The disclosure is realized by the following technical solution: a Fourier analysis method with a variable sampling frequency includes the following steps:
S100 preliminarily sampling a to-be-analyzed signal by comparing an initially set sampling frequency, and further analyzing its fundamental frequency and fundamental amplitude value by Fourier analysis;
S200 preliminarily judging the fundamental amplitude value obtained through analysis to determine a sampling frequency meeting integer period truncation, and resampling the signal;
S300 performing Fourier analysis on the sampled signal, and determining an optimum sampling frequency, i.e., an optimum sampling frequency both meeting the integer period truncation and meeting no spectrum leakage, by a three-spectral line method; and
S400 resampling the signal to obtain a frequency and amplitude value composition of each harmonic.
Further, in step S100, specifically, an operating frequency f is directly calculated through a formula f=np/60, and sampling is performed by using 2n times of a fundamental frequency estimated value {circumflex over (f)}1 as the initial sampling frequency.
Further, before step S100, the following step is further included:
step S000 obtaining a rotating speed signal fed back by a tested high-speed motor.
Further, the optimum sampling frequency has a plurality of values and is in periodic change.
Further, for the period change of the optimum sampling frequency, a period is gradually increased along with frequency increase.
The disclosure has the beneficial effects that the disclosure designs the Fourier analysis method with the variable sampling frequency. The initial sampling frequency can be fast estimated according to the rotating speed signal. Then, according to a calculation algorithm, the optimum sampling frequency can be fast determined. Acquisition errors and analysis errors of the high-frequency signals can be effectively reduced. The fence effect and the spectrum leakage can be reduced to 0. The analysis precision of the Fourier analysis algorithm can be effectively improved. At the same time, by using this algorithm, high-precision and high-frequency signal acquisition and analysis can be realized by using a smaller sampling frequency, the cost of an acquisition system is reduced, and an analysis speed is accelerated.
The technical solutions in embodiments of the disclosure will be described clearly and completely hereinafter in conjunction with the accompanying drawings in the embodiments of the disclosure, and obviously, the described embodiments are only a part of the embodiments of the disclosure, but not all of them. Based on the embodiments of the disclosure, all other embodiments obtained by those of ordinary skill in the art without creative labor are all within the protection scope of the disclosure.
As shown in
S100 A to-be-analyzed signal is preliminarily sampled by comparing an initially set sampling frequency, and its fundamental frequency and fundamental amplitude value are further analyzed by Fourier analysis.
S200 The fundamental amplitude value obtained through analysis is preliminarily judged to determine a sampling frequency meeting integer period truncation, and the signal is resampled.
S300 Fourier analysis is performed on the sampled signal, and an optimum sampling frequency, i.e., an optimum sampling frequency both meeting the integer period truncation and meeting no spectrum leakage, is determined by a three-spectral line method.
S400 The signal is resampled to obtain a frequency and amplitude value composition of each harmonic.
In partial preferred embodiments, in step S100, specifically, an operating frequency f is directly calculated through a formula f=np/60, and sampling is performed by using 2n times of a fundamental frequency estimated value {circumflex over (f)}1 as the initial sampling frequency.
In partial preferred embodiments, before step S100, the following step is further included:
Step S000 A rotating speed signal fed back by a tested high-speed motor is obtained.
In partial preferred embodiments, in step S300, the sampling frequency is subjected to optimization regulation according to a harmonic frequency and harmonic amplitude value object to be discriminated so as to determine the optimum sampling frequency.
In partial preferred embodiments, the optimum sampling frequency has a plurality of values and is in periodic change.
In partial preferred embodiments, for the period change of the optimum sampling frequency, a period is gradually increased along with frequency increase.
Specifically, referring to
In the formula, mmax is a highest number of times of the harmonic required to be analyzed.
If {circumflex over (f)}1=f1, the practical fundamental frequency should correspond to a (2p−n)th spectral line, spectral line, at the moment, the frequency of the m-time harmonic corresponds to a (m×2p−n)th spectral line, and the condition of coinciding with a frequency unit is met. However, because the estimated value fed back by the rotating speed has errors, a spectrum leakage result certainly occurs. In order to reduce the spectrum leakage to the maximum degree, the method obtains different signal sample sequences by changing the sampling frequency, performs Fourier analysis by using these sample sequences, and seeks a maximum value of amplitude value results. When the amplitude value reaches the maximum value, the spectrum leakage is almost totally eliminated. The continuously corrected sampling frequency gradually approaches to the optimum sampling frequency fsop=f1×2p−n.
Illustration is made in conjunction with
In the formula, [ ] is an integer and can be rounded. The amplitude value A*(k) of a corresponding spectral line is solved. Its calculation method is shown as follows:
However, a difference between {circumflex over (f)}1 and the practical fundamental frequency f1 may cause errors in results of a rounding function, so that an order difference of the corresponding spectral line positions is 1, i.e. the calculated spectral lines are spectral lines at two sides of a main lobe. In order to avoid such errors, the amplitude values of the left and right adjacent spectral lines of a kth spectral line are usually calculated, and a maximum value of the three values is found to be used as the estimated value of the fundamental frequency amplitude value. If the sampling frequency is changed again in a subsequent step, the calculated spectral line positions do not necessarily have errors. Therefore, the amplitude values of the three spectral lines do not need to be calculated each time. For this purpose, a start threshold for the three-spectral line analysis is set:
In the formula, At is the start threshold. When the amplitude value of the spectral line obtained through calculation is smaller than At, its sidelobes and main lobe are almost identical. It shows that the spectral line is certainly not the main lobe, and the main lobe needs to be found through other calculation. The frequency unit corresponding to the main lobe is a newly obtained fundamental frequency analysis result:
According to the newly obtained fundamental frequency value f*1, a new sampling frequency 2″×f*1 is defined according to the previous setting.
Referring to
Referring to
For the condition in
In the formula, f**1 is a corrected fundamental frequency value. A*(k), A*(k−1) and A*(k+1) are respectively amplitude values of spectral lines of the main lobe and the left and right side sidelobes. λ and λr are respectively correction coefficients when the sampling frequency is leftwards and rightwards corrected.
Referring to
The bisection method is a fast method suitable for searching in a large data volume interval. By the principle of a bisection region, the region can be reduced at an exponential speed. As shown in
A rounding function in the formula needs positive rounding. The sampling frequency after twice correction is very close to the optimum sampling frequency, so that the searching interval of the bisection method can be greatly compressed to reduce the number of searching times. As mentioned above, the corrected 2″×f**1 is certainly distributed between fsc and fsd, so that the searching interval can be calculated by subtracting fsc by fsd. However, the searching interval is too large, and a specific difference value between 2″×f**1 and the optimum sampling frequency fsop is unknown, it is not suitable to reduce the searching interval blindly, and the searching interval can be further reduced by using the principle of correcting the fundamental frequency by three spectral lines. A new sequence obtained through sampling by using 2″×f**1 is analyzed by using the three-spectral line method again. Identically, according to the above steps, a sampling frequency correction direction is judged according to amplitude values of the sidelobes. Then, a correction coefficient 2 of the fundamental frequency is solved by using an amplitude value proportion of the main lobe and the sidelobes:
In the formula, A**(k), A**(k−1) and A**(k−1) are amplitude value results of the main lobe and sidelobes of three-spectral line analysis performed again. In order to ensure that the optimum sampling frequency fsop is in the interval, a correction amount of the fundamental frequency is increased to twice.
End points of another corresponding interval are as follows:
After the searching interval is determined, one half of a sum of sampling frequencies of two end points is taken according to the principle of the bisection method to be used as a new sampling frequency. An amplitude value Aop of the main lobe of the new signal sequence based on the new sampling frequency is calculated by using the three-spectral line method. A different value between Aop and the amplitude value of the main lobe corresponding to the sampling frequencies of the two end points is as follows:
According to sizes of Δa and Δb, the end points of the interval can be determined and updated, and value reassignment is performed:
A calculation precision is set to be 10−4 grade. When the difference value Δab=|Aa−Ab| between the amplitude value results of the main lobe corresponding to the sampling frequencies of the two end points of a new interval is smaller than the precision, the searching of the bisection method can be considered completed.
Number | Date | Country | Kind |
---|---|---|---|
2019103496139 | Apr 2019 | CN | national |
Number | Date | Country | |
---|---|---|---|
Parent | PCT/CN2020/076318 | Feb 2020 | US |
Child | 17152844 | US |