Comprehensive cross-ambiguity function loop optimization method based on determined information

Abstract

The invention discloses a comprehensive cross-ambiguity function loop optimization method based on determined information. The method comprises the following steps: step 1, designing a discrete template of the determined information of a cross-ambiguity function; 2, determining a zero-sum game model of a radar party and an interference party; 3, improving a zero-sum game model based on the template model, and constructing a model independent of the template; and 4, degenerating the code elements into pulses, constructing a special model which is feasible in calculation, and designing a radar waveform according to the special model. According to the method, the optimal transmitted waveform and the matched waveform can be efficiently obtained through a reverse iteration method, so that the anti-interference capability of a radar system is improved, and the problem of calculation efficiency of attack and defense strategy selection in a large-scale network is effectively solved.

Description

Technical Field

[0001] The invention belongs to the technical field of radar waveform design, and in particular relates to a comprehensive mutual fuzzy function loop optimization method based on certain information. Background Art

[0002] In modern radar systems, waveform design is one of the key technologies related to its performance and effectiveness. Radar systems detect and identify targets by emitting specific electromagnetic waveforms, while dealing with various complex electromagnetic environments and hostile interference. Traditional waveform design methods often rely on fixed waveform patterns or empirical designs. Although these designs meet basic detection needs to a certain extent, their performance and effectiveness are significantly limited when faced with changing electromagnetic environments and strong electronic warfare interference.

[0003] The development of modern radar technology has driven the continuous evolution and optimization of waveform design methods. Based on advanced signal processing technology and mathematical modeling, the new generation of radar waveform design aims to improve the performance of radar systems in target detection, tracking and anti-interference capabilities by optimizing the energy distribution, spectrum characteristics and spatiotemporal characteristics of the waveform.

[0004] The background technology of waveform design covers many aspects: first, a comprehensive understanding and analysis of the radar working environment and usage scenarios. Different mission requirements may pose different challenges to the waveform design of the radar system, such as high-resolution target recognition, low-probability interception, and effective signal transmission in complex electromagnetic environments. Secondly, waveform design also needs to take into account the actual characteristics and limitations of radar hardware, such as the frequency response, power output capability, and signal processing speed of the transmitter and receiver.

[0005] In traditional waveform design, researchers usually rely on mathematical models and simulation tools to evaluate the performance of different waveform schemes. These models involve the propagation characteristics of radar signals, target feature extraction methods, and interference suppression techniques. With the improvement of computing power and algorithm optimization, modern waveform design increasingly adopts methods based on statistical learning and optimization theory, and maximizes the signal processing efficiency in dynamic scenarios through big data analysis and adaptive algorithms.

[0006] As an important part of modern radar system performance optimization, waveform design is constantly driven and influenced by technological innovation and scientific research. Future development directions include more intelligent waveform adaptive design, integrated multi-sensor collaboration, and adaptability to extreme conditions in complex electromagnetic environments. These technological advances will further enhance the wide application and actual effectiveness of radar systems in military, civilian and scientific research fields. Summary of the invention

[0007] In view of this, a comprehensive mutual ambiguity function loop optimization method based on deterministic information is disclosed. Based on the existing target and interference signal models, a series of waveform templates or optimization parameters are designed, so that the radar system can extract target features to the greatest extent when receiving target signals, and at the same time can suppress interference to the greatest extent when facing interference signals, thereby ensuring the stability and performance of the system.

[0008] The object of the present invention is achieved through the following technical solution, which is a comprehensive mutual fuzzy function loop optimization method based on determined information, comprising the following steps:

[0009] Step 1, designing a discrete template of mutual fuzzy functions to determine information;

[0010] Step 2: Determine the zero-sum game model between the radar party and the jammer;

[0011] Step 3, improve the zero-sum game model based on the template model and build a model that does not rely on the template;

[0012] Step 4, degenerate the code element into pulses, construct a computationally feasible special model, and design the radar waveform according to the special model;

[0013] The determined information at least includes the time delay t for the radar to accurately know the target target 、Doppler frequency shift f of the target d,target , the jammer's jamming pattern, the jammer's delay modulation based on the target delay t jammer , the Doppler frequency shift modulation of the jammer based on the target Doppler frequency shift f d,jammer one.

[0014] Preferably, the mutual fuzzy function is defined as:

[0015]

[0016] Among them, χ cross,S,W (τ, f) is the mutual ambiguity function, which is a function of time delay τ and frequency offset f, where τ represents the time delay parameter and f represents the frequency offset parameter. Indicates the double integral symbol. S(t) and W(t) are two different signals. t is time, S * (t) represents the complex conjugate of signal S(t), W(t) is another signal that performs mutual fuzzy operation with signal S(t), It is a complex exponential function containing frequency offset f and time delay τ, which is used to introduce the corresponding phase shift, dt differential integration, which means integrating the signal in time;

[0017] The mutual ambiguity function has the following properties:

[0018]

[0019] | cross,S,W (-τ,-f)|=|χ cross,S,W (τ,f)|;

[0020]

[0021] Among them, |χ cross,S,W (0,0)| represents the amplitude of the mutual ambiguity function at (τ,f)=(0,0), represents the square root of the energy of the signal S, represents the square root of the energy of signal W, | cross，S，W (τ, f)| and |χ cross,S,W (-τ,-f)| represents the magnitude of the mutual ambiguity function at (τ,f) and (-τ,-f), respectively. represents the double integral symbol, indicating that the integration is performed over the entire range of time delay τ and frequency offset f, |χ cross,S,W (τ,f)| 2 represents the square of the modulus of the mutual ambiguity function, represents the energy density, dτ represents the infinitesimal time delay τ, and df represents the infinitesimal frequency offset f, indicating that the mutual ambiguity function has symmetry about (τ, f).

[0022] Preferably, the zero-sum game model between the radar party and the interference party is:

[0023] Radar said:

[0024]

[0025] The interfering party is represented as follows:

[0026]

[0027] Among them, Q target and Q jammer Represents the characteristic set of target and interference source respectively, g k,p is a discrete template of the mutual ambiguity function, which is used to describe the characteristics of the signal in the radar system, where k and p are the discretized time delay and frequency offset indices, and d auto∨cross (τ,f) is the error between the discrete template of the mutual fuzzy function and its expected value, χ auto∨cross,S,W (τ, f) are the self-ambiguity and mutual-ambiguity functions used to describe the correlation between signals S and W under time delay τ and frequency offset f. and They represent the local signal-to-noise ratio and interference-to-noise ratio respectively, and their dynamic changes are determined by the internal state of the system and the external control input u radar (t) and v jammer (t) determined by and It is a function that describes the dynamic evolution and is used to calculate the rate of change of the local signal-to-noise ratio and interference-to-noise ratio;

[0028] The zero-sum game model is established as follows:

[0029]

[0030] Among them, min X,Y It indicates that the variables X and Y are minimized. represents the sum of the index set T×F of time delay and frequency offset, g kp Represents the discrete template of the mutual ambiguity function, describing the characteristics of the signal under time delay k and frequency offset p, represents the modulus of the expected value of the mutual fuzzy function template, X H represents the conjugate transpose of X, J kp represents the complex number representation of the mutual ambiguity function template, Y represents the received signal processing vector in the optimization model;

[0031] Equivalent to

[0032]

[0033] in, Represents the variables X, Y and phase angle φ k,p To minimize the operation, Represents a complex exponential function, introducing a phase angle φ kp The phase shift, φ k,p Represents the phase angle, which is used to describe the phase shift of a signal.

[0034] Preferably, the standard model form of the zero-sum game model is:

[0035]

[0036] Among them, min means minimizing the objective function, represents the sum of the index set T×F of time delay and frequency offset, |X k | represents the modulus of the kth element in X,

[0037] m∈[1,2,…,M],n∈[1,2,…,N] means to constrain all indexes k that meet the conditions, and st means the constraint conditions. represents the square root of MN, φ kp ∈[0,2π) represents the phase angle φ kp The value range constraints.

[0038] Preferably, the simplified model of the zero-sum game model is as follows:

[0039] Determine the symbol vector of the first pulse

[0040]

[0041] The following M-1 pulse symbol vectors are derived from the symbol vector of the first pulse.

[0042]

[0043] Then at this time

[0044]

[0045] Get the simplified model

[0046]

[0047] Among them, min means that the goal is to minimize the subsequent expression, g kp Represents the discrete template of the mutual ambiguity function, describing the characteristics of the signal under time delay k and frequency offset p, Represents a complex exponential function, introducing a phase angle φ kp The phase shift, φ kp Represents the phase angle, which is used to describe the phase shift of the signal. H represents the conjugate transpose of X, J kp represents the complex representation of the mutual ambiguity function template, Y represents the received signal processing vector in the optimization model, st represents the constraint condition, |X 1,n |=1 indicates a constraint condition, ensuring that the first element X in X 1,n The modulus is 1, and for all n, |ρ m |=1 indicates a constraint condition, ensuring ρ m The modulus is 1, for all m, X m,n =ρ m X 1,n Represents constraints, describing how the subsequent parts of X are constructed, where X m,n By m Multiply by X 1,n Get, X m,N+l =0 indicates a constraint condition, ensuring that some elements in X are zero. Represents a constraint condition, ensuring that the modulus of Y is φ kp ∈[0,2π) represents the constraint condition, limiting the phase angle φ kp The value range of .

[0048] Preferably, the template-independent model is:

[0049]

[0050] Among them, min means that the goal is to minimize the subsequent expression, X H represents the conjugate transpose of X, J kp The complex form describing the template of the mutual ambiguity function is related to the time delay k and the frequency offset p. Y is the received signal processing vector in the optimization model. st represents the constraint condition. |X k |=1 constraint, ensuring that the modulus of the kth row in X is 1, where k=(m,n), m represents the pulse index, n represents the symbol index, and X k = 0 constraint, ensuring that the elements at specific positions in X are zero, where k = (m, N + l), m represents the pulse index, l represents the additional index offset, Constraints ensure that the modulus of Y is Where M and N represent the size of the time and frequency dimensions, respectively, kp ∈[0,2π) constraint, limiting the phase angle φ kp The value range of , where (k, p) represents the combination of time delay and frequency offset;

[0051] The model is written in standard model form as:

[0052]

[0053] Preferably, in the computationally feasible special model, when N+L=1, the code element degenerates into a pulse, and one pulse repetition period is one code element. At this time, the phase coding degenerates into the initial phase modulation of the pulse:

[0054] X=(X1,X2,…,X m ,…,X M ) T is the initial phase modulation of M pulses, and the model becomes

[0055]

[0056] At this time, the index set is T = [-M, ..., 0, ..., M], min means that the goal is to minimize the subsequent expression, g kp Represents the discrete template of the mutual ambiguity function, describing the characteristics of the signal under time delay k and frequency offset p, Represents a complex exponential function, introducing a phase angle φ kp The phase shift, φ kp Represents the phase angle, which is used to describe the phase shift of the signal. H represents the conjugate transpose of X, J kprepresents the complex representation of the mutual ambiguity function template, Y represents the received signal processing vector in the optimization model, st represents the constraint condition, |X m |=1 represents the constraint condition, ensuring that the modulus of the mth row in X is 1. For all m, Represents a constraint condition, ensuring that the modulus of Y is φ kp ∈[0,2π) represents the constraint condition, limiting the phase angle φ kp The value range of .

[0057] Compared with the existing methods, the advantage of the method of the present invention is that network game is a research hotspot in recent years. However, traditional game solving algorithms are extremely difficult to solve in large-scale networks. The present invention provides a network game solving method based on a greedy algorithm. This method gradually approaches the global optimal solution through local optimal selection, effectively solving the computational efficiency problem of attack and defense strategy selection in large-scale networks. BRIEF DESCRIPTION OF THE DRAWINGS

[0058] Figure 1 A schematic diagram of a process flow of an embodiment of the present invention is shown;

[0059] Figure 2 The range Doppler diagram of the target echo and the interference before optimization in the embodiment of the present invention is shown;

[0060] Figure 3 A 2D mutual ambiguity function diagram of a transmission waveform before optimization according to an embodiment of the present invention is shown;

[0061] Figure 4 The optimized range Doppler diagram of the target echo and interference according to the embodiment of the present invention is shown. DETAILED DESCRIPTION

[0062] The present invention is further described below in conjunction with the accompanying drawings, but the present invention is not limited in any way. Any changes or substitutions made based on the teachings of the present invention belong to the protection scope of the present invention.

[0063] In this embodiment, linear frequency modulation is adopted to add noise and interference to the waveform, and the effectiveness of the method is verified by comparing the signal-to-interference-noise ratio before and after optimization.

[0064] like Figure 1 As shown, a comprehensive mutual fuzzy function loop optimization method based on determined information includes expressing

[0065] Step 1, designing a discrete template of mutual fuzzy functions to determine information;

[0066] Step 2: Determine the zero-sum game model between the radar party and the jammer;

[0067] Step 3, improve the zero-sum game model based on the template model, and build and design a model that does not rely on the template;

[0068] Step 4: degenerate the code elements into pulses and construct a computationally feasible special model.

[0069] The above-mentioned information is an ideal model, which means that the radar knows the target's delay t accurately. target 、Doppler frequency shift f of the target d,target , the jammer's jamming pattern (assumed to be forwarding deception jamming), the jammer's delay modulation based on the target delay t jammer , the Doppler frequency shift modulation of the jammer based on the target Doppler frequency shift f d,jammer And other information.

[0070] Specifically, the mutual fuzzy function is defined as

[0071]

[0072] Among them, χ cross,S,W (τ, f) is the mutual ambiguity function, which is a function of time delay (τ) and frequency offset (f), where τ represents the time delay parameter and f represents the frequency offset parameter. represents the double integral symbol, S(t) and W(t) are two signals, S * (t) represents the complex conjugate of signal S(t), W(t) is another signal that performs mutual fuzzy operation with signal S(t), It is a complex exponential function containing frequency offset f and time delay τ, which is used to introduce the corresponding phase shift, dt differential integration, which means integrating the signal over time.

[0073] It has properties similar to self-fuzzifying functions:

[0074]

[0075] | cross,S,W (-τ,-f)|=|χ cross,S,W (τ,f)|;

[0076]

[0077] Among them, |χ cross,S,W (0,0)| represents the amplitude of the mutual ambiguity function at (τ,f)=(0,0), represents the square root of the energy of the signal S, represents the square root of the energy of the signal W,

[0078] | cross,S,W (-τ,-f)| represents the magnitude of the mutual ambiguity function at (τ,f) and (-τ,-f), represents the double integral symbol, indicating that the integration is performed over the entire range of time delay τ and frequency offset f, |χ cross,S,W (τ,f)| 2 represents the square of the modulus of the mutual ambiguity function, represents the energy density, dτ represents the infinitesimal time delay τ, and df represents the infinitesimal frequency offset f, indicating that the mutual ambiguity function has symmetry about (τ, f).

[0079] Furthermore, the zero-sum game model between the radar party and the jammer is as follows:

[0080] The so-called deterministic information mode refers to the radar's target ,Q jammer is certain, so the discrete template g of the designed mutual fuzzy function k,p It is relatively accurate, then we can use the zero-sum game model of deterministic information mode + mismatch filter + self-fuzzy function

[0081] Radar said

[0082]

[0083] The interfering party said

[0084]

[0085] Among them, Q target and Q jammer Represents the characteristic set of target and interference source respectively, g k,p is a discrete template of the mutual ambiguity function, which is used to describe the characteristics of the signal in the radar system, where k and p are the discretized time delay and frequency offset indices, and d auto∨cross (τ,f) is the error between the discrete template of the mutual fuzzy function and its expected value, χ auto∨cross,S,W (τ, f) are the self-ambiguity and mutual-ambiguity functions used to describe the correlation between signals S and W under time delay τ and frequency offset f. and They represent the local signal-to-noise ratio and interference-to-noise ratio respectively, and their dynamic changes are determined by the internal state of the system and the external control input u radar (t) and v jammer (t) determined by and It is a function that describes the dynamic evolution and is used to calculate the rate of change of the local signal-to-noise ratio and interference-to-noise ratio;

[0086] Can build models

[0087]

[0088] Among them, min X,YIt indicates that the variables X and Y are minimized.

[0089] (Tex translation failed) means summing the index set T×F of time delay and frequency offset, g kp Represents the discrete template of the mutual ambiguity function, describing the characteristics of the signal under time delay k and frequency offset p, represents the modulus of the expected value of the mutual fuzzy function template, X H represents the conjugate transpose of X, J kp represents the complex number representation of the mutual ambiguity function template, Y represents the received signal processing vector in the optimization model;

[0090] Equivalent to

[0091]

[0092] in, Represents the variables X, Y and phase angle φ k,p To minimize the operation, Represents a complex exponential function, introducing a phase angle φ kp The phase shift, φ k,p Represents the phase angle, which is used to describe the phase shift of the signal;

[0093] Written in standard model form, we get

[0094]

[0095] Among them, min means minimizing the objective function, represents the sum of the index set T×F of time delay and frequency offset, |X k | represents the modulus of the kth element in X,

[0096] m∈[1,2,…,M],n∈[1,2,…,N] means to constrain all indexes k that meet the conditions, and st means the constraint conditions. represents the square root of MN, φ kp ∈[0,2π) represents the phase angle φ kp The value range constraints of

[0097] Many simplified models can be established. For example, we determine the code element vector of the first pulse

[0098]

[0099] The following M-1 pulse symbol vectors are derived from the symbol vector of the first pulse.

[0100]

[0101] Then at this time

[0102]

[0103] A simplified model can be obtained

[0104]

[0105] Among them, min means that the goal is to minimize the subsequent expression, g kp Represents the discrete template of the mutual ambiguity function, describing the characteristics of the signal under time delay k and frequency offset p, Represents a complex exponential function, introducing a phase angle φ kp The phase shift, φ kp Represents the phase angle, which is used to describe the phase shift of the signal. H represents the conjugate transpose of X, J kp represents the complex representation of the mutual ambiguity function template, Y represents the received signal processing vector in the optimization model, st represents the constraint condition, |X 1,n |=1 indicates a constraint condition, ensuring that the first element X in X 1,n The modulus is 1, and for all n, |ρ m |=1 indicates a constraint condition, ensuring ρ m The modulus is 1, for all m, X m,n =ρ m X 1,n Represents constraints, describing how the subsequent parts of X are constructed, where X m,n By m Multiply by X 1,n Get, X m,N+l =0 indicates a constraint condition, ensuring that some elements in X are zero. Represents a constraint condition, ensuring that the modulus of Y is φ kp ∈[0,2π) represents the constraint condition, limiting the phase angle φ kp The value range of .

[0106] In one embodiment, the code element vectors of the first and second pulses of the transmitted waveform are further optimized. Other pulse code element vectors are based on the code elements of the first and second pulses through the overall initial phase modulation. Through this idea, a large number of sub-models are established, and a special model that is computationally feasible can also be constructed. That is, when N+L=1, the code element degenerates into a pulse, and one pulse repetition period is one code element. At this time, the phase coding degenerates into the initial phase modulation of the pulse.

[0107] X=(X1,X2,…,X m ,…,XM ) T

[0108] is the initial phase modulation of M pulses, and the model becomes

[0109]

[0110] At this time, the index set is T = [-M, ..., 0, ..., M], min means that the goal is to minimize the subsequent expression, g kp Represents the discrete template of the mutual ambiguity function, describing the characteristics of the signal under time delay k and frequency offset p, Represents a complex exponential function, introducing a phase angle φ kp The phase shift, φ kp Represents the phase angle, which is used to describe the phase shift of the signal. H represents the conjugate transpose of X, J kp represents the complex representation of the mutual ambiguity function template, Y represents the received signal processing vector in the optimization model, st represents the constraint condition, |X m |=1 represents the constraint condition, ensuring that the modulus of the mth row in X is 1. For all m, Represents a constraint condition, ensuring that the modulus of Y is φ kp ∈[0,2π) represents the constraint condition, limiting the phase angle φ kp The value range of .

[0111] Furthermore, the template-independent model is represented by

[0112] The design of the template requires considerable accuracy to achieve good results, which requires the help of intelligent technology. We can improve the above model based on the template model and build and design a model that does not rely on the template.

[0113]

[0114] Among them, min means that the goal is to minimize the subsequent expression, X H represents the conjugate transpose of X, J kp The complex form describing the template of the mutual ambiguity function is related to the time delay k and the frequency offset p. Y is the received signal processing vector in the optimization model. st represents the constraint condition. |X k |=1 constraint, ensuring that the modulus of the kth row in X is 1, where k=(m,n), m represents the pulse index, n represents the symbol index, and X k = 0 constraint, ensuring that the elements at specific positions in X are zero, where k = (m, N + l), m represents the pulse index, l represents the additional index offset, Constraints ensure that the modulus of Y is Where M and N represent the size of the time and frequency dimensions, respectively, kp ∈[0,2π) constraint, limiting the phase angle φ kp The value range of , where (k, p) represents the combination of time delay and frequency offset;

[0115] It must be noted that the indicator set here

[0116]

[0117] is a proper subset of the original index, because if the proper subset is not selected, then according to the law of energy conservation, there is

[0118]

[0119] Then the optimization loses its meaning.

[0120] Model

[0121]

[0122] Written in standard model form, we get

[0123]

[0124] Many simplified models can be established. For example, we determine the code element vector of the first pulse

[0125]

[0126] The following M-1 pulse symbol vectors are derived from the symbol vector of the first pulse.

[0127]

[0128] Then at this time

[0129]

[0130] A simplified model can be obtained

[0131]

[0132] Among them, min means that the goal is to minimize the subsequent expression, X H represents the conjugate transpose of X, J kp The complex form describing the template of the mutual ambiguity function is related to the time delay k and the frequency offset p. Y is the received signal processing vector in the optimization model. st represents the constraint condition. |X 1,n |=1 constraint, ensuring that the first element X in X 1,nThe modulus is 1, and for all n, |ρ m |=1 constraint, ensuring ρ m The modulus is 1, for all m, ρ m is the complex scaling factor of the symbol vector, X m,n =ρ m X 1,n Constraints that describe how the subsequent parts of X are constructed, where X m,n By m Multiply by X 1,n Get, X m,N+l = 0 constraint, ensuring that some elements in X are zero, Constraints ensure that the modulus of Y is Where M and N represent the size of time and frequency dimensions respectively;

[0133] We can also optimize the code element vectors of the first and second pulses. The other pulse code element vectors are based on the code elements of the first and second pulses through the overall initial phase modulation. With this idea, we can establish a large number of sub-models.

[0134] Furthermore, the computationally feasible special model refers to the representation

[0135] When N+L=1, the code element degenerates into a pulse, and one pulse repetition period is one code element. At this time, the phase coding degenerates into the initial phase modulation of the pulse X=(X1,X2,…,X m ,…,X M ) T is the initial phase modulation of M pulses, and the model becomes

[0136]

[0137] Among them, min means that the goal is to minimize the subsequent expression, X H represents the conjugate transpose of X, J kp The complex form describing the template of the mutual ambiguity function is related to the time delay k and the frequency offset p. Y is the received signal processing vector in the optimization model. st represents the constraint condition. |X m |=1 constraint, ensuring that the modulus of the mth row in X is 1, where m represents the pulse index, Constraints ensure that the modulus of Y is Where M represents the number of pulses.

[0138] At this time, the indicator set is

[0139]

[0140] These sets of metrics describe the range of time delays and frequency offsets, ensuring that the model operates effectively under the specific conditions encountered in radar signal processing.

[0141] According to the above content, the range Doppler image before optimization in this embodiment is as follows: Figure 2 As shown, the mutual fuzzy function diagram before optimization is as follows Figure 3 As shown. To optimize the waveform, first, the radar initial waveform, filter and other necessary parameters are given. The two-dimensional color map and three-dimensional function map of the fuzzy function of the initial waveform and filter are calculated. According to the model of the target and interference, a range Doppler image is formed, and the local signal-to-interference-noise ratio of the target unit is calculated. The initial waveform and filter of the radar are modulated using the fuzzy function-based model and algorithm in the previous article, and the waveform and filter are output. The two-dimensional color map and three-dimensional function map of the fuzzy function of the modulated waveform and filter are calculated. Again, according to the signal model of the target and interference, a range Doppler image is formed, and the local signal-to-interference-noise ratio of the target unit is calculated. Compare the shapes of the fuzzy functions before and after, and compare the clarity and local signal-to-interference-noise ratio of the target units before and after. The optimized range Doppler map is shown as follows. Figure 4 shown.

[0142] Traditional game solving algorithms are extremely difficult to solve in large-scale networks. Compared with existing methods, the present invention can efficiently obtain the optimal transmission waveform and matching waveform through a reverse iteration method, thereby improving the anti-interference ability of the radar system.

[0143] As used herein, the word "preferred" is intended to be used as an example, instance, or illustration. Any aspect or design described herein as "preferred" is not necessarily to be construed as being more advantageous than other aspects or designs. On the contrary, the use of the word "preferred" is intended to present concepts in a specific way. The term "or" as used in this application is intended to mean an inclusive "or" rather than an exclusive "or". That is, unless otherwise specified or clear from the context, "X uses A or B" means any one of the naturally included permutations. That is, if X uses A; X uses B; or X uses both A and B, then "X uses A or B" is satisfied in any of the foregoing examples.

[0144] Moreover, although the present disclosure has been shown and described with respect to one or implementations, those skilled in the art will think of equivalent variations and modifications based on the reading and understanding of this specification and the accompanying drawings. The present disclosure includes all such modifications and variations, and is limited only by the scope of the appended claims. In particular, with respect to the various functions performed by the above-mentioned components (such as elements, etc.), the terms used to describe such components are intended to correspond to any component (unless otherwise indicated) that performs the specified function of the component (such as it is functionally equivalent), even if the structure is not equivalent to the disclosed structure of the function in the exemplary implementation of the present disclosure shown herein. In addition, although the specific features of the present disclosure have been disclosed with respect to only one of several implementations, such features can be combined with one or other features of other implementations that may be desired and advantageous for a given or specific application. Moreover, insofar as the terms "including", "having", "containing" or their variations are used in specific embodiments or claims, such terms are intended to be included in a manner similar to the term "comprising".

[0145] To sum up, the above embodiment is an implementation mode of the present invention, but the implementation mode of the present invention is not limited by the embodiment. Any other changes, modifications, substitutions, combinations, and simplifications that deviate from the spirit and principles of the present invention should be equivalent replacement methods and are included in the protection scope of the present invention.

Claims

1. A comprehensive mutual fuzzy function loop optimization method based on certain information, characterized in that: The following steps are involved: Step 1, designing a discrete template of mutual fuzzy functions to determine information; Step 2: Determine the zero-sum game model between the radar party and the jammer; Step 3, improve the zero-sum game model based on the template model and build a model that does not rely on the template; Step 4, degenerate the code element into pulses, construct a computationally feasible special model, and design the radar waveform according to the special model; Determine the information indicating that the radar knows the target's delay t accurately target 、Doppler frequency shift f of the target d,target , the jammer's jamming pattern, the jammer's delay modulation based on the target delay t jammer , the Doppler frequency shift modulation of the jammer based on the target Doppler frequency shift f d,jammer one.

2. The method for optimizing the comprehensive mutual fuzzy function loop based on the determined information according to claim 1, characterized in that: The mutual ambiguity function is defined as: Among them, χ cross,S,W (τ, f) is the mutual ambiguity function, which is a function of time delay τ and frequency offset f, where τ represents the time delay parameter and f represents the frequency offset parameter. represents the integral symbol, S(t) and W(t) are two different signals, t is time, S * (t) represents the complex conjugate of signal S(t), W(t) is another signal that performs mutual fuzzy operation with signal S(t), It is a complex exponential function containing frequency offset f and time delay τ, which is used to introduce the corresponding phase shift, dt differential integration, which means integrating the signal in time; The mutual ambiguity function has the following properties: Among them, |χ cross,S,W (0,0)| represents the amplitude of the mutual ambiguity function at (τ,f)=(0,0), represents the square root of the energy of the signal S, represents the square root of the energy of signal W, | cross，S，W (τ, f)| 2 and | cross,S,W (-τ,-f)| represents the magnitude of the mutual ambiguity function at (τ,f) and (-τ,-f), respectively. represents the double integral symbol, indicating that the integration is performed over the entire range of time delay τ and frequency offset f, |χ cross,S,W (τ,f)| 2 represents the square of the modulus of the mutual ambiguity function, represents the energy density, dτ represents the infinitesimal time delay τ, and df represents the infinitesimal frequency offset f, indicating that the mutual ambiguity function has symmetry about (τ, f).

3. The method for optimizing the comprehensive mutual fuzzy function loop based on the determined information according to claim 2, characterized in that: The zero-sum game model between the radar party and the jammer is: Radar said: The interfering party is represented as follows: Among them, Q target and Q jammer Represents the characteristic set of target and interference source respectively, g k,p is a discrete template of the mutual ambiguity function, which is used to describe the characteristics of the signal in the radar system, where k and p are the discretized time delay and frequency offset indices, and d auto∨cross (τ,f) is the error between the discrete template of the mutual fuzzy function and its expected value, χ auto∨cross,S,W (τ, f) are the self-ambiguity and mutual-ambiguity functions used to describe the correlation between signals S and W under time delay τ and frequency offset f. and They represent the local signal-to-noise ratio and interference-to-noise ratio respectively, and their dynamic changes are determined by the internal state of the system and the external control input u radar (t) and v jammer (t) determined by and It is a function that describes the dynamic evolution and is used to calculate the rate of change of the local signal-to-noise ratio and interference-to-noise ratio; The zero-sum game model is established as follows: Among them, min X,Y It indicates that the variables X and Y are minimized. represents the sum of the index set T×F of time delay and frequency offset, g kp Represents the discrete template of the mutual ambiguity function, describing the characteristics of the signal under time delay k and frequency offset p, represents the modulus of the expected value of the mutual fuzzy function template, X H represents the conjugate transpose of X, J kp represents the complex number representation of the mutual ambiguity function template, Y represents the received signal processing vector in the optimization model; Equivalent to in, Represents the variables X, Y and phase angle φ k,p To minimize the operation, Represents a complex exponential function, introducing a phase angle φ kp The phase shift, φ k,p Represents the phase angle, which is used to describe the phase shift of a signal.

4. The method for optimizing the comprehensive mutual fuzzy function loop based on the determined information according to claim 3, characterized in that: The standard model form of the zero-sum game model is: Among them, min means minimizing the objective function, represents the sum of the index set T×F of time delay and frequency offset, |X k | represents the modulus of the kth element in X, Indicates that all indexes k that meet the conditions are constrained, and st indicates the constraints. represents the square root of MN, φ kp ∈[0,2π) represents the phase angle φ kp The value range constraints.

5. The method for optimizing the comprehensive mutual fuzzy function loop based on the determined information according to claim 4, characterized in that: The simplified model of the zero-sum game model is as follows: Determine the symbol vector of the first pulse The following M-1 pulse symbol vectors are derived from the symbol vector of the first pulse. Then at this time Get the simplified model Among them, min means that the goal is to minimize the subsequent expression, g kp Represents the discrete template of the mutual ambiguity function, describing the characteristics of the signal under time delay k and frequency offset p, Represents a complex exponential function, introducing a phase angle φ kp The phase shift, φ kp Represents the phase angle, which is used to describe the phase shift of the signal. H represents the conjugate transpose of X, J kp represents the complex representation of the mutual ambiguity function template, Y represents the received signal processing vector in the optimization model, st represents the constraint condition, |X 1,n |=1 indicates a constraint condition, ensuring that the first element X in X 1,n The modulus is 1, and for all n, |ρ m |=1 indicates a constraint condition, ensuring ρ m The modulus is 1, for all m, X m,n =ρ m X 1,n Represents constraints, describing how the subsequent parts of X are constructed, where X m,n By m Multiply by X 1,n Get, X m,N+l =0 indicates a constraint condition, ensuring that some elements in X are zero. Represents a constraint condition, ensuring that the modulus of Y is φ kp ∈[0,2π) represents the constraint condition, limiting the phase angle φ kp The value range of .

6. The method for optimizing the comprehensive mutual fuzzy function loop based on the determined information according to claim 5, characterized in that: The template-independent model is: Among them, min means that the goal is to minimize the subsequent expression, X H represents the conjugate transpose of X, J kp The complex form describing the template of the mutual ambiguity function is related to the time delay k and the frequency offset p. Y is the received signal processing vector in the optimization model. st represents the constraint condition. |X k |=1 constraint, ensuring that the modulus of the kth row in X is 1, where k=(m,n), m represents the pulse index, n represents the symbol index, and X k = 0 constraint, ensuring that the elements at specific positions in X are zero, where k = (m, N + l), m represents the pulse index, l represents the additional index offset, Constraints ensure that the modulus of Y is Where M and N represent the size of the time and frequency dimensions, respectively, kp ∈[0,2π) constraint, limiting the phase angle φ kp The value range of , where (k, p) represents the combination of time delay and frequency offset; The model is written in standard model form as:

7. The method for optimizing the comprehensive mutual fuzzy function loop based on the determined information according to claim 6, characterized in that: The simplified model of the model that does not rely on templates is: Determine the symbol vector of the first pulse The following M-1 pulse symbol vectors are derived from the symbol vector of the first pulse. Then at this time Get the simplified model Among them, min means that the goal is to minimize the subsequent expression, X H represents the conjugate transpose of X, J kp The complex form describing the template of the mutual ambiguity function is related to the time delay k and the frequency offset p. Y is the received signal processing vector in the optimization model. st represents the constraint condition. |X 1,n |=1 constraint, ensuring that the first element X in X 1,n The modulus is 1, and for all n, |ρ m |=1 constraint, ensuring ρ m The modulus is 1, for all m, ρ m is the complex scaling factor of the symbol vector, X m,n =ρ m X 1,n Constraints that describe how the subsequent parts of X are constructed, where X m,n By m Multiply by X 1,n Get, X m,N+l = 0 constraint, ensuring that some elements in X are zero, Constraints ensure that the modulus of Y is Where M and N represent the sizes of time and frequency dimensions respectively.

8. The method for optimizing the comprehensive mutual fuzzy function loop based on the determined information according to claim 7, characterized in that: In the computationally feasible special model, when N+L=1, the code element degenerates into a pulse, and one pulse repetition period is one code element. At this time, the phase coding degenerates into the initial phase modulation of the pulse: X=(X1,X2,…,X m ,…,X M ) T is the initial phase modulation of M pulses, and the model becomes At this time, the indicator set is min means that the goal is to minimize the subsequent expression, g kp Represents the discrete template of the mutual ambiguity function, describing the characteristics of the signal under time delay k and frequency offset p, Represents a complex exponential function, introducing a phase angle φ kp The phase shift, φ kp Represents the phase angle, which is used to describe the phase shift of the signal. H represents the conjugate transpose of X, J kp represents the complex representation of the mutual ambiguity function template, Y represents the received signal processing vector in the optimization model, st represents the constraint condition, |X m |=1 represents the constraint condition, ensuring that the modulus of the mth row in X is 1. For all m, Represents a constraint condition, ensuring that the modulus of Y is φ kp ∈[0,2π) represents the constraint condition, limiting the phase angle φ kp The value range of .

Images