CN108008188A - A kind of inter-channel phase difference, the measuring method of amplitude difference - Google Patents
A kind of inter-channel phase difference, the measuring method of amplitude difference Download PDFInfo
- Publication number
- CN108008188A CN108008188A CN201711273036.7A CN201711273036A CN108008188A CN 108008188 A CN108008188 A CN 108008188A CN 201711273036 A CN201711273036 A CN 201711273036A CN 108008188 A CN108008188 A CN 108008188A
- Authority
- CN
- China
- Prior art keywords
- signal
- channel
- mrow
- phase difference
- passage
- 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 Current Or Voltage (AREA)
Abstract
The present invention is more particularly directed to a kind of inter-channel phase difference, the measuring method of amplitude difference, include the following steps:(A) multichannel input signal is denoted as r (k, n), selects any of which passage to be used as and refers to passage;(B) square operation is carried out to each channel signal and mean filter handles to obtain signal P (k):(C) signal Q (k) is obtained after reference channel being carried out multiplying and mean filter processing with other channel signals respectively;(D) the amplitude difference A between other passages and reference channel is calculated by signal P (k) and Q (k)Δ(k) and phase differenceBy being operated to the common digital function such as list entries square, multiplication, average, evolution, anticosine, the measurement result of amplitude difference and phase difference is provided in the form of formula, is dealt with very simply, conveniently and quickly;Program precision is adjustable at the same time, and the scope of application is more extensive, and need not carry out Digital Down Convert processing, more succinctly, quickly.
Description
Technical field
The present invention relates to communication reconnaissance and digital signal processing technique field, more particularly to a kind of inter-channel phase difference, width
Spend the measuring method of difference.
Background technology
The phase difference of traditional interchannel, amplitude difference measuring method are carried out using Fast Fourier Transform (FFT) (FFT)
Frequency domain measurement, there are many deficiencies for this measurement method:First, since spectral leakage caused by the FFT of non-integer can drop
Low amplitude and phase measurement accuracy, so that algorithm can introduce certain error in itself, cause measurement result not accurate enough;Second,
All frequency points can be calculated by being calculated by FFT, and most of frequency point is unwanted in practical application, therefore can increase computing
Measure, and operation time also greatly increases;Third, when being measured using Fast Fourier Transform (FFT), it is also necessary to which signal is carried out
Digital Down Convert is handled, and when computing needs to know the frequency of signal.
The content of the invention
It is an object of the invention to provide a kind of inter-channel phase difference, the measuring method of amplitude difference, arithmetic speed is fast and smart
Degree is high.
In order to achieve the above object, the technical solution adopted by the present invention is:A kind of inter-channel phase difference, the measurement side of amplitude difference
Method, includes the following steps:
(A) multichannel input signal is denoted as r (k, n), wherein channel number k=0,1,2 ..., K-1, sequential sampling moment n
=0,1,2 ..., N-1, selecting any of which passage, reference channel is denoted as r (m, n), m ∈ [0, K-1] as passage is referred to;
(B) square operation is carried out to each channel signal and mean filter handles to obtain signal P (k):
Wherein k=0,1,2 ..., K-1;
(C) by reference channel with obtaining signal Q after being handled respectively with other channel signal multiplyings and mean filter
(k):
Wherein k=0,1,2 ..., K-1 and k ≠ m, m ∈ [0, K-1];
(D) the amplitude difference A between other passages and reference channel is calculated by signal P (k) and Q (k)Δ(k) and phase
Difference
K=0 in two formula, 1,2 ..., K-1 and k ≠ m, m ∈ [0, K-1];
Mean filter processing in the step B and C is used to filter out the high fdrequency component in signal.
Compared with prior art, there are following technique effect by the present invention:By to list entries square, multiplication, it is average, open
The common digital function such as side, anticosine is operated, and the measurement result of amplitude difference and phase difference is provided in the form of formula,
Deal with very simply, conveniently and quickly;Can be by increasing the value of N, thus can be real to increase the precision of result of calculation
Existing precision is adjustable, and the scope of application is more extensive;Meanwhile this method is handled directly against multichannel input signal, without
Digital Down Convert processing is carried out, more succinctly, quickly.
Brief description of the drawings
Fig. 1 is the schematic diagram that signal P (k) is calculated according to each channel signal of input;
Fig. 2 is the schematic diagram that signal Q (k) is calculated according to each channel signal of input;
Fig. 3 is to calculate 0 passage and the schematic diagram of other channel amplitudes difference according to signal P (k);
Fig. 4 is to calculate 0 passage and the schematic diagram of other phase differences between channels according to signal P (k) and Q (k).
Embodiment
With reference to Fig. 1 to Fig. 4, the present invention is described in further detail.
Refering to Fig. 1 to Fig. 4, a kind of inter-channel phase difference, the measuring method of amplitude difference, include the following steps:A kind of passage
Between phase difference, the measuring method of amplitude difference, include the following steps:
(A) multichannel input signal is denoted as r (k, n), wherein channel number k=0,1,2 ..., K-1, sequential sampling moment n
=0,1,2 ..., N-1, selecting any of which passage, reference channel is denoted as r (m, n), m ∈ [0, K-1] as passage is referred to;
(B) square operation is carried out to each channel signal and mean filter handles to obtain signal P (k):
Wherein k=0,1,2 ..., K-1;
(C) by reference channel with obtaining signal Q after being handled respectively with other channel signal multiplyings and mean filter
(k):
Wherein k=0,1,2 ..., K-1 and k ≠ m, m ∈ [0, k-1];
(D) the amplitude difference A between other passages and reference channel is calculated by signal P (k) and Q (k)Δ(k) and phase
Difference
K=0 in two formula, 1,2 ..., K-1 and k ≠ m, m ∈ [0, k-1];
Mean filter processing in the step B and C is used to filter out the high fdrequency component in signal.
By being operated to the common digital function such as list entries square, multiplication, average, evolution, anticosine, by width
Degree difference and the measurement result of phase difference are provided in the form of formula, are dealt with very simply, conveniently and quickly;Increasing can be passed through
Add the value of N, to increase the precision of result of calculation, so can be achieved with that precision is adjustable, and the scope of application is more extensive;Meanwhile this method
Handled directly against multichannel input signal, without carrying out Digital Down Convert processing, more succinctly, quickly.Above-mentioned survey
Amount method can be realized by software, can also be realized by Field Programmable Gate Array (FPGA);Certainly, optimal mode
It is to be handled by the way of the FPGA combination softwares:P, Q are first obtained after N times of reduction of speed to the Data Stream Processing of input by FPGA
Data, are then uploaded to PC softwares, and final amplitude difference and phase difference are calculated finally by P, Q.Here by from time domain
Measurement overcomes spectrum leakage problem caused by when being realized using FFT, much smaller with hour operation quantity.
The above method is illustrated below by detailed derivation, wherein using the 0th passage as reference channel.
Assuming that original input signal expression formula is:
According to step B,
For high fdrequency component 2f, accumulator is equivalent to low-pass filter, therefore above formula is equivalent to:
According to above formula, the amplitude difference that we can obtain interchannel is measured as:
According to step C,
For high fdrequency component 2f, accumulator is equivalent to low-pass filter, therefore above formula is equivalent to:
Wherein,The average value of phase difference between kth passage and the 0th passage.
From above formula:
The phase of each interchannel we can confirm that, can accurately and rapidly be calculated by step A-D by above-mentioned calculating
Potential difference and amplitude difference.
As the preferred solution of the present invention, in the step A, the 0th passage r (0, n) is selected as referring to passage;Q
(k)、And AΔ(k) k=1 in, 2 ..., K-1.The selection of reference channel very arbitrarily can select the 0th passage to make
It is more convenient a bit for reference channel.
Claims (2)
1. a kind of inter-channel phase difference, the measuring method of amplitude difference, include the following steps:
(A) multichannel input signal is denoted as r (k, n), wherein channel number k=0,1,2 ..., K-1, sequential sampling moment n=0,
1,2 ..., N-1, selecting any of which passage, reference channel is denoted as r (m, n), m ∈ [0, K-1] as passage is referred to;
(B) square operation is carried out to each channel signal and mean filter handles to obtain signal P (k):
Wherein k=0,1,2 ..., K-1;
(C) by reference channel with obtaining signal Q (k) after being handled respectively with other channel signal multiplyings and mean filter:
<mrow>
<mi>Q</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>&lsqb;</mo>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
Wherein k=0,1,2 ..., K-1 and k ≠ m, m ∈ [0, K-1];
(D) the amplitude difference A between other passages and reference channel is calculated by signal P (k) and Q (k)Δ(k) and phase difference
K=0 in two formula, 1,2 ..., K-1 and k ≠ m, m ∈ [0, K-1];
Mean filter processing in the step B and C is used to filter out the high fdrequency component in signal.
2. inter-channel phase difference as claimed in claim 1, the measuring method of amplitude difference, it is characterised in that:The step A
In, select the 0th passage r (0, n) to be used as and refer to passage;Q(k)、And AΔ(k) k=1 in, 2 ..., K-1.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711273036.7A CN108008188A (en) | 2017-12-06 | 2017-12-06 | A kind of inter-channel phase difference, the measuring method of amplitude difference |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711273036.7A CN108008188A (en) | 2017-12-06 | 2017-12-06 | A kind of inter-channel phase difference, the measuring method of amplitude difference |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108008188A true CN108008188A (en) | 2018-05-08 |
Family
ID=62056675
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711273036.7A Pending CN108008188A (en) | 2017-12-06 | 2017-12-06 | A kind of inter-channel phase difference, the measuring method of amplitude difference |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108008188A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108710027A (en) * | 2018-05-23 | 2018-10-26 | 成都玖锦科技有限公司 | The high-precision measuring method of inter-channel phase difference, amplitude difference |
CN113702701A (en) * | 2021-08-09 | 2021-11-26 | 成都玖锦科技有限公司 | Amplitude-phase characteristic measuring method based on comb wave signal |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799704A (en) * | 2010-03-23 | 2010-08-11 | 电子科技大学 | Multichannel DDS signal generator with precise phase control function |
CN203275647U (en) * | 2013-06-17 | 2013-11-06 | 武汉大学 | Multichannel signal amplitude phase difference measuring system |
US20160337066A1 (en) * | 2015-05-11 | 2016-11-17 | Collision Communications, Inc. | Methods, Systems, And Computer Program Products For Calibrating Phase Hardware-Induced Distortion In A Long Term Evolution Communications System |
-
2017
- 2017-12-06 CN CN201711273036.7A patent/CN108008188A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799704A (en) * | 2010-03-23 | 2010-08-11 | 电子科技大学 | Multichannel DDS signal generator with precise phase control function |
CN203275647U (en) * | 2013-06-17 | 2013-11-06 | 武汉大学 | Multichannel signal amplitude phase difference measuring system |
US20160337066A1 (en) * | 2015-05-11 | 2016-11-17 | Collision Communications, Inc. | Methods, Systems, And Computer Program Products For Calibrating Phase Hardware-Induced Distortion In A Long Term Evolution Communications System |
Non-Patent Citations (2)
Title |
---|
师鹏宇 等: "通道间相位差求解算法分析", 《雷达与对抗》 * |
郑胜峰 等: "一种基于多重相关法的相位差测量方法", 《宇航计测技术》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108710027A (en) * | 2018-05-23 | 2018-10-26 | 成都玖锦科技有限公司 | The high-precision measuring method of inter-channel phase difference, amplitude difference |
CN108710027B (en) * | 2018-05-23 | 2020-05-15 | 成都玖锦科技有限公司 | High-precision measurement method for phase difference and amplitude difference between channels |
CN113702701A (en) * | 2021-08-09 | 2021-11-26 | 成都玖锦科技有限公司 | Amplitude-phase characteristic measuring method based on comb wave signal |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104914393A (en) | Apparatus and method for comb spectrum generator phase spectrum calibration | |
CN106610450B (en) | A kind of single-phase electric energy metering chip | |
CN105512431A (en) | Phase noise measuring method based on phase noise mathematical model | |
CN104122444B (en) | All-digital IF spectrum analyzer and frequency spectrum analysis method | |
CN102539944B (en) | Method for measuring phase noise based on phase discrimination method | |
CN104181577A (en) | Beam position and phase measurement system and method based on full digitalization technology | |
CN100461629C (en) | A digital phase-lock amplifier | |
CN104950173B (en) | A kind of passive device equivalent circuit structure and measurement method of parameters and device | |
US9960862B2 (en) | Method and device for detecting standing-wave ratio | |
CN105262503A (en) | Group delay calibration based multipath delay generation device and method | |
CN105510706B (en) | A kind of high-precision lack sampling frequency measuring method | |
CN103969508B (en) | A kind of Electric Power Harmonic Analysis method high-accuracy in real time and device | |
CN108008188A (en) | A kind of inter-channel phase difference, the measuring method of amplitude difference | |
CN103983849A (en) | Real-time high-accuracy power harmonic analysis method | |
CN108710027A (en) | The high-precision measuring method of inter-channel phase difference, amplitude difference | |
CN103823177B (en) | Performance of filter detection method based on THE DESIGN OF WINDOW FUNCTION and system | |
CN104407510A (en) | Method and device for time service | |
CN104483619A (en) | Frequency characteristic testing system based on virtual instrument | |
CN101308175A (en) | Phase spectrum analyzer | |
CN104062498B (en) | The mass center estimation method of the centre frequency of symmetrical spectrum signal | |
CN109088675B (en) | Channel calibration method and device for radio frequency signal source | |
CN106093905B (en) | Precision instrumentation radar frequency deviation classification frequency measuring method based on apFFT phase differences | |
US20220236308A1 (en) | Method for improving emc robustness of integrated capacitive sensors | |
CN103323063B (en) | Ultrasonic flow meter and measurement method for time difference thereof | |
CN106525180B (en) | Method of Calculating Compressibility Factors of Natural Gas |
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: 20180508 |
|
RJ01 | Rejection of invention patent application after publication |