CN105572650A - Broadband multiple correlation flow velocity measurement method - Google Patents
Broadband multiple correlation flow velocity measurement method Download PDFInfo
- Publication number
- CN105572650A CN105572650A CN201510938976.8A CN201510938976A CN105572650A CN 105572650 A CN105572650 A CN 105572650A CN 201510938976 A CN201510938976 A CN 201510938976A CN 105572650 A CN105572650 A CN 105572650A
- Authority
- CN
- China
- Prior art keywords
- sequence
- designated
- signal
- analog
- digital converter
- 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
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/415—Identification of targets based on measurements of movement associated with the target
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
- G01P5/24—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
- G01P5/26—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting optical wave
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Multimedia (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a broadband multiple correlation flow velocity measurement method. Echo signals are sampled by adopting a time-interleaved sampling technology so that a discrete sampling sequence is obtained; multiple correlation algorithm processing and cross spectrum autocorrelation algorithm processing are respectively performed on the discrete sampling sequence; the number of multiple correlation points is increased in the process of multiple correlation algorithm processing, and real-time performance of flow velocity measurement is guaranteed by adopting odd and even paths of parallel multiple correlation processing; and cross spectrum interpolation time delay estimation is performed on the time-interleaved sampled data by cross spectrum autocorrelation algorithm processing so that a time delay value is obtained, and the time delay value is adopted to roughly measure the flow velocity value to act as the judgment standard of a flow velocity size interval, assist multiple correlation algorithm flow velocity measurement and correct the range of flow velocity measurement. Advantages are that multiple correlation algorithm processing and cross spectrum autocorrelation algorithm processing are supplementary for each other, sampling rate of the signals is enhanced by adopting the time-interleaved sampling technology, and multiple correlation algorithm flow velocity measurement is assisted through combination of the roughly measured flow velocity value of the cross spectrum autocorrelation algorithm so that measurement precision is high and measurement range is wide.
Description
Technical field
The present invention relates to a kind of flow-speed measurement method, especially relate to a kind of broadband multiple correlation flow-speed measurement method.
Background technology
Along with resource is day by day deficient, exploitation maritime waters resource becomes the urgent demand of economic development.The flow velocity of maritime waters is the important parameter of maritime waters environment, is also the important evidence of maritime waters activity.Existing flow-speed measurement method is mainly based on broadband multiple correlation technology, and document " research of broadband doppler current meter speed-measuring method " proposes to study based on the flow velocity of broadband multiple correlation algorithm.In the document, the detailed process of fluid-velocity survey is: first utilize pulse signal generator to produce pure-tone polse signal, adopt M sequence code to carry out binary coding to this pure-tone polse signal and obtain subpulse signal, multiple subpulse signal is carried out splicing successively and obtains modulated pulse signal; Receive the signal that this modulated pulse signal is reflected back through maritime waters, and collect the echoed signal of this signal; Finally adopt broadband multiple correlation algorithm to process this echoed signal, obtain flow velocity.
In this flow-speed measurement method, echoed signal is the random signal of a Noise, the propagation loss etc. of the intensity of echoed signal, noise and echoed signal all affects the precision of fluid-velocity survey, multiple correlation algorithm can the certain white noise of filtering, there is good antijamming capability, this flow-speed measurement method precision is higher thus, but multiple correlation algorithm is in the process processed echoed signal, the phase information of the compound correlative function of the echoed signal obtained has periodically, the measurement range of flow velocity is restricted, and the method range is narrower thus.
Summary of the invention
It is high that technical matters to be solved by this invention is to provide a kind of measuring accuracy, and measure the broadband multiple correlation flow-speed measurement method of broad quantum.
The present invention solves the problems of the technologies described above adopted technical scheme: a kind of broadband multiple correlation flow-speed measurement method, comprises the following steps:
1. echoed signal is obtained
1.-1 utilize pulse signal generator generated frequency for f
0pure-tone polse signal, f
0=1.5MHz;
1.-2 adopt that code element numbers are L, symbol width is that the M sequence code of Δ τ carries out binary coding to pure-tone polse signal and obtains subpulse signal, code element number L=15, symbol width Δ τ is 10 pure-tone polse signal periods of filling,
the length of subpulse signal is designated as T, T=L × Δ τ, symbol × be multiplication symbol;
1.-3 by multiple subpulse signal sequential filming successively in the maritime waters for the treatment of velocity measurement, and multiple subpulse signals of sequential filming form a modulated pulse signal successively, and this modulated pulse signal is designated as s (t), and t is the time domain time;
1.-4 the signal that is reflected back through maritime waters of these modulated pulse signals s (t) is obtained, the signal this be reflected back is designated as s'(t), according to s'(t) obtain pending echoed signal rr (t), rr (t)=s'(t)+δ (t); Wherein δ (t) is white noise signal;
2. successively enlarge leadingly process, bandpass filtering treatment and low-pass filtering treatment are carried out to echoed signal rr (t), obtain pre-sampling signal, pre-sampling signal is designated as r (t);
3. adopt alternating sampling module to carry out four road time-interleaveds to pre-sampling signal r (t) and obtain discrete sampling sequence r (n), n is discrete series number, n=1,2,3, alternating sampling module comprises the first analog to digital converter, second analog to digital converter, 3rd analog to digital converter, 4th analog to digital converter and four selects a selector switch, first analog to digital converter, second analog to digital converter, 3rd analog to digital converter and the 4th analog to digital converter and four select four of a selector switch selecting sides to connect one to one, first analog to digital converter, second analog to digital converter, 3rd analog to digital converter and the 4th analog to digital converter are sampled to pre-sampling signal r (t), four select a selector switch to select the first analog to digital converter successively, second analog to digital converter, 3rd analog to digital converter and the 4th analog to digital converter, first analog to digital converter, second analog to digital converter, the sampled signal of the 3rd analog to digital converter and the 4th analog to digital converter is selected a selector switch to export successively through four and is spliced into discrete sampling sequence r (n) according to the order of sequence,
4. multiple correlation algorithm process and the process of cross-spectrum auto-correlation algorithm are done respectively to discrete sampling sequence r (n);
The detailed process of multiple correlation algorithm process is:
A. adopt DDS module to produce orthogonal demodulation signal, adopt this orthogonal demodulation signal to do digital modulation to discrete sampling sequence r (n), obtain two-way quadrature signal sequence, be designated as I (n) and Q (n) respectively; I (n) and Q (n) is adopted to construct the first sequence of complex numbers X (n), using the solid part signal of I (n) as the first sequence of complex numbers X (n), Q (n) is as the imaginary signals of the first sequence of complex numbers X (n), then the first sequence of complex numbers represents X (n)=I (n)+jQ (n), and wherein j represents the imaginary part unit of plural number;
B. the first sequence of complex numbers X (n) is divided into two sequence of complex numbers according to odd and even number discrete series number, the sequence of complex numbers being odd number by discrete series number is called the second sequence of complex numbers, be designated as X (2m+1), the sequence of complex numbers being even number by discrete series number is called the 3rd sequence of complex numbers, is designated as X (2m+2), m=0,1,2,3,
C. the compound correlative function of the first sequence of complex numbers X (n) is designated as R
1n (), is designated as R by the compound correlative function of the second sequence of complex numbers X (2m+1)
2n (), is designated as R by the compound correlative function of the 3rd sequence of complex numbers X (2m+2)
3(n),
Wherein N' represents the multiple correlation length of the second sequence of complex numbers X (2m+1), N " represent the multiple correlation length of the 3rd sequence of complex numbers X (2m+2), and
l is sample delay value and its value is the code length of 1 M sequence code, N=f
s× (t
1-t
0), f
sfor the sample frequency of alternating sampling module, t
0for the sampling start time of alternating sampling module, t
1the sampling finish time of alternating sampling module, X (2m+1+l)=I (2m+1+l)+jQ (2m+1+l), X (2m+2+l)=I (2m+2+l)+jQ (2m+2+l), I (2m+1+l) is the expression formula of I (2m+1) time delay l, I (2m+2+l) is the expression formula of I (2m+2) time delay l, Q (2m+1+l) is the expression formula of Q (2m+1) time delay l, and Q (2m+2+l) is the expression formula of Q (2m+2) time delay l; Symbol * is conjugate operation symbol, X
*(2m+1) be the conjugate complex number sequence of X (2m+1), X
*(2m+2) be the conjugate complex number sequence of X (2m+2);
D. by R
2(n) and R
3n () substitutes into formula R
1(n)=R
2(n)+R
3r is obtained in (n)
1n (), by R
1n the real part of () is designated as Re (R
1(n)), just R
1n the imaginary part of () is designated as Im (R
1(n));
E. by the compound correlative function R of the first sequence of complex numbers X (n)
1n the phasing degree of () is designated as θ, adopt arctan function to solve phasing degree θ, namely
F. Doppler frequency deviation is designated as f
d, adopt formula
wherein q is relevant repetition period, i.e. q=f
s× T;
G. the flow velocity that multiple correlation is measured is designated as v
c, adopt formula
calculate v
c, c is the velocity of sound in water, c=1500m/s;
The detailed process of cross-spectrum auto-correlation algorithm process is:
A. the frequency spectrum of discrete sampling sequence r (n) is designated as r'(k), the autocorrelation function frequency spectrum of discrete sampling sequence r (n) is designated as R (k), and k represents the sequence number of frequency domain, k=0, and 1,2,3,
B. FFT computing is adopted to calculate frequency spectrum r'(k), make R (k)=r'(k) × r'(k);
C. carry out IFFT computing again after interpolation processing being done to autocorrelation function frequency spectrum R (k), obtain the autocorrelation function of discrete sampling sequence r (n), the autocorrelation function of discrete sampling sequence r (n) is designated as R'(t);
D. by the autocorrelation function R'(t of discrete sampling sequence r (n)) determine the moment that the first secondary lobe occurs the moment that the first secondary lobe occurs to be designated as T
1; The moment that ideally the first secondary lobe occurs that frequency deviation does not occur is designated as T
2, wherein T
2=T
E. the difference that Doppler's time delay occurs is designated as Δ T, Δ T=T
1-T
2;
F. the flow velocity of autocorrelation measurement is designated as v
z, adopt formula
calculate the flow velocity v of autocorrelation measurement
z;
5. flow velocity corrected parameter is designated as α, α=[(v
z-v
c)/2v
max], v
maxfor fuzzy speed, v
max=c/4f
0t, symbol " [] " is rounding operation symbol;
6. formula v=v is adopted
c+ α 2v
maxcalculate and treat velocity measurement v.
Compared with prior art, the invention has the advantages that adopting time-interleaved technology to carry out sampling to echoed signal obtains discrete sampling sequence, respectively multiple correlation algorithm process and the process of cross-spectrum auto-correlation algorithm are done to discrete sampling sequence in addition, time-interleaved improves sampling rate under the prerequisite ensureing sampling precision, in multiple correlation algorithm process process, the corresponding multiple correlation that adds is counted, and adopt odd even two-way to walk abreast multiple correlation calculation process, on the basis of flow measurement precision of improving multiple correlation algorithm, improve computing velocity, ensure the real-time of fluid-velocity survey, the process of cross-spectrum auto-correlation algorithm is by making cross-spectrum interpolation Delay Estima-tion to time-interleaved data, improve the estimated accuracy of time delay value, reduce the error in judgement of the large minizone of flow velocity, realize wide-range fluid-velocity survey, adopt this delay value bigness scale flow speed value, there is not measurement fuzzy problem in this bigness scale flow speed value, can be used as the criterion of the large minizone of flow speed value, auxiliary multiple correlation algorithm fluid-velocity survey, overcome the confinement problems that multiple correlation algorithm fluid-velocity survey may exist fuzzy value and fluid-velocity survey range, revise fluid-velocity survey range, multiple correlation algorithm process and the process of cross-spectrum auto-correlation algorithm combine thus, by the sampling rate adopting time-interleaved technology to improve signal, flow speed value in conjunction with the bigness scale of cross-spectrum auto-correlation algorithm assists multiple correlation algorithm to carry out fluid-velocity survey, measuring accuracy is high, and measurement broad quantum.
Embodiment
Below in conjunction with embodiment, the present invention is described in further detail.
Embodiment: a kind of broadband multiple correlation flow-speed measurement method, comprises the following steps:
1. echoed signal is obtained
1.-1 utilize pulse signal generator generated frequency for f
0pure-tone polse signal, f
0=1.5MHz;
1.-2 adopt that code element numbers are L, symbol width is that the M sequence code of Δ τ carries out binary coding to pure-tone polse signal and obtains subpulse signal, code element number L=15, symbol width Δ τ is 10 pure-tone polse signal periods of filling,
the length of subpulse signal is designated as T, T=L × Δ τ, symbol × be multiplication symbol;
1.-3 by multiple (namely at least two) subpulse signal sequential filming successively in the maritime waters for the treatment of velocity measurement, multiple subpulse signals of sequential filming form a modulated pulse signal successively, this modulated pulse signal is designated as s (t), and t is the time domain time;
1.-4 the signal that is reflected back through maritime waters of these modulated pulse signals s (t) is obtained, the signal this be reflected back is designated as s'(t), according to s'(t) obtain pending echoed signal rr (t), rr (t)=s'(t)+δ (t); Wherein δ (t) is white noise signal;
2. successively enlarge leadingly process, bandpass filtering treatment and low-pass filtering treatment are carried out to echoed signal rr (t), obtain pre-sampling signal, pre-sampling signal is designated as r (t);
3. adopt alternating sampling module to carry out four road time-interleaveds to pre-sampling signal r (t) and obtain discrete sampling sequence r (n), n is discrete series number, n=1,2,3, alternating sampling module comprises the first analog to digital converter, second analog to digital converter, 3rd analog to digital converter, 4th analog to digital converter and four selects a selector switch, first analog to digital converter, second analog to digital converter, 3rd analog to digital converter and the 4th analog to digital converter and four select four of a selector switch selecting sides to connect one to one, first analog to digital converter, second analog to digital converter, 3rd analog to digital converter and the 4th analog to digital converter are sampled to pre-sampling signal r (t), four select a selector switch to select the first analog to digital converter successively, second analog to digital converter, 3rd analog to digital converter and the 4th analog to digital converter, first analog to digital converter, second analog to digital converter, the sampled signal of the 3rd analog to digital converter and the 4th analog to digital converter is selected a selector switch to export successively through four and is spliced into discrete sampling sequence r (n) according to the order of sequence,
4. multiple correlation algorithm process and the process of cross-spectrum auto-correlation algorithm are done respectively to discrete sampling sequence r (n);
The detailed process of multiple correlation algorithm process is:
A. DDS (DirectDigitalSynthesizer is adopted, Direct Digital Synthesizer) module generation orthogonal demodulation signal, this orthogonal demodulation signal is adopted to do digital modulation to discrete sampling sequence r (n), obtain two-way quadrature signal sequence, be designated as I (n) and Q (n) respectively; I (n) and Q (n) is adopted to construct the first sequence of complex numbers X (n), using the solid part signal of I (n) as the first sequence of complex numbers X (n), Q (n) is as the imaginary signals of the first sequence of complex numbers X (n), then the first sequence of complex numbers represents X (n)=I (n)+jQ (n), and wherein j represents the imaginary part unit of plural number;
B. the first sequence of complex numbers X (n) is divided into two sequence of complex numbers according to odd and even number discrete series number, the sequence of complex numbers being odd number by discrete series number is called the second sequence of complex numbers, be designated as X (2m+1), the sequence of complex numbers being even number by discrete series number is called the 3rd sequence of complex numbers, is designated as X (2m+2), m=0,1,2,3,
C. the compound correlative function of the first sequence of complex numbers X (n) is designated as R
1n (), is designated as R by the compound correlative function of the second sequence of complex numbers X (2m+1)
2n (), is designated as R by the compound correlative function of the 3rd sequence of complex numbers X (2m+2)
3(n),
Wherein N' represents the multiple correlation length of the second sequence of complex numbers X (2m+1), N " represent the multiple correlation length of the 3rd sequence of complex numbers X (2m+2), and
l is sample delay value and its value is the code length of 1 M sequence code, N=f
s× (t
1-t
0), f
sfor the sample frequency of alternating sampling module, t
0for the sampling start time of alternating sampling module, t
1the sampling finish time of alternating sampling module, f
s, t
0and t
1for adopting the parameter of existing conventional method to set up setting in alternating sampling module, X (2m+1+l)=I (2m+1+l)+jQ (2m+1+l), X (2m+2+l)=I (2m+2+l)+jQ (2m+2+l), I (2m+1+l) is the expression formula of I (2m+1) time delay l, I (2m+2+l) is the expression formula of I (2m+2) time delay l, Q (2m+1+l) is the expression formula of Q (2m+1) time delay l, and Q (2m+2+l) is the expression formula of Q (2m+2) time delay l; Symbol * is conjugate operation symbol, X
*(2m+1) be the conjugate complex number sequence of X (2m+1), X
*(2m+2) be the conjugate complex number sequence of X (2m+2);
D. by R
2(n) and R
3n () substitutes into formula R
1(n)=R
2(n)+R
3r is obtained in (n)
1n (), by R
1n the real part of () is designated as Re (R
1(n)), just R
1n the imaginary part of () is designated as Im (R
1(n));
E. by the compound correlative function R of the first sequence of complex numbers X (n)
1n the phasing degree of () is designated as θ, adopt arctan function to solve phasing degree θ, namely
F. Doppler frequency deviation is designated as f
d, adopt formula
wherein q is relevant repetition period, i.e. q=f
s× T;
G. the flow velocity that multiple correlation is measured is designated as v
c, adopt formula
calculate v
c, c is the velocity of sound in water, c=1500m/s;
The detailed process of cross-spectrum auto-correlation algorithm process is:
A. the frequency spectrum of discrete sampling sequence r (n) is designated as r'(k), the autocorrelation function frequency spectrum of discrete sampling sequence r (n) is designated as R (k), and k represents the sequence number of frequency domain, k=0, and 2,3,
B. FFT computing is adopted to calculate frequency spectrum r'(k), make R (k)=r'(k) × r'(k);
C. carry out IFFT computing again after interpolation processing being done to autocorrelation function frequency spectrum R (k), obtain the autocorrelation function of discrete sampling sequence r (n), the autocorrelation function of discrete sampling sequence r (n) is designated as R'(t);
D. by the autocorrelation function R'(t of discrete sampling sequence r (n)) determine the moment that the first secondary lobe occurs the moment that the first secondary lobe occurs to be designated as T
1; The moment that ideally the first secondary lobe occurs that frequency deviation does not occur is designated as T
2, wherein T
2=T
E. the difference that Doppler's time delay occurs is designated as Δ T, Δ T=T
1-T
2;
F. the flow velocity of autocorrelation measurement is designated as v
z, adopt formula
calculate the flow velocity v of autocorrelation measurement
z;
5. flow velocity corrected parameter is designated as α, α=[(v
z-v
c)/2v
max], v
maxfor fuzzy speed, v
max=c/4f
0t, symbol " [] " is rounding operation symbol;
6. formula v=v is adopted
c+ α 2v
maxcalculate and treat velocity measurement v.
In the present embodiment, alternating sampling module and DDS module all adopt the matured product of its technical field.
Claims (1)
1. a broadband multiple correlation flow-speed measurement method, is characterized in that comprising the following steps:
1. echoed signal is obtained
1.-1 utilize pulse signal generator generated frequency for f
0pure-tone polse signal, f
0=1.5MHz;
1.-2 adopt that code element numbers are L, symbol width is that the M sequence code of Δ τ carries out binary coding to pure-tone polse signal and obtains subpulse signal, code element number L=15, symbol width Δ τ is 10 pure-tone polse signal periods of filling,
the length of subpulse signal is designated as T, T=L × Δ τ, symbol × be multiplication symbol;
1.-3 by multiple subpulse signal sequential filming successively in the maritime waters for the treatment of velocity measurement, and multiple subpulse signals of sequential filming form a modulated pulse signal successively, and this modulated pulse signal is designated as s (t), and t is the time domain time;
1.-4 the signal that is reflected back through maritime waters of these modulated pulse signals s (t) is obtained, the signal this be reflected back is designated as s'(t), according to s'(t) obtain pending echoed signal rr (t), rr (t)=s'(t)+δ (t); Wherein δ (t) is white noise signal;
2. successively enlarge leadingly process, bandpass filtering treatment and low-pass filtering treatment are carried out to echoed signal rr (t), obtain pre-sampling signal, pre-sampling signal is designated as r (t);
3. adopt alternating sampling module to carry out four road time-interleaveds to pre-sampling signal r (t) and obtain discrete sampling sequence r (n), n is discrete series number, n=1,2,3, alternating sampling module comprises the first analog to digital converter, second analog to digital converter, 3rd analog to digital converter, 4th analog to digital converter and four selects a selector switch, first analog to digital converter, second analog to digital converter, 3rd analog to digital converter and the 4th analog to digital converter and four select four of a selector switch selecting sides to connect one to one, first analog to digital converter, second analog to digital converter, 3rd analog to digital converter and the 4th analog to digital converter are sampled to pre-sampling signal r (t), four select a selector switch to select the first analog to digital converter successively, second analog to digital converter, 3rd analog to digital converter and the 4th analog to digital converter, first analog to digital converter, second analog to digital converter, the sampled signal of the 3rd analog to digital converter and the 4th analog to digital converter is selected a selector switch to export successively through four and is spliced into discrete sampling sequence r (n) according to the order of sequence,
4. multiple correlation algorithm process and the process of cross-spectrum auto-correlation algorithm are done respectively to discrete series r (n);
The detailed process of multiple correlation algorithm process is:
A. adopt DDS module to produce orthogonal demodulation signal, adopt this orthogonal demodulation signal to do digital modulation to discrete sampling sequence r (n), obtain two-way quadrature signal sequence, be designated as I (n) and Q (n) respectively; I (n) and Q (n) is adopted to construct the first sequence of complex numbers X (n), using the solid part signal of I (n) as the first sequence of complex numbers X (n), Q (n) is as the imaginary signals of the first sequence of complex numbers X (n), then the first sequence of complex numbers represents X (n)=I (n)+jQ (n), and wherein j represents the imaginary part unit of plural number;
B. the first sequence of complex numbers X (n) is divided into two sequence of complex numbers according to odd and even number discrete series number, the sequence of complex numbers being odd number by discrete series number is called the second sequence of complex numbers, be designated as X (2m+1), the sequence of complex numbers being even number by discrete series number is called the 3rd sequence of complex numbers, is designated as X (2m+2), m=0,1,2,3,
C. the compound correlative function of the first sequence of complex numbers X (n) is designated as R
1n (), is designated as R by the compound correlative function of the second sequence of complex numbers X (2m+1)
2n (), is designated as R by the compound correlative function of the 3rd sequence of complex numbers X (2m+2)
3(n),
Wherein N' represents the multiple correlation length of the second sequence of complex numbers X (2m+1), N " represent the multiple correlation length of the 3rd sequence of complex numbers X (2m+2), and
l is sample delay value and its value is the code length of 1 M sequence code, N=f
s× (t
1-t
0), f
sfor the sample frequency of alternating sampling module, t
0for the sampling start time of alternating sampling module, t
1the sampling finish time of alternating sampling module, X (2m+1+l)=I (2m+1+l)+jQ (2m+1+l), X (2m+2+l)=I (2m+2+l)+jQ (2m+2+l), I (2m+1+l) is the expression formula of I (2m+1) time delay l, I (2m+2+l) is the expression formula of I (2m+2) time delay l, Q (2m+1+l) is the expression formula of Q (2m+1) time delay l, and Q (2m+2+l) is the expression formula of Q (2m+2) time delay l; Symbol * is conjugate operation symbol, X
*(2m+1) be the conjugate complex number sequence of X (2m+1), X
*(2m+2) be the conjugate complex number sequence of X (2m+2);
D. by R
2(n) and R
3n () substitutes into formula R
1(n)=R
2(n)+R
3r is obtained in (n)
1n (), by R
1n the real part of () is designated as Re (R
1(n)), just R
1n the imaginary part of () is designated as Im (R
1(n));
E. by the compound correlative function R of the first sequence of complex numbers X (n)
1n the phasing degree of () is designated as θ, adopt arctan function to solve phasing degree θ, namely
F. Doppler frequency deviation is designated as f
d, adopt formula
wherein q is relevant repetition period, i.e. q=f
s× T;
G. the flow velocity that multiple correlation is measured is designated as v
c, adopt formula
calculate v
c, c is the velocity of sound in water, c=1500m/s;
The detailed process of cross-spectrum auto-correlation algorithm process is:
A. the frequency spectrum of discrete sampling sequence r (n) is designated as r'(k), the autocorrelation function frequency spectrum of discrete sampling sequence r (n) is designated as R (k), and k represents the sequence number of frequency domain, k=0, and 1,2,3,
B. FFT computing is adopted to calculate frequency spectrum r'(k), make R (k)=r'(k) × r'(k);
C. carry out IFFT computing again after interpolation processing being done to autocorrelation function frequency spectrum R (k), obtain the autocorrelation function of discrete sampling sequence r (n), the autocorrelation function of discrete series r (n) is designated as R'(t);
D. by the autocorrelation function R'(t of discrete sampling sequence r (n)) determine the moment that the first secondary lobe occurs the moment that the first secondary lobe occurs to be designated as T
1; The moment that ideally the first secondary lobe occurs that frequency deviation does not occur is designated as T
2, wherein T
2=T
E. the difference that Doppler's time delay occurs is designated as Δ T, Δ T=T
1-T
2;
F. the flow velocity of autocorrelation measurement is designated as v
z, adopt formula
Δ T calculates the flow velocity v of autocorrelation measurement
z;
5. flow velocity corrected parameter is designated as α, α=[(v
z-v
c)/2v
max], v
maxfor fuzzy speed, v
max=c/4f
0t, symbol " [] " is rounding operation symbol;
6. formula v=v is adopted
c+ α 2v
maxcalculate and treat velocity measurement v.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510938976.8A CN105572650A (en) | 2015-12-15 | 2015-12-15 | Broadband multiple correlation flow velocity measurement method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510938976.8A CN105572650A (en) | 2015-12-15 | 2015-12-15 | Broadband multiple correlation flow velocity measurement method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105572650A true CN105572650A (en) | 2016-05-11 |
Family
ID=55882999
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510938976.8A Pending CN105572650A (en) | 2015-12-15 | 2015-12-15 | Broadband multiple correlation flow velocity measurement method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105572650A (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106248991A (en) * | 2016-07-15 | 2016-12-21 | 南京世海声学科技有限公司 | ADCP solution under a kind of low signal-to-noise ratio tests the speed blur method |
| CN107870034A (en) * | 2017-10-24 | 2018-04-03 | 宁波大学科学技术学院 | A kind of underwater sound sonic velocity measurement method based on phase difference |
| CN108737305A (en) * | 2018-05-29 | 2018-11-02 | 西安交通大学 | A kind of frequency difference UNB method for communication transmission of overlapping square phase-couple |
| CN108828602A (en) * | 2018-03-06 | 2018-11-16 | 北京大学 | A kind of pulsion phase dry method tests the speed the fuzzy signal processing method of middle release rate |
| CN109270295A (en) * | 2018-08-20 | 2019-01-25 | 南京世海声学科技有限公司 | A kind of underwater sound Doppler's flow-speed measurement method screened based on autocorrelation estimation and valid data |
| CN111352099A (en) * | 2018-12-20 | 2020-06-30 | 宁波大学科学技术学院 | Time delay estimation method based on cross-correlation signal phase decomposition |
| CN114324959A (en) * | 2021-12-03 | 2022-04-12 | 中国船舶重工集团公司七五0试验场 | FPGA and ARM-based complex correlation method frequency measurement method and system |
| CN114374407A (en) * | 2022-01-10 | 2022-04-19 | 哈尔滨工程大学 | Spatial channel characteristic prediction method and system based on m sequence and storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102928619A (en) * | 2012-11-05 | 2013-02-13 | 杭州电子科技大学 | Signal processing method of broadband beam-control acoustic Doppler flow testing system |
| CN105021843A (en) * | 2015-07-28 | 2015-11-04 | 江苏中海达海洋信息技术有限公司 | 600kHZ broadband acoustics Doppler current profiler and realization method |
-
2015
- 2015-12-15 CN CN201510938976.8A patent/CN105572650A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102928619A (en) * | 2012-11-05 | 2013-02-13 | 杭州电子科技大学 | Signal processing method of broadband beam-control acoustic Doppler flow testing system |
| CN105021843A (en) * | 2015-07-28 | 2015-11-04 | 江苏中海达海洋信息技术有限公司 | 600kHZ broadband acoustics Doppler current profiler and realization method |
Non-Patent Citations (4)
| Title |
|---|
| 叶忠辉 等: ""基于FPGA的自适应TIADC数据采集"", 《数学的实践与认识》 * |
| 韩礼波 等: ""宽带多普勒海流计测速方法研究"", 《声学与电子工程》 * |
| 韩礼波 等: ""宽带海流计速度解模糊方法研究"", 《声学与电子工程》 * |
| 饶中洋 等: ""宽带ADCP的复相关算法仿真"", 《海洋技术》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106248991B (en) * | 2016-07-15 | 2018-12-21 | 南京世海声学科技有限公司 | A kind of ADCP solution under low signal-to-noise ratio tests the speed blur method |
| CN106248991A (en) * | 2016-07-15 | 2016-12-21 | 南京世海声学科技有限公司 | ADCP solution under a kind of low signal-to-noise ratio tests the speed blur method |
| CN107870034A (en) * | 2017-10-24 | 2018-04-03 | 宁波大学科学技术学院 | A kind of underwater sound sonic velocity measurement method based on phase difference |
| CN107870034B (en) * | 2017-10-24 | 2019-12-24 | 宁波大学科学技术学院 | Underwater acoustic velocity measurement method based on phase difference |
| CN108828602A (en) * | 2018-03-06 | 2018-11-16 | 北京大学 | A kind of pulsion phase dry method tests the speed the fuzzy signal processing method of middle release rate |
| CN108737305A (en) * | 2018-05-29 | 2018-11-02 | 西安交通大学 | A kind of frequency difference UNB method for communication transmission of overlapping square phase-couple |
| CN109270295A (en) * | 2018-08-20 | 2019-01-25 | 南京世海声学科技有限公司 | A kind of underwater sound Doppler's flow-speed measurement method screened based on autocorrelation estimation and valid data |
| CN109270295B (en) * | 2018-08-20 | 2021-03-30 | 南京世海声学科技有限公司 | Underwater acoustic Doppler flow velocity measurement method based on autocorrelation estimation and effective data screening |
| CN111352099A (en) * | 2018-12-20 | 2020-06-30 | 宁波大学科学技术学院 | Time delay estimation method based on cross-correlation signal phase decomposition |
| CN114324959A (en) * | 2021-12-03 | 2022-04-12 | 中国船舶重工集团公司七五0试验场 | FPGA and ARM-based complex correlation method frequency measurement method and system |
| CN114324959B (en) * | 2021-12-03 | 2024-06-11 | 中国船舶重工集团公司七五0试验场 | Complex correlation method frequency measurement method and system based on FPGA and ARM |
| CN114374407A (en) * | 2022-01-10 | 2022-04-19 | 哈尔滨工程大学 | Spatial channel characteristic prediction method and system based on m sequence and storage medium |
| CN114374407B (en) * | 2022-01-10 | 2024-03-08 | 哈尔滨工程大学 | Spatial channel characteristic prediction method, system and storable medium based on m sequence |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105572650A (en) | Broadband multiple correlation flow velocity measurement method | |
| KR101779315B1 (en) | Radar operation method with increased doppler capability and system thereabout | |
| CN102628937B (en) | Radar detection method based on generalized keystone transformation and non-coherent accumulation | |
| CN105229431A (en) | The level gauging that the distance with improvement is determined | |
| EP1098206A3 (en) | Radar system and coherent integrating method therefor | |
| CN113640752B (en) | Waveform design method based on inter-pulse phase frequency spectrum double agility | |
| EP3771921B1 (en) | Method for processing non-interrupted phase synchronization signal of bistatic sar based on coded signal | |
| CN108562883B (en) | Maximum likelihood distance estimation algorithm of multi-carrier radar system | |
| JP4828295B2 (en) | Doppler measuring instrument and tide meter | |
| CN102226839B (en) | Estimation method for time delay of line scanning pulse with low sampling rate | |
| CN105738889A (en) | Frequency modulated continuous wave speed measurement and distance measurement method | |
| CN108828602B (en) | Signal processing method for eliminating velocity ambiguity in pulse phase dry method velocity measurement | |
| CN107390210B (en) | Digital processing method of beat signal in material level measurement | |
| CN115508820A (en) | Target detection method of linear frequency modulation pulse radar | |
| CN103728464A (en) | Method for combined pulse speed measurement for acoustic Doppler flow rate section plotter | |
| CN103336275B (en) | A kind of ambiguity solution method of step frequency pulsed radar signal fine motion detection | |
| CN113640753B (en) | LFM pulse train signal waveform design method based on pulse width agility | |
| CN112130140B (en) | Method for improving speed and distance measurement precision of pseudo code phase modulation fuze system | |
| JP4077092B2 (en) | Doppler frequency measurement method and Doppler sonar | |
| CN113627398B (en) | Signal characteristic detection method based on adaptive reconstruction filtering | |
| CN106248991B (en) | A kind of ADCP solution under low signal-to-noise ratio tests the speed blur method | |
| CN115586507A (en) | MIMO radar waveform design method for inhibiting speed ambiguity | |
| CN104020450B (en) | A kind of radar return reconstructing method being layered based on adjacent wave correlation and phase | |
| CN102998659A (en) | Doppler frequency spectrum shaping method and system based on interpulse modulation | |
| JP4249332B2 (en) | Frequency measuring method and Doppler sonar |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160511 |
|
| RJ01 | Rejection of invention patent application after publication |