CN107766293A - Signal spectral analysis method and system when fractional-sample data rule lacks - Google Patents
Signal spectral analysis method and system when fractional-sample data rule lacks Download PDFInfo
- Publication number
- CN107766293A CN107766293A CN201710900980.4A CN201710900980A CN107766293A CN 107766293 A CN107766293 A CN 107766293A CN 201710900980 A CN201710900980 A CN 201710900980A CN 107766293 A CN107766293 A CN 107766293A
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtr
- mtd
- msub
- msup
- 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
Classifications
-
- 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
-
- 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
- G06F17/141—Discrete Fourier transforms
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Discrete Mathematics (AREA)
- Complex Calculations (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
Disclose signal spectral analysis method and system during a kind of fractional-sample data rule missing.This method can include:For the measurement signal of fractional-sample data rule missing, using the Discrete Fourier transform operation of block sampling mode, the intermediate result of acquisition measurement signal frequency spectrum;The rule presented from orthogonality of trigonometric function in the case of missing data, influence matrix corresponding to calculating, and then obtain multiple degree reduction matrixes;Calculate the inverse matrix of degree reduction matrix;Based on inverse matrix, intermediate result is converted into final frequency spectrum data.Signal spectral analysis method is expanded to new scope by the present invention, has relatively low operand and the higher degree of accuracy.
Description
Technical field
The present invention relates to electromagnetic measurement, Fundamental Theory of Electrical Engineering, field of signal processing, more particularly, to a kind of fractional-sample number
Signal spectral analysis method and system when being lacked according to systematicness.
Background technology
Classical discrete fourier transform method (DFT) has a wide range of applications.The sampled data of measurement signal is carried out from
Fourier transform is dissipated, then can obtain the frequency spectrum of signal, i.e., the cosine coefficient and sinusoidal coefficients of each frequency component.Generally sampling is
Interval, the condition required for DFT is whole samplings in a cycle.But for various reasons, sampling
Data often produce missing.Wherein important one kind, belong to the missing of time rule.This missing is probably derived from artificial
Processing, such as delete some data intentionally to overcome the influence of transient process in ladder wave measurement, deposited in these data
In transient process, it is impossible to accurately reflect original signal, actual conditions can be preferably reflected after deletion;It may also derive from non-
Artificial process, such as loss of the signal during (or communication) is transmitted, this lose is often random, it is likely that depositing
In the situation of systematicness missing, or by processing, turn into the situation that systematicness lacks.
Due to shortage of data, it is impossible to meet the requirement of common classical discrete fourier transform method, tend not to carry out
The calculating of cosine coefficient and sinusoidal coefficients.A kind of therefore, it is necessary to signal frequency when developing fractional-sample data rule missing
Spectral analysis method and system.
The information for being disclosed in background of invention part is merely intended to deepen the reason of the general background technology to the present invention
Solution, and be not construed as recognizing or imply known to those skilled in the art existing of the information structure in any form
Technology.
The content of the invention
The present invention proposes signal spectral analysis method and system during a kind of fractional-sample data rule missing, will believe
Number spectrum analysis is expanded to new scope, has relatively low operand and the higher degree of accuracy.
A kind of according to an aspect of the invention, it is proposed that signal spectral analysis side during fractional-sample data rule missing
Method.Methods described can include:For fractional-sample data rule missing measurement signal, using block sampling mode from
Fourier transformation computation is dissipated, obtains the intermediate result of measurement signal frequency spectrum;From orthogonality of trigonometric function in missing data situation
The rule of lower presentation, influence matrix corresponding to calculating, and then obtain multiple degree reduction matrixes;Calculate the inverse matrix of degree reduction matrix;Base
In inverse matrix, intermediate result is converted into final frequency spectrum data.
Preferably, the measurement signal of fractional-sample data rule missing is:Appearance lacks and lacks quantity phase at equal intervals
Same measurement signal.
Preferably, the Discrete Fourier transform operation of the block sampling mode is:
Wherein,For intermediate result, aiAnd biThe respectively cosine coefficient and sinusoidal coefficients of ith harmonic wave, yn,mFor
The sampled data that measurement signal obtains, h are the sampling interval, and N is the sum of measurement signal even group-division in a cycle, and n is
The deck label of measurement signal, m are the sampled data label of every group of measurement signal equal intervals, and [0, s-1] and [t+1, M-1] are scarce
Data label range is lost, [s, t] is retention data label range, and M is the complete of the sampled data of every group of measurement signal equal intervals
Portion's quantity.
Preferably, influence matrix is:
Wherein, aa (i, j) is the element of cosine and cosine interaction matrix AA, and bb (i, j) is sinusoidal mutual with sine
Acting matrix BB element, ab (i, j) are the elements of cosine and sinusoidal interaction matrix AB.
Preferably, as t+s=M, influence matrix AB all elements are zero, and matrix A A is decomposed into three class depression of order squares
Battle array:Ak, A0.5N, and AN, influence matrix BB is decomposed into three class degree reduction matrixes:Bk, B0.5N, and BN。
Preferably, the cosine result of the final frequency spectrum data is:
Wherein, akIt is that frequency is kf1Frequency spectrum final result cosine coefficient, (Ak)-1,(A0.5N)-1, (AN)-1It is matrix A A
Three price reduction inverse of a matrix matrixes.
Preferably, the sinusoidal result of the final frequency spectrum data is:
Wherein, bkIt is that frequency is kf1Frequency spectrum final result sinusoidal coefficients, (Bk)-1,(B0.5N)-1, (BN)-1It is matrix B B
Three price reduction inverse of a matrix matrixes.
Preferably, signal spectral analysis method when fractional-sample data rule lacks is applied to the quantum of spectrum analysis
In measurement.
A kind of according to another aspect of the invention, it is proposed that signal spectral analysis during fractional-sample data rule missing
System, it can include:Memory, it is stored with computer executable instructions, and the signal sample data of excalation;Processing
Device, the computer executable instructions in processor run memory, perform following steps:Lacked for fractional-sample data rule
The measurement signal of mistake, classical Discrete Fourier transform operation is carried out using block sampling mode, obtained in measurement signal frequency spectrum
Between result;The rule presented from orthogonality of trigonometric function in the case of missing data, influence matrix corresponding to calculating, and then obtain
Obtain multiple degree reduction matrixes;Calculate the inverse matrix of degree reduction matrix;Based on inverse matrix, intermediate result is converted into final frequency spectrum data.
Preferably, the measurement signal of fractional-sample data rule missing is:Appearance lacks and lacks quantity phase at equal intervals
Same measurement signal.
Methods and apparatus of the present invention has other characteristics and advantage, and these characteristics and advantage are attached from what is be incorporated herein
It will be apparent in figure and subsequent embodiment, or by the accompanying drawing being incorporated herein and subsequent specific reality
Apply in mode and stated in detail, these the drawings and specific embodiments are provided commonly for explaining the certain principles of the present invention.
Brief description of the drawings
Exemplary embodiment of the present is described in more detail in conjunction with the accompanying drawings, of the invention is above-mentioned and other
Purpose, feature and advantage will be apparent, wherein, in exemplary embodiments of the present invention, identical reference number is usual
Represent same parts.
Fig. 1 shows the step of signal spectral analysis method when being lacked according to the fractional-sample data rule of the present invention
Flow chart.
Fig. 2 shows the step of signal spectral analysis method when being lacked according to the fractional-sample data rule of the present invention
Explanatory diagram.
Fig. 3 a, Fig. 3 b and Fig. 3 c respectively illustrate measurement signal curve according to an embodiment of the invention, all adopted
Sample data, the schematic diagram for lacking post-sampling data.
Fig. 4 shows the schematic diagram of staircase waveform aliasing effect according to an embodiment of the invention.
Fig. 5 shows the difference measurement of the non-sine periodic signal according to an embodiment of the invention based on staircase waveform
In, the alignment schematic diagram of two signals.
Embodiment
The present invention is more fully described below with reference to accompanying drawings.Although showing the preferred embodiments of the present invention in accompanying drawing,
However, it is to be appreciated that the present invention is may be realized in various forms without should be limited by embodiments set forth here.Conversely, there is provided
These embodiments are in order that the present invention is more thorough and complete, and can will fully convey the scope of the invention to ability
The technical staff in domain.
Fig. 1 shows the step of signal spectral analysis method when being lacked according to the fractional-sample data rule of the present invention
Flow chart.
Fig. 2 shows the step of signal spectral analysis method when being lacked according to the fractional-sample data rule of the present invention
Explanatory diagram.
In this embodiment, signal spectral analysis method when being lacked according to the fractional-sample data rule of the present invention can
With including:
Step 101, for the measurement signal of fractional-sample data rule missing, using discrete Fu of block sampling mode
Vertical leaf transformation computing, obtain the intermediate result of measurement signal frequency spectrum;
Step 102, the rule presented from orthogonality of trigonometric function in the case of missing data, corresponding influence square is calculated
Battle array, and then obtain multiple degree reduction matrixes;
Step 103, the inverse matrix of degree reduction matrix is calculated;
Step 104, based on inverse matrix, intermediate result is converted into final frequency spectrum data.
In one example, the measurement signal of fractional-sample data rule missing is:Appearance is lacked and lacked at equal intervals
Quantity identical measurement signal.
In one example, the Discrete Fourier transform operation of block sampling mode is:
Wherein,For the PRELIMINARY RESULTS of signal spectrum, aiAnd biThe respectively cosine coefficient and sine of ith harmonic wave
Coefficient, yn,mThe sampled data obtained for measurement signal, h are the sampling interval, and N is measurement signal even group-division in a cycle
Sum, n is the deck label of measurement signal, and m is the sampled data label of every group of measurement signal equal intervals, [0, s-1] and [t+
1, M-1] it is missing data label range, [s, t] is retention data label range, and M is adopting for every group of measurement signal equal intervals
The entire quantity of sample data.
In one example, influence matrix is:
Wherein, aa (i, j) is the element of cosine and cosine interaction matrix AA, and bb (i, j) is sinusoidal mutual with sine
Acting matrix BB element, ab (i, j) are the elements of cosine and sinusoidal interaction matrix AB.
In one example, as t+s=M, influence matrix AB all elements are zero, and matrix A A is decomposed into three classes drop
Rank matrix:Ak, A0.5N, and AN, influence matrix BB is decomposed into three class degree reduction matrixes:Bk, B0.5N, and BN。
In one example, the cosine result of final frequency spectrum data is:
Wherein, akIt is that frequency is kf1Frequency spectrum final result cosine coefficient, (Ak)-1,(A0.5N)-1, (AN)-1It is matrix A A
Three price reduction inverse of a matrix matrixes.
In one example, the sinusoidal result of final frequency spectrum data is:
Wherein, bkIt is that frequency is kf1Frequency spectrum final result sinusoidal coefficients, (Bk)-1,(B0.5N)-1, (BN)-1It is matrix B B
Three price reduction inverse of a matrix matrixes.
Specifically, the measurement signal based on fractional-sample data rule missing, using discrete Fu of block sampling mode
Vertical leaf transformation computing, obtain the intermediate result of the frequency spectrum of measurement signal;The systematicness missing of this data is probably derived from artificial
Processing, such as delete some data intentionally to overcome the influence of transient process in ladder wave measurement, deposited in these data
In transient process, it is impossible to accurately reflect original signal, actual conditions can be preferably reflected after deletion;It may also derive from non-
Artificial process, such as loss of the signal during (or communication) is transmitted, although this lose is often random, having can
The systematicness missing of description can be present, or by converting, turn into the situation of systematicness.Fractional-sample data rule
Missing can be further depicted as, and the sampled data of a cycle is divided into N groups, deck label n=0,1,2 ..., N-1;Every group of tool
There are M equally spaced sampled datas, every group of data are marked as m=0,1,2 ..., M-1;Identical missing is produced in every group of data,
The data of missing are marked as m=0,1,2 ..., s-1 and m=t+1, t+2 ..., and M-1, remaining data label is from s to t, i.e. m=
S, s+1 ..., t, shared (t-s+1) is individual, typically there is 0 < s < 0.5M < t < M-1.In order to describe conveniently, it is referred to as segmentation and adopts
Sample, (s, t-s+1, M-1-t) is designated as, includes DFT and only middle remaining (t-s+1) individual sampled data is entered
Capable implication.In the foregoing description, missing in the case of missing is not or not both sides, is not difficult to pass through positioned at the both sides of every group of data
The starting point in cycle is reselected, the both sides that missing is arranged into one group of data are gone, therefore, fractional-sample data rule missing
Measurement signal there is wide applicability.
Measurement signal is:
Wherein, y is measurement signal, aiAnd biThe respectively cosine coefficient and sinusoidal coefficients of ith harmonic wave, f1It is fundamental wave frequency
Rate, ω are circular frequency, use block sampling mode to carry out classical Discrete Fourier transform operation as formula (2), wherein, between sampling
Formula (6) is divided into, and then obtains intermediate result,
H=2 π/(NM) (6).
The rule presented from orthogonality of trigonometric function in the case of missing data, influence matrix corresponding to calculating are formula
(3).As t+s=M, influence matrix AB all elements are zero;As t+s=M, influence matrix AA is decomposed into three class depression of orders
Matrix:Ak, A0.5N, and AN;As t+s=M, influence matrix BB is decomposed into three class degree reduction matrixes:Bk, B0.5N, and BN。
Fig. 4 shows the schematic diagram of staircase waveform aliasing effect according to an embodiment of the invention.
Wherein, degree reduction matrix Ak(BkEqually) be original matrix AA (BB) a part, be that line label is k, N- in original matrix
K, N+k, 2N-k, 2N+k ... and row label is also k, N-k, N+k, 2N-k, 2N+k ... element composition, wherein final line
Label (row label) is by requiring that the highest order of analysis determines.If last line label (row label) is designated as wN-k (or wN+
K), 2 ω≤(t-s+1) < M should typically be met.Degree reduction matrix Ak(Bk) diagonal element be 1.Work as k=1,2 ..., 0.During 5N-1,
All Ak(BkEqually) there is identical form, in Fig. 4, these matrixes correspond to set out from k=1,2 ..., 0.5N-1 points
Dotted line.
Wherein, degree reduction matrix A0.5N(B0.5NEqually) be original matrix AA (BB) a part, be that line label is in original matrix
0.5N, 1.5N ... and row label is also 0.5N, 1.5N ... element composition, the same depression of order of requirement of final line label (row label)
Matrix Ak, in Fig. 4, this matrix corresponds to from point 0.5N, positioned at the line (not shown) of rightmost.
Wherein, degree reduction matrix AN(BNEqually) be original matrix AA (BB) a part, be that line label is N in original matrix,
2N ... and row label is also N, 2N ... element composition, the requirement of final line label (row label) is the same as degree reduction matrix Ak, in Fig. 4
In, this matrix corresponded to from 0 (the N)s, positioned at leftmost line (not shown) of point.
Respective inverse matrix is calculated by each degree reduction matrix:(Ak)-1, (A0.5N)-1, and (AN)-1;(Bk)-1, (B0.5N)-1, and
(BN)-1, those skilled in the art can calculate inverse matrix as the case may be.
Based on inverse matrix, intermediate result is converted into final frequency spectrum data.By taking cosine coefficient as an example, first by intermediate result
Press relation of plane and be divided into three class column vectors:And then obtain the cosine knot of final frequency spectrum data
Fruit is formula (4), and sinusoidal result is formula (5).
Signal spectral analysis method is expanded to new scope by this method, has relatively low operand and higher accurate
Degree.
For ease of understanding the scheme of the embodiment of the present invention and its effect, two concrete application examples given below.This area
It should be understood to the one skilled in the art that the example, only for the purposes of understanding the present invention, its any detail is not intended to be limited in any way
The system present invention.
Using example 1
Fig. 3 a, Fig. 3 b and Fig. 3 c respectively illustrate measurement signal curve according to an embodiment of the invention, all adopted
Sample data, the schematic diagram for lacking post-sampling data.
Table 1 gives the arranges value of a primary signal, including the amplitude of each frequency component and phase angle.The shape of its curve
Shape is as shown in Figure 3 a.From Fig. 3 a, for curve within ± 1V, this is easy to subsequent calculations error, and obtained absolute error is right
For 1V, therefore relative error can also be considered as.
Table 1
Primary signal arranges value | 1 | 3 | 4 | 15 | 24 |
Amplitude | 0.800 | -0.320 | 0.200 | -0.045 | 0.023 |
Phase angle | -π/5 | π/3 | 0 | -π/2.3 | π/10 |
This signal is done into uniform sampling, a cycle branch 960, is divided into 20 groups, every group of 48 samplings;Open in each group
12 to begin, and last 11 data are deleted and not had to, and middle 25 data retain, and do subsequent arithmetic.Each cycle is participated in
The sampling of computing has 500.Therefore 250 subharmonic can be at most analyzed, agreement highest harmonic wave is 60 times.
N=20, M=48, block sampling are designated as (12,25,11), s=12, t=36, meet s+t=M condition.According to
Formula (2) carries out being segmented DFT computings to sampling.The intermediate result of signal spectrum is obtained, as shown in table 2.
Fig. 4 shows the schematic diagram of staircase waveform aliasing effect according to an embodiment of the invention.When the frequency of analysis
During less than 0.5N, classical Fourier transformation can obtain accurate result without difficulty, when the frequency of analysis is more than or equal to
0.5N and when being less than 0.5N (t-s+1), should use the present invention to obtain accurate result.
Then influence matrix is calculated according to formula (3), they are:
AkIt is 6 × 6 matrixes, row (column) is marked as (k, 20-k, 20+k, 40-k, 40+k, 60-k), matrix
A10It is 3 × 3 matrixes, row (column) is marked as (10,30,50), matrix
A20It is 3 × 3 matrixes, row (column) is marked as (20,40,60), matrix
BkFor
B10For
B20For
Their inverse matrix is respectively:
(Ak)-1For
(A10)-1For
(A20)-1For
(Bk)-1For
(B10)-1For
(B20)-1For
Finally, using these inverse matrixs, recovery processing is done to the intermediate result obtained in the first step, obtains final spectrum number
According to as shown in table 2.Two row errors of rightmost are the differences of final frequency spectrum data and original arranges value in table.As before, due to measurement
Signal maximum amplitude is close to 1V, so the result can be considered as relative value.
Table 2
As seen from Table 2, the deviation of intermediate result and original arranges value is larger, and this is that block sampling carries out classical direct computation of DFT
The result of leaf transformation;But pass through after inverse matrix recovery, the difference of final result and arranges value is 10-9Magnitude, substantially count
Calculation error caused by the limited setting of the numerical digit of calculation machine in itself.(data are given to after decimal point 7 in table, but actual calculate does not have
Restricted numerical digit, so 10 can be obtained-9The error of magnitude).
From the principle, influence matrix of the invention is actually relevant with orthogonality of trigonometric function, block sampling
Deviation caused by classical discrete Fourier transform is exactly as caused by these influence matrixs.Its inverse matrix is then extensive by baseline results
Come again.Due to the systematicness of shortage of data in time, therefore as the basic function being related in classical discrete Fourier transform,
I.e. trigonometric function can directly write out, and in other words, this measurement pattern can design in advance, or determine afterwards
, therefore the present invention is with a wide range of applications.
Present embodiment shows that possessed ability of the present invention, compared with general classics discrete Fourier transform method
Compared with the case where sampled data lacks because of temporal regularity, when highest harmonic wave is primary demand as 60 times, it is necessary to ask
The influence matrix and its inverse matrix, workload of solution 120 × 120 are very huge.But using the present invention, it is only necessary to 26 × 6 matrixes,
The calculating of 43 × 3 matrixes and invert, compare under, amount of calculation is simultaneously uncomplicated, but the degree of accuracy obtained is very high.
Using example 2
In one example, signal spectral analysis method when fractional-sample data rule lacks is applied to spectrum analysis
Quantum measurement in.
The lecture experiment of the difference measurement of non-sine periodic signal of this example based on staircase waveform, actual presentation frequency analysis
Quantum measuring method.
The non-sine periodic signal of measurement is provided by signal source, the signal and quantum ladder ripple signal (in lecture experiment by
DAC exports in general stairstep signal) do difference sampling, sampled value and quantum ladder ripple signal into analog-digital converter (ADC)
Known step value is added (obtaining these step values by ADC measurements in lecture experiment), is then classical discrete Fu of block sampling
Vertical leaf transformation, its result are recovered by the inverse matrix of correlation, obtain the fourier coefficient of each harmonic of measurement signal.
Here what the step data for first having to solve staircase waveform come from, and this is related to the alignment issues of two waveforms.It is former
Then the step data of upper staircase waveform should follow the sampled data of measurement signal.It should specifically carry out as follows.
Measurement signal is sampled by ADC.If by description of the present invention to being lacked by temporal regularity, a cycle N groups
Sampling, every group of M sampled data, sampling interval should be taken as measurement signal and adopted by formula (6) expression, then the step value of the n-th step
That sampled data among sample n-th group, i.e.,:
zn=y (nMh+Mh/2) (7)
Fig. 5 shows the difference measurement of the non-sine periodic signal according to an embodiment of the invention based on staircase waveform
In, the alignment schematic diagram of two signals.
The alignment of two signal waveforms is as shown in Figure 5.So make the crosspoint of two waveforms among every group, when
When abandoning the data on every group of both sides, less difference value is positioned at centre, required for this is difference measurement.
The parameter of lecture experiment is as follows, and ladder wave frequency rate is 1.0416kHz, 20 steps, per the μ s of step width 48;It is anon-normal
The frequency component arranges value of string periodic signal is identical with table 1, and fundamental frequency is also 1.0416kHz;Sample frequency is 1MHz, because
This each step (or there be not group) has sampling M=48, and each cycle is sampled as 960;Block sampling selects (12,25,11);Synchronous letter
Non-sine periodic signal source number is derived from, controls DAC and ADC synchronous operations.
ADC has two input ports A and B, and staircase waveform and non-sine periodic signal are respectively connected into input port A and B.
Measured first under Mo ShiA &B, i.e., measure two waveforms respectively.Step is obtained by the sampled data of staircase waveform
Value, by the sampled data of non-sine periodic signal, the frequency spectrum of measurement signal is obtained by common classical DFT,
Data i.e. in a cycle all participate in calculating, do not abandon the situation of data, it is not required that the help of inverse matrix, this result
Calculating as error below refers to, and resulting error is exactly of the invention to do ratio with common classical DFT
Compared with.
Then measured under Mode A-B, obtain difference sampled value, become according to the classical discrete Fourier of block sampling
Swap-in row calculates, and obtains the intermediate result of signal spectrum, then the recovery by inverse matrix, obtains the final result of frequency spectrum.
In order to overcome the difference between common-mode rejection ratio present in difference measurement (CMRR), and two passages, in experiment,
Two signal exchanges are connected to input port A and B, and do multiple measurement.
Table 3 illustrates the result of experiment, and error is the final result of the present invention and the common discrete Fourier of classics in table
The difference of leaf transformation result, because measurement signal maximum amplitude is close to 1V, so the result can be considered as error relative value.
Table 3
As frequency fi/f1When=20,40,60, often there is larger error, in fi/f1When=10,30,50, sometimes
Larger error occurs.These errors are not included in table.Preliminary experiment shows, two in table block sampling pattern
Under, error is within ± 40 μ V/V.
Lecture experiment shows that error reaches 10-6Magnitude, much larger than the analog result 10 of application example 1-9Magnitude.This is real
Caused by the noise tested, this noise is typically 10-5Magnitude, the average effect of classical DFT make result in table
Reach 10-6Magnitude is rational.See simultaneously, in the case of N=20M=48, (12,25,11), error is up to ± 40.0
μ V/V, this is due to that the classical DFT of block sampling can be interpreted as roughly from the process for partly knowing entirety,
Similar to extrapolation, commonly greater than 1, this has been showed the design factor of presentation in above-mentioned inverse matrix.From such understanding,
We can increase sampling utilization rate, i.e. part and the ratio of entirety, reach the purpose for reducing this error.
Another block sampling is have selected for this, its pattern N=20M=100, (15,71,14), wherein s=15, t=
85, meet s+t=100=M condition, it is (71/100) now to sample utilization rate, more than (25/48) of former pattern.It is real
Test result to show, under this pattern, error is within ± 4.0 μ V/V on most of Frequency point, hence it is evident that smaller than former pattern.
It is additionally observed that in lecture experiment, frequency fi/f1When=20,40,60, often occur larger error, and frequency
fi/f1When=10,30,50, it sometimes appear that larger error, preliminary experiment shows, in above two block sampling pattern
Under, this error is up to ± 40 μ V/V.The reason for this error is due to frequency fi/f1=20,40,60 with step number (or point
Group number) N=20 is relevant.If similar noise be present on each step, it will produce error on these Frequency points;As for
fi/f1When=10,30,50, situation is substantially similar.But contribution of these noise profiles to other Frequency points, will be due to offsetting
Effect and error is diminished.
In actual measurement, this kind of error can be overcome by changing segments N.
In the quantum measurement of frequency analysis, in lecture experiment DAC produce that staircase waveform this link is changed to directly input can
Quantum voltage is programmed, for ADC also without the sampled measurements for being Mo ShiA &B, final result should be uncertainty, rather than error,
The level of uncertainty data and quantum voltage, the level with ADC are relevant with the method for measurement to be assessed, it is necessary to make a concrete analysis of
Go out.
In summary, signal spectral analysis method is expanded to new scope by the present invention, have relatively low operand and compared with
The high degree of accuracy.
It will be understood by those skilled in the art that the purpose of the description to embodiments of the invention is only for exemplarily saying above
The beneficial effect of bright embodiments of the invention, it is not intended to limit embodiments of the invention to given any example.
According to an embodiment of the invention, there is provided signal spectral analysis system during a kind of fractional-sample data rule missing
System, can include:Memory, it is stored with computer executable instructions, and the signal sample data of excalation;Processor,
Computer executable instructions in processor run memory, perform following steps:Lacked for fractional-sample data rule
Measurement signal, using the Discrete Fourier transform operation of block sampling mode, obtain the intermediate result of measurement signal frequency spectrum;From
The rule that orthogonality of trigonometric function is presented in the case of missing data, influence matrix corresponding to calculating, and then obtain multiple drops
Rank matrix;Calculate the inverse matrix of degree reduction matrix;Based on inverse matrix, intermediate result is converted into final frequency spectrum data.
In one example, the measurement signal of fractional-sample data rule missing is:Appearance is lacked and lacked at equal intervals
Quantity identical measurement signal.
Signal spectral analysis method is expanded to new scope by the present invention, has relatively low operand and higher accurate
Degree.
It will be understood by those skilled in the art that the purpose of the description to embodiments of the invention is only for exemplarily saying above
The beneficial effect of bright embodiments of the invention, it is not intended to limit embodiments of the invention to given any example.
It is described above various embodiments of the present invention, described above is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.In the case of without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes will be apparent from for the those of ordinary skill in art field.
Claims (10)
1. a kind of signal spectral analysis method during fractional-sample data rule missing, including:
For the measurement signal of fractional-sample data rule missing, transported using the DFT of block sampling mode
Calculate, obtain the intermediate result of the measurement signal frequency spectrum;
The rule presented from orthogonality of trigonometric function in the case of missing data, influence matrix corresponding to calculating, and then obtain
Multiple degree reduction matrixes;
Calculate the inverse matrix of the degree reduction matrix;
Based on the inverse matrix, the intermediate result is converted into final frequency spectrum data.
2. signal spectral analysis method during fractional-sample data rule missing according to claim 1, wherein, it is described
Fractional-sample data rule missing measurement signal be:Appearance lacks and lacks quantity identical measurement signal at equal intervals.
3. signal spectral analysis method during fractional-sample data rule missing according to claim 1, wherein, segmentation
The Discrete Fourier transform operation of sample mode is:
<mrow>
<mo>{</mo>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>a</mi>
<mi>i</mi>
<mi>P</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>s</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mi>s</mi>
</mrow>
<mi>t</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>y</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>m</mi>
</mrow>
</msub>
<mi>cos</mi>
<mi> </mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>,</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>b</mi>
<mi>i</mi>
<mi>P</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>s</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mi>s</mi>
</mrow>
<mi>t</mi>
</munderover>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>y</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>m</mi>
</mrow>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>,</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein,For intermediate result, aiAnd biThe respectively cosine coefficient and sinusoidal coefficients of ith harmonic wave, yn,mBelieve for measurement
Number obtained sampled data, h are the sampling interval, and N is the sum of measurement signal even group-division in a cycle, and n believes for measurement
Number deck label, m is the sampled data label of every group of measurement signal equal intervals, and [0, s-1] and [t+1, M-1] are missing data
Label range, [s, t] are retention data label range, and M is the entire quantity of the sampled data of every group of measurement signal equal intervals.
4. signal spectral analysis method during fractional-sample data rule missing according to claim 3, wherein, it is described
Influence matrix is:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>a</mi>
<mi>a</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>s</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mfrac>
<mfrac>
<mn>2</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mi>s</mi>
</mrow>
<mi>t</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mi>cos</mi>
<mi> </mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>,</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mi> </mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>,</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>b</mi>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>s</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mfrac>
<mfrac>
<mn>2</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mi>s</mi>
</mrow>
<mi>t</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mi>sin</mi>
<mi> </mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>,</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mi> </mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>,</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>a</mi>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>s</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mfrac>
<mfrac>
<mn>2</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mi>s</mi>
</mrow>
<mi>t</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mi>cos</mi>
<mi> </mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>,</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mi> </mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>,</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, aa (i, j) is the element of cosine and cosine interaction matrix AA, and bb (i, j) is sinusoidal and sinusoidal interaction
Matrix B B element, ab (i, j) are the elements of cosine and sinusoidal interaction matrix AB.
5. signal spectral analysis method during fractional-sample data rule missing according to claim 4, wherein, work as t+
During s=M, influence matrix AB all elements are zero, and matrix A A is decomposed into three class degree reduction matrixes:Ak, A0.5N, and AN, influence square
Battle array BB is decomposed into three class degree reduction matrixes:Bk, B0.5N, and BN。
6. signal spectral analysis method during fractional-sample data rule missing according to claim 1, wherein, it is described
Finally the cosine result of frequency spectrum data is:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfenced open = "|" close = "|">
<mtable>
<mtr>
<mtd>
<msub>
<mi>a</mi>
<mi>k</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>a</mi>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>a</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>a</mi>
<mrow>
<mn>2</mn>
<mi>N</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>a</mi>
<mrow>
<mn>2</mn>
<mi>N</mi>
<mo>+</mo>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>A</mi>
<mi>k</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfenced open = "|" close = "|">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>a</mi>
<mi>P</mi>
</msup>
<mi>k</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>a</mi>
<mi>P</mi>
</msup>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>a</mi>
<mi>P</mi>
</msup>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mi>k</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>a</mi>
<mi>P</mi>
</msup>
<mrow>
<mn>2</mn>
<mi>N</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>a</mi>
<mi>P</mi>
</msup>
<mrow>
<mn>2</mn>
<mi>N</mi>
<mo>+</mo>
<mi>k</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfenced open = "|" close = "|">
<mtable>
<mtr>
<mtd>
<msub>
<mi>a</mi>
<mrow>
<mn>0.5</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>a</mi>
<mrow>
<mn>1.5</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>A</mi>
<mrow>
<mn>0.5</mn>
<mi>N</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfenced open = "|" close = "|">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>a</mi>
<mi>P</mi>
</msup>
<mrow>
<mn>0.5</mn>
<mi>N</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>a</mi>
<mi>P</mi>
</msup>
<mrow>
<mn>1.5</mn>
<mi>N</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfenced open = "|" close = "|">
<mtable>
<mtr>
<mtd>
<msub>
<mi>a</mi>
<mi>N</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>a</mi>
<mrow>
<mn>2</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>A</mi>
<mi>N</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfenced open = "|" close = "|">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>a</mi>
<mi>P</mi>
</msup>
<mi>N</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>a</mi>
<mi>P</mi>
</msup>
<mrow>
<mn>2</mn>
<mi>N</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, akIt is that frequency is kf1Frequency spectrum final result cosine coefficient, (Ak)-1,(A0.5N)-1, (AN)-1It is the three of matrix A A
Individual price reduction inverse of a matrix matrix.
7. signal spectral analysis method during fractional-sample data rule missing according to claim 1, wherein, it is described
Finally the sinusoidal result of frequency spectrum data is:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfenced open = "|" close = "|">
<mtable>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mi>k</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mrow>
<mn>2</mn>
<mi>N</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mrow>
<mn>2</mn>
<mi>N</mi>
<mo>+</mo>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>B</mi>
<mi>k</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfenced open = "|" close = "|">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>b</mi>
<mi>P</mi>
</msup>
<mi>k</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>b</mi>
<mi>P</mi>
</msup>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>b</mi>
<mi>P</mi>
</msup>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mi>k</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>b</mi>
<mi>P</mi>
</msup>
<mrow>
<mn>2</mn>
<mi>N</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>b</mi>
<mi>P</mi>
</msup>
<mrow>
<mn>2</mn>
<mi>N</mi>
<mo>+</mo>
<mi>k</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfenced open = "|" close = "|">
<mtable>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mrow>
<mn>0.5</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mrow>
<mn>1.5</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>B</mi>
<mrow>
<mn>0.5</mn>
<mi>N</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfenced open = "|" close = "|">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>b</mi>
<mi>P</mi>
</msup>
<mrow>
<mn>0.5</mn>
<mi>N</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>b</mi>
<mi>P</mi>
</msup>
<mrow>
<mn>1.5</mn>
<mi>N</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfenced open = "|" close = "|">
<mtable>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mi>N</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mrow>
<mn>2</mn>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>B</mi>
<mi>N</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfenced open = "|" close = "|">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>b</mi>
<mi>P</mi>
</msup>
<mi>N</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>b</mi>
<mi>P</mi>
</msup>
<mrow>
<mn>2</mn>
<mi>N</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, bkIt is that frequency is kf1Frequency spectrum final result sinusoidal coefficients, (Bk)-1,(B0.5N)-1, (BN)-1It is the three of matrix B B
Individual price reduction inverse of a matrix matrix.
Signal spectral analysis method when 8. the fractional-sample data rule in claim 1-7 described in any one lacks,
Applied in the quantum measurement of spectrum analysis.
9. signal spectral analysis system during a kind of fractional-sample data rule missing, it is characterised in that the system includes:
Memory, it is stored with computer executable instructions, and the signal sample data of excalation;
Processor, the processor run the computer executable instructions in the memory, perform following steps:
For the measurement signal of fractional-sample data rule missing, classical discrete Fourier change is carried out using block sampling mode
Computing is changed, obtains the intermediate result of the measurement signal frequency spectrum;
The rule presented from orthogonality of trigonometric function in the case of missing data, influence matrix corresponding to calculating, and then obtain
Multiple degree reduction matrixes;
Calculate the inverse matrix of the degree reduction matrix;
Based on the inverse matrix, the intermediate result is converted into final frequency spectrum data.
10. signal spectral analysis system during fractional-sample data rule missing according to claim 9, wherein, institute
Stating the measurement signal that fractional-sample data rule lacks is:Appearance lacks and lacks quantity identical measurement signal at equal intervals.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710900980.4A CN107766293B (en) | 2017-09-28 | 2017-09-28 | Signal spectrum analysis method and system when part of sampling data is missing regularly |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710900980.4A CN107766293B (en) | 2017-09-28 | 2017-09-28 | Signal spectrum analysis method and system when part of sampling data is missing regularly |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107766293A true CN107766293A (en) | 2018-03-06 |
CN107766293B CN107766293B (en) | 2021-04-23 |
Family
ID=61266324
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710900980.4A Active CN107766293B (en) | 2017-09-28 | 2017-09-28 | Signal spectrum analysis method and system when part of sampling data is missing regularly |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107766293B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI718625B (en) * | 2019-08-16 | 2021-02-11 | 瑞昱半導體股份有限公司 | Computation circuit used in dct, dst, idct and idst |
CN112422979A (en) * | 2019-08-23 | 2021-02-26 | 瑞昱半导体股份有限公司 | Operation circuit applied to discrete and inverse discrete sine and cosine transform |
CN116087590A (en) * | 2022-06-27 | 2023-05-09 | 中国电力科学研究院有限公司 | Differential sampling control device and method for quantum voltage precise measurement |
CN116359605A (en) * | 2023-04-21 | 2023-06-30 | 中国计量科学研究院 | Harmonic signal analysis method based on secondary weighting |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101944917A (en) * | 2010-06-23 | 2011-01-12 | 王红星 | Method for generating waveforms of non-sinusoidal orthogonal bandpass signals in time domain |
CN104142424A (en) * | 2014-08-15 | 2014-11-12 | 中国计量科学研究院 | Method and system for accurately measuring non-sinusoidal voltage signals |
CN104407197A (en) * | 2014-11-27 | 2015-03-11 | 湖南大学 | Signal phasor measurement method based on trigonometric function iteration |
US20170089995A1 (en) * | 2014-05-19 | 2017-03-30 | The United States Of America, As Represented By The Secretary, Department Of Health And Human Serv | Magnetic resonance 2d relaxometry reconstruction using partial data |
CN106842321A (en) * | 2015-12-04 | 2017-06-13 | 中国石油天然气股份有限公司 | Seismic data reconstruction method and device |
-
2017
- 2017-09-28 CN CN201710900980.4A patent/CN107766293B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101944917A (en) * | 2010-06-23 | 2011-01-12 | 王红星 | Method for generating waveforms of non-sinusoidal orthogonal bandpass signals in time domain |
US20170089995A1 (en) * | 2014-05-19 | 2017-03-30 | The United States Of America, As Represented By The Secretary, Department Of Health And Human Serv | Magnetic resonance 2d relaxometry reconstruction using partial data |
CN104142424A (en) * | 2014-08-15 | 2014-11-12 | 中国计量科学研究院 | Method and system for accurately measuring non-sinusoidal voltage signals |
CN104407197A (en) * | 2014-11-27 | 2015-03-11 | 湖南大学 | Signal phasor measurement method based on trigonometric function iteration |
CN106842321A (en) * | 2015-12-04 | 2017-06-13 | 中国石油天然气股份有限公司 | Seismic data reconstruction method and device |
Non-Patent Citations (4)
Title |
---|
BEN ADCOCK: "Weighted frames of exponentials and stable recovery of multidimensional functions from nonuniform Fourier samples", 《APPLIED AND COMPUTATIONAL HARMONIC ANALYSIS》 * |
P. ESPEL: "Traceability of Voltage Measurements for Non-Sinusoidal Waveforms", 《MEASUREMENT SCIENCE REVIEW》 * |
ZULIANG LU: "Further Investigation for Piecewise Sampling to Overcome Transient Effect of Staircase Waveform", 《IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT》 * |
武艳华: "频谱分析理论在谐波及间谐波检测中的应用研究", 《中国优秀博硕士学位论文全文数据库 (硕士)工程科技Ⅱ辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI718625B (en) * | 2019-08-16 | 2021-02-11 | 瑞昱半導體股份有限公司 | Computation circuit used in dct, dst, idct and idst |
CN112422979A (en) * | 2019-08-23 | 2021-02-26 | 瑞昱半导体股份有限公司 | Operation circuit applied to discrete and inverse discrete sine and cosine transform |
CN116087590A (en) * | 2022-06-27 | 2023-05-09 | 中国电力科学研究院有限公司 | Differential sampling control device and method for quantum voltage precise measurement |
CN116087590B (en) * | 2022-06-27 | 2024-02-23 | 中国电力科学研究院有限公司 | Differential sampling control device and method for quantum voltage precise measurement |
CN116359605A (en) * | 2023-04-21 | 2023-06-30 | 中国计量科学研究院 | Harmonic signal analysis method based on secondary weighting |
CN116359605B (en) * | 2023-04-21 | 2023-09-12 | 中国计量科学研究院 | Harmonic signal analysis method based on secondary weighting |
Also Published As
Publication number | Publication date |
---|---|
CN107766293B (en) | 2021-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107766293A (en) | Signal spectral analysis method and system when fractional-sample data rule lacks | |
CN109633262A (en) | Three phase harmonic electric energy gauging method, device based on composite window multiline FFT | |
CN104237622B (en) | The method of sampling and wideband voltage/calibration of power device based on software based frequency tracking | |
CN103457603B (en) | A kind of method that ADC dynamic parameter is tested based on average frequency spectrum | |
CN102955068B (en) | A kind of harmonic detecting method based on compression sampling orthogonal matching pursuit | |
CN101900761B (en) | High-accuracy non-integer-period sampled harmonic analysis and measurement method | |
CN106556342A (en) | A kind of grating subdivision device and method based on FPGA | |
CN111865331B (en) | Phase calibration method for MWC extended system transfer matrix | |
CN109581018A (en) | Wide-band sampling oscillograph time base error synchroballistic method and device | |
CN109446552A (en) | Multiaxis correlation arbitrary excitation flowering structure fatigue life time-domain calculation method | |
CN116359605A (en) | Harmonic signal analysis method based on secondary weighting | |
CN109992748A (en) | A kind of river mouth tidal level harmonic analysis method | |
CN102735937A (en) | Method for measuring signal phase difference | |
CN104142424B (en) | A kind of method of accurate measurement non-sinusoidal voltage signal and system thereof | |
CN105067882B (en) | A kind of assay method of electricity amplitude | |
CN104731762B (en) | Cube phase modulated parameter estimating method based on cyclic shift | |
CN106649201A (en) | Data complementation method based on Vandermonde decomposition of exponential signal | |
CN105807128A (en) | Method and system for measuring alternating voltages by applying multicycle strategy to digital-to-analogue conversion | |
CN103576120B (en) | Third-harmonic component plesiochronous information transfer check and Self-healing Algorithm | |
CN110320400A (en) | Quasi-synchro sampling and the voltage flicker envelope parameters extracting method for improving energy operator | |
CN108982966B (en) | Harmonic phase angle analysis method based on linear correction algorithm | |
CN207147415U (en) | Signal subdividing device | |
CN105846822A (en) | Method and system for measuring mismatch error of time-interleaved analog-to-digital conversion system | |
Shao et al. | Nurbs model for chaotic time series | |
CN107809247A (en) | A kind of High Speed High Precision ADC dynamic input-output characteristic curve method for rapidly testing |
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 |