CN114696788B - Wave form and filter combined cognition design method for resisting multi-main lobe interference - Google Patents

Wave form and filter combined cognition design method for resisting multi-main lobe interference Download PDF

Info

Publication number
CN114696788B
CN114696788B CN202210379049.7A CN202210379049A CN114696788B CN 114696788 B CN114696788 B CN 114696788B CN 202210379049 A CN202210379049 A CN 202210379049A CN 114696788 B CN114696788 B CN 114696788B
Authority
CN
China
Prior art keywords
interference
code words
intercepting
expressed
section
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.)
Active
Application number
CN202210379049.7A
Other languages
Chinese (zh)
Other versions
CN114696788A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202210379049.7A priority Critical patent/CN114696788B/en
Publication of CN114696788A publication Critical patent/CN114696788A/en
Application granted granted Critical
Publication of CN114696788B publication Critical patent/CN114696788B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0201Wave digital filters
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0211Frequency selective networks using specific transformation algorithms, e.g. WALSH functions, Fermat transforms, Mersenne transforms, polynomial transforms, Hilbert transforms
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H2017/0072Theoretical filter design
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Abstract

The invention discloses a wave form and filter joint cognition design method for resisting multi-main lobe interference, which comprises the following steps: s1, constructing a C & I, direct intermittent sampling and repeated intermittent sampling interference signal model and a mismatch filter model based on a constant-mode phase coded signal; s2, under the condition that interference parameters are completely known or deviation exists, establishing a wave form correlation function peak value level, a correlation function main lobe template matching error and an interference energy weighting and minimizing criterion transmitting-receiving joint optimization problem under the constraint of signal to noise ratio; s3, converting the optimization problem into a minimization problem related to the independent variable function, and solving by using an iterative L-BFGS algorithm. The optimized signal designed by the invention has good Doppler tolerance, multi-main-lobe interference resistance and lower peak sidelobes after mismatch filtering.

Description

Wave form and filter combined cognition design method for resisting multi-main lobe interference
Technical Field
The invention belongs to the field of radar anti-interference surgery, and particularly relates to a wave form and filter combined cognition design method for resisting multi-main-lobe interference.
Background
In recent years, with the high-speed development of modern electronic technology, electronic countermeasure is increasingly strong, various interference patterns are layered, and the normal operation of a radar system is seriously hindered. For main lobe interference, the interference energy has absolute advantages, the main lobe interference is highly overlapped with a target in dimensions such as space time frequency, and the like, the effect of the existing interference suppression means is poor, and the main lobe interference is still one of the problems to be solved in the radar industry. The general main lobe interference resistance method comprises the following steps: signal processing and waveform design.
In terms of a signal processing main lobe interference resisting method, a common signal processing method includes: blind source separation and filtering. The blind source separation can separate the mixed signals without priori knowledge, and if the blind source separation main lobe interference resisting algorithm based on matrix combined diagonalization feature vectors is provided in the document [ G.Huang, L.Yang, G.Su.Blind source separation used for radar anti-jamming [ C ].2003International Conference on Neural Networks and Signal Processing,Nanjing,China,2003:1382-1385 ], the target signals and the interference signals are separated, so that interference suppression is realized, and the influence of main lobe interference on radar detection performance can be effectively reduced. However, after the mixed signal is separated, there is often a small amount of interference residual, and a single type of interference is aimed at. The filtering process mainly utilizes the difference of interference signals and target echoes in time, frequency, space, polarization and other domains, and designs filter parameters in different dimensions to achieve the effect of interference suppression, such as document [ S.Zhang, Y.Yang, G.Cui, et al range-velocity jamming suppression algorithm based on adaptive iterative filtering [ C ] 2016IEEE Radar Conference,Philadelphia,PA,USA,2016:1-6 ], and the distance-speed joint spoofing interference suppression is achieved by designing a mismatch filter for each distance-Doppler unit. However, this method requires that the waveform of the interfering signal be precisely known and is only applicable to a certain interference type.
The method for resisting main lobe interference by waveform design is characterized in that the difference of target echo and interference signals in time, frequency, space, polarization and other domains is amplified by designing the information of the phase, frequency and the like of intra-pulse and inter-pulse waveforms, so that the purpose of resisting interference is achieved, such as document [ Xu Leilei ], radar waveform design and a plurality of technical researches on resisting main lobe active interference [ D ]. Siemens, siemens electronic technology university, 2019 ] ], the interference signals are compressed by utilizing a waveform optimization method based on the amplitude difference between angular domain waveforms and the interference signals, and DRFM forward interference is effectively resisted. However, this method can only resist a single type of interference, and cannot effectively resist the composite interference where multiple interference types are superimposed.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide a waveform and filter combined cognition design method for resisting multi-main-lobe interference.
The aim of the invention is realized by the following technical scheme: a wave form and filter joint cognition design method for resisting multi-main lobe interference comprises the following steps:
s1, constructing a C & I, direct intermittent sampling and repeated intermittent sampling interference signal model and a mismatch filter model based on a constant-mode phase coded signal;
s2, under the condition that interference parameters are completely known or deviation exists, establishing a wave form correlation function peak value level, a correlation function main lobe template matching error and an interference energy weighting and minimizing criterion transmitting-receiving joint optimization problem under the constraint of signal to noise ratio;
s3, converting the optimization problem into a minimization problem related to the independent variable function, and solving by using an iterative L-BFGS algorithm.
Further, the step S1 specifically includes the following sub-steps:
s11, the constant-mode phase coding signal is expressed as:
Figure BDA0003591934730000021
in the formula [ (] T Representing the transpose s n For transmitting codewords of the waveform, n=1, 2, …, N s ,N s The number of codewords for the phase encoded waveform; the nth element s n Expressed as:
Figure BDA0003591934730000022
in the method, in the process of the invention,
Figure BDA0003591934730000023
representation->
Figure BDA0003591934730000024
Phase of the nth codeword, +.>
Figure BDA0003591934730000025
Encoding the phase of the signal s for a constant mode phase;
s12, establishing C&I interference signal model: for C&I interference, suppose that the DRFM jammer intercepts radar signals, from s i Starting, m code words are taken as one section, each section is forwarded for l times, k sections are intercepted in total, and the whole C&I interference signal codeword length N s =klm;C&The I interference signal is:
Figure BDA0003591934730000026
wherein H is t T is more than or equal to 1 and less than or equal to k in the t section intercepted by the jammer; h t The represented codeword is expressed as:
Η t =[s (t-1)ml+i ,s (t-1)ml+i+1 ,…,s (t-1)ml+i+m-1 ];
s13, directly intermittently sampling an interference signal model: for direct intermittent sampling interference, suppose that a DRFM (digital radio frequency modulation) jammer intercepts radar signals, and s is a reference number i Firstly, m code words are taken as a section, and are intercepted and forwarded once every ml code words, k sections are intercepted in total, and interference is directly and intermittently sampledThe signals are as follows:
J 2 =[Η′ 1 ,Η′ 0 ,Η′ 2 ,Η′ 0 ,…,Η′ k ,Η′ 0 ] T
in the formula, H' 0 =0 1×m(l-1) ,0 1×m(l-1) A 0 vector representing 1×m (l-1); h' t For the t-th segment intercepted by the jammer, the represented codeword is expressed as:
Η′ t =[s (t-1)ml+i ,s (t-1)ml+i+1 ,…s (t-1)ml+i+m-1 ];
s14, repeating the intermittent sampling interference signal model: for repeated intermittent sampling interference, assume that a DRFM (digital radio frequency modulation) jammer intercepts radar signals, and then s i Firstly, m code words are taken as a section, the code words are intercepted once every ml, the code words are repeated q times and then forwarded, k sections are intercepted in total, and the repeated intermittent sampling interference signals are as follows:
Figure BDA0003591934730000031
in the formula, H 0 =0 1×m(l-q) ,Η″ t For the t-th segment intercepted by the jammer, the represented codeword is expressed as:
Η″ t =[s (t-1)ml+i ,s (t-1)ml+i+1 ,…s (t-1)ml+i+m-1 ];
s15, a mismatch filter signal model: assume that the length of a corresponding mismatch filter h of a radar receiving end is N h ,N h ≥N s The mismatch filter h is denoted as:
Figure BDA0003591934730000032
the nth element h n Expressed as:
Figure BDA0003591934730000033
where a is the magnitude vector of the mismatch filter h, a n Is the nth element in a; θ is the phase vector of the mismatch filter h, θ n Is the n element in theta; h=a.exp (j θ), and a indicates the Hadamard product.
Further, the step S2 specifically includes:
s21, designing optimization criteria:
(1) Distance sidelobe: the output of the phase encoded signal s through the mismatch filter h is represented at a distance displacement j as:
Figure BDA0003591934730000034
wherein, represents conjugation, Γ 1 Represents a value interval, Γ 1 =[-(N s +N h )/2+1,(N s +N h )/2-1];
Let ω be j For the widened main lobe region, j=0, ±1, …, ±m, M is a main lobe width control parameter;
all the distance main lobe levels are stacked in vector omega main Expressed as:
ω main =[ω -M ,…,ω -101 ,…,ω M ] T
all range side lobe levels are stacked in vector ω side Expressed as:
Figure BDA0003591934730000041
let the sidelobe region be denoted as Γ 2 ,Γ 2 Expressed as:
Γ 2 =[-(N s +N h )/2+1,-M+1]∪[M-1,(N s +N h )/2-1]
the range side lobe is reduced by minimizing the peak side lobe, namely:
Figure BDA0003591934730000042
(2) Loss of signal-to-noise ratio:
|h H h-N s |≤η 1
0 -N s |≤η 2
in the formula, the constant eta 1 ≥0,η 2 More than or equal to 0, and the superscript H represents conjugate transposition;
(3) Main lobe control: assume that
Figure BDA0003591934730000043
Is the desired main lobe, i.e. q is a vector of dimension 2m+1; vector e is the error vector of the desired main lobe and the designed main lobe, expressed as:
e=ω main -q
the nth element e (n) of the error vector e is expressed as:
e(n)=ω n-M-1 -q(n),1≤n≤2M+1
the main lobe shape is maintained with a method that minimizes the maximum main lobe matching error, namely:
Figure BDA0003591934730000044
(4) Multi-main lobe interference resistance: interference signal J m The output result via the mismatch filter h is represented at a distance displacement j as:
Figure BDA0003591934730000045
J m (n) represents J m In (2) the n-th element, all distance levels in the above formula are stacked in vector ω' m Expressed as:
Figure BDA0003591934730000046
when the interference parameters are completely known: cognitive radar senses multi-main-lobe interference and generates interference parametersUnder completely known conditions, different interference signals J m Distance level stack vector ω 'generated by mismatch filter h' m ,m=1,2,…,N J All stacked on vector omega J ,ω J Expressed as:
Figure BDA0003591934730000047
N J indicating the number of interferences;
when the interference parameter has deviation: let it be assumed that actual C&The number of interception segments of the interference I is k, the forwarding frequency is l, the number of interception segments of the perceived parameter is k ', and the forwarding frequency is l'. Assuming that the error is small, there is k=k '+l, l=l' -L, L being a natural number greater than 0. In the presence of errors, C&I interference distance level stack vector
Figure BDA0003591934730000051
Expressed as:
Figure BDA0003591934730000052
in the formula omega' 1 C is generated by intercepting segment number k' and forwarding frequency l&Stacking vectors of all distance levels after filtering processing of the I interference and mismatch filter;
Figure BDA0003591934730000053
c is generated by intercepting segment number k '+i and forwarding frequency l' -i&Stacking vectors of all distance levels after filtering processing of the I interference and mismatch filter; />
Figure BDA0003591934730000054
C is generated by intercepting segment number k '-i and forwarding frequency l' +i&Stacking vectors of all distance levels after filtering processing of the I interference and mismatch filter, i=1, 2, …, L;
the actual direct intermittent sampling interference generation mode is assumed to be: taking m code words as a section, and every mIntercepting and forwarding l code words once, intercepting k sections altogether, wherein the perceived generation mode of direct intermittent sampling interference is as follows; taking m code words as a section, intercepting and forwarding every ml 'code words, intercepting k' sections altogether, and assuming that the error is smaller, k=k '+L, and l=l' -L are included; direct intermittent sampling of interference distance level stacking vectors in the presence of deviations
Figure BDA0003591934730000055
Expressed as:
Figure BDA0003591934730000056
ω′ 2 taking m code words as a section, intercepting and forwarding every ml code words, and intercepting the stacking vectors of direct intermittent sampling interference generated by the k section and all distance levels after filtering treatment of a mismatch filter;
Figure BDA0003591934730000057
taking m code words as a section, intercepting and forwarding every m (l '-i) code words, and intercepting direct intermittent sampling interference generated by k' +i sections and stacking vectors of all distance levels after filtering treatment of a mismatch filter; />
Figure BDA0003591934730000058
Taking m code words as a section, intercepting and forwarding every m (l '+i) code words, and intercepting the stacking vectors of direct intermittent sampling interference generated by the section k' -i and all distance levels after filtering treatment of a mismatch filter;
the actual repeated intermittent sampling interference generation mode is assumed to be: taking m code words as a section, intercepting and forwarding q times every ml code words, and intercepting k sections altogether; the perceived repeated intermittent sampling interference generation mode is as follows; taking m code words as a section, intercepting and forwarding q times every ml 'code words, intercepting k' sections altogether, and assuming that the error is smaller, k=k '+L, and l=l' -L are included; repeated intermittent sampling of interference distance level stacking vectors in the presence of deviations
Figure BDA0003591934730000059
Expressed as: />
Figure BDA00035919347300000510
ω′ 3 Taking m code words as a section, intercepting and forwarding q times every ml 'code words, and intercepting the stacking vectors of the repeated intermittent sampling interference generated by the k' section and all distance levels after filtering treatment of a mismatch filter altogether;
Figure BDA00035919347300000511
taking m code words as a section, intercepting and forwarding q times every m (l '-i) code words, and intercepting repeated intermittent sampling interference generated by k' +i sections and stacking vectors of all distance levels after filtering treatment of a mismatch filter; />
Figure BDA00035919347300000512
Taking m code words as a section, intercepting and forwarding q times every m (l '+i) code words, and intercepting the repeated intermittent sampling interference generated by the k' -i section and the stacking vectors of all distance levels after the filtering treatment of a mismatch filter;
all distance levels ω of all interference output via mismatched filters J The method comprises the following steps: if there is no deviation in several disturbance parameters, all distance levels omega J Is that
Figure BDA0003591934730000061
If there is a deviation in one or more of the interference parameters, then +.>
Figure BDA0003591934730000062
The corresponding vector in (a) is replaced by a distance level stack vector when errors exist;
all distance levels of all interference output via the mismatch filter are compressed using a min-max optimization criterion, namely:
Figure BDA0003591934730000063
s22, establishing an optimization problem: according to the optimization criterion of the S21 design, the optimization problem of the phase coding signal S and the mismatch filter h resisting the multi-main lobe interference is expressed as follows:
Figure BDA0003591934730000064
s.t.|h H h-N s |≤γ 1
wherein lambda is 1 ,λ 2 Is a preset weight coefficient.
Further, the specific implementation method of the S3 is as follows:
s31, optimizing problem conversion: converting the optimization problem into a minimization problem with respect to the independent variable function, namely:
Figure BDA0003591934730000065
wherein lambda is 3 For the weight coefficient to be set in advance, I Represents an infinite norm;
defining an objective function as:
f(x)=||ω side || p1 ||ω J || p2 ||e|| p3 ||h H h-N s || p
in the method, in the process of the invention, I p Represents the p-norm; vector quantity
Figure BDA0003591934730000066
Is composed of vector->
Figure BDA0003591934730000067
a and θ are sequentially formed column vectors expressed as:
Figure BDA0003591934730000068
s32, solving an optimization problem: and solving the problem by using an L-BFGS algorithm based on iteration, solving the minimum value of the objective function by using the L-BFGS algorithm, continuously iterating until the objective function f (x) reaches the minimum descent epsilon, stopping iterating, and outputting x.
The beneficial effects of the invention are as follows: firstly, constructing a C & I, direct intermittent sampling and repeated intermittent sampling interference signal model based on a constant-mode phase coding signal; then under the condition that the interference parameters are completely known or deviation exists, the problem of transmitting-receiving joint optimization of waveform correlation function peak level, correlation function main lobe template matching error and interference energy weighted sum minimization criterion under the constraint of signal to noise ratio is established; finally, the optimization problem of non-convex, non-smooth and constraint is converted into the minimization problem of the function of independent variables (composed of signal phase, filter amplitude and phase) by using methods such as penalty function method, norm conversion and the like, and the iterative L-BFGS algorithm is used for solving. The optimized signal designed by the invention has good Doppler tolerance, multi-main-lobe interference resistance and lower peak sidelobes after mismatch filtering.
Drawings
FIG. 1 is a flow chart of a method for joint cognitive design of waveforms and filters for multi-lobe interference resistance in accordance with the present invention;
FIG. 2 is a schematic diagram of C & I interference generation;
FIG. 3 is a schematic diagram of direct intermittent sampling interference generation;
FIG. 4 is a schematic diagram of a repeated intermittent sampling disturbance generation;
FIG. 5 is a flowchart of an iteration-based L-BFGS algorithm;
FIG. 6 is an optimized signal Doppler tolerance designed when the interference parameters are fully known and there is a deviation in the interference parameters;
FIG. 7 is a graph showing the result of the matched filtering of the LFM signal and its corresponding C & I interference, direct intermittent sampling interference and repeated intermittent sampling interference;
FIG. 8 is a result of mismatch filtering the optimized signal and its corresponding multi-lobe interference when the interference parameter is completely known and has a deviation from the interference parameter;
fig. 9 is a comparison diagram of detection results of an optimized signal designed when LFM signals and interference parameters are completely known and have deviations from the interference parameters.
Detailed Description
The technical scheme of the invention is further described below with reference to the accompanying drawings.
As shown in fig. 1, the method for designing the waveform and filter joint cognition for resisting multi-main lobe interference comprises the following steps:
s1, constructing a C & I, direct intermittent sampling and repeated intermittent sampling interference signal model and a mismatch filter model based on a constant-mode phase coded signal; the method specifically comprises the following sub-steps:
s11, the constant-mode phase coding signal is expressed as:
Figure BDA0003591934730000071
in the formula [ (] T Representing the transpose s n For transmitting codewords of the waveform, n=1, 2, …, N s ,N s The number of codewords for the phase encoded waveform; the nth element s n Expressed as:
Figure BDA0003591934730000072
in the method, in the process of the invention,
Figure BDA0003591934730000073
representation->
Figure BDA0003591934730000074
Phase of the nth codeword, +.>
Figure BDA0003591934730000075
Encoding the phase of the signal s for a constant mode phase;
s12, establishing C&I interference signal model: for C&I trunkScrambling, suppose that the DRFM jammer intercepts radar signals, from s i Starting, m code words are taken as one section, each section is forwarded for l times, k sections are intercepted in total, and the whole C&I interference signal codeword length N s =klm; as shown in FIG. 2, C&The I interference signal is:
Figure BDA0003591934730000081
wherein H is t T is more than or equal to 1 and less than or equal to k in the t section intercepted by the jammer; h t The represented codeword is expressed as:
Η t =[s (t-1)ml+i ,s (t-1)ml+i+1 ,…,s (t-1)ml+i+m-1 ]
A 1 as a matrix, it can be expressed as:
Figure BDA0003591934730000082
1. 1 represents that the elements at positions i to i+m-1 are all 1;
s13, directly intermittently sampling an interference signal model: for direct intermittent sampling interference, suppose that a DRFM (digital radio frequency modulation) jammer intercepts radar signals, and s is a reference number i Firstly, m code words are taken as a section, and are intercepted and forwarded once every ml code words, k sections are intercepted in total, as shown in fig. 3, and the direct intermittent sampling interference signals are as follows:
J 2 =[Η′ 1 ,Η′ 0 ,Η′ 2 ,Η′ 0 ,…,Η′ k ,Η′ 0 ] T
=A 2 s
in the formula, H' 0 =0 1×m(l-1) ,0 1×m(l-1) A 0 vector representing 1×m (l-1); h' t For the t-th segment intercepted by the jammer, the represented codeword is expressed as:
Η′ t =[s (t-1)ml+i ,s (t-1)ml+i+1 ,…s (t-1)ml+i+m-1 ]
A 2 as a matrix, it can be expressed as:
Figure BDA0003591934730000091
s14, repeating the intermittent sampling interference signal model: for repeated intermittent sampling interference, assume that a DRFM (digital radio frequency modulation) jammer intercepts radar signals, and then s i Firstly, m code words are taken as a section, the code words are intercepted once every ml, the code words are repeated q times and then forwarded, k sections are intercepted in total, and as shown in fig. 4, the repeated intermittent sampling interference signals are as follows:
Figure BDA0003591934730000092
in the formula, H 0 =0 1×m(l-q) ,Η″ t For the t-th segment intercepted by the jammer, the represented codeword is expressed as:
Η″ t =[s (t-1)ml+i ,s (t-1)ml+i+1 ,…s (t-1)ml+i+m-1 ]
A 3 as a matrix, it can be expressed as:
Figure BDA0003591934730000093
s15, a mismatch filter signal model: assume that the length of a corresponding mismatch filter h of a radar receiving end is N h ,N h ≥N s The mismatch filter h is denoted as:
Figure BDA0003591934730000094
the nth element h n Expressed as:
Figure BDA0003591934730000095
where a is the magnitude vector of the mismatch filter h, a n Is the nth element in a; θ is the phase vector of the mismatch filter h, θ n Is the n element in theta; h=a.exp (j θ), and a indicates the Hadamard product.
S2, under the condition that interference parameters are completely known or deviation exists, establishing a wave form correlation function peak value level, a correlation function main lobe template matching error and an interference energy weighting and minimizing criterion transmitting-receiving joint optimization problem under the constraint of signal to noise ratio;
the step S2 specifically comprises the following steps:
s21, designing optimization criteria:
(1) Distance sidelobe: in practical application, the distance sidelobe level output by the mismatch filter should be as low as possible, so as to ensure that a weak target is not submerged by a high sidelobe of a strong target echo signal;
the output of the phase encoded signal s through the mismatch filter h is represented at a distance displacement j as:
Figure BDA0003591934730000101
wherein, represents conjugation, Γ 1 Represents a value interval, Γ 1 =[-(N s +N h )/2+1,(N s +N h )/2-1];
Let ω be j For the widened main lobe region, j=0, ±1, …, ±m, M is a main lobe width control parameter;
all the distance main lobe levels are stacked in vector omega main Expressed as:
ω main =[ω -M ,…,ω -101 ,…,ω M ] T
all range side lobe levels are stacked in vector ω side Expressed as:
Figure BDA0003591934730000102
let the sidelobe region be denoted as Γ 2 ,Γ 2 Represented as:
Γ 2 =[-(N s +N h )/2+1,-M+1]∪[M-1,(N s +N h )/2-1]
The range side lobe is reduced by minimizing the peak side lobe, namely:
Figure BDA0003591934730000103
(2) Loss of signal-to-noise ratio, which can cause a certain loss in the process of designing a mismatch filter. There is therefore a need to control the signal-to-noise loss, namely:
|h H h-N s |≤η 1
0 -N s |≤η 2
in the formula, the constant eta 1 ≥0,η 2 More than or equal to 0, wherein the specific value can be defined by a user, and the superscript H represents conjugate transposition;
(3) Main lobe control: assume that
Figure BDA0003591934730000104
Is the desired main lobe, i.e. q is a vector of dimension 2m+1; vector e is the error vector of the desired main lobe and the designed main lobe, expressed as:
e=ω main -q
the nth element e (n) of the error vector e is expressed as:
e(n)=ω n-M-1 -q(n),1≤n≤2M+1
the main lobe shape is maintained with a method that minimizes the maximum main lobe matching error, namely:
Figure BDA0003591934730000111
(4) Multi-main lobe interference resistance: in practical applications, the distance level between the mismatched filter and the multi-main lobe interference output should be as low as possible, so as to reduce the influence of the interference signal on the real target detection. Interference signal J m Output result through mismatch filter hExpressed at distance displacement j as:
Figure BDA0003591934730000112
J m (n) represents J m In (2) the n-th element, all distance levels in the above formula are stacked in vector ω' m Expressed as:
Figure BDA0003591934730000113
when the interference parameters are completely known: the cognitive radar senses multi-main-lobe interference, and different interference signals J are transmitted under the condition that interference parameters are completely known m Distance level stack vector ω 'generated by mismatch filter h' m ,m=1,2,…,N J All stacked on vector omega J ,ω J Expressed as:
Figure BDA0003591934730000114
N J indicating the number of disturbances, in this example 3;
when the interference parameter has deviation: the cognitive radar senses multi-main-lobe interference, and under the condition that other types of interference signal parameters are known, C&The interference parameters of I have deviations. Let it be assumed that actual C&The number of interception segments of the interference I is k, the forwarding frequency is l, the number of interception segments of the perceived parameter is k ', and the forwarding frequency is l'. Assuming that the error is small, there is k=k '+l, l=l' -L, L being a natural number greater than 0. In the presence of errors, C&I interference distance level stack vector
Figure BDA0003591934730000115
Expressed as: />
Figure BDA0003591934730000116
Wherein omega is 1 'is C generated by intercepting segment number k' and forwarding frequency l&Stacking vectors of all distance levels after filtering processing of the I interference and mismatch filter;
Figure BDA0003591934730000117
c is generated by intercepting segment number k '+i and forwarding frequency l' -i&Stacking vectors of all distance levels after filtering processing of the I interference and mismatch filter; />
Figure BDA0003591934730000118
C is generated by intercepting segment number k '-i and forwarding frequency l' +i&Stacking vectors of all distance levels after filtering processing of the I interference and mismatch filter, i=1, 2, …, L;
similarly, under the condition that other types of interference signal parameters are known, the estimation of the direct intermittent sampling interference parameters has a certain error. The actual direct intermittent sampling interference generation mode is assumed to be: taking m code words as a section, intercepting and forwarding every ml code words, intercepting k sections altogether, and the perceived generation mode of direct intermittent sampling interference is as follows; taking m code words as a section, intercepting and forwarding every ml 'code words, intercepting k' sections altogether, and assuming that the error is smaller, k=k '+L, and l=l' -L are included; direct intermittent sampling of interference distance level stacking vectors in the presence of deviations
Figure BDA0003591934730000121
Expressed as:
Figure BDA0003591934730000122
ω′ 2 taking m code words as a section, intercepting and forwarding every ml code words, and intercepting the stacking vectors of direct intermittent sampling interference generated by the k section and all distance levels after filtering treatment of a mismatch filter;
Figure BDA0003591934730000123
m code words are used as a section, every otherm (l '-i) code words are intercepted once and forwarded, and stacking vectors of direct intermittent sampling interference generated by k' +i sections and all distance levels after filtering treatment of a mismatch filter are intercepted altogether; />
Figure BDA0003591934730000124
Taking m code words as a section, intercepting and forwarding every m (l '+i) code words, and intercepting the stacking vectors of direct intermittent sampling interference generated by the section k' -i and all distance levels after filtering treatment of a mismatch filter;
under the condition that other types of interference signal parameters are known, the estimation of the repeated intermittent sampling interference parameters has a certain error. The actual repeated intermittent sampling interference generation mode is assumed to be: taking m code words as a section, intercepting and forwarding q times every ml code words, and intercepting k sections altogether; the perceived repeated intermittent sampling interference generation mode is as follows; taking m code words as a section, intercepting and forwarding q times every ml 'code words, intercepting k' sections altogether, and assuming that the error is smaller, k=k '+L, and l=l' -L are included; repeated intermittent sampling of interference distance level stacking vectors in the presence of deviations
Figure BDA0003591934730000125
Expressed as:
Figure BDA0003591934730000126
ω′ 3 taking m code words as a section, intercepting and forwarding q times every ml 'code words, and intercepting the stacking vectors of the repeated intermittent sampling interference generated by the k' section and all distance levels after filtering treatment of a mismatch filter altogether;
Figure BDA0003591934730000127
taking m code words as a section, intercepting and forwarding q times every m (l '-i) code words, and intercepting repeated intermittent sampling interference generated by k' +i sections and stacking vectors of all distance levels after filtering treatment of a mismatch filter; />
Figure BDA0003591934730000128
Taking m code words as a section, intercepting and forwarding q times every m (l '+i) code words, and intercepting the repeated intermittent sampling interference generated by the k' -i section and the stacking vectors of all distance levels after the filtering treatment of a mismatch filter;
all distance levels ω of all interference output via mismatched filters J The method comprises the following steps: no deviation exists in a plurality of interference parameters (whether the deviation exists in the interference parameters or not is obtained by the interference sensing of the cognitive radar), and the total distance level omega J Is that
Figure BDA0003591934730000129
If there is a deviation in one or more of the interference parameters, then +.>
Figure BDA00035919347300001210
The corresponding vector in (a) is replaced by a distance level stack vector when errors exist; for example, when C&When the I interference parameter has errors, the obtained total distance level omega J The method comprises the following steps: />
Figure BDA00035919347300001211
All distance levels of all interference output via the mismatch filter are compressed using a min-max optimization criterion, namely:
Figure BDA0003591934730000131
s22, establishing an optimization problem: according to the optimization criterion of the S21 design, the optimization problem of the phase coding signal S and the mismatch filter h resisting the multi-main lobe interference is expressed as follows:
Figure BDA0003591934730000132
s.t.|h H h-N s |≤γ 1
wherein lambda is 1 ,λ 2 Is a preset weight coefficient.
S3, converting the optimization problem which is non-convex, non-smooth and provided with constraint into the minimization problem of the function of the independent variable (composed of signal phase, filter amplitude and phase) by using a penalty function method, a norm conversion method and the like, and solving by using an iterative L-BFGS algorithm.
The specific implementation method of S3 is as follows:
s31, optimizing problem conversion: the optimization problem of non-convex, non-smooth and constraint is converted into the minimization problem of the function of the independent variables (composed of signal phase, filter amplitude and phase) by using the methods of penalty function, norm conversion and the like, namely:
Figure BDA0003591934730000133
wherein lambda is 3 For the weight coefficient to be set in advance, I Represents an infinite norm;
defining an objective function as:
f(x)=||ω side || p1 ||ω J || p2 ||e|| p3 ||h H h-N s || p
in the method, in the process of the invention, I p Represents the p-norm; vector quantity
Figure BDA0003591934730000134
Is composed of vector->
Figure BDA0003591934730000135
a and θ are sequentially formed column vectors expressed as:
Figure BDA0003591934730000136
s32, solving an optimization problem: based on the analysis, the problem is solved by using an L-BFGS algorithm based on iteration, and the main idea is to solve the minimum value of the objective function by using the L-BFGS algorithm, iterate continuously until the objective function f (x) reaches the minimum descent epsilon, stop iterating, and output x. As shown in fig. 5, the specific steps are as follows:
step I: initializing variable x 0 Epsilon, p, mu: set x 0 As a starting value, epsilon is the minimum descent of the objective function, let iteration number w=1, parameter p=2, μ=2;
step II: calculate the current function value f (x w ): will x w-1 As an initial value, the minimization function f (x) is used to find x by using the L-BFGS algorithm w Let f w =f(x w );
Step III: if |f w-1 -f w If the I is less than or equal to epsilon, obtaining x w And stopping the iteration; if not, p=μp, w=w+1, and step II is skipped.
The step II, where the minimization of the function f (x) is required, can be implemented using the L-BFGS algorithm. Compared with the quasi-Newton method, the L-BFGS algorithm utilizes a small quantity of vectors to estimate the Hessian matrix, so that the calculated amount is greatly reduced, and the method is easy to realize. The L-BFGS algorithm flow is as follows:
first, an initial iteration number w=1 is set, and an initial vector x 0 And calculate the objective function value f 1 Sum of derivative values
Figure BDA0003591934730000141
To derive the sign and select the starting direction +.>
Figure BDA0003591934730000142
And the algorithm updates the variable n. The specific circulation process is as follows:
step 1: for the w-th iteration, calculate f w
Figure BDA0003591934730000143
And g w And obtaining the step length lambda by using a backtracking line searching algorithm w
Step 2: calculate q w =λ w g w ,x w+1 =x ww g w
Step 3: calculation of
Figure BDA0003591934730000144
And->
Figure BDA0003591934730000145
Step 4: order the
Figure BDA0003591934730000146
Step 5: fori=w, w-1, …, w-n+1
Figure BDA0003591934730000147
p=p-t i y i
end
Step 6:
Figure BDA0003591934730000148
step 7: fori=w, w-1, …, w-n+1
Figure BDA0003591934730000149
z=z+(t i -β)q i
end
Step 8: g w+1 =-z;
Step 9: w=w+1;
and finally, when the termination condition is met, obtaining an optimization result.
It should be explained that, in the step 5-7, the product of the Hessian matrix and the gradient p is mainly constructed by using the gradient information of the adjacent n iterative functions, and the search direction of the next step is obtained according to the obtained vector z. Step 1 requires obtaining a step length lambda by using a backtracking line search algorithm w The specific process of the backtracking search algorithm is as follows:
step a: selecting a lambda w >0,ρ,c 1 ∈(0,1);
Step b: if it is
Figure BDA00035919347300001410
Let lambda get w =ρλ w
Step c: repeating step b until the relation is satisfied
Figure BDA00035919347300001411
Output lambda w
Simulation verification and analysis
Simulation parameters:
the phase coded signal has a time width of 60 mu s, a bandwidth of 5MHz, an LFM signal has a time width of 60 mu s, a bandwidth of 5MHz, a mismatch filter has a time width of 120 mu s, a radar system carrier frequency of 2GHz, a number of transmitted pulses of 100, a pulse repetition period of 480 mu s, a signal to noise ratio of 0dB, a signal to interference ratio of 20dB, a target distance of 45km, a C & I interference distance of 46km, a direct intermittent sampling interference distance of 44km, a repeated intermittent sampling interference distance of 48km, a target speed of 30m/s, a C & I interference speed of 32m/s, a direct intermittent sampling speed of 28m/s, and a repeated intermittent sampling interference speed of 26m/s.
The interference parameters are completely known:
c & I interference generation mode: 10 segments are intercepted, each segment is forwarded 5 times, and each segment has 12 sampling points.
Direct intermittent sampling interference generation mode: 10 segments are intercepted, each segment is forwarded 1 time, and 8 sampling points are adopted in each segment.
Repeating intermittent sampling interference generation modes: 10 segments are intercepted, each segment is forwarded for 4 times, and each segment has 4 sampling points.
The interference parameters have deviation:
c & I interference generation mode: 9 segments are intercepted, each segment is forwarded 6 times, and each segment has 12 sampling points.
Direct intermittent sampling interference generation mode: 10 segments are intercepted, each segment is forwarded 1 time, and 8 sampling points are adopted in each segment.
Repeating intermittent sampling interference generation modes: 10 segments are intercepted, each segment is forwarded for 4 times, and each segment has 4 sampling points.
Weight coefficient lambda 1 =1,λ 2 =1,λ 3 Main lobe width control parameter m=3, iterate 100 times.
Simulation analysis: fig. 6 shows the doppler tolerance of an optimized signal designed under the condition that the interference parameters are completely known (fig. 6 (a)) and the interference parameters deviate (fig. 6 (b)). When the speed offset is 150m/s, the main lobe of the optimized waveform designed under the condition that the interference parameters are completely known is about 23dB higher than the peak side lobe level, and the main lobe of the optimized waveform designed under the condition that the interference parameters have deviation is about 20dB higher than the peak side lobe level. It can be seen that the cognitive waveform generated by the method has better Doppler tolerance, and the interference parameter is completely known to be better than the Doppler tolerance of the optimized waveform designed under the condition that the interference parameter has deviation. Fig. 7 shows the result of the matched filtering processing of the LFM signal and its corresponding C & I interference, direct intermittent sampling interference and repeated intermittent sampling interference. Fig. 8 shows the result of mismatch filtering processing of the optimized signal and the corresponding multi-main lobe interference under the conditions that the interference parameters are completely known and the interference parameters have deviation, wherein (a) is that the interference parameters are completely known, and (b) is that the interference parameters have deviation. As can be seen from fig. 7 and 8, the peak sidelobe level of the optimized signal after the mismatch processing designed by the method is at least 12dB lower than the peak sidelobe level of the LFM signal after pulse compression, and the designed optimized signal detection performance is superior to that of the LFM signal; the multi-main-lobe interference resistance of the optimized signal is also superior to that of the LFM signal, but a certain signal-to-noise ratio loss exists in the design process, and the interference parameter is completely known to have superior detection performance and multi-main-lobe interference resistance compared with the interference parameter under the deviation condition. Fig. 9 is a comparison chart of detection results of the designed optimized signal when the LFM signal and the interference parameter are completely known and have deviation from the interference parameter, as shown in fig. 9 (a), (b), and (c), respectively. As can be seen from fig. 9, under the simulation parameters, the LFM cannot effectively detect the target, and the optimized signal can effectively detect the distance and speed of the target, and the detection performance of the optimized signal designed under the condition that the interference parameters are completely known is slightly better than that of the optimized signal designed under the condition that the interference parameters have deviation.
Those of ordinary skill in the art will recognize that the embodiments described herein are for the purpose of aiding the reader in understanding the principles of the present invention and should be understood that the scope of the invention is not limited to such specific statements and embodiments. Those of ordinary skill in the art can make various other specific modifications and combinations from the teachings of the present disclosure without departing from the spirit thereof, and such modifications and combinations remain within the scope of the present disclosure.

Claims (2)

1. A wave form and filter joint cognition design method for resisting multi-main lobe interference is characterized by comprising the following steps:
s1, constructing a C & I, direct intermittent sampling and repeated intermittent sampling interference signal model and a mismatch filter model based on a constant-mode phase coded signal; the method specifically comprises the following sub-steps:
s11, the constant-mode phase coding signal is expressed as:
Figure FDA0004153852830000011
in the formula [ (] T Representing the transpose s n For transmitting codewords of the waveform, n=1, 2, …, N s ,N s The number of codewords for the phase encoded waveform;
the nth element s n Expressed as:
Figure FDA0004153852830000012
in the method, in the process of the invention,
Figure FDA0004153852830000013
representation->
Figure FDA0004153852830000014
Phase of the nth codeword, +.>
Figure FDA0004153852830000015
Encoding the phase of the signal s for a constant mode phase;
s12, establishing C&I interference signal model: for C&I interference, suppose DRFM jammer pairIntercepting radar signals from s i Starting, m code words are taken as one section, each section is forwarded for l times, k sections are intercepted in total, and the whole C&I interference signal codeword length N s =klm;C&The I interference signal is:
Figure FDA0004153852830000016
wherein H is t T is more than or equal to 1 and less than or equal to k in the t section intercepted by the jammer; h t The represented codeword is expressed as:
Η t =[s (t-1)ml+i ,s (t-1)ml+i+1 ,…,s (t-1)ml+i+m-1 ];
s13, directly intermittently sampling an interference signal model: for direct intermittent sampling interference, suppose that a DRFM (digital radio frequency modulation) jammer intercepts radar signals, and s is a reference number i Firstly, m code words are taken as a section, and are intercepted and forwarded once every ml code words, k sections are intercepted in total, and the direct intermittent sampling interference signal is as follows:
J 2 =[Η′ 1 ,Η′ 0 ,Η′ 2 ,Η′ 0 ,…,Η′ k ,Η′ 0 ] T
in the formula, H' 0 =0 1×m(l-1) ,0 1×m(l-1) A 0 vector representing 1×m (l-1); h' t For the t-th segment intercepted by the jammer, the represented codeword is expressed as:
Η′ t =[s (t-1)ml+i ,s (t-1)ml+i+1 ,…s (t-1)ml+i+m-1 ];
s14, repeating the intermittent sampling interference signal model: for repeated intermittent sampling interference, assume that a DRFM (digital radio frequency modulation) jammer intercepts radar signals, and then s i Firstly, m code words are taken as a section, the code words are intercepted once every ml, the code words are repeated q times and then forwarded, k sections are intercepted in total, and the repeated intermittent sampling interference signals are as follows:
Figure FDA0004153852830000021
in the formula, H 0 =0 1×m(l-q) ,Η″ t For the t-th segment intercepted by the jammer, the represented codeword is expressed as:
Η″ t =[s (t-1)ml+i ,s (t-1)ml+i+1 ,…s (t-1)ml+i+m-1 ];
s15, a mismatch filter signal model: assume that the length of a corresponding mismatch filter h of a radar receiving end is N h ,N h ≥N s The mismatch filter h is denoted as:
Figure FDA0004153852830000022
/>
the nth element h n Expressed as:
Figure FDA0004153852830000023
where a is the magnitude vector of the mismatch filter h, a n Is the nth element in a; θ is the phase vector of the mismatch filter h, θ n Is the n element in theta; h=a.exp (j θ), where a represents the Hadamard product;
s2, under the condition that interference parameters are completely known or deviation exists, establishing a wave form correlation function peak value level, a correlation function main lobe template matching error and an interference energy weighting and minimizing criterion transmitting-receiving joint optimization problem under the constraint of signal to noise ratio; the method comprises the following steps:
s21, designing optimization criteria:
(1) Distance sidelobe: the output of the phase encoded signal s through the mismatch filter h is represented at a distance displacement j as:
Figure FDA0004153852830000024
wherein, represents conjugation, Γ 1 Represents a value interval, Γ 1 =[-(N s +N h )/2+1,(N s +N h )/2-1];
Let ω be j For the widened main lobe region, j=0, ±1, …, ±m, M is a main lobe width control parameter;
all the distance main lobe levels are stacked in vector omega main Expressed as:
ω main =[ω -M ,…,ω -101 ,…,ω M ] T
all range side lobe levels are stacked in vector ω side Expressed as:
Figure FDA0004153852830000025
let the sidelobe region be denoted as Γ 2 ,Γ 2 Expressed as:
Γ 2 =[-(N s +N h )/2+1,-M+1]∪[M-1,(N s +N h )/2-1]
the range side lobe is reduced by minimizing the peak side lobe, namely:
Figure FDA0004153852830000031
(2) Loss of signal-to-noise ratio:
|h H h-N s |≤η 1
0 -N s |≤η 2
in the formula, the constant eta 1 ≥0,η 2 More than or equal to 0, and the superscript H represents conjugate transposition;
(3) Main lobe control: assume that
Figure FDA0004153852830000032
Is the desired main lobe, i.e. q is a vector of dimension 2m+1; vector e is the error vector of the desired main lobe and the designed main lobe, expressed as:
e=ω main -q
the nth element e (n) of the error vector e is expressed as:
e(n)=ω n-M-1 -q(n),1≤n≤2M+1
the main lobe shape is maintained with a method that minimizes the maximum main lobe matching error, namely:
Figure FDA0004153852830000033
(4) Multi-main lobe interference resistance: interference signal J m The output result via the mismatch filter h is represented at a distance displacement j as:
Figure FDA0004153852830000034
/>
J m (n) represents J m In (2) the n-th element, all distance levels in the above formula are stacked in vector ω' m Expressed as:
Figure FDA0004153852830000035
when the interference parameters are completely known: the cognitive radar senses multi-main-lobe interference, and different interference signals J are transmitted under the condition that interference parameters are completely known m Distance level stack vector ω 'generated by mismatch filter h' m ,m=1,2,…,N J All stacked on vector omega J ,ω J Expressed as:
Figure FDA0004153852830000036
N J indicating the number of interferences;
when the interference parameter has deviation: actual C&The number of interception segments of the interference I is k, the forwarding frequency is l, the number of interception segments of the perceived parameter is k ', and the forwarding frequency is l'; with k=k '+l, l=l' -L, L being a natural number greater than 0; in the presence of errors, C&I interference distance level stack vector
Figure FDA0004153852830000037
Expressed as:
Figure FDA0004153852830000038
in the formula omega' 1 C is generated by intercepting segment number k' and forwarding frequency l&Stacking vectors of all distance levels after filtering processing of the I interference and mismatch filter;
Figure FDA0004153852830000041
c is generated by intercepting segment number k '+i and forwarding frequency l' -i&Stacking vectors of all distance levels after filtering processing of the I interference and mismatch filter; omega' 1-i C is generated by intercepting segment number k '-i and forwarding frequency l' +i&Stacking vectors of all distance levels after filtering processing of the I interference and mismatch filter, i=1, 2, …, L;
the actual direct intermittent sampling interference generation mode is as follows: taking m code words as a section, intercepting and forwarding every ml code words, intercepting k sections altogether, and the perceived generation mode of direct intermittent sampling interference is as follows; taking m code words as a section, intercepting and forwarding every ml code words, and intercepting k sections altogether; direct intermittent sampling of interference distance level stacking vectors in the presence of deviations
Figure FDA0004153852830000042
Expressed as:
Figure FDA0004153852830000043
ω′ 2 taking m code words as a section, intercepting and forwarding every ml code words, and intercepting direct intermittent sampling interference generated by k section and stacking all distance levels after filtering treatment of a mismatch filterAn amount of;
Figure FDA0004153852830000044
taking m code words as a section, intercepting and forwarding every m (l '-i) code words, and intercepting direct intermittent sampling interference generated by k' +i sections and stacking vectors of all distance levels after filtering treatment of a mismatch filter; />
Figure FDA0004153852830000045
Taking m code words as a section, intercepting and forwarding every m (l '+i) code words, and intercepting the stacking vectors of direct intermittent sampling interference generated by the section k' -i and all distance levels after filtering treatment of a mismatch filter;
the actual repeated intermittent sampling interference generation mode is as follows: taking m code words as a section, intercepting and forwarding q times every ml code words, and intercepting k sections altogether; the perceived repeated intermittent sampling interference generation mode is as follows; taking m code words as a section, intercepting and forwarding q times every ml 'code words, and intercepting k' sections altogether; repeated intermittent sampling of interference distance level stacking vectors in the presence of deviations
Figure FDA0004153852830000046
Expressed as:
Figure FDA0004153852830000047
ω′ 3 taking m code words as a section, intercepting and forwarding q times every ml 'code words, and intercepting the stacking vectors of the repeated intermittent sampling interference generated by the k' section and all distance levels after filtering treatment of a mismatch filter altogether;
Figure FDA0004153852830000048
taking m code words as a section, intercepting and forwarding q times every m (l '-i) code words, and intercepting repeated intermittent sampling interference generated by k' +i sections and stacking vectors of all distance levels after filtering treatment of a mismatch filter;/>
Figure FDA0004153852830000049
taking m code words as a section, intercepting and forwarding q times every m (l '+i) code words, and intercepting the repeated intermittent sampling interference generated by the k' -i section and the stacking vectors of all distance levels after the filtering treatment of a mismatch filter;
all distance levels ω of all interference output via mismatched filters J The method comprises the following steps: if there is no deviation in several disturbance parameters, all distance levels omega J Is that
Figure FDA00041538528300000410
If there is a deviation in one or more interference parameters, then
Figure FDA00041538528300000411
The corresponding vector in (a) is replaced by a distance level stack vector when errors exist;
all distance levels of all interference output via the mismatch filter are compressed using a min-max optimization criterion, namely:
Figure FDA0004153852830000051
s22, establishing an optimization problem: according to the optimization criterion of the S21 design, the optimization problem of the phase coding signal S and the mismatch filter h resisting the multi-main lobe interference is expressed as follows:
Figure FDA0004153852830000052
s.t.|h H h-N s |≤γ 1
wherein lambda is 1 ,λ 2 The weight coefficient is preset;
s3, converting the optimization problem into a minimization problem related to the independent variable function, and solving by using an iterative L-BFGS algorithm.
2. The method for combined cognitive design of a waveform and a filter for resisting multi-main lobe interference according to claim 1, wherein the specific implementation method of S3 is as follows:
s31, optimizing problem conversion: converting the optimization problem into a minimization problem with respect to the independent variable function, namely:
Figure FDA0004153852830000053
wherein lambda is 3 For the weight coefficient to be set in advance, I Represents an infinite norm;
defining an objective function as:
f(x)=||ω side || p1 ||ω J || p2 ||e|| p3 ||h H h-N s || p
in the method, in the process of the invention, I p Represents the p-norm; vector quantity
Figure FDA0004153852830000054
Is composed of vector->
Figure FDA0004153852830000055
a and θ are sequentially formed column vectors expressed as:
Figure FDA0004153852830000056
s32, solving an optimization problem: and solving the problem by using an L-BFGS algorithm based on iteration, solving the minimum value of the objective function by using the L-BFGS algorithm, continuously iterating until the objective function f (x) reaches the minimum descent epsilon, stopping iterating, and outputting x.
CN202210379049.7A 2022-04-12 2022-04-12 Wave form and filter combined cognition design method for resisting multi-main lobe interference Active CN114696788B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210379049.7A CN114696788B (en) 2022-04-12 2022-04-12 Wave form and filter combined cognition design method for resisting multi-main lobe interference

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210379049.7A CN114696788B (en) 2022-04-12 2022-04-12 Wave form and filter combined cognition design method for resisting multi-main lobe interference

Publications (2)

Publication Number Publication Date
CN114696788A CN114696788A (en) 2022-07-01
CN114696788B true CN114696788B (en) 2023-05-05

Family

ID=82143701

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210379049.7A Active CN114696788B (en) 2022-04-12 2022-04-12 Wave form and filter combined cognition design method for resisting multi-main lobe interference

Country Status (1)

Country Link
CN (1) CN114696788B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116224246B (en) * 2022-12-14 2023-11-17 南京航空航天大学 Suppression interference suppression method for joint design of transmitting waveform and receiving filter

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106093877A (en) * 2016-07-19 2016-11-09 西安电子科技大学 Orthogonal wide main lobe phase coding signal and mismatched filter combined optimization method
CN106992813A (en) * 2017-02-16 2017-07-28 广东顺德中山大学卡内基梅隆大学国际联合研究院 A kind of inexpensive adaptive wireless information and energy multiplex transmission system
CN112464561A (en) * 2020-11-25 2021-03-09 南京航空航天大学 MIMO radar waveform design method based on improved Jaya algorithm
CN112540351A (en) * 2020-12-28 2021-03-23 哈尔滨工业大学(威海) Radar anti-shielding interference waveform design method
CN112633514A (en) * 2020-10-26 2021-04-09 华南师范大学 Multi-task function-to-function regression method

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020152250A1 (en) * 1999-02-24 2002-10-17 Robert B. Staszewski Finite-impulse-response (fir) filter employing a parallel architecture
US11321504B2 (en) * 2018-05-09 2022-05-03 Palo Alto Research Center Incorporated Learning constitutive equations of physical components with constraints discovery
WO2020037171A1 (en) * 2018-08-16 2020-02-20 Movano Inc. Calibration, classification and localization using channel templates

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106093877A (en) * 2016-07-19 2016-11-09 西安电子科技大学 Orthogonal wide main lobe phase coding signal and mismatched filter combined optimization method
CN106992813A (en) * 2017-02-16 2017-07-28 广东顺德中山大学卡内基梅隆大学国际联合研究院 A kind of inexpensive adaptive wireless information and energy multiplex transmission system
CN112633514A (en) * 2020-10-26 2021-04-09 华南师范大学 Multi-task function-to-function regression method
CN112464561A (en) * 2020-11-25 2021-03-09 南京航空航天大学 MIMO radar waveform design method based on improved Jaya algorithm
CN112540351A (en) * 2020-12-28 2021-03-23 哈尔滨工业大学(威海) Radar anti-shielding interference waveform design method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
一种相位编码信号和失配滤波器联合优化方法;徐磊磊;周生华;刘宏伟;马林;;西安电子科技大学学报(第05期);25-30 *
一种稳健的多主瓣干扰抑制方法;陈卓等;胸痛工程与电子技术;544-549 *
面向度主瓣干扰的波形与滤波器联合认知设计方法;林瑜等;系统工程与电子技术;1-15 *

Also Published As

Publication number Publication date
CN114696788A (en) 2022-07-01

Similar Documents

Publication Publication Date Title
CN109116337B (en) Sparse approximate minimum variance orientation estimation method based on matrix filtering
CN106662645B (en) Adaptive selection STAP (static adaptive station) method for beam-Doppler channel
CN114696788B (en) Wave form and filter combined cognition design method for resisting multi-main lobe interference
CN113239554B (en) Waveform optimization design method for resisting sidelobe shielding interference
CN112034434A (en) Radar radiation source identification method based on sparse time-frequency detection convolutional neural network
CN113376607B (en) Airborne distributed radar small sample space-time self-adaptive processing method
CN107102305A (en) Sane cognitive radar emission receives co-design method under a kind of clutter environment
CN110632573A (en) Airborne broadband radar space-time two-dimensional keystone transformation method
CN105223553B (en) Shift-frequency jamming recognition methods is realized in a kind of half frequency range matched filtering
CN107635181B (en) Multi-address sensing source feedback optimization method based on channel learning
CN114384477A (en) Detection and interference integrated waveform generation method based on intermittent sampling
CN108896963B (en) Airborne radar space-time self-adaptive dimension reduction processing method
CN114488055A (en) Coherent processing method suitable for carrier frequency and repetition frequency combined agile radar
CN108152796A (en) A kind of main lobe based on grey-kalman filtering moves interference elimination method
CN116953631A (en) Full-polarization radar emission waveform optimization method, computer device and storage medium
CN113640754B (en) Main lobe intermittent sampling forwarding interference resisting method based on low-rank sparse model
CN111025258B (en) Joint mismatch filter for radar waveform diversity and design method thereof
Li et al. A Radar Waveform Design of MCPC Method for Interrupted Sampling Repeater Jamming Suppression via Fractional Fourier Transform.
Oliinyk et al. Center weighted median filter application to time delay estimation in non-Gaussian noise environment
CN108008360B (en) Amplitude-weighted nonlinear frequency modulation waveform design method
CN108132460B (en) Pulse compression compensation algorithm based on frequency domain channel equalization
Li et al. Shared Waveform Design for Integrated Detection and Jamming Signal Based on Smart Noise Jamming of Convolution Modulation
CN114594425B (en) Short-time pulse train waveform design method for resisting clutter interference
CN110196412B (en) STAP method combining sparsity
CN113030880B (en) Signal estimation method of intermittent sampling forwarding type interference based on ADMM

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