CN109164298A - A kind of compressed sensing based ultra harmonics detection device and detection method - Google Patents
A kind of compressed sensing based ultra harmonics detection device and detection method Download PDFInfo
- Publication number
- CN109164298A CN109164298A CN201810826416.7A CN201810826416A CN109164298A CN 109164298 A CN109164298 A CN 109164298A CN 201810826416 A CN201810826416 A CN 201810826416A CN 109164298 A CN109164298 A CN 109164298A
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- low
- pass filter
- sampling
- 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
Classifications
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
A kind of compressed sensing based ultra harmonics detection device and detection method, sensing circuit is mounted on power distribution network measurement end, the signal of sensing circuit output first passes through high-pass filter, low-pass filter is sent into after being mixed secondly with the pseudo-random sequence that multifunctional data acquisition card D/A module exports, it is realized again by multifunctional data acquisition card A/D module and the low speed of signal is sampled, last signal is uploaded to the reconstruct and display of the complete pair signals of host computer by usb bus.High-pass filter is designed as adjustable gain, in addition to filtering out power frequency and the interference of other low frequency signals, can also be achieved the conditioning to signal.Low pass filter design is that cutoff frequency is adjustable, then Ratios detection can be achieved.The synchronization trigger module of multifunctional data acquisition card can realize that D/A is synchronous with A/D, construct observing matrix when phase deviation is 0, improve reconstruction accuracy.What the invention had is clearly advantageous that: in the case where ensuring ultra harmonics and its correct situation of frequency spectrum reconfiguration, the sample frequency and sampling number of measuring device is greatly reduced.
Description
Technical field
The present invention relates to compressed sensing and harmonic detectings, and in particular to a kind of compressed sensing based ultra harmonics inspection
Survey device and detection method.
Background technique
The important infrastructure that power distribution network develops as economy and society fights to realization smart grid and energy internet
Slightly target plays key effect.Ultra harmonics is the new electric energy quality problems that China's New Generation of Intelligent power distribution network derives,
It can carry out many detrimental effects to electric distribution network electrical equipment and communication band, cause the great attention of industry and academia.
In June, 2017, International Conference on Electricity Distribution CIRED specially set up innovation forum and technology lecture for ultra harmonics problem.
Currently, the harmonic measure instrument harmonic measure number of the market mainstream is not more than 50 times, the anti-spectral aliasing electricity of instrument
Road makes it be unable to measure the harmonic component greater than 50 times or more.By Nyquist sampling thheorem it is found that meet ultra harmonics
Measurement, it is desirable that the pendulous frequency of harmonic measure instrument be greater than 100 times or more, this just to the ADC sample rate of measuring instrument, adopt
More stringent requirements are proposed for sample digit, frequency overlapped-resistable filter order and cutoff frequency, and it is real to increase harmonic measure instrument hardware
Ready-made.Meanwhile because in real time, the sampling of high-speed easily lead to the generation of mass data, further initiation data it is computationally intensive,
A series of problems, such as storage is difficult, transmission efficiency is low.
The compressive sensing theory proposed by Candes and Donoho et al. in 2004, breaches Nyquist sampling thheorem
Limitation.Its thought are as follows: signal x ∈ RN×1Pass through basic function Ψ ∈ RN×NIt is transformed to θ=ΨTX, if θ is K- sparse in the domain Ψ
(nonzero coefficient K < < N).Then using one and the incoherent observing matrix Φ ∈ R of ΨM×N(M < < N) executes compression to x and adopts
Compressed signal y=Φ x ∈ R then can be obtained in sampleM×1, which contains enough information when reconstruction signal.But it is directly extensive from y
There are infinite multiple solutions by multiple x, therefore y can be regarded as to y=Φ x=Φ Ψ θ=A θ, wherein A=Φ Ψ ∈ RM×NReferred to as sense square
Battle array.Since θ is K- sparse, then θ is first recovered from y becomes the committed step of reconstruct x, and restructuring procedure is under 0- norm meaning
Optimization problem model solution.
Summary of the invention
The object of the present invention is to provide a kind of compressed sensing based ultra harmonics detection device and detection method, needles
The problems of detection method and technology to aforementioned existing harmonic measure instrument, and compressed in conjunction with compressive sensing theory
Sampling and the progress that has obtained of reconstruct aspect, make up conventional harmonic measuring instrument ultra harmonics signal detection field not
Be enough and compressed sensing harmonic detecting application in blank.Ultra harmonics has the characteristics that frequency-domain sparse signal, therefore
It can realize that Sampling Compression and reconstruct, Sampling Compression frequency depend no longer on the most high frequency of signal by compressive sensing theory
Rate ingredient and bandwidth, but depend on the degree of rarefication of ultra harmonics.The detection method for using compressed sensing simultaneously, can will count
The restructuring procedure that complexity is transferred to background computer data from collection process is calculated, is effectively prevented because of high-speed sampling generation
A series of problems.
To achieve the goals above, the technical solution adopted by the present invention is that:
A kind of compressed sensing based ultra harmonics detection device, including sensing circuit, the high-pass filtering of adjustable gain
Device, mixing filter, multifunctional data acquisition card, host computer, sensing circuit connect adjustable gain by voltage and current signal
The high-pass filter of high-pass filter, adjustable gain connects mixing filter, frequency mixer, the low-pass filter of mixing filter
Successively signal is connect with the A/D module of multifunctional data acquisition card, D/A module, the synchronous triggering of multifunctional data acquisition card
Module is connect with A/D module, D/A module by signal, and multifunctional data acquisition card connects host computer by USB；
The sensing circuit is arranged in power distribution network measurement end, including voltage sensor and current sensor；
The high-pass filter of the adjustable gain is used for signal condition, on the one hand can filter out power frequency in signal and its
On the other hand its low-frequency component plays the role of signal amplitude conditioning, is the simulation of gain, cutoff frequency, regulable center frequency
Filter；
The mixing filter is made of frequency mixer and low-pass filter, low pass filter cutoff frequency, centre frequency
It is adjustable, Ratios detection can be played the role of；
The multifunctional data acquisition card is generated for random sequence, low speed samples, synchronous triggering and usb data pass
It is defeated；
The A/D module is sampled for low speed, D/A module for generating pseudo-random sequence, synchronous triggering for realizing
D/A, A/D two-way press rising edge and failing edge triggering respectively, and USB module is used for data-transformation facility.
The host computer is for realizing the reconstruct and display to ultra harmonics time frequency signal, it can be achieved that pseudorandom sequence
The multiple tasks such as column generation, data acquisition, signal reconstruction and display.
A kind of compressed sensing based ultra harmonics detection method, comprising the following steps:
Step 1) sensing circuit is mounted on power distribution network measurement end, includes voltage transformer and current transformer, can be by distribution
Big voltage, current signal in net are converted to small voltage, current signal and play the role of being isolated with power grid；
The signal x (t) of step 2) sensing circuit output send the high-pass filter to adjustable gain to obtain filtered signal
x1(t), the effect of high-pass filter can play power frequency and the other low-frequency components filtered out in signal, only retain ultra harmonics
The signal of bandwidth ingredient；
Step 3) uses the synchronization trigger source on multifunctional data acquiring board, pseudorandom in guarantee step 4), step 5)
Sequence output channel AO is synchronous with low speed sampling input channel AI, it is ensured that is not in phase deviation situation；
Step 4) is by signal x1(t) it sending into mixing filter, mixing filter includes multiplier and low-pass filter,
Wherein pseudo-random sequence pc(t) it is exported by the AO mouth of multifunctional data acquisition card, x1(t) and pc(t) it is first multiplied in multiplier
Signal y (t) is obtained, according to the Convolution Properties of continuous fourier transform:
FT [y (t)]=FT [x1(t)pc(t)]=X1(j Ω) * P (j Ω)=Y (j Ω)
It is found that the frequency spectrum of y (t) is x1(t)、pc(t) convolution of frequency spectrum, since the spectrum distribution of pseudo-random sequence is entire
In frequency domain, therefore can will be at the frequency spectrum shift tremendously low frequency of ultra harmonics by convolutional calculation；
Then, y (t) obtains signal y by low-pass filter again1(t), it can be indicated by the process of low-pass filter
Are as follows:
H (t) is the unit impulse response of low-pass filter；
After passband maximum attenuation is given, the cut-off frequecy of passband f of the low-pass filter0It can be according to required Sampling Compression
Than being manually adjusted, Ratios filtering is realized；
Step 5) is by signal y1(t) it send to the channel AI of multifunctional data acquisition card 4 and carries out low speed A/D sampling, sampled
Journey is as follows:
Sampling time interval T is by sample frequency fsIt determines.Enable fs≥2f0, as the cut-off frequecy of passband f of low-pass filter0
When lower, sample frequency fsLower, sampling time interval T is bigger, and sampling number N is fewer, it can realizes to ultra harmonics
Ratios compression sampling；
Step 6) is by y1(n) it is sent into host computer, weight is carried out to compression sampling signal using orthogonal matching pursuit OMP algorithm
Structure, recover in the frequency spectrum θ, θ of signal first for 0 coefficient be exactly ultra harmonics spectrum component, secondly θ is passed through
Reconstruction signal x ' can be obtained after FFT inverse transformation1(n)；
Ultra harmonics time-frequency waveform after step 7) display reconstruct.
The signal of the described step 6) output is reconstructed that specific step is as follows:
Step A: according to analog information converter AIC actual circuit, constructing random sequence matrix P, low-pass filtering link
Unit impulse response matrix H, the sampling matrix B of low speed sampling element；
Step B: observing matrix Φ=BHP is obtained；
Step C: construction discrete Fourier transform basic matrix Ψ；
Step D: sensing matrix A=Φ Ψ is obtained；
Step E: the frequency spectrum θ of ultra harmonics is reconstructed using orthogonal matching pursuit algorithm, and inverse by discrete fourier
Transformation recovers ultra harmonics signal.
Physically realizable observing matrix is analog information converter AIC, and the low pass in AIC in the step 3)
Filter be gain, cutoff frequency, regulable center frequency analog filter.
The advantages and positive effects of the present invention are as follows:
1) the method for the present invention has the characteristics that high frequency, sparse using ultra harmonics signal, using compressive sensing theory and
Technology is realized to the compression sampling of ultra harmonics, reconstruct and analysis, can not only realize the correct inspection to ultra harmonics
It surveys, is also greatly reduced the sample frequency and sampling number of detection device.
2) high-pass filter in the present invention program is designed as adjustable gain, can not only effectively reduce power frequency and other
The interference of low frequency signal, moreover it is possible to realize the conditioning to ultra harmonics signal.Synchronization trigger module in the present invention program can be real
Existing D/A is synchronous with A/D, constructs observing matrix when phase deviation is 0, improves reconstruction accuracy.
3) detection method of the invention and scheme can also extend in the high-frequency signal detection in other fields.If will inspection
High-pass filter in survey method and in implementation removes, and could be applicable to the inspection of transient power quality signal in power distribution network
It surveys.
Detailed description of the invention
Fig. 1 is the implementation flow chart of detection method；
Fig. 2 is the structural schematic diagram of implementation of the present invention,
Wherein, 1 is sensing circuit；2 be the high-pass filter of adjustable gain；3 be mixing filter；4 be multi-functional data
Capture card；5 be host computer；
Fig. 3 is the timing diagram that phase deviation is 0 and phase deviation is 1；
Fig. 4 is the reconstruction accuracy under out of phase deviation, different cutoff frequencies；
Fig. 5 is the waveform by surveying after high-pass filter；
Fig. 6 is the waveform by surveying after frequency mixer.
Fig. 7 is the waveform diagram after low-pass filtering.
Fig. 8 is the LabVIEW software display figure that the waveform inputted and its amplitude spectrum are shown in PC machine.
Fig. 9 is the waveform and its amplitude spectrum after low speed sampling.
Figure 10 is the waveform and its amplitude spectrum reconstructed.
Specific embodiment
The present invention is further discussed below below in conjunction with attached drawing and implementation benefit.
As shown in Fig. 2, a kind of compressed sensing based ultra harmonics detection device, including sensing circuit 1, gain can
The high-pass filter 2 of tune, mixing filter 3, multifunctional data acquisition card 4, host computer 5, sensing circuit 1 pass through voltage and current
Signal connects the high-pass filter 2 of adjustable gain, and the high-pass filter 2 of adjustable gain connects mixing filter 3, mixing and filtering
Successively signal connects for the frequency mixer of device 3, the A/D module of low-pass filter and multifunctional data acquisition card 4, D/A module, more function
The module of setting out that synchronizes of energy data collecting card 4 is connect with A/D module, D/A module by signal, and multifunctional data acquisition card 4 passes through
USB connection host computer 5；
The sensing circuit 1 is arranged in power distribution network measurement end, including voltage sensor and current sensor；
The high-pass filter 2 of the adjustable gain is used for signal condition, on the one hand can filter out power frequency in signal and its
On the other hand its low-frequency component plays the role of signal amplitude conditioning, is the simulation of gain, cutoff frequency, regulable center frequency
Filter；
The mixing filter 3 is made of frequency mixer and low-pass filter, low pass filter cutoff frequency, center frequency
Rate is adjustable, can play the role of Ratios detection；
The multifunctional data acquisition card 4 is generated for random sequence, low speed samples, synchronous triggering and usb data pass
It is defeated；
The A/D module is sampled for low speed, D/A module for generating pseudo-random sequence, synchronous triggering for realizing
D/A, A/D two-way press rising edge and failing edge triggering respectively, and USB module is used for data-transformation facility.
The host computer 5 is for realizing the reconstruct and display to ultra harmonics time frequency signal, it can be achieved that pseudorandom sequence
The multiple tasks such as column generation, data acquisition, signal reconstruction and display.
As shown in Figure 1, 2, a kind of compressed sensing based ultra harmonics detection method, comprising the following steps:
Step 1), sensing circuit 1 are mounted on power distribution network measurement end, include voltage transformer and current transformer, can will match
Big voltage, current signal in power grid are converted to small voltage, current signal and play the role of being isolated with power grid；
The signal x (t) that step 2) sensing circuit 1 exports send the high-pass filter 2 to adjustable gain to obtain filtered letter
Number x1(t), the effect of high-pass filter can play power frequency and the other low-frequency components filtered out in signal, and it is humorous only to retain superelevation time
The signal of the wide ingredient of wavestrip；
Step 3) using the synchronization trigger source on multifunctional data acquiring board 4, guarantee step 4), in step 5) puppet with
Machine sequence output channel AO is synchronous with low speed sampling input channel AI, it is ensured that is not in phase deviation situation；
Step 4) is by signal x1(t) it send into mixing filter 3, mixing filter 3 includes multiplier and low-pass filtering
Device, wherein pseudo-random sequence pc(t) it is exported by the AO mouth of multifunctional data acquisition card 4, x1(t) and pc(t) first in multiplier
Multiplication obtains signal y (t), according to the Convolution Properties of continuous fourier transform:
FT [y (t)]=FT [x1(t)pc(t)]=X1(j Ω) * P (j Ω)=Y (j Ω)
It is found that the frequency spectrum of y (t) is x1(t)、pc(t) convolution of frequency spectrum, since the spectrum distribution of pseudo-random sequence is entire
In frequency domain, therefore can will be at the frequency spectrum shift tremendously low frequency of ultra harmonics by convolutional calculation；
Then, y (t) obtains signal y by low-pass filter again1(t), it can be indicated by the process of low-pass filter
Are as follows:
H (t) is the unit impulse response of low-pass filter；
After passband maximum attenuation is given, the cut-off frequecy of passband f of the low-pass filter0It can be according to required Sampling Compression
Than being manually adjusted, Ratios filtering is realized；
Step 5) is by signal y1(t) it send to the channel AI of multifunctional data acquisition card 4 and carries out low speed A/D sampling, sampled
Journey is as follows:
Sampling time interval T is by sample frequency fsIt determines.Enable fs≥2f0, as the cut-off frequecy of passband f of low-pass filter0
When lower, sample frequency fsLower, sampling time interval T is bigger, and sampling number N is fewer, it can realizes to ultra harmonics
Ratios compression sampling；
Step 6) is by y1(n) it is sent into host computer, weight is carried out to compression sampling signal using orthogonal matching pursuit OMP algorithm
Structure, recover in the frequency spectrum θ, θ of signal first for 0 coefficient be exactly ultra harmonics spectrum component, secondly θ is passed through
Reconstruction signal x can be obtained after FFT inverse transformation1(n)；
Ultra harmonics time-frequency waveform after step 7) display reconstruct.
The signal of the described step 6) output is reconstructed that specific step is as follows:
Step A: it is analog information converter AIC actual circuit according to physically realizable observing matrix, constructs stochastic ordering
Column matrix P, the unit impulse response matrix H of low-pass filtering link, the sampling matrix B of low speed sampling element；
Step B: observing matrix Φ=BHP is obtained；
Step C: construction Fourier transformation basic matrix Ψ；
Step D: sensing matrix A=Φ Ψ is obtained；
Step E: the frequency spectrum θ of ultra harmonics is reconstructed using orthogonal matching pursuit algorithm, and inverse by discrete fourier
Transformation recovers ultra harmonics signal.
Physically realizable observing matrix is analog information converter AIC, and the low pass in AIC in the step 3)
Filter be gain, cutoff frequency, regulable center frequency analog filter.
Embodiment one:
The embodiment influences to clearly illustrate phase deviation caused by signal reconstruction precision and low-pass filter
Influence of the cut-off frequecy of passband to signal reconstruction precision.It is the driving clock for generating pseudo-random sequence, p referring to Fig. 3, CLK_mc
It (t) is the pseudo-random sequence generated, CLK_AD is the low speed sampling clock that phase deviation is 0, and change frequency is driving clock
The 1/5 of frequency, first value sampled at this time is first value of pseudo-random sequence initial time, and is successively spaced 5
The driving clock cycle is sampled, as shown in red curve in Fig. 3.When there is no phase deviation, when reconstruct constructed it is pseudo- with
Machine sequence matrix P, low speed sampling matrix B are as follows:
P=diag (pc(0), pc(1) ..., pc(N-1))
CLD_AD1 is the low speed sampling clock that phase deviation is 1, and first of sampling value is the of m-sequence initial time
Two value, and successively be spaced 5 driving clock cycle sampled, as shown in blue curve.When phase deviation is 1, participate in
Corresponding change can occur for the pseudo-random sequence matrix P and low speed sampling matrix B of reconstruct:
P=diag (pc(1) ..., pc(N-1), pc(0))
This variation can cause observing matrix Φ also to change.If observing matrix Φ when signal reconstruction is according to nothing
It is constructed when phase deviation, when there are deviation, will lead to observing matrix corresponding to actual hardware and reconstruction model
Used in observing matrix it is inconsistent, reconstruction accuracy is had an impact.Referring to shown in Fig. 4, phase deviation is bigger, after reconstruct
Signal and original signal between residual error it is bigger, and when there is no phase deviation, reduce the passband cut-off of low-pass filter
Frequency does not influence reconstruction accuracy.The embodiment shows specific implementation step of the present invention 3) use Ratios filtering side
Method will not influence reconstruction accuracy, and the sample frequency of low speed sampling can be greatly reduced；Specific implementation step of the present invention (5) uses
Synchronization trigger method can ensure generate random sequence driving clock and low speed sampling when sampling clock there is not phase
Position offset issue.
Embodiment two:
The embodiment simulates the signal generated in practical power distribution network with a signal generator, contains power frequency in the signal
50Hz and 1kHz harmonic signal.Using detection method proposed by the present invention and implementation shown in Fig. 1,2, a set of base has been built
In the test platform of virtual instrument.The platform multiplier selects AD835, and high-pass filter and low-pass filter selection UAF42 have
Source filter, multifunctional data acquisition card select NI ELVIS II, and PC machine selects Intel (R) Core (TM) i7-4500U.
High Pass Filter Cutoff Frequency is set as 800Hz, and the cutoff frequency of low-pass filter is set as 1kHZ, the driving of random sequence
Clock setting is 10kHz, and low speed sample frequency is set as 2kHz, and acquisition time length is 0.1s, sampling number 200, compression
Sampling is more effective than being that the channel 20%, AO trigger source is set as rising edge, and it is effective that the channel AI trigger source is set as failing edge, triggering
Clock is the crt0 that NI ELVIS II is carried, and trigger output end is NI ELVIS II plate upper port PFI12.
According to specific implementation step 2), signal feeding high-pass filter is filtered, filtered waveform passes through NI
ELVIS II included oscillograph function is shown as shown in Figure 5.
According to specific implementation step 3), synchronous triggering mode is set.
According to specific implementation step 4), the signal after high-pass filtering is sent into next stage mixing filter, trigger source rises
Along the pulse it is flushed to when coming, the mixed waveform of output and waveform after low-pass filtering are as shown in Figure 6,7.
According to specific implementation step 5), the signal after low-pass filtering is sent into the channel AI of next stage NI ELVIS II, touching
When the failing edge pulse that rises arrives, start low speed sampling, the LabVIEW that the waveform and its amplitude spectrum of input are shown in PC machine is soft
On the front panel of part as shown in Figure 8.
According to specific implementation step 6), using OMP algorithm, complete signal reconstruction.
According to specific implementation step 7), show that the waveform reconstructed and its amplitude spectrum are as shown in Figure 10.
Compression ratio is adjusted, it is 50 that compression ratio, which is reduced to the sampling number in 5%, 0.1s, then low speed sample frequency is
500Hz, the anti-low pass filter cutoff frequency postponed are not more than 250Hz.The cutoff frequency for readjusting low-pass filter arrives
200Hz repeats specific implementation step of the present invention 2)~7), the waveform and its amplitude spectrum after low speed sampling are as shown in figure 9, weight
The waveform and its amplitude spectrum that structure goes out are as shown in Figure 10.
The above results show that the method for the present invention and implementation are in the premise for ensuring that ultra harmonics signal accurately detects
Under, greatly reduce the sample frequency and sampling number of ultra harmonics.
The superiority of method and scheme in ultra harmonics detection application in order to further illustrate the present invention, is adopted simultaneously
10 power frequency periods (0.2s) are carried out to identical harmonic wave with conventional harmonic measuring instrument with the method for the present invention and scheme to adopt
Sample, the results are shown in Table 1.
Table 1
In conclusion compressed sensing based ultra harmonics detection method and implementation of the invention meet it is humorous
Under the premise of wave accurately detects, the ADC sample rate of measuring instrument, section of sampling number and low-pass filter can be substantially reduced
Only frequency.This method can also apply in the monitoring and analysis of power distribution network transient power quality after slightly adjusting simultaneously.
The content that description in the present invention is not described in detail belongs to existing skill well known to professional and technical personnel in the field
Art.
Finally bright: above-described embodiments are merely to illustrate the technical scheme, rather than to it
Limitation；Although the present invention is described in detail referring to the foregoing embodiments, those skilled in the art can manage
Solution: without departing from the spirit and scope of the invention and the appended claims, various substitutions, changes and modifications are all that have can
Can.Therefore, protection scope of the present invention should be subject to the protection scope in claims.
Claims (4)
1. a kind of compressed sensing based ultra harmonics detection device, including sensing circuit (1), the high-pass filtering of adjustable gain
Device (2), mixing filter (3), multifunctional data acquisition card (4), host computer (5), which is characterized in that sensing circuit (1) passes through
Voltage and current signal connects the high-pass filter (2) of adjustable gain, and the high-pass filter (2) of adjustable gain connects mixing filter
(3), the frequency mixer of mixing filter (3), the A/D module of low-pass filter and multifunctional data acquisition card (4), D/A module according to
Secondary signal connection, the synchronization trigger module of multifunctional data acquisition card (4) is connect with A/D module, D/A module by signal, multi-functional
Data collecting card (4) passes through USB connection host computer (5)；
The sensing circuit (1) is arranged in power distribution network measurement end, including voltage sensor and current sensor；
The high-pass filter (2) of the adjustable gain is used for signal condition, on the one hand can filter out power frequency in signal and other
On the other hand low-frequency component plays the role of signal amplitude conditioning, be the simulation filter of gain, cutoff frequency, regulable center frequency
Wave device；
The mixing filter (3) is made of frequency mixer and low-pass filter, wherein section of low pass filtered device
Only frequency, regulable center frequency can play the role of Ratios detection；
The multifunctional data acquisition card (4) generates for random sequence, low speed sampling, synchronizes triggering and USB data transmission；
The A/D module is sampled for low speed, D/A module for generating pseudo-random sequence, synchronous trigger module for realizing
D/A, A/D two-way press rising edge and failing edge triggering respectively, and USB module is used for data-transformation facility.
The host computer (5) is for realizing the reconstruct and display to ultra harmonics time frequency signal, it can be achieved that pseudo-random sequence
The multiple tasks such as generation, data acquisition, signal reconstruction and display.
2. a kind of compressed sensing based ultra harmonics detection method, which comprises the following steps:
Step 1) sensing circuit (1) is mounted on power distribution network measurement end, includes voltage transformer and current transformer, can be by power distribution network
In big voltage, current signal is converted to small voltage, current signal and plays the role of being isolated with power grid；
The signal x (t) of step 2) sensing circuit (1) output, which is sent to the high-pass filter of adjustable gain (2), obtains filtered letter
Number x1(t), the effect of high-pass filter can play power frequency and the other low-frequency components filtered out in signal, only retain ultra harmonics
The signal of bandwidth ingredient；
Step 3) uses the synchronization trigger source on multifunctional data acquiring board (4), pseudorandom sequence in guarantee step 4), step 5)
Column output channel AO is synchronous with low speed sampling input channel AI, it is ensured that is not in phase deviation situation；
Step 4) is by signal x1(t) it sending into mixing filter (3), mixing filter (3) includes multiplier and low-pass filter,
Wherein pseudo-random sequence pc(t) it is exported by the AO mouth of multifunctional data acquisition card (4), x1(t) and pc(t) the first phase in multiplier
It is multiplied to arrive signal y (t), according to the Convolution Properties of continuous fourier transform:
FT [y (t)]=FT [x1(t)pc(t)]=X1(j Ω) * P (j Ω)=Y (j Ω)
It is found that the frequency spectrum of y (t) is x1(t)、pc(t) convolution of frequency spectrum, since the spectrum distribution of pseudo-random sequence is in entire frequency domain
It is interior, therefore can will be at the frequency spectrum shift tremendously low frequency of ultra harmonics by convolutional calculation；
Then, y (t) obtains signal y by low-pass filter again1(t), it may be expressed as: by the process of low-pass filter
H (t) is the unit impulse response of low-pass filter；
After passband maximum attenuation is given, the cut-off frequecy of passband f of the low-pass filter0It can be according to required Sampling Compression than carrying out
It manually adjusts, realizes Ratios filtering；
Step 5) is by signal y1(t) it send to the channel AI of multifunctional data acquisition card 4 and carries out low speed A/D sampling, sampling process is such as
Under:
Sampling time interval T is by sample frequency fsIt determines.Enable fs≥2f0, as the cut-off frequecy of passband f of low-pass filter0When lower,
Sample frequency fsLower, sampling time interval T is bigger, and sampling number N is fewer, it can realizes the variable compression to ultra harmonics
Compare compression sampling；
Step 6) is by y1(n) it is sent into host computer, compression sampling signal is reconstructed using orthogonal matching pursuit OMP algorithm, first
Recover in the frequency spectrum θ, θ of signal for 0 coefficient be exactly ultra harmonics spectrum component, FFT inverse transformation secondly is passed through to θ
After can obtain reconstruction signal x '1(n)；
Ultra harmonics time-frequency waveform after step 7) display reconstruct.
3. a kind of compressed sensing based ultra harmonics detection method according to claim 2, which is characterized in that described
The signal of step 6) output be reconstructed that specific step is as follows:
Step A: according to analog information converter AIC actual circuit, random sequence matrix P, the unit of low-pass filtering link are constructed
Impulse response matrix H, the sampling matrix B of low speed sampling element；
Step B: observing matrix Φ=BHP is obtained；
Step C: construction discrete Fourier transform basic matrix Ψ；
Step D: sensing matrix A=Φ Ψ is obtained；
Step E: the frequency spectrum θ of ultra harmonics is reconstructed using OMP algorithm, and superelevation is recovered by inverse discrete Fourier transform
Rd harmonic signal.
4. a kind of compressed sensing based ultra harmonics detection method according to claim 3, which is characterized in that described
Step 3) in physically realizable observing matrix be analog information converter AIC, and the low-pass filter in AIC be gain,
The analog filter of cutoff frequency, regulable center frequency.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810826416.7A CN109164298A (en) | 2018-07-25 | 2018-07-25 | A kind of compressed sensing based ultra harmonics detection device and detection method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810826416.7A CN109164298A (en) | 2018-07-25 | 2018-07-25 | A kind of compressed sensing based ultra harmonics detection device and detection method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109164298A true CN109164298A (en) | 2019-01-08 |
Family
ID=64898074
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810826416.7A Pending CN109164298A (en) | 2018-07-25 | 2018-07-25 | A kind of compressed sensing based ultra harmonics detection device and detection method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109164298A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110045184A (en) * | 2019-04-02 | 2019-07-23 | 国网上海市电力公司 | A kind of sub-harmonic wave measurement method based on compressed sensing MACSMP |
| CN110082592A (en) * | 2019-05-29 | 2019-08-02 | 国网上海市电力公司 | A kind of ultra harmonics measurement method based on flexible atom filtering |
| CN110716088A (en) * | 2019-10-30 | 2020-01-21 | 国网上海市电力公司 | Ultrahigh harmonic measurement method based on compressive sensing MACSMP |
| CN111355493A (en) * | 2020-04-03 | 2020-06-30 | 哈尔滨工业大学 | Support set screening and reconstructing method for modulation broadband converter |
| CN112450941A (en) * | 2020-11-11 | 2021-03-09 | 南昌大学 | Electrocardiosignal compression sampling device and method based on random demodulation structure |
| CN115902397A (en) * | 2023-01-06 | 2023-04-04 | 四川大学 | Rapid super-harmonic estimation method based on simplified compressed sensing model |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103983850A (en) * | 2014-05-13 | 2014-08-13 | 天津大学 | Power system harmonious wave compressed signal reconstruction and detection method based on compressed sensing |
| CN107192878A (en) * | 2017-04-07 | 2017-09-22 | 中国农业大学 | A kind of trend of harmonic detection method of power and device based on compressed sensing |
-
2018
- 2018-07-25 CN CN201810826416.7A patent/CN109164298A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103983850A (en) * | 2014-05-13 | 2014-08-13 | 天津大学 | Power system harmonious wave compressed signal reconstruction and detection method based on compressed sensing |
| CN107192878A (en) * | 2017-04-07 | 2017-09-22 | 中国农业大学 | A kind of trend of harmonic detection method of power and device based on compressed sensing |
Non-Patent Citations (5)
| Title |
|---|
| XIFENG GUO等: "Harmonic Detection Method Based on Compressed Sensing Theory in the Application of APF", 《2017 29TH CHINESE CONTROL AND DECISION CONFERENCE》 * |
| 刘嫣等: "基于CS的电能质量扰动数据稀疏分析与改进重构算法", 《电工技术学报》 * |
| 庄双勇等: "基于压缩感知OMP的超谐波测量新算法", 《仪器仪表学报》 * |
| 朱国锋等: "基于压缩感知的配电网谐波记录分析系统", 《计算机应用》 * |
| 陈国强: "电能质量压缩感知采样及重构算法研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110045184A (en) * | 2019-04-02 | 2019-07-23 | 国网上海市电力公司 | A kind of sub-harmonic wave measurement method based on compressed sensing MACSMP |
| CN110082592A (en) * | 2019-05-29 | 2019-08-02 | 国网上海市电力公司 | A kind of ultra harmonics measurement method based on flexible atom filtering |
| CN110716088A (en) * | 2019-10-30 | 2020-01-21 | 国网上海市电力公司 | Ultrahigh harmonic measurement method based on compressive sensing MACSMP |
| CN111355493A (en) * | 2020-04-03 | 2020-06-30 | 哈尔滨工业大学 | Support set screening and reconstructing method for modulation broadband converter |
| CN111355493B (en) * | 2020-04-03 | 2023-05-23 | 哈尔滨工业大学 | Support set screening and reconstructing method for modulation broadband converter |
| CN112450941A (en) * | 2020-11-11 | 2021-03-09 | 南昌大学 | Electrocardiosignal compression sampling device and method based on random demodulation structure |
| CN115902397A (en) * | 2023-01-06 | 2023-04-04 | 四川大学 | Rapid super-harmonic estimation method based on simplified compressed sensing model |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109164298A (en) | A kind of compressed sensing based ultra harmonics detection device and detection method | |
| CN104749498B (en) | A kind of detection of Portable PD On-Line and diagnostic device | |
| CN103257271A (en) | Device and method for detecting micro grid harmonic wave and inter-harmonics based on STM32F107VCT6 | |
| CN203287435U (en) | A micro electrical network harmonic wave and inter-harmonic wave test apparatus based on an STM32F107VCT6 | |
| CN103454497A (en) | Phase difference measuring method based on improved windowing discrete Fourier transform | |
| CN109946512B (en) | A kind of dynamic power analysis method for improving frequency domain interpolation | |
| CN203133168U (en) | Power harmonic detector | |
| CN102520245A (en) | Micro-grid harmonic and inter-harmonic analysis method based on cubic spline interpolation waveform reconstruction | |
| CN110020452B (en) | Electromagnetic compatibility test simulation analysis method and system thereof | |
| CN108318852A (en) | A kind of intelligent electric energy meter square wave influence test method | |
| CN106680585B (en) | Harmonic wave/m-Acetyl chlorophosphonazo detection method | |
| CN103969508B (en) | A kind of Electric Power Harmonic Analysis method high-accuracy in real time and device | |
| CN102749505A (en) | Method and device for voltage proportion measurement | |
| CN107093357A (en) | Multifunction teaching integrated signal analysis system | |
| CN105445541A (en) | Method for adaptively calculating power under arbitrary frequencies | |
| CN108152584A (en) | A kind of high ferro tractive power supply system harmonic wave Multi-path synchronous rapid detection method | |
| CN101806835B (en) | Interharmonics measuring meter based on envelope decomposition | |
| CN105866576A (en) | Simulation detection system of intelligent transformer station secondary-side electric energy metering error influence, and detection analysis method | |
| CN104991119A (en) | Co-prime spectrum analysis method and apparatus for eliminating pseudo peak and spectrum leakage effects | |
| CN203069643U (en) | Multiple variable bandwidth digital storage oscilloscope | |
| CN102006129A (en) | FPGA (Field Programmable Gate Array)-based method and device for simulating power line channel | |
| CN106872820B (en) | Method for measuring eye pattern based on vector network analyzer | |
| CN107167757A (en) | A kind of method for checking electronic transducer and system using improvement digital filtering algorithm | |
| CN104483546B (en) | Spectrum analysis method of FPGA (Field Programmable Gate Array) digital logic signal | |
| CN113009201A (en) | Electric signal sampling 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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190108 |
|
| RJ01 | Rejection of invention patent application after publication |