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 PDFInfo
- 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
Links
Images
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/28—Details of pulse systems
- G01S7/282—Transmitters
-
- 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/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting 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
Description
技术领域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:
上式中,Tr为子脉冲的平均脉冲周期,ΔTr为最小脉位跳变间隔,tp为脉冲宽度,ts为相位编码码元宽度,为第n个子脉冲的载波,μ为调频斜率,为第n个子脉冲的第m个码元,从发射信号的数学模型可知,波形优化参数为用向量表示记为Ω=[ξ1,ξ2,L,ξN]、 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, is the carrier of the nth sub-pulse, μ is the frequency modulation slope, is the m-th symbol of the n-th sub-pulse. From the mathematical model of the transmitted signal, the waveform optimization parameters are represented by a vector as Ω=[ξ 1 ,ξ 2 ,L,ξ N ],
波形参数中,每个脉冲内的相位编码位数为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:
上式中,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:
该模型的优化参数为离散的相位、频率和脉位编码,可采用遗传算法对模型进行优化,得到满足条件的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,目标到达方向为 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
定义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:
式中ξ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:
上式τk只与天线的布局有关,与目标的运动状态无关,固定时延τk可以用来表示,因此目标处合成信号为: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 to represent, so the synthesized signal at the target is:
上式中为发射信号导向矢量,为目标处接收信号向量,可表示为其中τ0=2R0/(c-v);In the above formula Steering vector for the transmitted signal, is the received signal vector at the target, which can be expressed as 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:
式中,为第q个天线相位差,令式中ξ1ξ2=ξ,将式(1)代入式(8),则第q个接收天线的回波为:In the formula, 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:
接收信号向量为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(θr)αT(θt)s[ω(t-τ0)](10)x(t)=ξb(θ r )α T (θ t )s[ω(t-τ 0 )](10)
式(10)中的b(θr)为接收天线导向矢量,现考虑第q个接收天线的接收回波的宽带匹配处理,由式(10)可得第q个接收通道回波信号为:b(θ r ) in equation (10) is the steering vector of the receiving antenna, 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:
现对第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:
第q个接收通道接收到第p个发射通道信号的匹配输出为:The matching output of the qth receiving channel receiving the signal of the pth transmitting channel is:
上式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:
由上式中,当ω=ηm,τ0=ξm时,匹配滤波器有最大的输出,匹配滤波的最大输出记为 From the above formula, when ω=η m , τ 0 =ξ m , the matched filter has the maximum output, and the maximum output of the matched filter is recorded as
回波信号通过宽带匹配网络后可以得到目标的距离速度信息,再通过波束形成可得到目标的方位角,波束形成的输出为: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:
式(15)中,为宽带匹配滤波的输出,τ0与ω为目标的时延和多普勒频率,当时,上式有最大值,此时目标方向为:In formula (15), is the output of the broadband matched filter, τ 0 and ω are the target delay and Doppler frequency, when When , the above formula has a maximum value, and the target direction is:
优选的,所述产生雷达发射信号的同时产生宽带回波匹配信号,正交波形的输出是由DDS波形发生器根据频率、相位编码序列产生每个对应的相位处的波形的幅度值,这些幅度值产生了一系列的子脉冲,每个的持续时间是由雷达的带宽所决定的;脉冲时序控制模块改变脉冲之间的脉冲周期,实现波形脉位捷变调制;通过多路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. The amplitude values of the waveform at the corresponding phases that produce a series of sub-pulses, each 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:
上式中,Tr为子脉冲的平均脉冲周期,ΔTr为最小脉位跳变间隔,tp为脉冲宽度,ts为相位编码码元宽度,为第n个子脉冲的载波,μ为调频斜率,为第n个子脉冲的第m个码元,从发射信号的数学模型可知,波形优化参数为用向量表示记为Ω=[ξ1,ξ2,L,ξN]、 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, is the carrier of the nth sub-pulse, μ is the frequency modulation slope, is the m-th symbol of the n-th sub-pulse. From the mathematical model of the transmitted signal, the waveform optimization parameters are represented by a vector as Ω=[ξ 1 ,ξ 2 ,L,ξ N ],
波形参数中,每个脉冲内的相位编码位数为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:
上式中,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:
该模型的优化参数为离散的相位、频率和脉位编码,可采用遗传算法对模型进行优化,得到满足条件的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,目标到达方向为 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
定义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:
式中ξ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:
上式τk只与天线的布局有关,与目标的运动状态无关,固定时延τk可以用来表示,因此目标处合成信号为: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 to represent, so the synthesized signal at the target is:
上式中为发射信号导向矢量,为目标处接收信号向量,可表示为其中τ0=2R0/(c-v);In the above formula Steering vector for the transmitted signal, is the received signal vector at the target, which can be expressed as 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:
式中,为第q个天线相位差,令式中ξ1ξ2=ξ,将式(1)代入式(8),则第q个接收天线的回波为:In the formula, 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:
接收信号向量为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(θr)αT(θt)s[ω(t-τ0)] (10)x(t)=ξb(θ r )α T (θ t )s[ω(t-τ 0 )] (10)
式(10)中的b(θr)为接收天线导向矢量,现考虑第q个接收天线的接收回波的宽带匹配处理,由式(10)可得第q个接收通道回波信号为:b(θ r ) in equation (10) is the steering vector of the receiving antenna, 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:
现对第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:
第q个接收通道接收到第p个发射通道信号的匹配输出为:The matching output of the qth receiving channel receiving the signal of the pth transmitting channel is:
上式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:
由上式中,当ω=ηm,τ0=ξm时,匹配滤波器有最大的输出,匹配滤波的最大输出记为 From the above formula, when ω=η m , τ 0 =ξ m , the matched filter has the maximum output, and the maximum output of the matched filter is recorded as
回波信号通过宽带匹配网络后可以得到目标的距离速度信息,再通过波束形成可得到目标的方位角,波束形成的输出为: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:
式(15)中,为宽带匹配滤波的输出,τ0与ω为目标的时延和多普勒频率,当时,上式有最大值,此时目标方向为:In formula (15), is the output of the broadband matched filter, τ 0 and ω are the target delay and Doppler frequency, when When , the above formula has a maximum value, and the target direction is:
产生雷达发射信号的同时产生宽带回波匹配信号,正交波形的输出是由DDS波形发生器根据频率、相位编码序列产生每个对应的相位处的波形的幅度值,这些幅度值产生了一系列的子脉冲,每个的持续时间是由雷达的带宽所决定的;脉冲时序控制模块改变脉冲之间的脉冲周期,实现波形脉位捷变调制;通过多路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. The amplitude values of the waveform at the corresponding phases that produce a series of sub-pulses, each 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)
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)
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)
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 |
-
2021
- 2021-01-11 CN CN202110027839.4A patent/CN114764136B/en active Active
Patent Citations (4)
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)
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 |