[go: up one dir, main page]

CN114764136A - Radar anti-interference waveform generation method based on multi-time scale coupling network - Google Patents

Radar anti-interference waveform generation method based on multi-time scale coupling network Download PDF

Info

Publication number
CN114764136A
CN114764136A CN202110027839.4A CN202110027839A CN114764136A CN 114764136 A CN114764136 A CN 114764136A CN 202110027839 A CN202110027839 A CN 202110027839A CN 114764136 A CN114764136 A CN 114764136A
Authority
CN
China
Prior art keywords
signal
waveform
pulse
radar
broadband
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
Application number
CN202110027839.4A
Other languages
Chinese (zh)
Other versions
CN114764136B (en
Inventor
熊清
胡文
赵月
董浩
姚翼荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu Yunhefeng Intelligent Technology Co ltd
Original Assignee
Jiangsu Yunhefeng Intelligent Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jiangsu Yunhefeng Intelligent Technology Co ltd filed Critical Jiangsu Yunhefeng Intelligent Technology Co ltd
Priority to CN202110027839.4A priority Critical patent/CN114764136B/en
Publication of CN114764136A publication Critical patent/CN114764136A/en
Application granted granted Critical
Publication of CN114764136B publication Critical patent/CN114764136B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems
    • G01S7/282Transmitters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems
    • G01S7/285Receivers
    • G01S7/292Extracting wanted echo-signals

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention relates to the technical field of radar signal processing, and discloses a radar anti-interference waveform generation method in a multi-time scale coupling network. The invention aims to design an MIMO radar waveform with better orthogonality and anti-jamming capability, a reference signal of broadband matched filtering is generated by a multi-time scale coupling network, and a processing method of an echo is given; the transmit waveforms between each subarray of the MIMO radar are required to satisfy orthogonality, so that the transmit waveforms between each subarray can be distinguished at a receiving end through matched filtering.

Description

基于多时间尺度耦合网络的雷达抗干扰波形产生方法Radar anti-jamming waveform generation method based on multi-time scale coupling network

技术领域technical field

本发明涉及雷达信号处理技术领域,具体为基于多时间尺度耦合网络的雷达抗干扰波形产生方法。The invention relates to the technical field of radar signal processing, in particular to a method for generating radar anti-jamming waveforms based on a multi-time scale coupling network.

背景技术Background technique

波形设计是雷达研制过程中非常基础性的工作,既包括信号波形的选择和相关参数的优化,也包含匹配处理权值的优化。不同类型的波形设计方式存在较大的差距,但对于MIMO雷达,旁瓣抑制与波形的正交性始终是波形设计的焦点。Waveform design is a very basic work in the process of radar development, including the selection of signal waveforms and the optimization of related parameters, as well as the optimization of matching processing weights. There is a big gap between different types of waveform design methods, but for MIMO radar, sidelobe suppression and waveform orthogonality are always the focus of waveform design.

MIMO雷达的正交波形分为频分类正交波形和编码类正交波形,若只考虑使用频分类正交波形或者编码类正交波形,当MIMO雷达的子阵较多时,子阵之间的正交性会降低,影响MIMO雷达的若目标探测。同时,频分类正交波形若要获取较低的自相关旁瓣,需要的脉冲个数较多,因此相干处理周期比较长,不利于高速目标的测量。Orthogonal waveforms of MIMO radars are divided into frequency-classified orthogonal waveforms and coded orthogonal waveforms. If only the frequency-classified orthogonal waveforms or coded orthogonal waveforms are used, when there are many subarrays of the MIMO radar, the difference between the subarrays will be reduced. Orthogonality will be reduced, which affects the detection of MIMO radars. At the same time, in order to obtain lower autocorrelation side lobes in frequency-classified orthogonal waveforms, a large number of pulses are required, so the coherent processing period is relatively long, which is not conducive to the measurement of high-speed targets.

根据转发式干扰机的干扰原理,截获接收机需要准确提取雷达发射信号的频率、相位、幅度以及脉位信息才能够对雷达起到干扰作用。本发明的MIMO雷达发射波形的参数自由度高,信号空间辐射功率低,增加干扰机精确提取雷达信号频率、相位和脉位信息带来了较大的难度。本发明的波形设计联合了考虑了波形的正交性和抗干扰能力。在线性调频脉冲串的基础上,单脉冲内进行了编码长度为M的相位编码,脉冲间进行了频率和脉位的随机捷变,通过遗传算法对波形优化后,使波形的频率、相位和脉位具有随机性,在确保波形的正交性的同时使波形具有抗转发式干扰的能力。According to the jamming principle of the repeater jammer, the intercepting receiver needs to accurately extract the frequency, phase, amplitude and pulse position information of the radar transmitted signal to interfere with the radar. The MIMO radar transmit waveform of the present invention has high parameter freedom and low signal space radiation power, and increasing the jammer to accurately extract the radar signal frequency, phase and pulse position information brings great difficulty. The waveform design of the present invention combines the orthogonality of the waveform and the anti-interference ability. On the basis of the chirp pulse train, a phase encoding with a coding length of M is carried out within a single pulse, and random agility of frequency and pulse position is carried out between pulses. After the waveform is optimized by genetic algorithm, the frequency, phase and pulse position of the waveform are The pulse position is random, which makes the waveform have the ability to resist transponder interference while ensuring the orthogonality of the waveform.

在回波信号处理上,传统信号处理是将回波信号作为窄带信号,对同一时刻各子阵回波脉冲进行匹配滤波,再通过相干积累和恒虚警处理得到目标的距离速度信息。传统窄带信号处理方法对宽带调制信号的回波处理需要对脉冲间进行复杂的相位补偿。此外,当目标处于高速运动状态,在相干处理周期内位移超过一个距离分辨单元,则按传统窄带回波处理误差较大。针对窄带处理算法误差较大的问题,本发明采用基于多时间尺度的宽带信号模糊函数匹配算法,提高目标距离速度测量精度。In the echo signal processing, the traditional signal processing is to use the echo signal as a narrowband signal, perform matched filtering on the echo pulses of each subarray at the same time, and then obtain the distance and speed information of the target through coherent accumulation and constant false alarm processing. The echo processing of broadband modulated signals by traditional narrowband signal processing methods requires complex phase compensation between pulses. In addition, when the target is in a high-speed motion state and the displacement exceeds one range resolution unit in the coherent processing period, the traditional narrowband echo processing error is large. Aiming at the problem that the error of the narrowband processing algorithm is relatively large, the present invention adopts a broadband signal fuzzy function matching algorithm based on multiple time scales to improve the measurement accuracy of the target distance and speed.

发明内容SUMMARY OF THE INVENTION

本发明提供了基于多时间尺度耦合网络的雷达抗干扰波形产生方法,针对上述背景技术中的缺点,综合考虑波形的正交性和抗干扰性,给出了基于多时间尺度耦合网络的雷达抗干扰波形产生方法,本发明用遗传算法对波形的参数进行优化,得到了较好的正交性和抗干扰能力。The invention provides a method for generating radar anti-jamming waveforms based on a multi-time scale coupling network. In view of the shortcomings of the above-mentioned background technology, the orthogonality and anti-jamming properties of waveforms are comprehensively considered, and a radar anti-jamming method based on a multi-time scale coupling network is provided. In the method for generating interference waveform, the present invention uses genetic algorithm to optimize the parameters of the waveform, and obtains better orthogonality and anti-interference ability.

本发明提供如下技术方案:基于多时间尺度耦合网络的雷达抗干扰波形产生方法,其主要内容包括发射波形的产生与回波信号的处理,包括以下步骤:The present invention provides the following technical solutions: a method for generating radar anti-jamming waveforms based on a multi-time scale coupling network, the main contents of which include the generation of transmit waveforms and the processing of echo signals, including the following steps:

步骤一:构建MIMO雷达发射波形数学模型,给出阵列发射信号的时频关系图,得出波形设计需要优化的编码参数;Step 1: Build a mathematical model of the MIMO radar transmit waveform, give the time-frequency relationship diagram of the array transmit signal, and obtain the coding parameters that need to be optimized for waveform design;

步骤二:根据MIMO雷达波形设计正交性原则构建波形优化代价函数,用遗传算法对编码参数进行优化,构造MIMO雷达发射波形;Step 2: Construct the waveform optimization cost function according to the orthogonality principle of MIMO radar waveform design, optimize the coding parameters with genetic algorithm, and construct the MIMO radar transmit waveform;

步骤三:结合宽带信号匹配滤波原理给出本发明MIMO雷达的波束形成和目标距离速度信息提取处理过程;Step 3: Combining with the principle of broadband signal matching filtering, the beamforming and target distance and speed information extraction processing procedures of the MIMO radar of the present invention are given;

步骤四:根据宽带信号模糊函数理论,在产生MIMO雷达发射波形的同时产生一组匹配滤波基准信号,设计多时间尺度耦合网络硬件系统,在产生雷达发射信号的同时产生相应的一组匹配滤波基准信号。Step 4: According to the broadband signal ambiguity function theory, generate a set of matched filter reference signals while generating the MIMO radar transmit waveform, design a multi-time-scale coupling network hardware system, and generate a corresponding set of matched filter benchmarks while generating the radar transmit signal Signal.

优选的,所述MIMO雷达抗干扰复合调制信号,其主要建模过程如下:Preferably, the main modeling process of the MIMO radar anti-jamming composite modulation signal is as follows:

设MIMO雷达有L个子阵,一个CPI内的脉冲个数为N,每个脉冲内采用码长为M的相位编码,现给出第个子阵发射信号的数学模型:Assuming that the MIMO radar has L sub-arrays, the number of pulses in one CPI is N, and the phase encoding with code length M is used in each pulse. Now the mathematical model of the signal transmitted by the th sub-array is given:

Figure BDA0002890978560000031
Figure BDA0002890978560000031

上式中,Tr为子脉冲的平均脉冲周期,ΔTr为最小脉位跳变间隔,tp为脉冲宽度,ts为相位编码码元宽度,

Figure BDA0002890978560000032
为第n个子脉冲的载波,
Figure BDA0002890978560000033
μ为调频斜率,
Figure BDA0002890978560000034
为第n个子脉冲的第m个码元,从发射信号的数学模型可知,波形优化参数为
Figure BDA0002890978560000035
用向量表示记为
Figure BDA0002890978560000036
Ω=[ξ12,L,ξN]、
Figure BDA0002890978560000037
In the above formula, T r is the average pulse period of the sub-pulse, ΔT r is the minimum pulse jump interval, t p is the pulse width, ts is the phase encoding symbol width,
Figure BDA0002890978560000032
is the carrier of the nth sub-pulse,
Figure BDA0002890978560000033
μ is the frequency modulation slope,
Figure BDA0002890978560000034
is the m-th symbol of the n-th sub-pulse. From the mathematical model of the transmitted signal, the waveform optimization parameters are
Figure BDA0002890978560000035
represented by a vector as
Figure BDA0002890978560000036
Ω=[ξ 12 ,L,ξ N ],
Figure BDA0002890978560000037

波形参数中,每个脉冲内的相位编码位数为M位,每个码元为四相码;脉间频率跳变的最小范围为Δf=1/tp,确保脉冲间的正交性;每个脉冲的脉间间隔也随机的在一个范围内跳变,增加信号的抗干扰性能。In the waveform parameters, the number of phase coding bits in each pulse is M bits, and each symbol is a quaternary code; the minimum range of frequency jump between pulses is Δf=1/t p to ensure the orthogonality between pulses; The inter-pulse interval of each pulse also randomly jumps within a range, which increases the anti-interference performance of the signal.

优选的,所述MIMO雷达抗干扰波形优化方法,构造波形优化代价函数的方法具体为:Preferably, in the MIMO radar anti-jamming waveform optimization method, the method for constructing the waveform optimization cost function is specifically:

MIMO雷达的正交性体现在子阵间发射信号的自相关和互相关性能上,本发明中要求发射的信号具有较低的自相关旁瓣峰值和互相关峰值,因此采用峰值旁瓣电平准则(PSL)进行正交波形设计:The orthogonality of the MIMO radar is reflected in the auto-correlation and cross-correlation performance of the transmitted signals between the sub-arrays. In the present invention, the transmitted signals are required to have lower auto-correlation side lobe peaks and cross-correlation peaks, so the peak side lobe level is adopted. Criterion (PSL) for quadrature waveform design:

Figure BDA0002890978560000038
Figure BDA0002890978560000038

Figure BDA0002890978560000039
Figure BDA0002890978560000039

上式中,sp(t)和sq(t)为不同子阵的发射信号,当p=q时,cpq(τ)为发射信号的自相关函数;当p≠q时,cpq(τ)为发射信号的互相关函数,本发明以发射波形的相位编码矩阵、载波的频率编码序列以及脉位捷变序列作为优化变量,将发射信号的峰值旁瓣比作为代价函数,建立如下优化模型:In the above formula, s p (t) and s q (t) are the transmitted signals of different subarrays. When p=q, c pq (τ) is the autocorrelation function of the transmitted signal; when p≠q, c pq (τ) is the cross-correlation function of the transmission signal, the present invention takes the phase encoding matrix of the transmission waveform, the frequency encoding sequence of the carrier wave and the pulse position agility sequence as the optimization variable, and the peak sidelobe ratio of the transmission signal is used as the cost function, and is established as follows Optimizing the model:

Figure BDA0002890978560000041
Figure BDA0002890978560000041

该模型的优化参数为离散的相位、频率和脉位编码,可采用遗传算法对模型进行优化,得到满足条件的MIMO正交波形。The optimized parameters of the model are discrete phase, frequency and pulse position coding. Genetic algorithm can be used to optimize the model to obtain MIMO orthogonal waveforms that meet the conditions.

优选的,所述MIMO雷达宽带模糊函数波形分离与波束形成结果:Preferably, the MIMO radar broadband ambiguity function waveform separation and beamforming results:

考虑在收发共置MIMO雷达系统中,发射天线阵元和接收天线阵元个数分别为P和Q,收发均为均匀线阵,发射天线和接收天线的阵元间距分别为dt和dr,目标到达方向为

Figure BDA0002890978560000042
Considering that in the transceiver co-located MIMO radar system, the number of transmitting and receiving antenna elements is P and Q, respectively, the transmitting and receiving are uniform linear arrays, and the distances between transmitting and receiving antennas are d t and d r respectively. , the target arrival direction is
Figure BDA0002890978560000042

定义sk(t)表示第m个发射天线阵元的发射信号,s=[s1(t),s2(t),L,sP(t)]表示发射信号向量,M个发射信号经空间传播,到达位于斜距R(对应第r个距离单元)、方位向θ的目标的合成信号为:Definition sk (t) represents the transmit signal of the mth transmit antenna element, s=[s 1 (t), s 2 (t), L, s P (t)] represents the transmit signal vector, M transmit signals After space propagation, the synthetic signal reaching the target at the slant range R (corresponding to the rth distance unit) and the azimuth θ is:

Figure BDA0002890978560000043
Figure BDA0002890978560000043

式中ξ1为传播衰减因子,假设对各子阵发射信号相同,ω=(c-v)/(c+v)为尺度参数,τk为第k(k=1,2,L,M)个阵元发射的信号相对于参考阵元的时延差,可表示为:In the formula, ξ 1 is the propagation attenuation factor, assuming that the transmitted signals for each sub-array are the same, ω=(cv)/(c+v) is the scale parameter, and τ k is the kth (k=1, 2, L, M) The delay difference of the signal transmitted by the array element relative to the reference array element can be expressed as:

Figure BDA0002890978560000044
Figure BDA0002890978560000044

上式τk只与天线的布局有关,与目标的运动状态无关,固定时延τk可以用

Figure BDA0002890978560000045
来表示,因此目标处合成信号为:The above formula τ k is only related to the layout of the antenna and has nothing to do with the motion state of the target. The fixed time delay τ k can be used
Figure BDA0002890978560000045
to represent, so the synthesized signal at the target is:

Figure BDA0002890978560000051
Figure BDA0002890978560000051

上式中

Figure BDA0002890978560000052
为发射信号导向矢量,
Figure BDA0002890978560000053
为目标处接收信号向量,可表示为
Figure BDA0002890978560000054
其中τ0=2R0/(c-v);In the above formula
Figure BDA0002890978560000052
Steering vector for the transmitted signal,
Figure BDA0002890978560000053
is the received signal vector at the target, which can be expressed as
Figure BDA0002890978560000054
where τ 0 =2R 0 /(cv);

信号r(t)经一定RCS的目标反射,第q个接收天线接收到的信号为:The signal r(t) is reflected by the target of a certain RCS, and the signal received by the qth receiving antenna is:

Figure BDA0002890978560000055
Figure BDA0002890978560000055

式中,

Figure BDA0002890978560000056
为第q个天线相位差,令式中ξ1ξ2=ξ,将式(1)代入式(8),则第q个接收天线的回波为:In the formula,
Figure BDA0002890978560000056
is the phase difference of the qth antenna, let ξ 1 ξ 2 =ξ in the formula, and substitute formula (1) into formula (8), then the echo of the qth receiving antenna is:

Figure BDA0002890978560000057
Figure BDA0002890978560000057

接收信号向量为x(t)=[x1(t),x2(t),L,xQ(t)],式(9)的矩阵形式如下:The received signal vector is x(t)=[x 1 (t), x 2 (t), L, x Q (t)], the matrix form of equation (9) is as follows:

x(t)=ξb(θrTt)s[ω(t-τ0)](10)x(t)=ξb(θ rTt )s[ω(t-τ 0 )](10)

式(10)中的b(θr)为接收天线导向矢量,

Figure BDA0002890978560000058
现考虑第q个接收天线的接收回波的宽带匹配处理,由式(10)可得第q个接收通道回波信号为:b(θ r ) in equation (10) is the steering vector of the receiving antenna,
Figure BDA0002890978560000058
Now consider the broadband matching processing of the received echo of the qth receiving antenna, the echo signal of the qth receiving channel can be obtained from equation (10) as:

Figure BDA0002890978560000059
Figure BDA0002890978560000059

现对第q个接收天线回波信号进行宽带匹配滤波,由宽带模糊函数理论可知第p个发射通道的参考信号为:Now perform broadband matched filtering on the echo signal of the qth receiving antenna. According to the theory of broadband ambiguity function, it can be known that the reference signal of the pth transmit channel is:

Figure BDA00028909785600000510
Figure BDA00028909785600000510

第q个接收通道接收到第p个发射通道信号的匹配输出为:The matching output of the qth receiving channel receiving the signal of the pth transmitting channel is:

Figure BDA0002890978560000061
Figure BDA0002890978560000061

上式T为脉冲串持续时间,由发射信号的正交性可知,发射信号矢量中只有第p个发射通道的信号才有输出,匹配输出可简化为:The above formula T is the duration of the burst. According to the orthogonality of the transmitted signal, only the signal of the pth transmit channel in the transmitted signal vector has an output, and the matching output can be simplified as:

Figure BDA0002890978560000062
Figure BDA0002890978560000062

由上式中,当ω=ηm,τ0=ξm时,匹配滤波器有最大的输出,匹配滤波的最大输出记为

Figure BDA0002890978560000063
From the above formula, when ω=η m , τ 0m , the matched filter has the maximum output, and the maximum output of the matched filter is recorded as
Figure BDA0002890978560000063

回波信号通过宽带匹配网络后可以得到目标的距离速度信息,再通过波束形成可得到目标的方位角,波束形成的输出为:After the echo signal passes through the broadband matching network, the distance and speed information of the target can be obtained, and then the azimuth angle of the target can be obtained through beamforming. The output of the beamforming is:

Figure BDA0002890978560000064
Figure BDA0002890978560000064

式(15)中,

Figure BDA0002890978560000065
为宽带匹配滤波的输出,τ0与ω为目标的时延和多普勒频率,当
Figure BDA0002890978560000066
时,上式有最大值,此时目标方向为:In formula (15),
Figure BDA0002890978560000065
is the output of the broadband matched filter, τ 0 and ω are the target delay and Doppler frequency, when
Figure BDA0002890978560000066
When , the above formula has a maximum value, and the target direction is:

Figure BDA0002890978560000067
Figure BDA0002890978560000067

优选的,所述产生雷达发射信号的同时产生宽带回波匹配信号,正交波形的输出是由DDS波形发生器根据频率、相位编码序列产生每个

Figure BDA0002890978560000068
对应的相位处的波形的幅度值,这些幅度值产生了一系列的子脉冲,每个
Figure BDA0002890978560000069
的持续时间是由雷达的带宽所决定的;脉冲时序控制模块改变脉冲之间的脉冲周期,实现波形脉位捷变调制;通过多路DDS电路后可以产生MIMO雷达正交波形,产生MIMO雷达正交波形后,一路信号给MIMO雷达发射天线,一路作为多时间尺度耦合网络电路的输入,产生回波信号的宽带匹配滤波的参考信号;Preferably, the broadband echo matching signal is generated while the radar transmit signal is generated, and the output of the quadrature waveform is generated by the DDS waveform generator according to the frequency and phase coding sequence.
Figure BDA0002890978560000068
The amplitude values of the waveform at the corresponding phases that produce a series of sub-pulses, each
Figure BDA0002890978560000069
The duration is determined by the bandwidth of the radar; the pulse timing control module changes the pulse period between the pulses to realize the waveform pulse position agile modulation; after multiple DDS circuits, the MIMO radar quadrature waveform can be generated, and the MIMO radar positive waveform can be generated. After the waveform is crossed, one signal is sent to the MIMO radar transmitting antenna, and the other is used as the input of the multi-time scale coupling network circuit to generate a broadband matched filter reference signal of the echo signal;

宽带匹配网络实质是寻找回波信号二维模糊函数的最大值,在发射端产生发射信号的同时也产生了宽带匹配网络的参考信号,多时间尺度耦合网络的目的是产生一组每个子阵发射信号在时间和频率上联合移位的信号组,用于回波信号的宽带匹配滤波;在多时间尺度耦合电路部分,设每个子阵的发射信号延时τk为第k个距离门的接收信号,然后在第k个距离门内对所有多普勒频率进行匹配滤波,实现对第k个距离门回波信号的宽带脉压处理。The essence of the broadband matching network is to find the maximum value of the two-dimensional ambiguity function of the echo signal. When the transmitting signal is generated at the transmitting end, the reference signal of the broadband matching network is also generated. The purpose of the multi-time scale coupling network is to generate a set of transmissions for each sub-array. The signal group in which the signal is jointly shifted in time and frequency is used for the broadband matched filtering of the echo signal; in the multi-time scale coupling circuit part, the transmission signal delay τ k of each sub-array is set as the reception of the kth distance gate Then, all Doppler frequencies are matched and filtered in the kth range gate to realize the broadband pulse pressure processing of the echo signal of the kth range gate.

本发明具备以下有益效果:The present invention has the following beneficial effects:

本发明目的在于设计一种具有较好正交性、抗干扰能力的MIMO雷达波形,用多时间尺度耦合网络产生宽带匹配滤波的参考信号,给出回波的处理方法;MIMO雷达各子阵之间的发射波形要满足正交性,以便在接收端能够通过匹配滤波区分各子阵之间的发射波形。The purpose of the invention is to design a MIMO radar waveform with good orthogonality and anti-interference ability, use a multi-time scale coupling network to generate a reference signal of broadband matched filtering, and provide a processing method for the echo; The transmit waveforms between the sub-arrays should satisfy the orthogonality, so that the transmit waveforms between the sub-arrays can be distinguished by matched filtering at the receiving end.

本发明的波形在线性调频的基础上,对脉冲串进行脉内相位编码,脉间频率、脉位联合捷变,使波形具有较好的正交性与抗干扰性能,回波信号采用宽带信号匹配滤波理论,对回波信号进行匹配滤波区分各个发射通道的回波信号,最后结合数字波束形成得到目标的距离、速度和角度信息。On the basis of linear frequency modulation, the waveform of the invention performs intrapulse phase encoding on the pulse train, and the inter-pulse frequency and pulse position are jointly agile, so that the waveform has better orthogonality and anti-interference performance, and the echo signal adopts a broadband signal. According to the matched filtering theory, matched filtering is performed on the echo signal to distinguish the echo signals of each transmitting channel, and finally the distance, speed and angle information of the target are obtained by combining with digital beamforming.

本发明的特点在于发射机产生复合调制发射信号的同时,通过多时间尺度耦合网络产生了宽带信号匹配滤波的参考信号组,可对每一个接收通道的信号进行宽带匹配滤波得到目标的速度和距离,与传统回波信号处理方法相比,大大降低回波信号的处理复杂度,并且目标检测精度较高。The invention is characterized in that when the transmitter generates a composite modulated transmission signal, a reference signal group of broadband signal matching filtering is generated through a multi-time scale coupling network, and the signal of each receiving channel can be subjected to broadband matching filtering to obtain the speed and distance of the target. , compared with the traditional echo signal processing method, the processing complexity of the echo signal is greatly reduced, and the target detection accuracy is high.

附图说明Description of drawings

图1为本发明第l个子阵发射信号示意图;Fig. 1 is the schematic diagram of the lth subarray transmission signal of the present invention;

图2为本发明MIMO雷达发射波形时频关系图;Fig. 2 is the time-frequency relation diagram of the MIMO radar transmit waveform of the present invention;

图3为本发明回波信号处理流程图;Fig. 3 is the echo signal processing flow chart of the present invention;

图4为本发明多时间尺度耦合网络波形产生系统图。FIG. 4 is a diagram of a multi-time scale coupling network waveform generation system diagram of the present invention.

具体实施方式Detailed ways

下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

请参阅图1,基于多时间尺度耦合网络的雷达抗干扰波形产生方法,基于多时间尺度耦合网络的雷达抗干扰波形产生方法,其主要内容包括发射波形的产生与回波信号的处理,包括以下步骤:Please refer to Figure 1. The generation method of radar anti-jamming waveform based on multi-time scale coupling network, the main content of radar anti-jamming waveform generation method based on multi-time scale coupling network, including the generation of transmit waveform and the processing of echo signal, including the following step:

步骤一:构建MIMO雷达发射波形数学模型,给出阵列发射信号的时频关系图,得出波形设计需要优化的编码参数;Step 1: Build a mathematical model of the MIMO radar transmit waveform, give the time-frequency relationship diagram of the array transmit signal, and obtain the coding parameters that need to be optimized for waveform design;

步骤二:根据MIMO雷达波形设计正交性原则构建波形优化代价函数,用遗传算法对编码参数进行优化,构造MIMO雷达发射波形;Step 2: Construct the waveform optimization cost function according to the orthogonality principle of MIMO radar waveform design, optimize the coding parameters with genetic algorithm, and construct the MIMO radar transmit waveform;

步骤三:结合宽带信号匹配滤波原理给出本发明MIMO雷达的波束形成和目标距离速度信息提取处理过程;Step 3: Combining with the principle of broadband signal matching filtering, the beamforming and target distance and speed information extraction processing procedures of the MIMO radar of the present invention are given;

步骤四:根据宽带信号模糊函数理论,在产生MIMO雷达发射波形的同时产生一组匹配滤波基准信号,设计多时间尺度耦合网络硬件系统,在产生雷达发射信号的同时产生相应的一组匹配滤波基准信号。Step 4: According to the broadband signal ambiguity function theory, generate a set of matched filter reference signals while generating the MIMO radar transmit waveform, design a multi-time-scale coupling network hardware system, and generate a corresponding set of matched filter benchmarks while generating the radar transmit signal Signal.

MIMO雷达抗干扰复合调制信号,其主要建模过程如下:The main modeling process of MIMO radar anti-jamming composite modulation signal is as follows:

设MIMO雷达有L个子阵,一个CPI内的脉冲个数为N,每个脉冲内采用码长为M的相位编码,现给出第个子阵发射信号的数学模型:Assuming that the MIMO radar has L sub-arrays, the number of pulses in one CPI is N, and the phase encoding with code length M is used in each pulse. Now the mathematical model of the signal transmitted by the th sub-array is given:

Figure BDA0002890978560000081
Figure BDA0002890978560000081

上式中,Tr为子脉冲的平均脉冲周期,ΔTr为最小脉位跳变间隔,tp为脉冲宽度,ts为相位编码码元宽度,

Figure BDA0002890978560000091
为第n个子脉冲的载波,
Figure BDA0002890978560000092
μ为调频斜率,
Figure BDA0002890978560000093
为第n个子脉冲的第m个码元,从发射信号的数学模型可知,波形优化参数为
Figure BDA0002890978560000094
用向量表示记为
Figure BDA0002890978560000095
Ω=[ξ12,L,ξN]、
Figure BDA0002890978560000096
In the above formula, T r is the average pulse period of the sub-pulse, ΔT r is the minimum pulse jump interval, t p is the pulse width, ts is the phase encoding symbol width,
Figure BDA0002890978560000091
is the carrier of the nth sub-pulse,
Figure BDA0002890978560000092
μ is the frequency modulation slope,
Figure BDA0002890978560000093
is the m-th symbol of the n-th sub-pulse. From the mathematical model of the transmitted signal, the waveform optimization parameters are
Figure BDA0002890978560000094
represented by a vector as
Figure BDA0002890978560000095
Ω=[ξ 12 ,L,ξ N ],
Figure BDA0002890978560000096

波形参数中,每个脉冲内的相位编码位数为M位,每个码元为四相码;脉间频率跳变的最小范围为Δf=1/tp,确保脉冲间的正交性;每个脉冲的脉间间隔也随机的在一个范围内跳变,增加信号的抗干扰性能。In the waveform parameters, the number of phase coding bits in each pulse is M bits, and each symbol is a quaternary code; the minimum range of frequency jump between pulses is Δf=1/t p to ensure the orthogonality between pulses; The inter-pulse interval of each pulse also randomly jumps within a range, which increases the anti-interference performance of the signal.

MIMO雷达抗干扰波形优化方法,构造波形优化代价函数的方法具体为:The waveform optimization method of MIMO radar anti-jamming, the method of constructing the waveform optimization cost function is as follows:

MIMO雷达的正交性体现在子阵间发射信号的自相关和互相关性能上,本发明中要求发射的信号具有较低的自相关旁瓣峰值和互相关峰值,因此采用峰值旁瓣电平准则(PSL)进行正交波形设计:The orthogonality of the MIMO radar is reflected in the auto-correlation and cross-correlation performance of the transmitted signals between the sub-arrays. In the present invention, the transmitted signals are required to have lower auto-correlation side lobe peaks and cross-correlation peaks, so the peak side lobe level is adopted. Criterion (PSL) for quadrature waveform design:

Figure BDA0002890978560000097
Figure BDA0002890978560000097

Figure BDA0002890978560000098
Figure BDA0002890978560000098

上式中,sp(t)和sq(t)为不同子阵的发射信号,当p=q时,cpq(τ)为发射信号的自相关函数;当p≠q时,cpq(τ)为发射信号的互相关函数,本发明以发射波形的相位编码矩阵、载波的频率编码序列以及脉位捷变序列作为优化变量,将发射信号的峰值旁瓣比作为代价函数,建立如下优化模型:In the above formula, s p (t) and s q (t) are the transmitted signals of different subarrays. When p=q, c pq (τ) is the autocorrelation function of the transmitted signal; when p≠q, c pq (τ) is the cross-correlation function of the transmission signal, the present invention takes the phase encoding matrix of the transmission waveform, the frequency encoding sequence of the carrier wave and the pulse position agility sequence as the optimization variable, and the peak sidelobe ratio of the transmission signal is used as the cost function, and is established as follows Optimizing the model:

Figure BDA0002890978560000101
Figure BDA0002890978560000101

该模型的优化参数为离散的相位、频率和脉位编码,可采用遗传算法对模型进行优化,得到满足条件的MIMO正交波形。The optimized parameters of the model are discrete phase, frequency and pulse position coding. Genetic algorithm can be used to optimize the model to obtain MIMO orthogonal waveforms that meet the conditions.

所述MIMO雷达宽带模糊函数波形分离与波束形成结果:The MIMO radar broadband ambiguity function waveform separation and beamforming results:

考虑在收发共置MIMO雷达系统中,发射天线阵元和接收天线阵元个数分别为P和Q,收发均为均匀线阵,发射天线和接收天线的阵元间距分别为dt和dr,目标到达方向为

Figure BDA0002890978560000102
Considering that in the transceiver co-located MIMO radar system, the number of transmitting and receiving antenna elements is P and Q, respectively, the transmitting and receiving are uniform linear arrays, and the distances between transmitting and receiving antennas are d t and d r respectively. , the target arrival direction is
Figure BDA0002890978560000102

定义sk(t)表示第m个发射天线阵元的发射信号,s=[s1(t),s2(t),L,sP(t)]表示发射信号向量,M个发射信号经空间传播,到达位于斜距R(对应第r个距离单元)、方位向θ的目标的合成信号为:Definition sk (t) represents the transmit signal of the mth transmit antenna element, s=[s 1 (t), s 2 (t), L, s P (t)] represents the transmit signal vector, M transmit signals After space propagation, the synthetic signal reaching the target at the slant range R (corresponding to the rth distance unit) and the azimuth θ is:

Figure BDA0002890978560000103
Figure BDA0002890978560000103

式中ξ1为传播衰减因子,假设对各子阵发射信号相同,ω=(c-v)/(c+v)为尺度参数,τk为第k(k=1,2,L,M)个阵元发射的信号相对于参考阵元的时延差,可表示为:In the formula, ξ 1 is the propagation attenuation factor, assuming that the transmitted signals for each sub-array are the same, ω=(cv)/(c+v) is the scale parameter, and τ k is the kth (k=1, 2, L, M) The delay difference of the signal transmitted by the array element relative to the reference array element can be expressed as:

Figure BDA0002890978560000104
Figure BDA0002890978560000104

上式τk只与天线的布局有关,与目标的运动状态无关,固定时延τk可以用

Figure BDA0002890978560000105
来表示,因此目标处合成信号为:The above formula τ k is only related to the layout of the antenna and has nothing to do with the motion state of the target. The fixed time delay τ k can be used
Figure BDA0002890978560000105
to represent, so the synthesized signal at the target is:

Figure BDA0002890978560000106
Figure BDA0002890978560000106

上式中

Figure BDA0002890978560000111
为发射信号导向矢量,
Figure BDA0002890978560000112
为目标处接收信号向量,可表示为
Figure BDA0002890978560000113
其中τ0=2R0/(c-v);In the above formula
Figure BDA0002890978560000111
Steering vector for the transmitted signal,
Figure BDA0002890978560000112
is the received signal vector at the target, which can be expressed as
Figure BDA0002890978560000113
where τ 0 =2R 0 /(cv);

信号r(t)经一定RCS的目标反射,第q个接收天线接收到的信号为:The signal r(t) is reflected by the target of a certain RCS, and the signal received by the qth receiving antenna is:

Figure BDA0002890978560000114
Figure BDA0002890978560000114

式中,

Figure BDA0002890978560000115
为第q个天线相位差,令式中ξ1ξ2=ξ,将式(1)代入式(8),则第q个接收天线的回波为:In the formula,
Figure BDA0002890978560000115
is the phase difference of the qth antenna, let ξ 1 ξ 2 =ξ in the formula, and substitute formula (1) into formula (8), then the echo of the qth receiving antenna is:

Figure BDA0002890978560000116
Figure BDA0002890978560000116

接收信号向量为x(t)=[x1(t),x2(t),L,xQ(t)],式(9)的矩阵形式如下:The received signal vector is x(t)=[x 1 (t), x 2 (t), L, x Q (t)], the matrix form of equation (9) is as follows:

x(t)=ξb(θrTt)s[ω(t-τ0)] (10)x(t)=ξb(θ rTt )s[ω(t-τ 0 )] (10)

式(10)中的b(θr)为接收天线导向矢量,

Figure BDA0002890978560000117
现考虑第q个接收天线的接收回波的宽带匹配处理,由式(10)可得第q个接收通道回波信号为:b(θ r ) in equation (10) is the steering vector of the receiving antenna,
Figure BDA0002890978560000117
Now consider the broadband matching processing of the received echo of the qth receiving antenna, the echo signal of the qth receiving channel can be obtained from equation (10) as:

Figure BDA0002890978560000118
Figure BDA0002890978560000118

现对第q个接收天线回波信号进行宽带匹配滤波,由宽带模糊函数理论可知第p个发射通道的参考信号为:Now perform broadband matched filtering on the echo signal of the qth receiving antenna. According to the theory of broadband ambiguity function, it can be known that the reference signal of the pth transmit channel is:

Figure BDA0002890978560000119
Figure BDA0002890978560000119

第q个接收通道接收到第p个发射通道信号的匹配输出为:The matching output of the qth receiving channel receiving the signal of the pth transmitting channel is:

Figure BDA00028909785600001110
Figure BDA00028909785600001110

上式T为脉冲串持续时间,由发射信号的正交性可知,发射信号矢量中只有第p个发射通道的信号才有输出,匹配输出可简化为:The above formula T is the duration of the burst. According to the orthogonality of the transmitted signal, only the signal of the pth transmit channel in the transmitted signal vector has an output, and the matching output can be simplified as:

Figure BDA0002890978560000121
Figure BDA0002890978560000121

由上式中,当ω=ηm,τ0=ξm时,匹配滤波器有最大的输出,匹配滤波的最大输出记为

Figure BDA0002890978560000122
From the above formula, when ω=η m , τ 0m , the matched filter has the maximum output, and the maximum output of the matched filter is recorded as
Figure BDA0002890978560000122

回波信号通过宽带匹配网络后可以得到目标的距离速度信息,再通过波束形成可得到目标的方位角,波束形成的输出为:After the echo signal passes through the broadband matching network, the distance and speed information of the target can be obtained, and then the azimuth angle of the target can be obtained through beamforming. The output of the beamforming is:

Figure BDA0002890978560000123
Figure BDA0002890978560000123

式(15)中,

Figure BDA0002890978560000124
为宽带匹配滤波的输出,τ0与ω为目标的时延和多普勒频率,当
Figure BDA0002890978560000125
时,上式有最大值,此时目标方向为:In formula (15),
Figure BDA0002890978560000124
is the output of the broadband matched filter, τ 0 and ω are the target delay and Doppler frequency, when
Figure BDA0002890978560000125
When , the above formula has a maximum value, and the target direction is:

Figure BDA0002890978560000126
Figure BDA0002890978560000126

产生雷达发射信号的同时产生宽带回波匹配信号,正交波形的输出是由DDS波形发生器根据频率、相位编码序列产生每个

Figure BDA0002890978560000127
对应的相位处的波形的幅度值,这些幅度值产生了一系列的子脉冲,每个
Figure BDA0002890978560000128
的持续时间是由雷达的带宽所决定的;脉冲时序控制模块改变脉冲之间的脉冲周期,实现波形脉位捷变调制;通过多路DDS电路后可以产生MIMO雷达正交波形,产生MIMO雷达正交波形后,一路信号给MIMO雷达发射天线,一路作为多时间尺度耦合网络电路的输入,产生回波信号的宽带匹配滤波的参考信号;While generating the radar transmission signal, a broadband echo matching signal is generated. The output of the quadrature waveform is generated by the DDS waveform generator according to the frequency and phase coding sequence.
Figure BDA0002890978560000127
The amplitude values of the waveform at the corresponding phases that produce a series of sub-pulses, each
Figure BDA0002890978560000128
The duration is determined by the bandwidth of the radar; the pulse timing control module changes the pulse period between the pulses to realize the waveform pulse position agile modulation; after multiple DDS circuits, the MIMO radar quadrature waveform can be generated, and the MIMO radar positive waveform can be generated. After the waveform is crossed, one signal is sent to the MIMO radar transmitting antenna, and the other is used as the input of the multi-time scale coupling network circuit to generate a broadband matched filter reference signal of the echo signal;

宽带匹配网络实质是寻找回波信号二维模糊函数的最大值,在发射端产生发射信号的同时也产生了宽带匹配网络的参考信号,多时间尺度耦合网络的目的是产生一组每个子阵发射信号在时间和频率上联合移位的信号组,用于回波信号的宽带匹配滤波;在多时间尺度耦合电路部分,设每个子阵的发射信号延时τk为第k个距离门的接收信号,然后在第k个距离门内对所有多普勒频率进行匹配滤波,实现对第k个距离门回波信号的宽带脉压处理。The essence of the broadband matching network is to find the maximum value of the two-dimensional ambiguity function of the echo signal. When the transmitting signal is generated at the transmitting end, the reference signal of the broadband matching network is also generated. The purpose of the multi-time scale coupling network is to generate a set of transmissions for each sub-array. The signal group in which the signal is jointly shifted in time and frequency is used for the broadband matched filtering of the echo signal; in the multi-time scale coupling circuit part, the transmission signal delay τ k of each sub-array is set as the reception of the kth distance gate Then, all Doppler frequencies are matched and filtered in the kth range gate to realize the broadband pulse pressure processing of the echo signal of the kth range gate.

需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。It should be noted that, in this document, relational terms such as first and second are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any relationship between these entities or operations. any such actual relationship or sequence exists. Moreover, the terms "comprising", "comprising" or any other variation thereof are intended to encompass a non-exclusive inclusion such that a process, method, article or device that includes a list of elements includes not only those elements, but also includes not explicitly listed or other elements inherent to such a process, method, article or apparatus.

尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。Although embodiments of the present invention have been shown and described, it will be understood by those skilled in the art that various changes, modifications, and substitutions can be made in these embodiments without departing from the principle and spirit of the invention and modifications, the scope of the present invention is defined by the appended claims and their equivalents.

Claims (5)

1.基于多时间尺度耦合网络的雷达抗干扰波形产生方法,其主要内容包括发射波形的产生与回波信号的处理,其特征在于:包括以下步骤:1. based on the radar anti-jamming waveform generation method of multi-time scale coupling network, its main content comprises the generation of transmission waveform and the processing of echo signal, it is characterized in that: comprise the following steps: 步骤一:构建MIMO雷达发射波形数学模型,给出阵列发射信号的时频关系图,得出波形设计需要优化的编码参数;Step 1: Build a mathematical model of the MIMO radar transmit waveform, give the time-frequency relationship diagram of the array transmit signal, and obtain the coding parameters that need to be optimized for waveform design; 步骤二:根据MIMO雷达波形设计正交性原则构建波形优化代价函数,用遗传算法对编码参数进行优化,构造MIMO雷达发射波形;Step 2: Construct the waveform optimization cost function according to the orthogonality principle of MIMO radar waveform design, optimize the coding parameters with genetic algorithm, and construct the MIMO radar transmit waveform; 步骤三:结合宽带信号匹配滤波原理给出本发明MIMO雷达的波束形成和目标距离速度信息提取处理过程;Step 3: Combining with the principle of broadband signal matching filtering, the beamforming and target distance and speed information extraction processing procedures of the MIMO radar of the present invention are given; 步骤四:根据宽带信号模糊函数理论,在产生MIMO雷达发射波形的同时产生一组匹配滤波基准信号,设计多时间尺度耦合网络硬件系统,在产生雷达发射信号的同时产生相应的一组匹配滤波基准信号。Step 4: According to the broadband signal ambiguity function theory, generate a set of matched filter reference signals while generating the MIMO radar transmit waveform, design a multi-time-scale coupling network hardware system, and generate a corresponding set of matched filter benchmarks while generating the radar transmit signal Signal. 2.根据权利要求1所述的基于多时间尺度耦合网络的雷达抗干扰波形产生方法,其特征在于:所述MIMO雷达抗干扰复合调制信号,其主要建模过程如下:2. the radar anti-jamming waveform generation method based on multi-time scale coupling network according to claim 1, is characterized in that: described MIMO radar anti-jamming composite modulation signal, its main modeling process is as follows: 设MIMO雷达有L个子阵,一个CPI内的脉冲个数为N,每个脉冲内采用码长为M的相位编码,现给出第个子阵发射信号的数学模型:Assuming that the MIMO radar has L sub-arrays, the number of pulses in one CPI is N, and the phase encoding with code length M is used in each pulse. Now the mathematical model of the signal transmitted by the th sub-array is given:
Figure FDA0002890978550000011
Figure FDA0002890978550000011
上式中,Tr为子脉冲的平均脉冲周期,ΔTr为最小脉位跳变间隔,tp为脉冲宽度,ts为相位编码码元宽度,
Figure FDA0002890978550000012
为第n个子脉冲的载波,
Figure FDA0002890978550000013
μ为调频斜率,
Figure FDA0002890978550000014
为第n个子脉冲的第m个码元,从发射信号的数学模型可知,波形优化参数为
Figure FDA0002890978550000015
用向量表示记为
Figure FDA0002890978550000016
Ω=[ξ12,L,ξN]、
Figure FDA0002890978550000017
In the above formula, T r is the average pulse period of the sub-pulse, ΔT r is the minimum pulse jump interval, t p is the pulse width, ts is the phase encoding symbol width,
Figure FDA0002890978550000012
is the carrier of the nth sub-pulse,
Figure FDA0002890978550000013
μ is the frequency modulation slope,
Figure FDA0002890978550000014
is the m-th symbol of the n-th sub-pulse. From the mathematical model of the transmitted signal, the waveform optimization parameters are
Figure FDA0002890978550000015
represented by a vector as
Figure FDA0002890978550000016
Ω=[ξ 12 ,L,ξ N ],
Figure FDA0002890978550000017
波形参数中,每个脉冲内的相位编码位数为M位,每个码元为四相码;脉间频率跳变的最小范围为Δf=1/tp,确保脉冲间的正交性;每个脉冲的脉间间隔也随机的在一个范围内跳变,增加信号的抗干扰性能。In the waveform parameters, the number of phase coding bits in each pulse is M bits, and each symbol is a quaternary code; the minimum range of frequency jump between pulses is Δf=1/t p to ensure the orthogonality between pulses; The inter-pulse interval of each pulse also randomly jumps within a range, which increases the anti-interference performance of the signal.
3.根据权利要求2所述的基于多时间尺度耦合网络的雷达抗干扰波形产生方法,其特征在于:所述MIMO雷达抗干扰波形优化方法,构造波形优化代价函数的方法具体为:3. the radar anti-jamming waveform generation method based on multi-time scale coupling network according to claim 2, is characterized in that: described MIMO radar anti-jamming waveform optimization method, the method for constructing waveform optimization cost function is specifically: MIMO雷达的正交性体现在子阵间发射信号的自相关和互相关性能上,本发明中要求发射的信号具有较低的自相关旁瓣峰值和互相关峰值,因此采用峰值旁瓣电平准则(PSL)进行正交波形设计:The orthogonality of the MIMO radar is reflected in the auto-correlation and cross-correlation performance of the transmitted signals between the sub-arrays. In the present invention, the transmitted signals are required to have lower auto-correlation side lobe peaks and cross-correlation peaks, so the peak side lobe level is adopted. Criterion (PSL) for quadrature waveform design:
Figure FDA0002890978550000021
Figure FDA0002890978550000021
Figure FDA0002890978550000022
Figure FDA0002890978550000022
上式中,sp(t)和sq(t)为不同子阵的发射信号,当p=q时,cpq(τ)为发射信号的自相关函数;当p≠q时,cpq(τ)为发射信号的互相关函数,本发明以发射波形的相位编码矩阵、载波的频率编码序列以及脉位捷变序列作为优化变量,将发射信号的峰值旁瓣比作为代价函数,建立如下优化模型:In the above formula, s p (t) and s q (t) are the transmitted signals of different subarrays. When p=q, c pq (τ) is the autocorrelation function of the transmitted signal; when p≠q, c pq (τ) is the cross-correlation function of the transmission signal, the present invention takes the phase encoding matrix of the transmission waveform, the frequency encoding sequence of the carrier wave and the pulse position agility sequence as the optimization variable, and the peak sidelobe ratio of the transmission signal is used as the cost function, and is established as follows Optimizing the model:
Figure FDA0002890978550000023
Figure FDA0002890978550000023
该模型的优化参数为离散的相位、频率和脉位编码,可采用遗传算法对模型进行优化,得到满足条件的MIMO正交波形。The optimized parameters of the model are discrete phase, frequency and pulse position coding. Genetic algorithm can be used to optimize the model to obtain MIMO orthogonal waveforms that meet the conditions.
4.根据权利要求3所述的基于多时间尺度耦合网络的雷达抗干扰波形产生方法,其特征在于:所述MIMO雷达宽带模糊函数波形分离与波束形成结果:4. the radar anti-jamming waveform generation method based on multi-time scale coupling network according to claim 3, is characterized in that: described MIMO radar wideband ambiguity function waveform separation and beamforming result: 考虑在收发共置MIMO雷达系统中,发射天线阵元和接收天线阵元个数分别为P和Q,收发均为均匀线阵,发射天线和接收天线的阵元间距分别为dt和dr,目标到达方向为
Figure FDA0002890978550000031
Considering that in the transceiver co-located MIMO radar system, the number of transmitting and receiving antenna elements is P and Q, respectively, the transmitting and receiving are uniform linear arrays, and the distances between transmitting and receiving antennas are d t and d r respectively. , the target arrival direction is
Figure FDA0002890978550000031
定义sk(t)表示第m个发射天线阵元的发射信号,s=[s1(t),s2(t),L,sP(t)]表示发射信号向量,M个发射信号经空间传播,到达位于斜距R(对应第r个距离单元)、方位向θ的目标的合成信号为:Definition sk (t) represents the transmit signal of the mth transmit antenna element, s=[s 1 (t), s 2 (t), L, s P (t)] represents the transmit signal vector, M transmit signals After space propagation, the synthetic signal reaching the target at the slant range R (corresponding to the rth distance unit) and the azimuth θ is:
Figure FDA0002890978550000032
Figure FDA0002890978550000032
式中ξ1为传播衰减因子,假设对各子阵发射信号相同,ω=(c-v)/(c+v)为尺度参数,τk为第k(k=1,2,L,M)个阵元发射的信号相对于参考阵元的时延差,可表示为:In the formula, ξ 1 is the propagation attenuation factor, assuming that the transmitted signals for each sub-array are the same, ω=(cv)/(c+v) is the scale parameter, and τ k is the kth (k=1, 2, L, M) The delay difference of the signal transmitted by the array element relative to the reference array element can be expressed as:
Figure FDA0002890978550000033
Figure FDA0002890978550000033
上式τk只与天线的布局有关,与目标的运动状态无关,固定时延τk可以用
Figure FDA0002890978550000034
来表示,因此目标处合成信号为:
The above formula τ k is only related to the layout of the antenna and has nothing to do with the motion state of the target. The fixed time delay τ k can be used
Figure FDA0002890978550000034
to represent, so the synthesized signal at the target is:
Figure FDA0002890978550000035
Figure FDA0002890978550000035
上式中
Figure FDA0002890978550000036
为发射信号导向矢量,
Figure FDA0002890978550000037
为目标处接收信号向量,可表示为
Figure FDA0002890978550000038
其中τ0=2R0/(c-v);
In the above formula
Figure FDA0002890978550000036
Steering vector for the transmitted signal,
Figure FDA0002890978550000037
is the received signal vector at the target, which can be expressed as
Figure FDA0002890978550000038
where τ 0 =2R 0 /(cv);
信号r(t)经一定RCS的目标反射,第q个接收天线接收到的信号为:The signal r(t) is reflected by the target of a certain RCS, and the signal received by the qth receiving antenna is:
Figure FDA0002890978550000041
Figure FDA0002890978550000041
式中,
Figure FDA0002890978550000042
为第q个天线相位差,令式中ξ1ξ2=ξ,将式(1)代入式(8),则第q个接收天线的回波为:
In the formula,
Figure FDA0002890978550000042
is the phase difference of the qth antenna, let ξ 1 ξ 2 =ξ in the formula, and substitute formula (1) into formula (8), then the echo of the qth receiving antenna is:
Figure FDA0002890978550000043
Figure FDA0002890978550000043
接收信号向量为x(t)=[x1(t),x2(t),L,xQ(t)],式(9)的矩阵形式如下:The received signal vector is x(t)=[x 1 (t), x 2 (t), L, x Q (t)], the matrix form of equation (9) is as follows: x(t)=ξb(θrTt)s[ω(t-τ0)] (10)x(t)=ξb(θ rTt )s[ω(t-τ 0 )] (10) 式(10)中的b(θr)为接收天线导向矢量,
Figure FDA0002890978550000044
现考虑第q个接收天线的接收回波的宽带匹配处理,由式(10)可得第q个接收通道回波信号为:
b(θ r ) in equation (10) is the steering vector of the receiving antenna,
Figure FDA0002890978550000044
Now consider the broadband matching processing of the received echo of the qth receiving antenna, the echo signal of the qth receiving channel can be obtained from equation (10) as:
Figure FDA0002890978550000045
Figure FDA0002890978550000045
现对第q个接收天线回波信号进行宽带匹配滤波,由宽带模糊函数理论可知第p个发射通道的参考信号为:Now perform broadband matched filtering on the echo signal of the qth receiving antenna. According to the theory of broadband ambiguity function, it can be known that the reference signal of the pth transmit channel is:
Figure FDA0002890978550000046
Figure FDA0002890978550000046
第q个接收通道接收到第p个发射通道信号的匹配输出为:The matching output of the qth receiving channel receiving the signal of the pth transmitting channel is:
Figure FDA0002890978550000047
Figure FDA0002890978550000047
上式T为脉冲串持续时间,由发射信号的正交性可知,发射信号矢量中只有第p个发射通道的信号才有输出,匹配输出可简化为:The above formula T is the duration of the burst. According to the orthogonality of the transmitted signal, only the signal of the pth transmit channel in the transmitted signal vector has an output, and the matching output can be simplified as:
Figure FDA0002890978550000051
Figure FDA0002890978550000051
由上式中,当ω=ηm,τ0=ξm时,匹配滤波器有最大的输出,匹配滤波的最大输出记为
Figure FDA0002890978550000052
From the above formula, when ω=η m , τ 0m , the matched filter has the maximum output, and the maximum output of the matched filter is recorded as
Figure FDA0002890978550000052
回波信号通过宽带匹配网络后可以得到目标的距离速度信息,再通过波束形成可得到目标的方位角,波束形成的输出为:After the echo signal passes through the broadband matching network, the distance and speed information of the target can be obtained, and then the azimuth angle of the target can be obtained through beamforming. The output of the beamforming is:
Figure FDA0002890978550000053
Figure FDA0002890978550000053
式(15)中,
Figure FDA0002890978550000054
为宽带匹配滤波的输出,τ0与ω为目标的时延和多普勒频率,当
Figure FDA0002890978550000055
时,上式有最大值,此时目标方向为:
In formula (15),
Figure FDA0002890978550000054
is the output of the broadband matched filter, τ 0 and ω are the target delay and Doppler frequency, when
Figure FDA0002890978550000055
When , the above formula has a maximum value, and the target direction is:
Figure FDA0002890978550000056
Figure FDA0002890978550000056
5.根据权利要求4所述的基于多时间尺度耦合网络的雷达抗干扰波形产生方法,其特征在于:所述产生雷达发射信号的同时产生宽带回波匹配信号,正交波形的输出是由DDS波形发生器根据频率、相位编码序列产生每个
Figure FDA0002890978550000057
对应的相位处的波形的幅度值,这些幅度值产生了一系列的子脉冲,每个
Figure FDA0002890978550000058
的持续时间是由雷达的带宽所决定的;脉冲时序控制模块改变脉冲之间的脉冲周期,实现波形脉位捷变调制;通过多路DDS电路后可以产生MIMO雷达正交波形,产生MIMO雷达正交波形后,一路信号给MIMO雷达发射天线,一路作为多时间尺度耦合网络电路的输入,产生回波信号的宽带匹配滤波的参考信号;
5. The method for generating radar anti-jamming waveforms based on a multi-time scale coupling network according to claim 4, characterized in that: said generating the radar transmit signal simultaneously generates a broadband echo matching signal, and the output of the quadrature waveform is determined by the DDS The waveform generator generates each
Figure FDA0002890978550000057
The amplitude values of the waveform at the corresponding phases that produce a series of sub-pulses, each
Figure FDA0002890978550000058
The duration is determined by the bandwidth of the radar; the pulse timing control module changes the pulse period between pulses to realize the waveform pulse position agile modulation; after multiple DDS circuits, the MIMO radar quadrature waveform can be generated, and the MIMO radar positive waveform can be generated. After the waveform is crossed, one signal is sent to the MIMO radar transmitting antenna, and the other is used as the input of the multi-time scale coupling network circuit to generate the reference signal of the broadband matched filter of the echo signal;
宽带匹配网络实质是寻找回波信号二维模糊函数的最大值,在发射端产生发射信号的同时也产生了宽带匹配网络的参考信号,多时间尺度耦合网络的目的是产生一组每个子阵发射信号在时间和频率上联合移位的信号组,用于回波信号的宽带匹配滤波;在多时间尺度耦合电路部分,设每个子阵的发射信号延时τk为第k个距离门的接收信号,然后在第k个距离门内对所有多普勒频率进行匹配滤波,实现对第k个距离门回波信号的宽带脉压处理。The essence of the broadband matching network is to find the maximum value of the two-dimensional ambiguity function of the echo signal. When the transmitting signal is generated at the transmitting end, the reference signal of the broadband matching network is also generated. The purpose of the multi-time scale coupling network is to generate a set of transmissions for each sub-array. The signal group in which the signal is jointly shifted in time and frequency is used for the broadband matched filtering of the echo signal; in the multi-time scale coupling circuit part, the transmission signal delay τ k of each sub-array is set as the reception of the kth distance gate Then, all Doppler frequencies are matched and filtered in the kth range gate to realize the broadband pulse pressure processing of the echo signal of the kth range gate.
CN202110027839.4A 2021-01-11 2021-01-11 Radar anti-interference waveform generation method based on multi-time scale coupling network Active CN114764136B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110027839.4A CN114764136B (en) 2021-01-11 2021-01-11 Radar anti-interference waveform generation method based on multi-time scale coupling network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110027839.4A CN114764136B (en) 2021-01-11 2021-01-11 Radar anti-interference waveform generation method based on multi-time scale coupling network

Publications (2)

Publication Number Publication Date
CN114764136A true CN114764136A (en) 2022-07-19
CN114764136B CN114764136B (en) 2025-06-10

Family

ID=82364325

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110027839.4A Active CN114764136B (en) 2021-01-11 2021-01-11 Radar anti-interference waveform generation method based on multi-time scale coupling network

Country Status (1)

Country Link
CN (1) CN114764136B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115166658A (en) * 2022-07-27 2022-10-11 中国船舶重工集团公司第七二四研究所 Online reconfigurable array radar target simulation method
CN115857027A (en) * 2022-10-20 2023-03-28 东莞市森威电子有限公司 Pipeline blockage measuring instrument
CN116027280A (en) * 2023-03-30 2023-04-28 西安电子科技大学 A Waveform Design Method for Low Peak Sidelobe Frequency Coded Radar
CN117890898A (en) * 2024-03-01 2024-04-16 清华大学 Bistatic radar encryption target detection method based on phase center agile array

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180325456A1 (en) * 2017-05-10 2018-11-15 Boston Scientific Scimed Inc. Multi-time scale waveform for display of sensor measurements
US20190142291A1 (en) * 2015-03-23 2019-05-16 Temple University-Of The Commonwealth System Of Higher Education System and Method for Automatic Interpretation of EEG Signals Using a Deep Learning Statistical Model
CN110632559A (en) * 2019-09-17 2019-12-31 南京航空航天大学 A Waveform Design and Optimization Method for MIMO Pulse Radar Based on Application Characteristics
CN111693964A (en) * 2020-06-05 2020-09-22 西安电子科技大学 Frequency agile signal forwarding type interference suppression method based on MIMO radar

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190142291A1 (en) * 2015-03-23 2019-05-16 Temple University-Of The Commonwealth System Of Higher Education System and Method for Automatic Interpretation of EEG Signals Using a Deep Learning Statistical Model
US20180325456A1 (en) * 2017-05-10 2018-11-15 Boston Scientific Scimed Inc. Multi-time scale waveform for display of sensor measurements
CN110632559A (en) * 2019-09-17 2019-12-31 南京航空航天大学 A Waveform Design and Optimization Method for MIMO Pulse Radar Based on Application Characteristics
CN111693964A (en) * 2020-06-05 2020-09-22 西安电子科技大学 Frequency agile signal forwarding type interference suppression method based on MIMO radar

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115166658A (en) * 2022-07-27 2022-10-11 中国船舶重工集团公司第七二四研究所 Online reconfigurable array radar target simulation method
CN115857027A (en) * 2022-10-20 2023-03-28 东莞市森威电子有限公司 Pipeline blockage measuring instrument
CN116027280A (en) * 2023-03-30 2023-04-28 西安电子科技大学 A Waveform Design Method for Low Peak Sidelobe Frequency Coded Radar
CN117890898A (en) * 2024-03-01 2024-04-16 清华大学 Bistatic radar encryption target detection method based on phase center agile array
CN117890898B (en) * 2024-03-01 2024-05-14 清华大学 Bistatic radar encryption target detection method based on phase center agile array

Also Published As

Publication number Publication date
CN114764136B (en) 2025-06-10

Similar Documents

Publication Publication Date Title
EP3589970B1 (en) Method and system for obtaining an adaptive angle-doppler ambiguity function in mimo radars
CN114764136A (en) Radar anti-interference waveform generation method based on multi-time scale coupling network
CN114152916B (en) A MIMO radar waveform orthogonal method based on pilot information
CN111693983B (en) Complementary waveform construction method and module, cognitive radar system and waveform emission method
Rabaste et al. Signal waveforms and range/angle coupling in coherent colocated MIMO radar
CN112596033B (en) Array element-pulse coding MIMO radar deception jamming inhibition method
Wang et al. High resolution MIMO-HFSWR using sparse frequency waveform
US20050237236A1 (en) Method and apparatus for performing bistatic radar functions
CN104777460A (en) Method of self-adaptive clutter canceling through double-waveform phase encoding in PD radar
EP2828683A1 (en) Detection techniques
CN116930907B (en) A target coherent accumulation detection method for DDMA-MIMO radar
CN114019462B (en) Method for improving low interception performance and anti-interference performance of radar
CN115494469B (en) Slow-time MIMO radar distance ambiguity suppression method based on waveform agile phase coding
Hickman et al. Non-recurrent wideband continuous active sonar
CN109212489A (en) A kind of fuzzy clutter suppression method of the FDA-MIMO radar based on false impulse
CN110221283A (en) A kind of moving target detecting method based on FDA-MIMO
CN110988811A (en) A Grating Lobe Deblurring Method Applied to Sparse Broadband Phased Arrays
Hague et al. The generalized sinusoidal frequency modulated waveform for continuous active sonar
CN116804737A (en) MIMO radar device and MIMO radar method
Yang et al. DDMA MIMO radar system for low, slow, and small target detection
CN114895261A (en) Clutter suppression method based on multi-frequency sub-pulse coding array
CN104914430B (en) Array radar system adaptively selecting transmitting array elements according to target distance
CN104280720A (en) Method for designing transmitting directional diagram of foresight airborne radar
CN101266294A (en) A Compensation Method for Restraining the Movement of Radar Echo Signal Envelope
Landi et al. Time-orthogonal-waveform-space-time adaptive processing for distributed aperture radars

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