CN106872779A - A kind of harmonic source determination methods based on harmonic phase jittering characteristic - Google Patents
A kind of harmonic source determination methods based on harmonic phase jittering characteristic Download PDFInfo
- Publication number
- CN106872779A CN106872779A CN201710121164.3A CN201710121164A CN106872779A CN 106872779 A CN106872779 A CN 106872779A CN 201710121164 A CN201710121164 A CN 201710121164A CN 106872779 A CN106872779 A CN 106872779A
- Authority
- CN
- China
- Prior art keywords
- phase
- signal
- frequency
- harmonic
- devices
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- 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
Abstract
The present invention relates to a kind of harmonic source determination methods based on harmonic phase jittering characteristic, the output port to the Devices to test of single port output carries out collecting Devices to test output signal for the first time；Devices to test output signal to getting carries out Fourier transformation；Second sampling and data processing：Second collection is carried out with data collecting card device under test output port, Devices to test output signal is obtained, the Devices to test output signal to second collection carries out Fourier transformation, the phase of signal is calculated according to phase definition；The phase that frequency multiplication is referred to as to harmonic component carries out solution winding, when the ratio between phase difference value of frequency-doubled signal is equal with the ratio between overtone order in two phase frequency figures, and frequency-doubled signal is when the ratio between the phase change rate in sampling interval is equal with the ratio between overtone order, think that these have consistent jittering characteristic into the frequency component of multiple proportion, so as to judge that these frequency-doubled signals are derived from same harmonic source, otherwise it is judged as from different independent sources.
Description
Technical field
The invention belongs to electromagnetic compatibility test field, it is related to a kind of harmonic source judgement side based on harmonic phase jittering characteristic
Method.
Background technology
In the Conduction Interference transmitting test of electromagnetic compatibility, there are many independent interference sources inside Device under test, these are done
The interference signal for disturbing source generation is coupled to port to be measured by coupling channel, together form the interference signal of port to be measured.When
In the port output of Device under test at a certain frequency signal transmitting it is exceeded, it is necessary to the corresponding emission source of the frequency is carried out positioning with
Rectification.But there is a kind of situation, the internal interference source of Device under test has frequency multiplication relation in the signal that port to be measured is launched
Several frequencies when, it is necessary to judge whether these are derived from the harmonic wave of certain signal source into the signal of frequency multiplication relation, or
Some independent signal sources, and then take corresponding braking measure to radiating exceeded frequency.
When solving the problems, such as to judge harmonic source, existing solution is studied respectively by the mutual information in information theory
Mutual information between harmonic, this method is complicated, computationally intensive, and needs the engineering staff to have professional knowledge higher, number
Gain knowledge, implement complexity.
The content of the invention
The technology of the present invention solve problem：The method for overcoming prior art is complicated, computationally intensive, there is provided one kind is based on harmonic wave phase
The harmonic source determination methods of position jittering characteristic, principle are it can be readily appreciated that operation is simple with data processing.
The technology of the present invention solution：A kind of harmonic source determination methods based on harmonic phase jittering characteristic, for equipment
Internal signal sources are stationary signal source, when frequency-doubled signal occurs in single port device, using the phase jitter of these signals
Characteristic judges whether these frequency-doubled signals are derived from the harmonic wave of a certain signal, or from different independent sources.
Realize that step is as follows：
Step one, signal acquisition, the output port of the Devices to test exported to single port with data collecting card carry out first
Secondary collection, obtains Devices to test output signal；
Step 2, data processing, the Devices to test output signal to being got in step one carry out Fourier transformation, obtain
The signal data of plural form, according to the definition of phaseIm is the imaginary part of Fourier transformation, and Re becomes for Fourier
The real part for changing,It is signal phase, calculates the frequency-phase relation of signal, draws frequency plot graph of a relation, i.e. phase frequency figure；
Step 3, second sampling and data processing, repeat step one and step 2, use data collecting card device under test
Output port carries out second collection, obtains Devices to test output signal.Devices to test output signal to second collection is entered
Row Fourier transformation, the phase of signal is calculated according to phase definition, draws phase frequency figure.The front and rear sampling interval twice will be as far as possible
Less than a quarter signal period, with before and after ensureing twice the phase difference of signal respective frequencies in 2 π, to make double sampling between
Every the requirement less than a quarter signal period is met, double sampling first can be carried out to equipment output end mouth, then to collection
Data do Fourier transformation；
Step 4, harmonic phase jittering characteristic is judged, into multiple proportion in the phase frequency figure drawn to step 2 and step 3
The phase of frequency component (i.e. harmonic wave, abbreviation frequency multiplication) carry out solution winding, when the phase of corresponding frequency-doubled signal in two phase frequency figures
The ratio between potential difference value is equal with the ratio between overtone order, and the ratio between the phase change rate of frequency-doubled signal it is equal with the ratio between overtone order when,
Think that these have consistent jittering characteristic into the frequency component of multiple proportion, it may be determined that these frequency-doubled signals are from same
One harmonic source, be otherwise judged as from different independent source (in the present invention, harmonic source is defined as producing the signal source of harmonic wave,
Independent source is defined as producing the signal source of simple signal).
Above-mentioned operation principle is illustrated below：
When the circuit system of Devices to test is nonlinear system, come from the signal amplitude of each harmonic wave of a certain signal source
May be affected, it is impossible to show some particular kind of relationships, it is difficult to determine whether that these frequencies are by the amplitude of harmonic wave
It is no to come from same signal source, but when harmonic wave comes from same signal source, the phase of each harmonic can show identical
Jittering characteristic, the judgement of harmonic source can be carried out with this.The signal source that this hair is related to is periodic signal.
Because the transmission equation of almost all of electronic component or circuit system is all nonlinear equation, but smaller
During range of motion, these relations are seen linear, very big error will not be produced, at this moment can be with Taylor series come table
Show operating point equation.
Here diode is taken as elaboration object.According in microwave engineering theory, typical RF diodes V-I characteristics are：
I (V)=Is(eαV-1) (1)
Wherein α=q/nkT, q are electron charges, and k is Boltzmann constant, and T is temperature, and n is the idealization factor, IsIt is full
And electric current, e is natural constant, output current I when I (V) represents that diode input voltage is V.When diode voltage is：
V=V0+v (2)
V0It is direct current biasing, v is small AC signal voltage.(2) formula is to V0Do Taylor series expansion and take three ranks, write as
DC bias current I0It is with alternating current i sum forms：
Wherein GdRepresent I to V in V0First derivative.Because preferable continuous wave signal is：
V (t)=A0sin(w0t) (4)
And true continuous wave signal is
A0It is signal amplitude, w0It is signal frequency, E (t) shakes for random magnitude,It is random phase fluctuation；
When diode voltage is (5) formula,
Above formula can be described as：
Wherein, D0、A1、A2、A3The respectively amplitude of direct current signal, fundamental wave, second harmonic, triple-frequency harmonics.
When diode output signal I (V) is exported by nonlinear circuit, can there is additional amplitude decay and additional phase
Move, it is assumed that I (V) is output as through nonlinear circuit：
Wherein, a0、a1、a2、a3The respectively additional amplitude decay of direct current signal, fundamental wave, second harmonic, triple-frequency harmonics,It is additional phase shift.
Assuming that the sampling time of data collecting card is very short, the random phase fluctuation of continuous signal within the sampling timeFor
Constant.
Sampled when data acquisition is stuck in the t=0 moment, the signal of collection can be expressed as：
Sampled when data acquisition is stuck in t=Δ ts, the signal of collection can be expressed as：
So frequency w0、2w0、3w0Phase change rate be respectively：
It can be seen that, sinusoidal signal can produce harmonic wave when by RF diodes and nonlinear circuit, and harmonic wave phase
Position and phase jitter characteristic are with fundamental frequency into multiple, and multiple is relevant with overtone order.
Different independent sources is derived from for non-harmonic signals, i.e. frequency-doubled signal.Because the phase jitter of each independent source is
Random, so the phase at each frequency is uncorrelated, now circuit output is：
Wherein, E0、B1、B2、B3The respectively amplitude of direct current signal, fundamental wave, second harmonic, triple-frequency harmonics,Respectively frequency w0、2w0、3w0Random phase fluctuation,It is additional phase shift.
Sampled when data acquisition is stuck in t=0, the signal of collection can be expressed as：
Sampled when data acquisition is stuck in t=Δs t, the signal of collection can be expressed as：
So frequency w0、2w0、3w0Phase change rate be respectively：
Due toIt is each independent, it can be seen that when these frequency signals are independent signal source,
Frequency w0、2w0、3w0Phase without contact, and without consistent phase jitter characteristic.
Present invention advantage compared with prior art is：What is run into blind source separating is solved the problems, such as judges harmonic source
During problem, existing solution studies the mutual information between each harmonic, this side by the mutual information in information theory
Method is complicated, computationally intensive, and needs the engineering staff to have professional knowledge higher, a math quality, and the present invention proposes research times
The method of the phase jitter between frequency component is then simple to operate, and data processing is simple, and principle should be readily appreciated that.Believed by frequency multiplication in principle
Phase jitter characteristic between number falls within the category of mutual information, but this patent employs the mode of research signal phase, thinking
Novelty, it is simple to operate.
Brief description of the drawings
Fig. 1 is the flow chart of harmonic source determination methods of the present invention based on harmonic phase jittering characteristic.
Specific embodiment
Below in conjunction with drawings and Examples, the present invention is described in further detail.
As shown in figure 1, the present invention proposes a kind of harmonic source determination methods based on harmonic phase jittering characteristic, including it is following
Process step：
Step one, signal acquisition, the output port of the Devices to test exported to single port with data collecting card in t=0
First time collection is carried out, Devices to test output signal is obtained, the output signal of collection is represented with formula (16) for the first time,
Step 2, data processing, the Devices to test output signal to being got in step one carry out Fourier transformation, obtain
The signal data of plural form, according to the definition of phaseIm is the imaginary part of Fourier transformation, and Re becomes for Fourier
The real part for changing,It is signal phase, the frequency-phase relation of signal is calculated, due to the corresponding phase of each frequency component in phase frequency figure
Position is initial phase of the frequency to induction signal, so w0、2w0、3w0Corresponding phase is respectively in phase frequency figure
Step 3, second sampling and data processing, repeat step one and step 2, the data collecting card in t=Δ t
Device under test output port carries out second collection, obtains Devices to test output signal, and second output signal of collection is used
Formula (17) expression,
Devices to test output signal to second collection carries out Fourier transformation, and the phase of signal is calculated according to phase definition
Position, obtains the w in phase frequency figure0、2w0、3w0Corresponding phase is respectively
Step 4, harmonic phase jittering characteristic is judged, into multiple proportion in the phase frequency figure drawn to step 2 and step 3
The phase of frequency component (i.e. harmonic wave) carry out solution winding.
The data result sampled for the first time：In phase frequency figure, w0、2w0、3w0Corresponding phase is respectively
Second data result of sampling：In phase frequency figure, w0、2w0、3w0Corresponding phase is respectively
The data of double sampling are in w0、2w0、3w0Difference at frequency is respectively
The harmonic source is in w0、2w0、3w0Phase change rate at frequency is：
According to the above results, frequency-doubled signal is in w in two phase frequency figures0、2w0、3w0The ratio between phase difference value at frequency is 1:
2:3, the ratio between phase change rate is 1:2:3, the ratio between overtone order is 1:2:3, it can be seen that frequency-doubled signal is in w0、2w0、3w0Frequency
The ratio between the phase difference value at place is equal to the ratio between phase change rate, equal to the ratio between overtone order, it is believed that these are into multiple proportion
Frequency component has consistent jittering characteristic, judges that these frequency-doubled signals are derived from same harmonic source.
Above example is provided just for the sake of the description purpose of the present invention, and is not intended to limit the scope of the present invention.This
The scope of invention is defined by the following claims.The various equivalents that do not depart from spirit and principles of the present invention and make and repair
Change, all should cover within the scope of the present invention.
Claims (2)
1. a kind of harmonic source determination methods based on harmonic phase jittering characteristic, it is characterised in that comprise the following steps：
Step one, signal acquisition：The output port of the Devices to test exported to single port with data collecting card adopt for the first time
Collection, obtains Devices to test output signal；
Step 2, data processing：Devices to test output signal to being got in step one carries out Fourier transformation, obtains plural number
The signal data of form, according to the definition of phaseThe frequency-phase relation of signal is calculated, wherein Im is Fourier
The imaginary part of conversion, Re is the real part of Fourier transformation,It is signal phase；Then frequency plot graph of a relation, i.e. phase frequency figure are drawn；
Step 3, second sampling and data processing：Repeat step one and step 2, are exported with data collecting card device under test
Port carries out second collection, obtains Devices to test output signal, and the Devices to test output signal to second collection carries out Fu
In leaf transformation, according to phase definition calculate signal phase, draw phase frequency figure；The front and rear sampling interval twice is less than four/
One signal period, to ensure the phase difference of front and rear signal respective frequencies twice in 2 π；
Step 4, harmonic phase jittering characteristic is judged, into the frequency of multiple proportion in the phase frequency figure drawn to step 2 and step 3
The phase that rate component, i.e. harmonic component are referred to as frequency multiplication carries out solution winding, when the phase difference value of frequency-doubled signal in two phase frequency figures
The ratio between it is equal with the ratio between overtone order, and frequency-doubled signal the sampling interval the ratio between phase change rate it is equal with the ratio between overtone order
When, it is believed that these have consistent jittering characteristic into the frequency component of multiple proportion, so as to judge that these frequency-doubled signals come
Same harmonic source is come from, is otherwise judged as from different independent sources, the harmonic source is defined as producing the signal source of harmonic wave,
The independent source is defined as producing the signal source of simple signal.
2. harmonic source determination methods based on harmonic phase jittering characteristic according to claim 1, it is characterised in that：It is described
In step 2, to make double sampling be spaced the requirement met less than a quarter signal period, must first device under test output end
Mouth carries out double sampling, then the data to collection do Fourier transformation respectively.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710121164.3A CN106872779B (en) | 2017-03-02 | 2017-03-02 | A kind of harmonic source judgment method based on harmonic phase jittering characteristic |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710121164.3A CN106872779B (en) | 2017-03-02 | 2017-03-02 | A kind of harmonic source judgment method based on harmonic phase jittering characteristic |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106872779A true CN106872779A (en) | 2017-06-20 |
| CN106872779B CN106872779B (en) | 2019-03-29 |
Family
ID=59169251
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710121164.3A Active CN106872779B (en) | 2017-03-02 | 2017-03-02 | A kind of harmonic source judgment method based on harmonic phase jittering characteristic |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106872779B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107192877A (en) * | 2017-05-19 | 2017-09-22 | 上海航天控制技术研究所 | A kind of phase interpretation method based on Fast Fourier Transform (FFT) |
| CN110927425A (en) * | 2018-09-20 | 2020-03-27 | 中国石油化工股份有限公司 | Harmonic source section positioning method of power distribution network based on harmonic current monitoring |
| CN110927426A (en) * | 2018-09-20 | 2020-03-27 | 中国石油化工股份有限公司 | Harmonic source section positioning method for power distribution network |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4412299A (en) * | 1981-02-02 | 1983-10-25 | Teltone Corporation | Phase jitter detector |
| US20030001557A1 (en) * | 2001-06-28 | 2003-01-02 | Conexant Systems, Inc. | System and method for detecting phase offset in a phase-locked loop |
| CN102128982A (en) * | 2010-12-23 | 2011-07-20 | 东莞市创锐电子技术有限公司 | Harmonic analysis method based on windowed interpolation FFT (Fast Fourier Transform) base frequency tracking technology |
| CN201974478U (en) * | 2010-12-23 | 2011-09-14 | 大连艾科科技开发有限公司 | Detection apparatus of harmonic signal of photoelectric conversion capable of reducing noise and raising sensitivity |
-
2017
- 2017-03-02 CN CN201710121164.3A patent/CN106872779B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4412299A (en) * | 1981-02-02 | 1983-10-25 | Teltone Corporation | Phase jitter detector |
| US20030001557A1 (en) * | 2001-06-28 | 2003-01-02 | Conexant Systems, Inc. | System and method for detecting phase offset in a phase-locked loop |
| CN102128982A (en) * | 2010-12-23 | 2011-07-20 | 东莞市创锐电子技术有限公司 | Harmonic analysis method based on windowed interpolation FFT (Fast Fourier Transform) base frequency tracking technology |
| CN201974478U (en) * | 2010-12-23 | 2011-09-14 | 大连艾科科技开发有限公司 | Detection apparatus of harmonic signal of photoelectric conversion capable of reducing noise and raising sensitivity |
Non-Patent Citations (2)
| Title |
|---|
| 艾永乐等: ""基于谐波有功功率贡献量的主谐波源定位"", 《电力系统保护与控制》 * |
| 陈厦门等: ""复杂电网中谐波源搜索定位的仿真分析"", 《电力与电工》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107192877A (en) * | 2017-05-19 | 2017-09-22 | 上海航天控制技术研究所 | A kind of phase interpretation method based on Fast Fourier Transform (FFT) |
| CN107192877B (en) * | 2017-05-19 | 2019-08-02 | 上海航天控制技术研究所 | A kind of phase interpretation method based on Fast Fourier Transform (FFT) |
| CN110927425A (en) * | 2018-09-20 | 2020-03-27 | 中国石油化工股份有限公司 | Harmonic source section positioning method of power distribution network based on harmonic current monitoring |
| CN110927426A (en) * | 2018-09-20 | 2020-03-27 | 中国石油化工股份有限公司 | Harmonic source section positioning method for power distribution network |
| CN110927426B (en) * | 2018-09-20 | 2021-09-07 | 中国石油化工股份有限公司 | Harmonic source section positioning method for power distribution network |
| CN110927425B (en) * | 2018-09-20 | 2021-09-07 | 中国石油化工股份有限公司 | Harmonic source section positioning method of power distribution network based on harmonic current monitoring |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106872779B (en) | 2019-03-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Senroy et al. | An improved Hilbert–Huang method for analysis of time-varying waveforms in power quality | |
| Wang et al. | A machine learning approach to detection of geomagnetically induced currents in power grids | |
| CN106872779B (en) | A kind of harmonic source judgment method based on harmonic phase jittering characteristic | |
| CN108414833B (en) | Accurate estimation method for signal component frequency | |
| Flowers-Jacobs et al. | A Boltzmann constant determination based on Johnson noise thermometry | |
| Ravindran et al. | Interharmonics in PV systems: a review of analysis and estimation methods; considerations for selection of an apt method | |
| Kuznetsov et al. | Characterization of the cyclostationary emissions in the near-field of electronic device | |
| Kuznetsov et al. | Cyclostationary characterization of radiated emissions in digital electronic devices | |
| Remley et al. | Measurement bandwidth extension using multisine signals: Propagation of error | |
| Mariscotti et al. | Variability caused by setup and operating conditions for conducted EMI of switched mode power supplies over the 2–1000 kHz interval | |
| Musznicki et al. | Adaptive estimation of the transformer stray capacitances for DC–DC converter modelling | |
| CN101806835B (en) | Interharmonics measuring meter based on envelope decomposition | |
| Wu et al. | Optimal pmu design based on sampling model and sensitivity analysis | |
| Zamani et al. | A novel improved Hilbert-Huang transform technique for implementation of power system local oscillation monitoring | |
| CN109581518A (en) | One kind is without synchronizing current acquisition transmission and the acquisition of frequency domain measurement SIP data and processing method | |
| Gunasekara | Contactless estimation of breathing rate using UWB radar | |
| Malkhandi et al. | A Fourier‐based single phase PLL algorithm: Design, analysis, and implementation in FPGA controller | |
| Lin | Development of interharmonics identification using enhanced‐fft algorithm | |
| CN107273620B (en) | Electromagnetic emission modeling method based on joint estimation | |
| Pramanik et al. | A two-stage adaptive symmetric-strong-tracking square-root cubature Kalman filter for harmonics and interharmonics estimation | |
| Kuznetsov et al. | Cyclostationary Source Separation Based on Electrornagnetic Measurements in the Near-field of PCB | |
| Yi et al. | Novel methods for in-situ direct magnetic loss measurement in a DC-DC converter | |
| Srivastava et al. | Harmonic/Inter-harmonic Estimation: Key Issues and Challenges | |
| Imtiaz et al. | Frequency domain feature extraction for power quality disturbance classification | |
| Chen et al. | Development of a high function portable PD measurement system for periodic inspection of power cable circuits |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |