CN101251596A - Wind profile radar non-phase parameter radio interference inhibition method - Google Patents

Abstract

The invention discloses a non-coherent radio interference suppression method of a wind profile radar, belonging to the signal processing and detection technical field. The method comprises the following steps: an interpulse pseudo random sequence and an intrapulse pseudo random sequence are generated at a transmitting end; coding processing is completed; then, modulation transmitting is completed and both pseudo random sequences are stored in a received signal processor; finally, receiving is completed at a receiving end and two times of decoding processing of the received signal is completed; thus, the processing procedure of the method is completed. The primary decoding operation of the received signal of a wind profile radar processed by the method is carried out through the interpulse pseudo random sequence synchronical to a transmitting modulation pulse, thereby restoring coherent wind measurement signal and converting non-coherent radio interference into white noise so as to reach the effect of reducing the power spectrum of radio interference signal; then, according to an intrapulse pseudo random codebook sequence, secondary decoding operation is carried out to reach the effect of increasing the power of wind signal; moreover, the interference suppression effect of the method is about (10log L)dB and the power of wind signal is increased about (20log P)dB.

Description

A kind of wind profile radar non-phase parameter radio interference inhibition method

Technical field

The present invention is a kind of radio interference inhibition method, and particularly a kind of wind profile radar non-phase parameter radio interference inhibition method belongs to signal Processing and detection technique field.

Background technology

Wind profile radar is to the Doppler remote sensing equipment of atmosphere wind direction and wind velocity with a kind of weather proof type of height profile Continuous Observation, and it realizes surveying the purpose of wind by the feeble signal that detects the atmospheric turbulence reflection.This Study on Technology starts from the U.S.'s 20th century early eighties, because it has higher survey wind real-time and precision, it is one of gordian technique of Small and Medium Sized atmospheric exploration, in recent years, wind profile radar more and more is subjected to the attention of relevant departments such as various countries' meteorology, environmental protection, military affairs, also is about to carry out relevant application and technical research in fields such as China's manned space flight, Eleventh Five-Year Plan hazard weather country's monitoring net and national defense construction.The working frequency range of wind profile radar is subjected to the influence of desired height coverage and resolution, divides according to International Telecommunications Union (ITU) and determines near 3 frequency range: 50MHz, 400MHz and 1000MHz.Because more approaching with the working frequency range of equipment such as broadcast communication, in addition wind profile radar echoed signal power weak (even less than-140dBm), so be subjected to the interference of radiowave easily.Wind profile radar mainly adopts classical signal processing technology at present, and for example pulse compression encodes, goes direct current, spectrum to accumulate, spectrum offsets etc., though Active Jamming is had certain inhibiting effect, signal interference ratio improves not quite, DeGrain.

Summary of the invention

The object of the present invention is to provide a kind of wind profile radar non-phase parameter radio interference inhibition method, this method can effectively solve the non-phase parameter radio perturbation technique problem of wind profile radar practical application.

Technology contents of the present invention is:

Because the pulse number that the wind profile radar system will send in a detect cycle when surveying usually is up to hundreds thousand of, the present invention proposes the method for associating pseudorandomcode inhibition non-phase parameter radio wave interference in the full coherent intercadence based on this work characteristics.Concrete steps are:

Step 1: at the wind profile radar transmitting terminal, produce between arteries and veins synchronously, two kinds of pseudo-random code sequences in the arteries and veins; Wherein pseudo-random sequence adopts the M sequence between arteries and veins, and sequence length is selected according to the real work parameter, supposes that its code length is L, and coded sequence adopts Barker code in the arteries and veins, supposes that its code length is P,

If the code book of pseudorandomcode is R between arteries and veins:

R＝[R ₁?R ₂?…?R _i?…?R _L] ^T (1)

R in the formula _iBe M sequence pseudo-random code code element, i=1,2 ... L, T representing matrix transposition;

The coding code book is a Barker code in the arteries and veins, and its figure place is P,

BK＝[bk(1)?bk(2)?…?bk(i)?…?bk(P)] ^T (2)

Code element is bk (i) in the formula, i=1,2 ... P, T representing matrix transposition;

Step 2: at the wind profile radar transmitting terminal, above-mentioned formula (1) and formula (2) code book sequence through the inline yard sequence of compiling in collaboration with of intercadence are:

$\mathrm{RB}=\left[\begin{array}{c}{\mathrm{RB}}_1\\ {\mathrm{RB}}_2\\ .\\ .\\ .\\ {\mathrm{RB}}_i\\ .\\ .\\ .\\ {\mathrm{RB}}_L\end{array}\right]---\left(3\right)$

RB in the formula _iBe the code book signal in i the pulse in back of encoding, i=1,2 ... L,

${\mathrm{RB}}_i=(\mathrm{BK}+R_i)\mathrm{mod}\left(2\right)=\left[\begin{array}{c}\mathrm{bk}\left(1\right)\mathrm{xor}{R}_i\\ \mathrm{bk}\left(2\right)\mathrm{xor}{R}_i\\ .\\ .\\ .\\ \mathrm{bk}\left(P\right)\mathrm{xor}{R}_i\end{array}\right]---\left(4\right)$

Mod () is the complementation computing in the formula, and xor is an xor operator;

Step 3: the code book sequence that formula (3) is produced is sent to wind profile radar transmit modulator and emission, and the code book sequence with formula (1) and formula (2) generation is kept in the received signal processor simultaneously;

Step 4: transmitting in the step 3 received at the wind profile radar receiving end;

If time domain accumulation is that FFT counts and is N J time, the range gate number is G, and then the signal processor sample of signal space I that need handle is:

$I={\left[\begin{array}{cccc}{A}_{11}& A_{12}& ...& A_{1J}\\ A_{21}& A_{22}& ...& A_{2J}\\ .& .& & .\\ .& .& A_{\mathrm{ij}}& .\\ .& .& & .\\ A_{N1}& A_{N2}& ...& A_{\mathrm{NJ}}\end{array}\right]}_{N\×J}---\left(5\right)$

A in the formula _IjRepresent i * j sampled signal a respectively _Ij(1), a _Ij(2) ... a _Ij(G+P-1) sample space of Gou Chenging, i=1 wherein, 2 ... N, j=1,2 ... J respectively with signal Processing in FFT count relevant with the coherent accumulation number of times, L=N * J;

A _ij＝[a _ij(1)?a _ij(2)?…?a _ij(G)?a _ij(G+1)?…?a _ij(G+P-1)] ^T (6)

Step 5: in the receiving end signal processor, with sample formula (5) carry out between arteries and veins, twice decoding processing in the arteries and veins, can obtain to suppress the result after the non-phase parameter radio interference method is handled;

Being decoded as for the first time between the arteries and veins of formula (5) signal and formula (1) generation pseudorandom code book sequence carries out decoding processing and gets:

$\mathrm{RI}=\left[\begin{array}{c}\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{A}_{1j}{R}_j\\ \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{A}_{2j}{R}_{J+j}\\ .\\ .\\ .\\ \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{A}_{\mathrm{Nj}}{R}_{(N-1)J+j}\end{array}\right]={\left[\begin{array}{cccc}\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{1j}\left(1\right)R_j& \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{1j}\left(2\right)R_j& ...& \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{1j}(G+P-1)R_j\\ \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{2j}\left(1\right)R_{J+j}& \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{2j}\left(2\right)R_{J+j}& ...& \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{2j}(G+P-1)R_{J+j}\\ .& .& & .\\ .& .& ...& .\\ .& .& & .\\ \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{\mathrm{Nj}}\left(1\right)R_{(N-1)J+j}& \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{\mathrm{Nj}}\left(2\right)R_{(N-1)J+j}& ...& \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{\mathrm{Nj}}(G+P-1)R_{(N-1)J+j}\end{array}\right]}_{N\×G}---\left(7\right)$

The interior code book sequence of arteries and veins that the result of calculation that is decoded as for the second time formula (7) and formula (2) produce is carried out decoding processing and is got:

$\mathrm{BRI}{=\left[\begin{array}{cccc}\underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{1j}\left(p\right)R_j]& \underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{1j}(p+1)R_j]& ...& \underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{1j}(G+p-1)R_j]\\ \underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{2j}\left(p\right)R_{J+j}]& \underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{2j}(p+1)R_{J+j}]& ...& \underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{2j}(G+p-1)R_{J+j}]\\ .& .& & .\\ .& .& ...& .\\ .& .& & .\\ \underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{\mathrm{Nj}}\left(p\right)R_{(N-1)J+j}& \underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{\mathrm{Nj}}(p+1)R_{(N-1)J+j}]& ...& \underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{\mathrm{Nj}}(G+p-1)R_{(N-1)J+j}]\end{array}\right]}_{N\×G}---\left(8\right)$

So far, this wind profile radar non-phase parameter radio interference inhibition method disposes.

Beneficial effect:

Wind profile radar received signal through the method processing, pseudo-random code sequence carries out the computing of decoding the first time between process and the complete synchronous arteries and veins of emission modulating pulse, with the survey wind signal recovery of coherent and with incoherent unlimited electrical interference whiten, thereby reach the effect that reduces the radio interference power spectrum signal, the computing of carrying out decoding for the second time reaches the effect of the power that improves the wind signal according to pseudorandom code book sequence in the arteries and veins again, its interference suppressioning effect is (10log L) dB approximately, and wind signal power improves (20log P) dB approximately.

Description of drawings

Fig. 1 is a kind of wind profile radar non-phase parameter radio interference inhibition method treatment scheme of the present invention

Fig. 2 is that a kind of wind profile radar non-phase parameter radio interference inhibition method of the present invention is realized theory diagram.

Embodiment

Fig. 1 unites pseudorandomcode inhibition non-phase parameter radio interference method treatment scheme in the full coherent intercadence, Fig. 2 is that associating pseudorandomcode bi-phase modulated suppresses non-phase parameter radio interference realization theory diagram in the full coherent intercadence.Pseudo-random code generator in the wind profile radar monitoring combination (adopting the linear logic shift register to produce) output pseudo-random code sequence is given radio-frequency modulator, transmitter and signal processor then.Pseudo-random code is carried out the decoding processing first time between received signal and the complete synchronous arteries and veins of emission modulating pulse, the pseudorandom code book carries out the decoding processing second time in its output result and the arteries and veins, can obtain the result that the non-phase parameter radio wave interference suppresses after above twice decoding computing, this result can be used as the input signal that the wind profile radar follow-up signal is handled.

Illustrate, if the atmospheric signal frequency is 10Hz, radio interference is 30Hz, both power is identical, pseudo-random code code selection length is 255 8 level structure M sequences between arteries and veins, and the inter-stage feedback is connected to 8,6,5,4, and coding is selected 5 Barker codes in the arteries and veins, then suppress about 20dB through the unlimited electrical interference after the associating pseudorandomcode decoding processing in the full coherent intercadence, wind signal power improves about 13dB.

Claims (1)

1, a kind of wind profile radar non-phase parameter radio interference inhibition method is characterized in that these method concrete steps are:

Step 1: at the wind profile radar transmitting terminal, produce between arteries and veins synchronously, two kinds of pseudo-random code sequences in the arteries and veins; Wherein pseudo-random sequence adopts the M sequence between arteries and veins, and sequence length is selected according to the real work parameter, supposes that its code length is L, and coded sequence adopts Barker code in the arteries and veins, supposes that its code length is P,

If the code book of pseudorandomcode is R between arteries and veins:

R＝[R ₁?R ₂?…?R _i?…?R _L] ^T (1)

R in the formula _iBe M sequence pseudo-random code code element, i=1,2 ... L, T representing matrix transposition;

The coding code book is a Barker code in the arteries and veins, and its figure place is P,

BK＝[bk(1)?bk(2)?…?bk(i)?…?bk(P)] ^T (2)

Code element is bk (i) in the formula, i=1,2 ... P, T representing matrix transposition;

Step 2: at the wind profile radar transmitting terminal, above-mentioned formula (1) and formula (2) code book sequence through the inline yard sequence of compiling in collaboration with of intercadence are:

$\mathrm{RB}=\left[\begin{array}{c}{\mathrm{RB}}_1\\ {\mathrm{RB}}_2\\ .\\ .\\ .\\ {\mathrm{RB}}_i\\ .\\ .\\ .\\ {\mathrm{RB}}_L\end{array}\right]---\left(3\right)$

RB in the formula _iBe the code book signal in i the pulse in back of encoding, i=1,2 ... L,

${\mathrm{RB}}_i=(\mathrm{BK}+R_i)\mathrm{mod}\left(2\right)=\left[\begin{array}{c}\mathrm{bk}\left(1\right)\mathrm{xor}{R}_i\\ \mathrm{bk}\left(2\right)\mathrm{xor}{R}_i\\ .\\ .\\ .\\ \mathrm{bk}\left(P\right)\mathrm{xor}{R}_i\end{array}\right]---\left(4\right)$

Mod () is the complementation computing in the formula, and xor is an xor operator;

Step 3: the code book sequence that formula (3) is produced is sent to wind profile radar transmit modulator and emission, and the code book sequence with formula (1) and formula (2) generation is kept in the received signal processor simultaneously;

Step 4: transmitting in the step 3 received at the wind profile radar receiving end;

If time domain accumulation is that FFT counts and is N J time, the range gate number is G, and then the signal processor sample of signal space I that need handle is:

$I={\left[\begin{array}{cccc}{A}_{11}& A_{12}& ...& A_{1J}\\ A_{21}& A_{22}& ...& A_{2J}\\ .& .& & .\\ .& .& A_{\mathrm{ij}}& .\\ .& .& & .\\ A_{N1}& A_{N2}& ...& A_{\mathrm{NJ}}\end{array}\right]}_{N\×J}---\left(5\right)$

A in the formula _IjRepresent i * j sampled signal a respectively _Ij(1), a _Ij(2) ... a _Ij(G+P-1) sample space of Gou Chenging, i=1 wherein, 2 ... N, j=1,2 ... J respectively with signal Processing in FFT count relevant with the coherent accumulation number of times, L=N * J;

A _ij＝[a _ij(1)?a _ij(2)?…?a _ij(G)?a _ij(G+1)?…?a _ij(G+P-1)] ^T (6)

Step 5: in the receiving end signal processor, with sample formula (5) carry out between arteries and veins, twice decoding processing in the arteries and veins, can obtain to suppress the result after the non-phase parameter radio interference method is handled;

Being decoded as for the first time between the arteries and veins of formula (5) signal and formula (1) generation pseudorandom code book sequence carries out decoding processing and gets:

$\mathrm{RI}=\left[\begin{array}{c}\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{A}_{1j}{R}_j\\ \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{A}_{2j}{R}_{J+j}\\ .\\ .\\ .\\ \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{A}_{\mathrm{Nj}}{R}_{(N-1)J+j}\end{array}\right]={\left[\begin{array}{cccc}\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{1j}\left(1\right)R_j& \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{1j}\left(2\right)R_j& ...& \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{1j}(G+P-1)R_j\\ \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{2j}\left(1\right)R_{J+j}& \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{2j}\left(2\right)R_{J+j}& ...& \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{2j}(G+P-1)R_{J+j}\\ .& .& & .\\ .& .& ...& .\\ .& .& & .\\ \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{\mathrm{Nj}}\left(1\right)R_{(N-1)J+j}& \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{\mathrm{Nj}}\left(2\right)R_{(N-1)J+j}& ...& \underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{\mathrm{Nj}}(G+P-1)R_{(N-1)J+j}\end{array}\right]}_{N\×G}---\left(7\right)$

The interior code book sequence of arteries and veins that the result of calculation that is decoded as for the second time formula (7) and formula (2) produce is carried out decoding processing and is got:

$\mathrm{BRI}{=\left[\begin{array}{cccc}\underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{1j}\left(p\right)R_j]& \underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{1j}(p+1)R_j]& ...& \underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{1j}(G+p-1)R_j]\\ \underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{2j}\left(p\right)R_{J+j}]& \underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{2j}(p+1)R_{J+j}]& ...& \underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{2j}(G+p-1)R_{J+j}]\\ .& .& & .\\ .& .& ...& .\\ .& .& & .\\ \underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{\mathrm{Nj}}\left(p\right)R_{(N-1)J+j}& \underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{\mathrm{Nj}}(p+1)R_{(N-1)J+j}]& ...& \underset{p=1}{\overset{P}{\mathrm{\Σ}}}[\mathrm{bk}\left(p\right)\·\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{a}_{\mathrm{Nj}}(G+p-1)R_{(N-1)J+j}]\end{array}\right]}_{N\×G}---\left(8\right)$

So far, this wind profile radar non-phase parameter radio interference inhibition method disposes.

Images