CN107017636A - A kind of charging pile stress_responsive genes and filtering method - Google Patents
A kind of charging pile stress_responsive genes and filtering method Download PDFInfo
- Publication number
- CN107017636A CN107017636A CN201710138741.XA CN201710138741A CN107017636A CN 107017636 A CN107017636 A CN 107017636A CN 201710138741 A CN201710138741 A CN 201710138741A CN 107017636 A CN107017636 A CN 107017636A
- Authority
- CN
- China
- Prior art keywords
- subsegment
- harmonic
- frequency
- formula
- denoising
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/01—Arrangements for reducing harmonics or ripples
-
- 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
- G01R23/165—Spectrum analysis; Fourier analysis using filters
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/40—Arrangements for reducing harmonics
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Image Processing (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
The invention discloses a kind of charging pile stress_responsive genes and filtering method, before the present invention is tested harmonic signal, denoising first is carried out to harmonic signal, the waveform after denoising is more smooth, while reducing the decay of signal；The windows and interpolated FFT Harmonic Detecting Algorithm of the present invention, on the basis of Fourier transformation, is calculated using a spectral line of amplitude maximum, can suitably be made up the inaccurate defect of peak point measurement that short scope leakage is caused, be reduced detection error, improve Detection results；The present invention quickly follows harmonic signal using PWM converter, can quickly eliminate each harmonic detected.
Description
Art
The present invention relates to a kind of charging pile stress_responsive genes and filtering method.
Background technology
With the continuous intensification of global energy crisis, in exhausted and the atmosphere pollution, Global Temperature increasingly of petroleum resources
The harm aggravation risen, energy-conservation and emission reduction are the main directions of future automobile technology development.Electric automobile is used as traffic of new generation
Instrument, possesses the incomparable advantage of orthodox car in energy-saving and emission-reduction, the reduction mankind to the relying party face of traditional fossil energy.
When electric automobile rechargeable battery energy resource consumption to a certain extent when, it is necessary to using energy resource supply device to the battery
Charged.Current charging device mainly has two kinds of forms, and one kind is direct-current charging post, and the charging pile power is larger, and 100kW is left
The right side, the charging interval is short, and volume ratio is larger, therefore is typically mounted at fixed place；Another is alternating-current charging pile, is directly utilized
AC network, exports AC energy, and AC energy is converted to direct current energy and is by the vehicle-mounted charge stake carried by electric automobile
Rechargeable battery is charged.This kind of charged form power is smaller, and generally 10kW or so, the charging interval is long, small volume, therefore can
To make full use of each corner in city to be charged as electric automobile.
By the load of the charging pile and institute's band of charging pile use is non-linear equipment, therefore it can operationally give power network
The quality of power supply brings deleterious effect, is mainly reflected in terms of power factor of electric network declines and bring harmonic pollution to power network.
Especially, in the market starts a wide range of promote and fills stake soon at present, and it is newly-built it is substantial amounts of fill station soon, stake big work(a kind of is filled soon
Rate power electronic devices, will bring a large amount of harmonic waves to power network when in use.Theory analysis shows to fill stake production soon with actual test
Raw overtone order is mainly m=6k ± 1, and (k=1,2,3 ..., n), i.e., mainly produce 5,7,11,13,17,19 inferior harmonic waves.And
Harmonic amplitude is inversely proportional with overtone order, i.e., overtone order is higher, and harmonic content is lower.
To improve the quality of power supply, it is necessary to accurately detect the harmonic component of charging pile first, being possible to compensation, these are humorous
Wave component.Therefore, how harmonic component is quickly and easily detected, is always important study a question.Existing harmonic wave inspection
Survey scheme mainly has following four：(1) use after analog bandpass or bandstop filter, the amplified circuit of signal of input by each
Subharmonic wave filter, the centre frequency of each wave filter is the integral multiple of power frequency, and its advantage is that theory structure is simple, output impedance
Low, it is easy to control, realize, it has the disadvantage that centre frequency is vulnerable to the influence of extraneous factor, and unstable, precision is low；(2) based on small
The harmonic detecting of wave conversion, obtains optimal time resolution and frequency domain resolution, effectively from letter in the different parts of signal
Useful information is extracted in number, some indeterminable problems of Fourier transformation are can solve the problem that, small echo to sudden signal very
Sensitivity, it is possible to for detecting dynamic harmonic, and harmonic wave of the Fourier transformation commonly used to detection stable state；(3) instantaneous nothing
The harmonic detecting of work(Power Theory, when only detecting reactive current, completely no-delay can draw testing result, but examining
When surveying harmonic current, because the composition of harmonic current is different with wave filter, the delay at this moment caused is not more than a power supply
Cycle, so its real-time is very well, but its hardware is more, and cost is high；(4) harmonic detecting based on Fourier transformation, should
The method for planting harmonic detecting is current using a kind of most wide method, because its computational accuracy is higher, principle is simple, it is easy to connect
By, but the harmonic detecting of Fourier transformation needs certain time to complete, speed is slower, because its is computationally intensive, is
Overcome this point, it is necessary to consider the cooperation of amount of calculation and hardware, be certain section of finite digital signal is entered certainly during due to sampling
Row blocks sampling, so spectrum leakage will necessarily be so brought, then just can be directly by Fourier transformation if synchronized sampling
Obtain the amplitude of harmonic wave, frequency, phase, but often due to delay, reaction time of hardware device of AD samplings etc., synchronized sampling
It is difficult to accomplish to be, so will result in fence effect, influences the result of frequency analysis, particularly phase by influenceed can very
Greatly.
The content of the invention
To solve the above problems, the present invention provides a kind of charging pile stress_responsive genes and filtering method, with it, can
Charging pile net side each harmonic content is fast and accurately tested out, and can quickly filter out harmonic current in time, so as to effectively solve
Determine the problems of the prior art.
To achieve these goals, the present invention provides a kind of charging pile stress_responsive genes and filtering method, and this method is included such as
Lower step：
S1. under different operating modes, charging pile is sampled with power distribution network junction current waveform；
S2. current waveform is carried out, except processing of making an uproar, obtaining the waveform after denoising；
S3. utilize windows and interpolated FFT Harmonic Detecting Algorithm, to after denoising waveform carry out stress_responsive genes, obtain each time it is humorous
Wave number；
S4. each harmonic value is filtered out into each harmonic automatically as the input signal of filter.
It is preferred that, in step sl, under different operating modes, charging pile is adopted with power distribution network junction current waveform
Sample, definition gained sample is x1,x2,...,xn, sampling time interval is Δ t, then this n sampled point is divided into m son
Section, each subsegment includes 128 sampled points.
It is preferred that, in step s 2, according to the definition of fractal dimension, the Short Time Fractal Numbers of each subsegment are calculated, according to institute
Each subsegment Short Time Fractal Numbers obtained, calculate the wavelet threshold parameter of k-th of subsegment, and wherein constant term is needed according to actual noise
Variance chosen；Wavelet decomposition is carried out to each subsegment respectively, denoising is carried out using the threshold value tried to achieve, then carries out small
Reconstructed wave, obtains the waveform after denoising.
It is preferred that, in the step S1, definition gained sample is x1,x2,...,xn, sampling time interval is Δ
T, then this n sampled point is divided into m subsegment, each subsegment includes 128 sampled points, then m=n/128, for k-th of son
128 sampled points of section, are designated as sequence xi (k)(i=1,2 ..., 128；K=1,2 ..., m), definition
It is preferred that, in step s 2, the Short Time Fractal Numbers of each subsegment are calculated in the following way
Order
And
In formula (2), N(k)(Δ) and N(k)(2 Δ) expression covers this with the square net that width is Δ and 2 Δs respectively
Lattice number required for the functional image of segment signals, according to the definition of fractal dimension, calculates the Short-time fractal dimension of each subsegment
NumberFor：
It is preferred that, in step S2, according to each subsegment Short Time Fractal Numbers of gained, according to formulaMeter
Calculate the wavelet threshold parameter T of k-th of subsegment(k), wherein constant c needs to be chosen according to the variance of actual noise, respectively to every
Individual subsegment carries out wavelet decomposition, carries out denoising using tried to achieve threshold value, then carries out wavelet reconstruction, obtain the waveform after denoising.
It is preferred that, the step S3 specifically includes following steps：
S31. single-frequency signals x (t) frequency is set as f0, amplitude is A, and initial phase is θ, then uses fsSample frequency pair
Waveform after denoising is sampled, then obtains following discrete signal：
xn=Asin (2 π f0n/fs)+θ) (4)；
S32. adding window detection is carried out for above-mentioned discrete signal, if the time domain of institute's adding window is ω (n), its continuous frequency spectrum is W (2 π
F), then the continuous Fourier transform of the discrete signal is after adding window：
Ignore negative frequency-f0Locate the secondary lobe influence at frequency peak, then in positive frequency f0Neighbouring continuous frequency spectrum is：
Discrete sampling is carried out to above formula, DFT expression formula is obtained：
Discrete frequency intervals are Δ f=f in above formulas/ N, N are sampling number；
S33. the two piece spectral lines nearest apart from peak point are set as kth1And kth2Bar spectral line, it is clear that this two spectral lines are peak values
Point spectral line greatest around and time maximum, and one in k0Left side, one in k0Right side；Found in discrete spectrum
This two spectral lines, so as to can determine that k1And k2；
Make kth1And kth2The amplitude of bar spectral line is y respectively1=| X (k1Δ f) |, y2=| X (k2Δf)|；If f0=(k1+
λm) Δ f, | λ m |≤0.5, λ m are asynchronous degree；K1 is realized by finding spectral peak, in order to determine asynchronous degree λ m, makes k1With k2Spectrum
The ratio between line value is：
Synchronism deviation λ m are calculated by above formula, the updating formula of revised frequency, amplitude and phase is respectively：
Frequency correction formula：f0=(k1+λm)Δf；
It is worth updating formula：A=2 | X (ki)|/|W(2π·|λm|)/N)|；
Phasing formula：θ=arg | X (ki)|+π/2-λmπ；
I in amplitude rectification formula and phasing formula selects 1 or 2.
It is preferred that, in step s 4, specifically adopt and realize harmonic filtration with the following method：
S41. by the individual harmonic current size detected, as command signal, carried out with the output current of PWM converter
Compare, error size obtains one group of PWM ripple compared with the ring width of hysteresis comparator；
S42.PWM ripples are sent to the switch of the control end control power device of power device, follow individual harmonic current；
S43.PWM current transformers are equal in magnitude by individual harmonic current, and electric current in opposite direction is injected into distribution net side, with
The individual harmonic current included in grid side is cancelled out each other, so as to reach the purpose for eliminating distribution net side individual harmonic current.
The invention has the advantages that：(1) before testing harmonic signal, denoising first is carried out to harmonic signal, gone
Waveform after making an uproar is more smooth, while reducing the decay of signal；(2) windows and interpolated FFT Harmonic Detecting Algorithm of the invention,
On the basis of Fourier transformation, calculated using a spectral line of amplitude maximum, can suitably make up short scope leakage and cause
The inaccurate defect of peak point measurement, reduce detection error, improve Detection results；(3) quickly harmonic wave is followed to believe using PWM converter
Number, it can quickly eliminate each harmonic detected.
Brief description of the drawings
Fig. 1 shows a kind of flow chart of the charging pile stress_responsive genes and filtering method of the present invention.
Embodiment
Fig. 1 shows a kind of charging pile stress_responsive genes and filtering method, and this method comprises the following steps：
S1. under different operating modes, charging pile is sampled with power distribution network junction current waveform；
S2. current waveform is carried out, except processing of making an uproar, obtaining the waveform after denoising；
S3. utilize windows and interpolated FFT Harmonic Detecting Algorithm, to after denoising waveform carry out stress_responsive genes, obtain each time it is humorous
Wave number；
S4. each harmonic value is filtered out into each harmonic automatically as the input signal of filter.
It is preferred that, in step sl, under different operating modes, charging pile is adopted with power distribution network junction current waveform
Sample, definition gained sample is x1,x2,...,xn, sampling time interval is Δ t, then this n sampled point is divided into m son
Section, each subsegment includes 128 sampled points.
It is preferred that, in step s 2, according to the definition of fractal dimension, the Short Time Fractal Numbers of each subsegment are calculated, according to institute
Each subsegment Short Time Fractal Numbers obtained, calculate the wavelet threshold parameter of k-th of subsegment, and wherein constant term is needed according to actual noise
Variance chosen；Wavelet decomposition is carried out to each subsegment respectively, denoising is carried out using the threshold value tried to achieve, then carries out small
Reconstructed wave, obtains the waveform after denoising.
It is preferred that, in the step S1, definition gained sample is x1,x2,...,xn, sampling time interval is Δ
T, then this n sampled point is divided into m subsegment, each subsegment includes 128 sampled points, then m=n/128, for k-th of son
128 sampled points of section, are designated as sequence xi (k)(i=1,2 ..., 128；K=1,2 ..., m), definition
It is preferred that, in step s 2, the Short Time Fractal Numbers of each subsegment are calculated in the following way
Order
And
In formula (2), N(k)(Δ) and N(k)(2 Δ) expression covers this with the square net that width is Δ and 2 Δs respectively
Lattice number required for the functional image of segment signals, according to the definition of fractal dimension, calculates the Short-time fractal dimension of each subsegment
NumberFor：
It is preferred that, in step S2, according to each subsegment Short Time Fractal Numbers of gained, according to formulaMeter
Calculate the wavelet threshold parameter T of k-th of subsegment(k), wherein constant c needs to be chosen according to the variance of actual noise, respectively to every
Individual subsegment carries out wavelet decomposition, carries out denoising using tried to achieve threshold value, then carries out wavelet reconstruction, obtain the waveform after denoising.
It is preferred that, the step S3 specifically includes following steps：
S31. single-frequency signals x (t) frequency is set as f0, amplitude is A, and initial phase is θ, then uses fsSample frequency pair
Waveform after denoising is sampled, then obtains following discrete signal：
xn=Asin (2 π f0n/fs)+θ) (4)；
S32. adding window detection is carried out for above-mentioned discrete signal, if the time domain of institute's adding window is ω (n), its continuous frequency spectrum is W (2 π
F), then the continuous Fourier transform of the discrete signal is after adding window：
Ignore negative frequency-f0Locate the secondary lobe influence at frequency peak, then in positive frequency f0Neighbouring continuous frequency spectrum is：
Discrete sampling is carried out to above formula, DFT expression formula is obtained：
Discrete frequency intervals are Δ f=f in above formulas/ N, N are sampling number；
S33. the two piece spectral lines nearest apart from peak point are set as kth1And kth2Bar spectral line, it is clear that this two spectral lines are peak values
Point spectral line greatest around and time maximum, and one in k0Left side, one in k0Right side；Found in discrete spectrum
This two spectral lines, so as to can determine that k1And k2；
Make kth1And kth2The amplitude of bar spectral line is y respectively1=| X (k1Δ f) |, y2=| X (k2Δf)|；If f0=(k1+
λm) Δ f, | λ m |≤0.5, λ m are asynchronous degree；K1 is realized by finding spectral peak, in order to determine asynchronous degree λ m, makes k1With k2Spectrum
The ratio between line value is：
Synchronism deviation λ m are calculated by above formula, the updating formula of revised frequency, amplitude and phase is respectively：
Frequency correction formula：f0=(k1+λm)Δf；
It is worth updating formula：A=2 | X (ki)|/|W(2π·|λm|)/N)|；
Phasing formula：θ=arg | X (ki)|+π/2-λmπ；
I in amplitude rectification formula and phasing formula selects 1 or 2.
It is preferred that, in step s 4, specifically adopt and realize harmonic filtration with the following method：
S41. by the individual harmonic current size detected, as command signal, carried out with the output current of PWM converter
Compare, error size obtains one group of PWM ripple compared with the ring width of hysteresis comparator；
S42.PWM ripples are sent to the switch of the control end control power device of power device, follow individual harmonic current；
S43.PWM current transformers are equal in magnitude by individual harmonic current, and electric current in opposite direction is injected into distribution net side, with
The individual harmonic current included in grid side is cancelled out each other, so as to reach the purpose for eliminating distribution net side individual harmonic current.
Using the method for the present invention, can control the distribution network voltage total harmonic distortion factor threshold values is 5.0%, odd harmonic
Voltage containing ratio threshold values is 4.0%, and even harmonic voltages containing ratio threshold values is 2.0%.The determination of these numerical value is establishing criteria
Number be GB_T_14549-1993 national standard, title be quality of power supply utility network harmonic wave.
Above content is to combine specific preferred embodiment further description made for the present invention, it is impossible to assert
The specific implementation of the present invention is confined to these explanations.For general technical staff of the technical field of the invention,
On the premise of not departing from present inventive concept, some equivalent substitutes or obvious modification are made, and performance or purposes are identical, all should
It is considered as belonging to protection scope of the present invention.
Claims (8)
1. a kind of charging pile stress_responsive genes and filtering method, this method comprise the following steps：
S1. under different operating modes, charging pile is sampled with power distribution network junction current waveform；
S2. current waveform is carried out, except processing of making an uproar, obtaining the waveform after denoising；
S3. windows and interpolated FFT Harmonic Detecting Algorithm is utilized, stress_responsive genes are carried out to the waveform after denoising, each harmonic value is obtained；
S4. each harmonic value is filtered out into each harmonic automatically as the input signal of filter.
2. the method as described in claim 1, it is characterised in that in step sl, under different operating modes, to charging pile and distribution
Net junction current waveform is sampled, and definition gained sample is x1,x2,...,xn, sampling time interval is Δ t, then will
This n sampled point is divided into m subsegment, and each subsegment includes 128 sampled points.
3. the method as described in claim 1, it is characterised in that in step s 2, according to the definition of fractal dimension, calculates each son
The Short Time Fractal Numbers of section, according to each subsegment Short Time Fractal Numbers of gained, calculate the wavelet threshold parameter of k-th of subsegment, its
Middle constant term needs to be chosen according to the variance of actual noise；Wavelet decomposition is carried out to each subsegment respectively, using being tried to achieve
Threshold value carry out denoising, then carry out wavelet reconstruction, obtain the waveform after denoising.
4. method as claimed in claim 2 or claim 3, it is characterised in that in the step S1, definition gained sample is x1,
x2,...,xn, sampling time interval is Δ t, then this n sampled point is divided into m subsegment, and each subsegment includes 128 samplings
Point, then m=n/128, for 128 sampled points of k-th of subsegment, is designated as sequence xi (k)(i=1,2 ..., 128；K=
1,2 ..., m), definition
5. method as claimed in claim 4, it is characterised in that in step s 2, calculates each subsegment in the following way
Short Time Fractal Numbers
Order
And
In formula (2), N(k)(Δ) and N(k)(2 Δ) expression covers this section with the square net that width is Δ and 2 Δs respectively
Lattice number required for the functional image of signal, according to the definition of fractal dimension, calculates the Short Time Fractal Numbers of each subsegment
For：
6. method as claimed in claim 5, it is characterised in that in step S2, according to each subsegment Short Time Fractal Numbers of gained,
According to formulaCalculate the wavelet threshold parameter T of k-th of subsegment(k), wherein constant c needs are according to actual noise
Variance chosen, respectively to each subsegment carry out wavelet decomposition, using tried to achieve threshold value carry out denoising, then carry out small echo
Reconstruct, obtains the waveform after denoising.
7. method as claimed in claim 6, it is characterised in that the step S3 specifically includes following steps：
S31. single-frequency signals x (t) frequency is set as f0, amplitude is A, and initial phase is θ, then uses fsSample frequency to denoising
Waveform afterwards is sampled, then obtains following discrete signal：
xn=Asin (2 π f0n/fs)+θ) (4)；
S32. adding window detection is carried out for above-mentioned discrete signal, if the time domain of institute's adding window is ω (n), its continuous frequency spectrum is W (2 π f),
Then the continuous Fourier transform of the discrete signal is after adding window：
Ignore negative frequency-f0Locate the secondary lobe influence at frequency peak, then in positive frequency f0Neighbouring continuous frequency spectrum is：
Discrete sampling is carried out to above formula, DFT expression formula is obtained：
Discrete frequency intervals are Δ f=f in above formulas/ N, N are sampling number；
S33. the two piece spectral lines nearest apart from peak point are set as kth1And kth2Bar spectral line, it is clear that this two spectral lines are that peak point is attached
Near maximum and time maximum spectral line, and one in k0Left side, one in k0Right side；Found in discrete spectrum this two
Bar spectral line, so as to can determine that k1And k2；
Make kth1And kth2The amplitude of bar spectral line is y respectively1=| X (k1Δ f) |, y2=| X (k2Δf)|；If f0=(k1+λm)Δ
F, | λ m |≤0.5, λ m are asynchronous degree；K1 is realized by finding spectral peak, in order to determine asynchronous degree λ m, makes k1With k2Spectral line value
The ratio between be：
Synchronism deviation λ m are calculated by above formula, the updating formula of revised frequency, amplitude and phase is respectively：
Frequency correction formula：f0=(k1+λm)Δf；
It is worth updating formula：A=2 | X (ki)|/|W(2π·|λm|)/N)|；
Phasing formula：θ=arg | X (ki)|+π/2-λmπ；
I in amplitude rectification formula and phasing formula selects 1 or 2.
8. the method as described in claim 1-7 is any, it is characterised in that in step s 4, specifically adopts and realizes with the following method
Harmonic filtration：
S41. by the individual harmonic current size detected, as command signal, compared with the output current of PWM converter
Compared with error size obtains one group of PWM ripple compared with the ring width of hysteresis comparator；
S42.PWM ripples are sent to the switch of the control end control power device of power device, follow individual harmonic current；
S43.PWM current transformers are equal in magnitude by individual harmonic current, and electric current in opposite direction is injected into distribution net side, with power distribution network
The individual harmonic current included in side is cancelled out each other, so as to reach the purpose for eliminating distribution net side individual harmonic current.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710138741.XA CN107017636B (en) | 2017-03-09 | 2017-03-09 | Charging pile harmonic testing and filtering method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710138741.XA CN107017636B (en) | 2017-03-09 | 2017-03-09 | Charging pile harmonic testing and filtering method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107017636A true CN107017636A (en) | 2017-08-04 |
| CN107017636B CN107017636B (en) | 2021-07-20 |
Family
ID=59440237
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710138741.XA Active CN107017636B (en) | 2017-03-09 | 2017-03-09 | Charging pile harmonic testing and filtering method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107017636B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113675852A (en) * | 2021-08-25 | 2021-11-19 | 安徽信息工程学院 | Harmonic suppression method and system in electric vehicle charging process in grid-connected state |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103207319A (en) * | 2013-03-12 | 2013-07-17 | 广东电网公司电力科学研究院 | Harmonic wave measurement method of electricity signal of digital substation under non-synchronous sampling condition |
| CN105512469A (en) * | 2015-12-01 | 2016-04-20 | 江苏省电力公司淮安供电公司 | Charging pile harmonic wave detection algorithm based on windowing interpolation FFT and wavelet packet |
| CN106230081A (en) * | 2016-07-01 | 2016-12-14 | 李红彪 | The harmonic detecting of a kind of charging pile and filtering method |
-
2017
- 2017-03-09 CN CN201710138741.XA patent/CN107017636B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103207319A (en) * | 2013-03-12 | 2013-07-17 | 广东电网公司电力科学研究院 | Harmonic wave measurement method of electricity signal of digital substation under non-synchronous sampling condition |
| CN105512469A (en) * | 2015-12-01 | 2016-04-20 | 江苏省电力公司淮安供电公司 | Charging pile harmonic wave detection algorithm based on windowing interpolation FFT and wavelet packet |
| CN106230081A (en) * | 2016-07-01 | 2016-12-14 | 李红彪 | The harmonic detecting of a kind of charging pile and filtering method |
Non-Patent Citations (1)
| Title |
|---|
| 姚强: "基于FPGA的电力系统谐波检测系统的研究与设计", 《万方硕士学位论文》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113675852A (en) * | 2021-08-25 | 2021-11-19 | 安徽信息工程学院 | Harmonic suppression method and system in electric vehicle charging process in grid-connected state |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107017636B (en) | 2021-07-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101806832B (en) | Measuring method for frequencies of low-frequency signals | |
| CN102662106B (en) | Method for measuring electric energy of harmonic power network | |
| CN102288807A (en) | Method for measuring electric network voltage flicker | |
| CN104236702A (en) | System and method for judging interior looseness of power transformer | |
| CN101587146A (en) | Parameter detection method for harmonic wave and indirect harmonic wave | |
| CN103116064A (en) | Method and device for detecting voltage fluctuation and flicker based on energy operator and spectrum correction | |
| CN108318852B (en) | Square wave influence test method for intelligent electric energy meter | |
| CN101701984A (en) | Fundamental wave and harmonic wave detecting method based on three-coefficient Nuttall windowed interpolation FFT | |
| CN103575984A (en) | Harmonic analysis method based on Kaiser window double-spectral-line interpolation FFT | |
| CN102590618B (en) | Detection method of positive sequence voltage phase of fundamental wave for power grid | |
| CN201773140U (en) | Alternating current and direct current hybrid electric field intensity measuring device | |
| CN107390022B (en) | Electric energy metering method based on discrete spectrum correction | |
| CN109061345B (en) | Effective value measuring method and system suitable for power system | |
| CN206725655U (en) | A kind of electric car BMS Insulation Inspection Systems | |
| CN103399204A (en) | Rife-Vincent (II) window interpolation FFT (Fast Fourier Transform)-based harmonic and inter-harmonic detection method | |
| CN105137181A (en) | Double-spectrum-line interpolation harmonic analysis algorithm based on Nuttall-Kaiser composite window | |
| CN105486921A (en) | Kaiser third-order mutual convolution window triple-spectrum-line interpolation harmonic wave and inter-harmonic wave detection method | |
| CN101718816B (en) | Fundamental wave and harmonic wave detection method based on four-item coefficient Nuttall window interpolation FFT | |
| CN114966194A (en) | ADC-based low-frequency metering three-phase electric energy meter | |
| CN108414838A (en) | A kind of inverter parallel system line impedance measurement method | |
| CN109581045B (en) | Inter-harmonic power metering method meeting IEC standard framework | |
| CN107017636A (en) | A kind of charging pile stress_responsive genes and filtering method | |
| LU101400B1 (en) | High-precision estimation method and apparatus for carrier frequency of amplitude modulation signal | |
| CN203672975U (en) | Harmonic-impedance measuring device for power system | |
| CN112485522B (en) | Electric energy data perception-based flat-top window function synchronous phasor measurement method and device |
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 | ||
| CB02 | Change of applicant information |
Address after: No. 12002-2021-1, 20th floor, unit 1, building 1, block C, city gate, Tangyan South Road, Zhangba Street office, high tech Zone, Xi'an, Shaanxi 710000 Applicant after: Chuang Rui equipment intelligent manufacturing (Xi'an) Co.,Ltd. Address before: 523170 room 06, 11th floor, Goldman Sachs technology building, No.5, Longxi Road, Zhouxi, Nancheng District, Dongguan City, Guangdong Province Applicant before: GUANGDONG TOPFER ELECTRONIC TECHNOLOGY Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20210622 Address after: 542400 No.5, East Industrial Avenue, Ertang Comprehensive Industrial Park, industrial concentration area, Ertang Town, Pingle County, Guilin City, Guangxi Zhuang Autonomous Region Applicant after: Guilin chuangrei Equipment Co.,Ltd. Address before: No. 12002-2021-1, 20th floor, unit 1, building 1, block C, city gate, Tangyan South Road, Zhangba Street office, high tech Zone, Xi'an, Shaanxi 710000 Applicant before: Chuang Rui equipment intelligent manufacturing (Xi'an) Co.,Ltd. |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220919 Address after: Room 5072, Unit 1, Building 1, No. 2, Nancheng Section, Tiyu Road, Nancheng Street, Dongguan City, Guangdong 523000 Patentee after: DONGGUAN CHUANGRUI NEW ENERGY Co.,Ltd. Address before: 542400 No.5, East Industrial Avenue, Ertang Comprehensive Industrial Park, industrial concentration area, Ertang Town, Pingle County, Guilin City, Guangxi Zhuang Autonomous Region Patentee before: Guilin chuangrei Equipment Co.,Ltd. |
|
| TR01 | Transfer of patent right |