CN107945147A - One kind is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT - Google Patents
One kind is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT Download PDFInfo
- Publication number
- CN107945147A CN107945147A CN201711259381.5A CN201711259381A CN107945147A CN 107945147 A CN107945147 A CN 107945147A CN 201711259381 A CN201711259381 A CN 201711259381A CN 107945147 A CN107945147 A CN 107945147A
- Authority
- CN
- China
- Prior art keywords
- amplitude
- peak
- frequency
- formula
- kaiser
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000012937 correction Methods 0.000 claims abstract description 68
- 230000003595 spectral effect Effects 0.000 claims abstract description 45
- 238000001228 spectrum Methods 0.000 claims abstract description 18
- 238000006073 displacement reaction Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000001514 detection method Methods 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 claims description 4
- 210000001747 pupil Anatomy 0.000 claims description 4
- 238000013459 approach Methods 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 7
- 238000004088 simulation Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000009413 insulation Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 206010037660 Pyrexia Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
The invention discloses one kind to be based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT, including based on tetra- peak interpolation correction formulas of Kaiser windows FFT, the frequency displacement of FFT and the frequency spectrum of window function Kaiser fundamental wave, harmonic wave and the symmetrical spectral line of left and right four of m-Acetyl chlorophosphonazo crest frequency point amplitude, using fitting polynomial formulas, general amplitude, the phase and frequency interpolation correction formula of signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo are obtained；The tested fundamental wave of picture signal, the amplitude of harmonic wave and m-Acetyl chlorophosphonazo, phase and frequency parameter are calculated according to general amplitude, phase and frequency interpolation correction formula；And then calculate fundamental wave, harmonic wave, m-Acetyl chlorophosphonazo.Present invention tool utilizes four peak interpolation correction algorithms, takes full advantage of the symmetry of four spectral lines of amplitude Frequency point or so and excellent master, the side lobe performance of Kaiser window functions, has higher accuracy in computation.
Description
Technical field
The present invention relates to electric power system electrical energy measurement field and digitalized electric energy metering field, and in particular to one kind is based on
The tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT.
Background technology
Insulation performance deterioration high-tension apparatus insulate it is entirely ineffective before often with early stage abnormal heating or
Discharge sign signal, and fever is accompanied by the radiation of infrared ray, and discharges and be accompanied by ultraviolet radiation, therefore passes through thermal infrared imager
The insulation status of operational outfit and timely discovering device defect can be assessed indirectly with ultraviolet imager.
But in engineering in practice, some equipment deficiencies can not be detected by ultraviolet imager and infrared thermoviewer at the same time,
Easily there is missing inspection in a kind of alone instrument, lacks the mechanism of infrared and ultraviolet joint-detection analysis.Charged electric power apparatus inspection at the same time
It is time-consuming and laborious, maintenance task weight, cycle length, it is difficult to which discovering device whether there is hidden danger in time, with the hair of on-line monitoring technique
Exhibition, substation equipment detection work, which will be also realized to automation, unmanned, on-line analysis direction, to be developed.
Three wave band optical imagery on-Line Monitor Devices are researched and developed, realize the infrared temperature and SF6 of substation's whole station power equipment
Leak detection, the monitoring in real time of ultraviolet, high definition visible ray joint, innovatively combine infrared, ultraviolet and high definition visible light image information into
Row analysis, with reference to power industry standard, carries out feature extraction to image deflects information using adaptive neural network, draws rapidly
Defect information simultaneously carries out analysis judgement, and related foundation is proposed for maintenance.
The content of the invention
It is an object of the invention to propose that one kind is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT, fill
Divide excellent master, the side lobe performance of the symmetry and Kaiser window functions using four spectral lines of crest frequency point or so, to reach
The requirement of high precision image fusion.
One kind of the invention is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT, comprises the following steps：S1：
Based on tetra- peak interpolation correction formulas of Kaiser windows FFT, the frequency displacement of FFT and the frequency spectrum of window function Kaiser fundamental wave, harmonic wave and
The amplitude of the symmetrical spectral line of m-Acetyl chlorophosphonazo crest frequency point arranged on left and right sides four obtains the general width of signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo
Value, phase and frequency interpolation correction formula；S2：Calculated according to the general amplitude, phase and frequency correction formula by altimetric image
The fundamental wave of signal, the amplitude of harmonic wave and m-Acetyl chlorophosphonazo, phase and frequency parameter；S3：According to the fundamental wave of the tested picture signal, humorous
The amplitude of ripple and m-Acetyl chlorophosphonazo, phase and frequency parameter calculate fundamental wave, harmonic wave, m-Acetyl chlorophosphonazo respectively.
One kind of the invention is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT, and Kaiser window functions are one
The adjustable window function of group, the proportion of main lobe width and side lobe height is adjusted by unrestricted choice form parameter β, β value is bigger, window
The main lobe width of spectrum of function is wider, and side lobe peak is smaller.This patent uses the Kaiser window functions of form parameter β=22, its
Main lobe width is 28 π/N, and side lobe peak is -172dB, side lobe attenuation speed 14dB/otc, has preferable master, side lobe performance.
In addition, the present invention is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT, can also have following attached
The technical characteristic added：
Further, step S1 further comprises：S101：Based on Kaiser windows FFT tetra- peak interpolation correction formulas, FFT
The frequency spectrum of frequency displacement and window function Kaiser are in fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo crest frequency point four symmetrical spectral lines of arranged on left and right sides
Amplitude give expression to the formula of four peak spectral line amplitude ratio parameter γ；S102：Contrast value parameter γ carries out fitting of a polynomial and approaches,
Fit the formula of the universal frequency correction factor δ based on tetra- peak spectral lines of Kaiser windows FFT；S103：According to the width of measured signal
Value formula fitting goes out the formula of general amplitude correction coefficient g (δ)；S104：According to the public affairs of general amplitude correction coefficient g (δ)
The formula of formula and universal frequency correction factor δ obtain general amplitude correction formula, general phase only pupil filter formula and universal frequency and repair
Positive formula.
Further, step S2 further comprises：S201：Apply Kaiser window functions respectively to picture signal, and carry out
FFT transform is to obtain the frequency spectrum that the signal applies Kaiser window functions；S202：Using peak detection technology, detect respectively
It is tested the amplitude of four symmetrical spectral lines of voltage, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo peak point arranged on left and right sides；
S203：Using four peak spectral line amplitude ratio parameter ratio formula, the fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo of picture signal are calculated respectively
The value of four peak spectral line amplitude ratio parameter γ at peak point；S204：According to the fundamental wave of picture signal, harmonic wave and m-Acetyl chlorophosphonazo peak
The value of four peak spectral line amplitude ratio parameter γ at value point, and calculated according to the tetra- peak interpolation correction formulas of Kaiser windows FFT
The value of frequency correction factor δ at measured signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo peak point；S205：According to the value of correction factor δ
Obtain the value of the amplitude correction coefficient g (δ) at measured signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo peak point；S206：It is according to correcting
The value of number δ and the value interpolation of correction factor g (δ) correct out the width of tested voltage, the fundamental wave of current signal, harmonic wave and m-Acetyl chlorophosphonazo
Value, phase and frequency.
Further, the calculation formula of universal frequency correction factor δ is：
δ=+ 0.7046599 γ of 3.1041273 γ3+0.35296224·γ5+0.22441945·γ7+
0.16913521·γ9。
Further, the calculation formula of general amplitude correction coefficient g (δ) is：
G (δ)=1.50096168+0.22288389 δ2+0.01752579·δ4+0.00098765·δ6。
Further, the correction formula of general amplitude is：
A=2 (y1+2y2+2y3+y4)·g(δ)/N；
Wherein, y1、y2、y3And y4Four symmetrical spectral lines of respectively spectrum curve peak point arranged on left and right sides are respectively kth1
Bar, kth2Bar, kth3Bar and kth4The amplitude of bar spectral line, N are sampled data length.
The additional aspect and advantage of the present invention will be set forth in part in the description, and will partly become from the following description
Obtain substantially, or recognized by the practice of the present invention.
The beneficial effects of the invention are as follows the flexibility using Kaiser window function form parameters, by setting suitable parameter
The special efficacy for excellent master, the Kaiser window functions of side lobe performance, restraining by Kaiser window functions and preferentially extracting is obtained,
Realize the purpose shown in high precision that picture signal is dynamically tested under four peak interpolation algorithms.
Brief description of the drawings
Fig. 1 is a kind of flow chart based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT of the present invention；
Fig. 2 a are four symmetrical spectral lines of this electricity net fundamental frequency more than peak point arranged on left and right sides in the case of 50Hz
Schematic diagram；
Fig. 2 b are four symmetrical spectral lines of this electricity net fundamental frequency less than peak point arranged on left and right sides in the case of 50Hz
Schematic diagram；
Fig. 3 is spectral characteristic of the Kaiser window functions when form parameter β is respectively 0,4,8,22.It is other with the increase of β
Valve peak value is smaller, but main lobe width broadens.As β=22, there is relatively low side lobe peak and faster side lobe attenuation speed, together
When with moderate main lobe width.
Embodiment
Such as Fig. 1-3, one kind is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT, comprises the following steps：
S1：Existed based on tetra- peak interpolation correction formulas of Kaiser windows FFT, the frequency displacement of FFT and the frequency spectrum of window function Kaiser
The amplitude of fundamental wave, harmonic wave and the symmetrical spectral line of m-Acetyl chlorophosphonazo crest frequency point arranged on left and right sides four obtains signal fundamental wave, harmonic wave and humorous
General amplitude, the phase and frequency interpolation correction formula of ripple.
Step S1 further comprises：
S101：Based on tetra- peak interpolation correction formulas of Kaiser windows FFT, the frequency displacement of FFT and the frequency spectrum of window function Kaiser
Four peak spectral line Amplitude Rations are given expression in the amplitude of fundamental wave, harmonic wave and the symmetrical spectral line of m-Acetyl chlorophosphonazo crest frequency point arranged on left and right sides four
The formula of value parameter γ.
S102：Contrast value parameter γ carries out fitting of a polynomial and approaches, and fits logical based on tetra- peak spectral lines of Kaiser windows FFT
With the formula of frequency correction factor δ.
S103：Go out the formula of general amplitude correction coefficient g (δ) according to the amplitude formula fitting of measured signal.
S104：Led to according to the formula of the formula of general amplitude correction coefficient g (δ) and universal frequency correction factor δ
With amplitude correction formula, general phase only pupil filter formula and universal frequency correction formula.
S2：The base of tested voltage signal, current signal is calculated according to the general amplitude, phase and frequency correction formula
Ripple, the amplitude of harmonic wave and m-Acetyl chlorophosphonazo, phase and frequency parameter.
Step S2 further comprises：
S201：Apply Kaiser window functions respectively to picture signal, and carry out FFT transform to obtain the Kaiser windows letter
Several frequency spectrums；
S202：Using peak detection technology, detect to be tested fundamental wave, harmonic wave and the m-Acetyl chlorophosphonazo peak value of picture signal respectively
The amplitude of the point symmetrical spectral line of arranged on left and right sides four；
S203：Using four peak spectral line amplitude ratio parameter ratio formula, calculate respectively fundamental wave, the harmonic wave of picture signal with
And the value of four peak spectral line amplitude ratio parameter γ at m-Acetyl chlorophosphonazo peak point；
S204：According to four peak spectral line amplitude ratio parameters at the fundamental wave of picture signal, harmonic wave and m-Acetyl chlorophosphonazo peak point
The value of γ, and measured signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo peak are calculated according to the tetra- peak interpolation correction formulas of Kaiser windows FFT
The value of frequency correction factor δ at value point；
S205：Repaiied according to the amplitude being worth at measured signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo peak point of correction factor δ
The value of positive coefficient g (δ)；
S206：According to the value interpolation of the value of correction factor δ and correction factor g (δ) correct out tested picture signal fundamental wave,
The amplitude of harmonic wave and m-Acetyl chlorophosphonazo, phase and frequency.
S3：Counted respectively according to the fundamental wave of the tested picture signal, the amplitude of harmonic wave and m-Acetyl chlorophosphonazo, phase and frequency parameter
Calculate fundamental energy, harmonic electric energy, m-Acetyl chlorophosphonazo electric energy and total electric energy.
In one embodiment of the invention, the calculation formula of universal frequency correction factor δ is：
δ=+ 0.7046599 γ of 3.1041273 γ3+0.35296224·γ5+0.22441945·γ7+
0.16913521·γ9。
In one embodiment of the invention, the calculation formula of general amplitude correction coefficient g (δ) is：
G (δ)=1.50096168+0.22288389 δ2+0.01752579·δ4+0.00098765·δ6。
The correction formula of general amplitude is：
A=2 (y1+2y2+2y3+y4)·g(δ)/N；
Wherein, y1、y2、y3And y4The respectively symmetrical spectral line of spectrum curve peak point arranged on left and right sides four is respectively kth1Bar,
Kth2Bar, kth3Bar and kth4The amplitude of bar spectral line, N are sampled data length.
To make it is further understood that the present invention, will be described further by following embodiments.
(1) tetra- peak interpolation correction formulas of Kaiser windows FFT are based on
The time-domain expression of Kaiser window functions：
In formula, I0(β) is the 1st class deformation zero Bessel function；β is the form parameter of Kaiser functions；
In formula, α is the difference between the main lobe value and side lobe levels of Kaiser windows, by setting different β value, thus it is possible to vary
The main lobe and side lobe performance of Kaiser window functions.
Universal frequency correction factors formula δ,
δ=+ 0.7046599 γ of 3.1041273 γ3+0.35296224·γ5+0.22441945·γ7+
0.16913521·γ9 (1)
General amplitude correction coefficient formula g (δ),
G (δ)=1.50096168+0.22288389 δ2+0.01752579·δ4+0.00098765·δ6 (2)
General amplitude correction formula：
A=2 (y1+2y2+2y3+y4)·g(δ)/N (3)
General phase only pupil filter formula：
Universal frequency correction formula：
f0=k0Δ f=(δ+k2-0.5)fs/N (5)
(2) derivation
Assuming that the signal containing multiple harmonic (m-Acetyl chlorophosphonazo), with sample frequency fsAfter being sampled, discrete series x is obtained
(n)：
In formula, m represents overtone order；AmAnd θmThe amplitude and phase of m subharmonic are represented respectively；f0Represent fundamental frequency.
Formula (6) is used instead Euler's formula to be expressed as：
The forms of time and space of Kaiser window functions is w (n), its discrete spectrum isThe then FFT tables after x (n) adding windows
It is up to formula:
If ignore negative frequency point-mf0Locating the secondary lobe of spectral peak influences, in positive frequency point mf0Neighbouring continuous frequency spectrum function is：
Formula (9) is sampled, the expression formula that can obtain its discrete Fourier transform is：
In formula, Δ f=fs/N is discrete frequency intervals；N is sampled data length.
With fundamental frequency f0=k0Exemplified by Δ f.Since picture frequency is there are fluctuation, and it is generally not accurate 50Hz, such as
This, samples measured signal with fixed frequency, and fence effect will occur, i.e. the Frequency point of its highest discrete spectral line is very
The frequency of hardly possible exactly spectrum curve peak value, that is, k0It is generally not integer.If spectrum curve peak point arranged on left and right sides four is right
It is respectively kth to claim spectral line1Bar, kth2Bar, kth3Bar and kth4Bar spectral line (k1≤k2≤k0≤k3≤k4), corresponding amplitude is respectivelyAs shown in Figure 2.Introduce parameter
δ=k0-k1- 0.5, obtain δ and obtain [- 0.5,0.5], make four peak spectral line amplitude ratio parameter ratio γ=((2y3+y4)-(2y2+
y1))/(y1+2y2+2y3+y4), order：
R=2 | W [2 π (- δ+0.5)/N] |+| W [2 π (- δ+1.5)/N] | (11)
S=2 | W [2 π (- δ -0.5)/N] |+| W [2 π (- δ -1.5)/N] | (12)
According to fitting of a polynomial regression formulas polyfit (γ, δ, 9), inverse function δ=f is obtained-1(γ), that is, obtain general
The formula of frequency correction factor δ, is shown in formula (1).
Due to the kth near peak point2Bar and kth3The amplitude information amount of bar spectral line is maximum, when then calculating, gives this two
Spectral line is with the weight of bigger, and specifically, in turn, the weight of four spectral lines is respectively 1,2,2,1, and the correction formula of amplitude is：
G δ=2N/ (R+S) are made, when n is large, approximation by polynomi-als is carried out by polyfit (δ, g (δ), 7) program, is asked
Go out the formula of amplitude correction coefficient g (δ), see formula (2).Bring formula (2) into formula (14), then obtain general amplitude correction formula,
See formula (3).
It can be seen that by formula (3), the actual fundamental wave peak value of measured signal is to utilize four symmetrical spectral lines of peak point arranged on left and right sides
Amplitude it is weighted correct what is obtained, therefore the algorithm is known as four peak interpolation correction algorithms.
(3) calculation formula
There are the actual fundametal compoment amplitude of measured signal, the amendment type of phase and frequency, so that it may sufficiently accurately calculate
The tested fundamental wave of actual signal, the amplitude of each harmonic and m-Acetyl chlorophosphonazo component, phase and frequency, and then can sufficiently accurately count
Calculate fundamental wave, each harmonic and m-Acetyl chlorophosphonazo.From orthogonality of trigonometric function, different integral frequency harmonizing wave components do not produce active
Electric energy, therefore the electric energy calculating formula of fundamental wave, harmonic wave, m-Acetyl chlorophosphonazo is：
In formula, WnFor n-th harmonic or m-Acetyl chlorophosphonazo, when n takes 1, W1For fundamental wave；When n takes the integer more than 1, WnTo be humorous
Ripple electric energy；When the negated integers of n, WnFor m-Acetyl chlorophosphonazo；Un、InThe respectively n-th harmonic of picture signal or the amplitude of m-Acetyl chlorophosphonazo,
It is fundamental voltage amplitude when n takes 1；αn、βnThe respectively n-th harmonic of picture signal or the initial phase of m-Acetyl chlorophosphonazo, are base when n takes 1
Wave phase；T is the primitive period of network voltage electric current；K is the number of time window, according to IEC61000-4-7 standards, is generally taken
10 primitive periods.
Meter and fundamental wave, harmonic wave, total electric energy calculation formula of m-Acetyl chlorophosphonazo are：
W=CjWj-CposxWposx+CoppxWoppx-CposjxWposjx+CoppjxWoppjx(16)
Wherein, WjFor fundamental energy；WposxFor positive harmonic electric energy；WoppxFor reverse harmonic electric energy；WposjxPositive m-Acetyl chlorophosphonazo
Electric energy；WoppjxReverse m-Acetyl chlorophosphonazo electric energy；CjFor fundamental energy weighting coefficient；CposxFor positive harmonic electric energy weighting coefficient；CoppxFor
Reverse harmonic electric energy weighting coefficient；CposjxFor positive m-Acetyl chlorophosphonazo electric energy weighting coefficient；CoppjxWeighted for reverse m-Acetyl chlorophosphonazo electric energy and be
Number.It is recommended that weighting coefficient value：Cj=1,0<Cposx、Cposjx<1, Coppx、Coppjx>1。
(4) emulate
Image simulation model of powering is as shown in table 1, wherein containing multiple harmonic component, is specially
In formula, f0For fundamental frequency, usually fluctuated between 49.7Hz~50.3Hz, sample frequency fs=4000Hz, according to
IEC61000-4-7 standards, the spectrum analysis time window of harmonic measure are 200ms by Uniform provisions.
The fundamental wave and harmonic parameters of 1 image of table
Hanning windows, Hamming windows, Blackman windows, 45 ranks of Nuttall windows and Kaiser windows (β is respectively adopted
=22) electrical energy measurement simulation calculation is carried out based on four peak interpolation algorithms, the length of window function elects as at 800 points, and fundamental frequency is
50.1Hz, simulation result are given in table 2.
The simulation result of 2 image of table
The calculating based on the modified electrical energy measurement algorithm of tetra- peak interpolations of Kaiser windows FFT is accurate it can be seen from simulation result
Exactness is compared to high one of other four kinds of algorithms even two orders of magnitude, while the accuracy of even-order harmonic electric energy is also carried significantly
Height, illustrates proposed by the present invention there is higher meter really based on the tetra- modified Image Fusions of peak interpolation of Kaiser windows FFT
Calculate exactness accurately.
Claims (6)
1. one kind is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT, it is characterised in that comprises the following steps：
S1：Based on tetra- peak interpolation correction formulas of Kaiser windows FFT, the frequency displacement of FFT and the frequency spectrum of window function Kaiser fundamental wave,
The amplitude of harmonic wave and the symmetrical spectral line of m-Acetyl chlorophosphonazo crest frequency point arranged on left and right sides four, using fitting of a polynomial, obtains signal base
Ripple, the general amplitude of harmonic wave and m-Acetyl chlorophosphonazo, phase and frequency interpolation correction formula；
S2：Fundamental wave, harmonic wave and the m-Acetyl chlorophosphonazo of tested picture signal are calculated according to the general amplitude, phase and frequency correction formula
Amplitude, phase and frequency parameter；
S3：Calculated respectively according to the fundamental wave of the tested picture signal, the amplitude of harmonic wave and m-Acetyl chlorophosphonazo, phase and frequency parameter
Fundamental wave, harmonic wave, m-Acetyl chlorophosphonazo.
2. one kind according to claim 1 is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT, its feature
It is, step S1 further comprises：
S101：Frequency spectrum based on tetra- peak interpolation correction formulas of Kaiser windows FFT, the frequency displacement of FFT and window function Kaiser is in base
The amplitude of ripple, harmonic wave and the symmetrical spectral line of m-Acetyl chlorophosphonazo crest frequency point arranged on left and right sides four gives expression to four peak spectral line amplitude ratio ginsengs
The formula of number γ；
S102：Contrast value parameter γ carries out fitting of a polynomial and approaches, and fits the general frequency based on tetra- peak spectral lines of Kaiser windows FFT
The formula of rate correction factor δ；
S103：Go out the formula of general amplitude correction coefficient g (δ) according to the amplitude formula fitting of measured signal；
S104：General width is obtained according to the formula of the formula of general amplitude correction coefficient g (δ) and universal frequency correction factor δ
It is worth correction formula, general phase only pupil filter formula and universal frequency correction formula.
3. one kind according to claim 2 is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT, its feature
It is, step S2 further comprises：
S201：Apply Kaiser window functions respectively to picture signal, and carry out FFT transform to obtain the Kaiser window functions
Frequency spectrum；
S202：Using peak detection technology, detect to be tested fundamental wave, harmonic wave and the m-Acetyl chlorophosphonazo peak point of picture signal respectively
The amplitude of the symmetrical spectral line of four of arranged on left and right sides；
S203：Using four peak spectral line amplitude ratio parameter ratio formula, calculate respectively the fundamental wave of picture signal, harmonic wave and
The value of four peak spectral line amplitude ratio parameter γ at harmonic spike point；
S204：According to four peak spectral line amplitude ratio parameter γ's at the fundamental wave of picture signal, harmonic wave and m-Acetyl chlorophosphonazo peak point
Value, and measured signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo peak point are calculated according to the tetra- peak interpolation correction formulas of Kaiser windows FFT
The value of the frequency correction factor δ at place；
S205：According to the amplitude correction system being worth at measured signal fundamental wave, harmonic wave and m-Acetyl chlorophosphonazo peak point of correction factor δ
The value of number g (δ)；
S206：The fundamental wave of tested picture signal, harmonic wave are corrected out according to the value interpolation of the value of correction factor δ and correction factor g (δ)
Amplitude, phase and frequency parameter with m-Acetyl chlorophosphonazo.
4. one kind according to claim 2 is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT, its feature
It is, the calculation formula of universal frequency correction factor δ is：
δ=+ 0.7046599 γ of 3.1041273 γ3+0.35296224·γ5+0.22441945·γ7+
0.16913521·γ9。
5. one kind according to claim 4 is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT, its feature
It is, the calculation formula of general amplitude correction coefficient g (δ) is：
G (δ)=1.50096168+0.22288389 δ2+0.01752579·δ4+0.00098765·δ6。
6. one kind according to claim 5 is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT, its feature
It is, the correction formula of general amplitude is：
A=2 (y1+2y2+2y3+y4)·g(δ)/N；
Wherein, y1、y2、y3And y4Respectively the symmetrical spectral line of spectrum curve peak point arranged on left and right sides four is respectively kth1Article,
k2Bar, kth3Bar and kth4The amplitude of bar spectral line, N are sampled data length.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711259381.5A CN107945147A (en) | 2017-12-04 | 2017-12-04 | One kind is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711259381.5A CN107945147A (en) | 2017-12-04 | 2017-12-04 | One kind is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107945147A true CN107945147A (en) | 2018-04-20 |
Family
ID=61948506
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711259381.5A Pending CN107945147A (en) | 2017-12-04 | 2017-12-04 | One kind is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107945147A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111476761A (en) * | 2020-03-18 | 2020-07-31 | 河北科技大学 | Visibility measuring method and system based on system identification |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106526312A (en) * | 2016-10-10 | 2017-03-22 | 清华大学 | Electric energy metering method based on R-V(III) window FFT double-peak interpolation |
| CN106814246A (en) * | 2016-12-15 | 2017-06-09 | 贵州电网有限责任公司贵阳供电局 | Electric energy gauging method based on the Kaiser windows unimodal interpolation of FFT |
| CN106814230A (en) * | 2016-12-15 | 2017-06-09 | 贵州电网有限责任公司贵阳供电局 | Electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra- |
-
2017
- 2017-12-04 CN CN201711259381.5A patent/CN107945147A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106526312A (en) * | 2016-10-10 | 2017-03-22 | 清华大学 | Electric energy metering method based on R-V(III) window FFT double-peak interpolation |
| CN106814246A (en) * | 2016-12-15 | 2017-06-09 | 贵州电网有限责任公司贵阳供电局 | Electric energy gauging method based on the Kaiser windows unimodal interpolation of FFT |
| CN106814230A (en) * | 2016-12-15 | 2017-06-09 | 贵州电网有限责任公司贵阳供电局 | Electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra- |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111476761A (en) * | 2020-03-18 | 2020-07-31 | 河北科技大学 | Visibility measuring method and system based on system identification |
| CN111476761B (en) * | 2020-03-18 | 2023-07-21 | 河北科技大学 | Visibility measurement method and system based on system identification |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107991591A (en) | One kind is based on the modified image interfusion method of the unimodal interpolation of Kaiser windows FFT | |
| US8244489B2 (en) | Parametric multi-cycle averaging in an intelligent electronic device | |
| CN101807795B (en) | Method for forming electric energy metering simulation system and device thereof | |
| CN103926456B (en) | Flicker value calculation method based on improved FFT and inverter | |
| CN102288804B (en) | Method for calculating resistive current of arrester based on orthogonal transformation | |
| CN109061302A (en) | A kind of wind power generator incorporated in power network group harmonic measure system converted based on EEMD and Hilbert | |
| CN106814230A (en) | Electric energy gauging method based on the peak interpolations of Kaiser windows FFT tetra- | |
| Gallo et al. | Low cost smart power metering | |
| CN108169540A (en) | A kind of measuring method of wind power generating set voltage flicker | |
| CN108490257A (en) | A kind of short time-window inter-harmonic wave measuring method based on spectrum curve fitting | |
| CN106526312A (en) | Electric energy metering method based on R-V(III) window FFT double-peak interpolation | |
| CN105068035B (en) | A kind of voltage transformer error horizontal dynamic detection method and system | |
| CN103018709A (en) | Method and system for monitoring error of gateway metering device | |
| CN108896944B (en) | Laboratory calibrator of synchronous measuring device and synchronous phasor measuring method thereof | |
| CN109375065A (en) | Travelling wave identification method and positioning device based on three-dimensional grey absolute correlation degree | |
| CN106814246A (en) | Electric energy gauging method based on the Kaiser windows unimodal interpolation of FFT | |
| CN107945147A (en) | One kind is based on the tetra- modified image interfusion methods of peak interpolation of Kaiser windows FFT | |
| CN105021980B (en) | AC filter and breaker double-fracture voltage's distribiuting characteristic evaluation system and method | |
| Hamid et al. | Wavelet packet transform for RMS and power measurements | |
| CN104297604B (en) | Pilot system and method for a kind of harmonic wave on shunt capacitor temperature rise and loss influence | |
| Wu et al. | Effect of frequency offset on power measurement error in digital input electricity meters | |
| KR101664010B1 (en) | An Estimation Method of Line Parameter based on Synchrophasor Measurements in Power System | |
| Zhang et al. | A novelty digital algorithm for online measurement of dielectric loss factor of electronic transformers | |
| CN103247996B (en) | Compensation method for secondary current distortion caused by current transformer saturation | |
| CN109342816A (en) | The detection method of spectral leakage in electric energy quality monitoring |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180420 |
|
| WD01 | Invention patent application deemed withdrawn after publication |