METHOD AND DEVICE FOR DESIGNING AND OPTIMIZING MULTI-DEGREE-OF-FREEDOM FREQUENCY-MODULATION SIGNAL

Abstract

Provided are a method and device for designing and optimizing a multi-degree-of-freedom frequency-modulation (FM) signal, and a computer storage medium. The method includes that: a time domain function of the multi-degree-of-freedom FM signal is determined; the constraint condition and the objective function of multi-degree-of-freedom FM signal optimization are established; an algorithm model of augmented Lagrangian genetic algorithm is determined based on the constraint condition and the objective function; characteristic parameters of the first iteration of the algorithm model are initialized, where the characteristic parameters of the first iteration at least include the Lagrange multiplier of the first iteration and the offset of the first iteration; initialization signal of the multi-degree-of-freedom FM signal is acquired, and initialization frequency control points of the initialization signal are determined based on the initialization signal.

Claims

1. A method for designing and optimizing a multi-degree-of-freedom frequency-modulation (FM) signal, comprising: determining a time domain function of the multi-degree-of-freedom FM signal; establishing a constraint condition and an objective function for optimizing the multi-degree-of-freedom FM signal; determining an algorithm model of an augmented Lagrangian genetic algorithm based on the constraint condition and the objective function; initializing characteristic parameters of a first iteration of the algorithm model, wherein the characteristic parameters of the first iteration at least comprise a Lagrange multiplier of the first iteration and an offset of the first iteration; acquiring an initialization signal of the multi-degree-of-freedom FM signal, and determining initialization frequency control points of the initialization signal based on the initialization signal; calculating fitness of the initialization frequency control points, a Lagrange multiplier of a second iteration and an offset of the second iteration based on the characteristic parameters of the first iteration and the initialization frequency control points by using the algorithm model; determining frequency control points of the second iteration from the initialization frequency control points based on the fitness by using a selection strategy; determining frequency control points of a third iteration based on the frequency control points of the second iteration, the Lagrange multiplier of the second iteration and the offset of the second iteration by using the algorithm model and the selection strategy; continuously performing iterations of the Lagrange multiplier and the offset based on the algorithm model and the selection strategy until the augmented Lagrangian genetic algorithm converges to obtain an optimized multi-degree-of-freedom FM signal; and performing signal transmission by adopting the optimized multi-degree-of-freedom FM signal.

2. The method for designing and optimizing the multi-degree-of-freedom FM signal according to claim 1, wherein the operation of determining the time domain function of the multi-degree-of-freedom FM signal comprises: establishing the time-frequency function of the multi-degree-of-freedom FM signal in the form of a piecewise linear function, and determining a time domain function of the multi-degree-of-freedom FM signal based on the time-frequency function, wherein the operation of establishing the time-frequency function of the multi-degree-of-freedom FM signal in the form of the piecewise linear function comprises: in a Cartesian coordinate system, defining a coordinate of a time-frequency relation of the multi-degree-of-freedom FM signal to be (t, f), defining a pulse width of the multi-degree-of-freedom FM signal to be T, defining a bandwidth to be B, and constructing the time-frequency function of the multi-degree-of-freedom FM signal by 2n+2 stages of piecewise linear functions, wherein in a time-frequency relation coordinate plane, time piecewise points are uniformly distributed, and there are 2n+3 time piecewise point as follows: $t_s={\left[-\frac{T}{2},T_{2.Math..Math.n},\mathrm{.Math.}.Math.,T_{21},0,T_{11},\mathrm{.Math.}.Math.,T_{1.Math..Math.n},\frac{T}{2}\right]}^T,.Math.T_{1.Math..Math.i}=\frac{T}{2.Math..Math.n+2}.Math.i,T_{2.Math..Math.i}=-\frac{T}{2.Math..Math.n+2}.Math.i,$ (i=1, 2, . . . , n) being known quantities; in the time-frequency relation coordinate plane, 2n frequency control points are given and expressed as follows: B.sub.c=[B.sub.2n, . . . , B.sub.21, B.sub.11, . . . , B.sub.1n].sup.T, corresponding 2n+3 frequency piecewise points are expressed as follows: $B_s={\left[-\frac{B}{2},B_{2.Math.n},\mathrm{.Math.}.Math.,B_{21},0,B_{11},\mathrm{.Math.}.Math.,B_{1.Math.n},\frac{B}{2}\right]}^T,$ and the piecewise linear function of the time-frequency function is expressed as follows: $f\left(t\right)=\{\begin{array}{cc}{B}_{2.Math.n}+k_{2.Math.n}\left(t-T_{2.Math.n}\right)& -T/2t<T_{2.Math.n}\\ \mathrm{.Math.}& \mathrm{.Math.}\\ B_{21}+k_{21}\left(t-T_{21}\right)& T_{22}t<T_{21}\\ k_{20}.Math.t& T_{21}t<0\\ k_{10}.Math.t& 0t<T_{11}\\ B_{11}+k_{11}\left(t-T_{11}\right)& T_{11}t<T_{12}\\ \mathrm{.Math.}& \mathrm{.Math.}\\ B_{1.Math.n}+k_{1.Math.n}\left(t-T_{1.Math.n}\right)& T_{1.Math.n}t<T/2\end{array},$ where k.sub.1i and k.sub.2i are frequency modulation rates of stages of the piecewise linear function and are expressed as follows: $k_{10}=\frac{B_{11}}{T_{11}},\mathrm{.Math.}.Math.,\underset{1i<n}{k_{1.Math.i}}=\frac{B_{1.Math.\left(i+1\right)}-B_{1.Math.i}}{T_{1.Math.\left(i+1\right)}-T_{1.Math.i}},\mathrm{.Math.}.Math.,.Math.k_{1.Math.n}=\frac{\frac{B}{2}-B_{1.Math.n}}{\frac{T}{2}-T_{1.Math.n}},.Math.k_{20}=\frac{B_{21}}{T_{21}},\mathrm{.Math.}.Math.,\underset{1i<n}{k_{2.Math.i}}=\frac{B_{2.Math.\left(i+1\right)}-B_{2.Math.i}}{T_{2.Math.\left(i+1\right)}-T_{2.Math.i}},\mathrm{.Math.}.Math.,.Math.k_{2.Math.n}=\frac{\frac{B}{2}+B_{2.Math.n}}{\frac{T}{2}+T_{2.Math.n}}.$

3. The method for designing and optimizing the multi-degree-of-freedom FM signal according to claim 1, wherein the operation of establishing the constraint condition and the objective function for optimizing the multi-degree-of-freedom FM signal comprises: establishing the constraint condition and the objective function for optimizing the multi-degree-of-freedom FM signal based on a main lobe width (MW) and a peak sidelobe ratio (PSLR) in a performance of an autocorrelation function of the multi-degree-of-freedom FM signal.

4. The method for designing and optimizing the multi-degree-of-freedom FM signal according to claim 1, wherein the selection strategy comprises a roulette wheel selection (RWS) rule and a crossover and mutation rule.

5. The method for designing and optimizing the multi-degree-of-freedom FM signal according to claim 1, wherein the time-frequency function of the multi-degree-of-freedom FM signal is expressed as follows:
s(t)=A exp(j2.sub.0.sup.if(t)dt), where A is a signal amplitude, and f(t) is the time-frequency function.

6. The method for designing and optimizing the multi-degree-of-freedom FM signal according to claim 3, wherein the operation of establishing the constraint condition and the objective function for optimizing the multi-degree-of-freedom FM signal based on the MW and the PSLR in the performance of the autocorrelation function of the multi-degree-of-freedom FM signal comprises: establishing the constraint condition and the objective function for optimizing the multi-degree-of-freedom FM signal according to an ideal performance of the autocorrelation function of the multi-degree-of-freedom FM signal based on a principle of sidelobe suppression without widening a main lobe: $\underset{B_c}{\min }.Math.f.Math..Math.\left(\mathrm{PSLR}\right),\mathrm{such}.Math..Math.\mathrm{that}.Math.$ $c\left(\mathrm{MW}\right)0,-B/2B_{c}B/2,$ where c(MW) is a nonlinear constraint of the MW.

7. The method for designing and optimizing the multi-degree-of-freedom FM signal according to claim 1, wherein the operation of acquiring an initialization signal of the multi-degree-of-freedom FM signal and determining the initialization frequency control points of the initialization signal based on the initialization signal comprises: selecting a window function P(f) as a power spectrum of the multi-degree-of-freedom FM signal, and calculating a group delay function as follows: $T\left(f\right)={}_0^f.Math.\frac{P\left(f\right)}{C}.Math.\mathrm{df},.Math.\mathrm{where}.Math..Math.C=\frac{1}{T}.Math.{}_{-B/2}^{B/2}.Math.P\left(f\right).Math.\mathrm{df},$ the time-frequency function of the initialization signal is expressed as follows:
f(t)=T.sup.1(f), based on the above inverse function, solving the initialization signal by using a numerical method, obtaining B.sub.s according to the time piecewise point t.sub.s, and obtaining frequency piecewise points B.sub.c, the frequency piecewise points B.sub.c being the initialization frequency control points.

8. The method for designing and optimizing the multi-degree-of-freedom FM signal according to claim 1, wherein the operation of determining the algorithm model of the augmented Lagrangian genetic algorithm based on the constraint condition and the objective function comprises: converting the constraint condition and the objective function to:
(B.sub.c,,s)=f(B.sub.c)s log(sc(B.sub.c)), where f(B.sub.c)=PSLR(B), PSLR(B.sub.c) being the PSLR of the multi-degree-of-freedom FM signal obtained according to the frequency control points B.sub.c; c(B.sub.c)=MW(B.sub.c), MW(B.sub.c) being the MW of the multi-degree-of-freedom FM signal obtained according to the frequency control points B.sub.c; and is a constraint value of a main lobe.

9. A device for designing and optimizing a multi-degree-of-freedom frequency-modulation (FM) signal, comprising: a processor; and a memory storing instructions, which, when executed by the processor, cause the processor to execute operations comprising: determining a time domain function of the multi-degree-of-freedom FM signal; establishing a constraint condition and an objective function for optimizing the multi-degree-of-freedom FM signal; and determining an algorithm model of an augmented Lagrangian genetic algorithm based on the constraint condition and the objective function; initializing characteristic parameters of a first iteration of the algorithm model, wherein the characteristic parameters of the first iteration at least comprise a Lagrange multiplier of the first iteration and an offset of the first iteration; acquiring an initialization signal of the multi-degree-of-freedom FM signal, and determine initialization frequency control points of the initialization signal based on the initialization signal; calculating fitness of the initialization frequency control points, a Lagrange multiplier of a second iteration and an offset of the second iteration based on the characteristic parameters of the first iteration and the initialization frequency control points by using the algorithm model; determining frequency control points of the second iteration from the initialization frequency control points based on the fitness by using selection strategy; determining frequency control points of a third iteration based on the frequency control points of the second iteration, the Lagrange multiplier of the second iteration and the offset of the second iteration by using the algorithm model and the selection strategy; continuously performing iterations of the Lagrange multiplier and the offset based on the algorithm model and the selection strategy, until the augmented Lagrangian genetic algorithm converges to obtain an optimized multi-degree-of-freedom FM signal; and performing signal transmission by adopting the optimized multi-degree-of-freedom FM signal.

10. The device for designing and optimizing a multi-degree-of-freedom FM signal according to claim 9, wherein the operation of determining the time domain function of the multi-degree-of-freedom FM signal comprises: establishing the time-frequency function of the multi-degree-of-freedom FM signal in the form of a piecewise linear function, and determining a time domain function of the multi-degree-of-freedom FM signal based on the time-frequency function, wherein the operation of establishing the time-frequency function of the multi-degree-of-freedom FM signal in the form of the piecewise linear function comprises: in a Cartesian coordinate system, defining a coordinate of a time-frequency relation of the multi-degree-of-freedom FM signal to be (t, f), defining a pulse width of the multi-degree-of-freedom FM signal to be T, defining a bandwidth to be B, and constructing the time-frequency function of the multi-degree-of-freedom FM signal by 2n+2 stages of piecewise linear functions, wherein in a time-frequency relation coordinate plane, time piecewise points are uniformly distributed, and there are 2n+3 time piecewise point as follows: $t_s={\left[-\frac{T}{2},T_{2.Math.n},\mathrm{.Math.}.Math.,T_{21},0,T_{11},\mathrm{.Math.}.Math.,T_{1.Math.n},\frac{T}{2}\right]}^T,.Math.T_{1.Math.i}=\frac{T}{2.Math.n+2}.Math.i,T_{2.Math.i}=-\frac{T}{2.Math.n+2}.Math.i,$ (i=1 2, . . . , n) being known quantities; in a time-frequency relation coordinate plane, 2n frequency control points are given and expressed as follows: B.sub.c=[B.sub.2n, . . . , B.sub.21, B.sub.11, . . . , B.sub.1n].sup.T, corresponding 2n+3 frequency piecewise points are expressed as follows: $B_s={\left[-\frac{B}{2},B_{2.Math.n},\mathrm{.Math.}.Math.,B_{21},0,B_{11},\mathrm{.Math.}.Math.,B_{1.Math.n},\frac{B}{2}\right]}^T,$ and the piecewise linear function of the time-frequency function is expressed as follows: $f\left(t\right)=\{\begin{array}{cc}{B}_{2.Math.n}+k_{2.Math.n}\left(t-T_{2.Math.n}\right)& -T/2t<T_{2.Math.n}\\ \mathrm{.Math.}& \mathrm{.Math.}\\ B_{21}+k_{21}\left(t-T_{21}\right)& T_{22}t<T_{21}\\ k_{20}.Math.t& T_{21}t<0\\ k_{10}.Math.t& 0t<T_{11}\\ B_{11}+k_{11}\left(t-T_{11}\right)& T_{11}t<T_{12}\\ \mathrm{.Math.}& \mathrm{.Math.}\\ B_{1.Math.n}+k_{1.Math.n}\left(t-T_{1.Math.n}\right)& T_{1.Math.n}t<T/2\end{array},$ where k.sub.1i and k.sub.2i are frequency modulation rates of stages of the piecewise linear function and are expressed as follows: $k_{10}=\frac{B_{11}}{T_{11}},\mathrm{.Math.}.Math.,\underset{1i<n}{k_{1.Math.i}}=\frac{B_{1.Math.\left(i+1\right)}-B_{1.Math.i}}{T_{1.Math.\left(i+1\right)}-T_{1.Math.i}},\mathrm{.Math.}.Math.,.Math.k_{1.Math.n}=\frac{\frac{B}{2}-B_{1.Math.n}}{\frac{T}{2}-T_{1.Math.n}},.Math.k_{20}=\frac{B_{21}}{T_{21}},\mathrm{.Math.}.Math.,\underset{1i<n}{k_{2.Math.i}}=\frac{B_{2.Math.\left(i+1\right)}-B_{2.Math.i}}{T_{2.Math.\left(i+1\right)}-T_{2.Math.i}},\mathrm{.Math.}.Math.,.Math.k_{2.Math.n}=\frac{\frac{B}{2}+B_{2.Math.n}}{\frac{T}{2}+T_{2.Math.n}}.$

11. The device for designing and optimizing the multi-degree-of-freedom FM signal according to claim 9, wherein the operation of establishing the constraint condition and the objective function for optimizing the multi-degree-of-freedom FM signal comprises: establishing the constraint condition and the objective function for optimizing the multi-degree-of-freedom FM signal based on a main lobe width (MW) and a peak sidelobe ratio (PSLR) in a performance of an autocorrelation function of the multi-degree-of-freedom FM signal.

12. The device for designing and optimizing the multi-degree-of-freedom FM signal according to claim 9, wherein the selection strategy comprises a roulette wheel selection (RWS) rule and a crossover and mutation rule.

13. The device for designing and optimizing the multi-degree-of-freedom FM signal according to claim 9, wherein the time-frequency function of the multi-degree-of-freedom FM signal is expressed as follows:
s(t)=A exp(j2.sub.0.sup.if(t)dt), where A is a signal amplitude, and f(t) is the time-frequency function.

14. The device for designing and optimizing the multi-degree-of-freedom FM signal according to claim 11, wherein the operation of establishing the constraint condition and the objective function for optimizing the multi-degree-of-freedom FM signal based on the MW and the PSLR in the performance of the autocorrelation function of the multi-degree-of-freedom FM signal comprises: establishing the constraint condition and the objective function for optimizing the multi-degree-of-freedom FM signal according to an ideal performance of the autocorrelation function of the multi-degree-of-freedom FM signal based on a principle of sidelobe suppression without widening a main lobe: $\underset{B_c}{\min }.Math.f.Math..Math.\left(\mathrm{PSLR}\right),\mathrm{such}.Math..Math.\mathrm{that}.Math.$ $c\left(\mathrm{MW}\right)0,-B/2B_{c}B/2,$ where c(MW) is a nonlinear constraint of the MW.

15. The device for designing and optimizing the multi-degree-of-freedom FM signal according to claim 9, wherein the operation of acquiring an initialization signal of the multi-degree-of-freedom FM signal and determining the initialization frequency control points of the initialization signal based on the initialization signal comprises: selecting a window function P(f) as a power spectrum of the multi-degree-of-freedom FM signal, and calculating a group delay function as follows: $T\left(f\right)={}_0^f.Math.\frac{P\left(f\right)}{C}.Math.\mathrm{df},.Math.\mathrm{where}.Math..Math.C=\frac{1}{T}.Math.{}_{-B/2}^{B/2}.Math.P\left(f\right).Math.\mathrm{df},$ the time-frequency function of the initialization signal is expressed as follows:
f(t)=T.sup.1(f), based on the above inverse function, solving the initialization signal by using a numerical method, obtaining B.sub.s according to the time piecewise point t.sub.s, and obtaining frequency piecewise points B.sub.c, the frequency piecewise points B.sub.c being the initialization frequency control points.

16. The device for designing and optimizing the multi-degree-of-freedom FM signal according to claim 9, wherein the operation of determining the algorithm model of the augmented Lagrangian genetic algorithm based on the constraint condition and the objective function comprises: converting the constraint condition and the objective function to:
(B.sub.c,,s)=f(B.sub.c)s log(sc(B.sub.c)), where f(B.sub.c)=PSLR(B.sub.c), PSLR(B.sub.c) being the PSLR of the multi-degree-of-freedom FM signal obtained according to the frequency control points B.sub.c; c(B.sub.c)=MW(B.sub.c), MW(B.sub.c) being the MW of the multi-degree-of-freedom FM signal obtained according to the frequency control points B.sub.c; and is a constraint value of a main lobe.

17. A non-transitory computer-readable storage medium having stored thereon computer-executable instructions, which when being executed by a processor, to execute operations comprising: determining a time domain function of the multi-degree-of-freedom frequency-modulation (FM) signal; establishing a constraint condition and an objective function for optimizing the multi-degree-of-freedom FM signal; determining an algorithm model of an augmented Lagrangian genetic algorithm based on the constraint condition and the objective function; initializing characteristic parameters of a first iteration of the algorithm model, wherein the characteristic parameters of the first iteration at least comprise a Lagrange multiplier of the first iteration and an offset of the first iteration; acquiring an initialization signal of the multi-degree-of-freedom FM signal, and determining initialization frequency control points of the initialization signal based on the initialization signal; calculating fitness of the initialization frequency control points, a Lagrange multiplier of second iteration and an offset of the second iteration based on the characteristic parameters of the first iteration and the initialization frequency control points by using the algorithm model; determining frequency control points of the second iteration from the initialization frequency control points based on the fitness by using a selection strategy; determining frequency control points of a third iteration based on the frequency control points of the second iteration, the Lagrange multiplier of the second iteration and the offset of the second iteration by using the algorithm model and the selection strategy; continuously performing iterations of the Lagrange multiplier and the offset based on the algorithm model and the selection strategy until the augmented Lagrangian genetic algorithm converges to obtain an optimized multi-degree-of-freedom FM signal; and performing signal transmission by adopting the optimized multi-degree-of-freedom FM signal.

18. The non-transitory computer-readable storage medium according to claim 17, wherein the operation of determining the time domain function of the multi-degree-of-freedom FM signal comprises: establishing the time-frequency function of the multi-degree-of-freedom FM signal in the form of a piecewise linear function, and determining a time domain function of the multi-degree-of-freedom FM signal based on the time-frequency function, wherein the operation of establishing the time-frequency function of the multi-degree-of-freedom FM signal in the form of the piecewise linear function comprises: in a Cartesian coordinate system, defining a coordinate of a time-frequency relation of the multi-degree-of-freedom FM signal to be (t, f), defining a pulse width of the multi-degree-of-freedom FM signal to be T, defining a bandwidth to be B, and constructing the time-frequency function of the multi-degree-of-freedom FM signal by 2n+2 stages of piecewise linear functions, wherein in a time-frequency relation coordinate plane, time piecewise points are uniformly distributed, and there are 2n+3 time piecewise point as follows: $t_s={\left[-\frac{T}{2},T_{2.Math.n},\mathrm{.Math.}.Math.,T_{21},0,T_{11},\mathrm{.Math.}.Math.,T_{1.Math.n},\frac{T}{2}\right]}^T,.Math.T_{1.Math.i}=\frac{T}{2.Math.n+2}.Math.i,T_{2.Math.i}=-\frac{T}{2.Math.n+2}.Math.i,$ (i=1, 2, . . . , n) being known quantities; in the time-frequency relation coordinate plane, 2n frequency control points are given and expressed as follows: B.sub.c=[B.sub.2n, . . . , B.sub.21, B.sub.11, . . . , B.sub.1n].sup.T, corresponding 2n+3 frequency piecewise points are $B_s={\left[-\frac{B}{2},B_{2.Math.n},\mathrm{.Math.}.Math.,B_{21},0,B_{11},\mathrm{.Math.}.Math.,B_{1.Math.n},\frac{B}{2}\right]}^T,$ and the piecewise linear function of the time-frequency function is expressed as follows: $f\left(t\right)=\{\begin{array}{cc}{B}_{2.Math.n}+k_{2.Math.n}\left(t-T_{2.Math.n}\right)& -T/2t<T_{2.Math.n}\\ \mathrm{.Math.}& \mathrm{.Math.}\\ B_{21}+k_{21}\left(t-T_{21}\right)& T_{22}t<T_{21}\\ k_{20}.Math.t& T_{21}t<0\\ k_{10}.Math.t& 0t<T_{11}\\ B_{11}+k_{11}\left(t-T_{11}\right)& T_{11}t<T_{12}\\ \mathrm{.Math.}& \mathrm{.Math.}\\ B_{1.Math.n}+k_{1.Math.n}\left(t-T_{1.Math.n}\right)& T_{1.Math.n}t<T/2\end{array},$ where k.sub.1i and k.sub.2i are frequency modulation rates of stages of the piecewise linear function and are expressed as follows: $k_{10}=\frac{B_{11}}{T_{11}},\mathrm{.Math.}.Math.,\underset{1i<n}{k_{1.Math.i}}=\frac{B_{1.Math.\left(i+1\right)}-B_{1.Math.i}}{T_{1.Math.\left(i+1\right)}-T_{1.Math.i}},\mathrm{.Math.}.Math.,.Math.k_{1.Math.n}=\frac{\frac{B}{2}-B_{1.Math.n}}{\frac{T}{2}-T_{1.Math.n}},.Math.k_{20}=\frac{B_{21}}{T_{21}},\mathrm{.Math.}.Math.,\underset{1i<n}{k_{2.Math.i}}=\frac{B_{2.Math.\left(i+1\right)}-B_{2.Math.i}}{T_{2.Math.\left(i+1\right)}-T_{2.Math.i}},\mathrm{.Math.}.Math.,.Math.k_{2.Math.n}=\frac{\frac{B}{2}+B_{2.Math.n}}{\frac{T}{2}+T_{2.Math.n}}.$

19. The non-transitory computer-readable storage medium according to claim 17, wherein the operation of establishing the constraint condition and the objective function for optimizing the multi-degree-of-freedom FM signal comprises: establishing the constraint condition and the objective function for optimizing the multi-degree-of-freedom FM signal based on a main lobe width (MW) and a peak sidelobe ratio (PSLR) in a performance of an autocorrelation function of the multi-degree-of-freedom FM signal.

20. The non-transitory computer-readable storage medium according to claim 17, wherein the time-frequency function of the multi-degree-of-freedom FM signal is expressed as follows:
s(t)=A exp(j2.sub.0.sup.if(t)dt), where A is a signal amplitude, and f(t) is the time-frequency function.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0033] FIG. 1 is a schematic diagram of time-frequency relation of a multi-degree-of-freedom FM signal according to an embodiment of the disclosure;

[0034] FIG. 2 is a schematic flowchart of a method for designing and optimizing a multi-degree-of-freedom FM signal according to an embodiment of the disclosure;

[0035] FIG. 3 is a schematic diagram of time-frequency relation of the optimized multi-degree-of-freedom FM signal according to an embodiment of the disclosure;

[0036] FIG. 4 is a time domain waveform diagram of the optimized multi-degree-of-freedom FM signal according to an embodiment of the disclosure;

[0037] FIG. 5 is an autocorrelation function diagram of the optimized multi-degree-of-freedom FM signal according to an embodiment of the disclosure; and

[0038] FIG. 6 is a schematic structural diagram of a device for designing and optimizing a multi-degree-of-freedom FM signal according to an embodiment of the disclosure.

DETAILED DESCRIPTION

[0039] The disclosure is further explained in detail in conjunction with the drawings.

[0040] As illustrated in FIG. 2, the method for designing and optimizing a multi-degree-of-freedom FM signal in the disclosure includes the following operations.

[0041] At block 101, a time domain function of the multi-degree-of-freedom FM signal is determined.

[0042] Specifically, a time-frequency function of the multi-degree-of-freedom FM signal is established in the form of piecewise linear function, and then the time domain function of the multi-degree-of-freedom FM signal is determined based on the time-frequency function.

[0043] Specifically, in Cartesian coordinate system, the coordinate of time-frequency relation of the multi-degree-of-freedom FM signal is defined to be (t, f). It is assumed that the pulse width of the multi-degree-of-freedom FM signal is defined to be T, the bandwidth of the multi-degree-of-freedom FM signal is defined to be B, and the time-frequency function of the multi-degree-of-freedom FM signal is constructed by 2n+2 stages of piecewise linear functions, as illustrated in FIG. 1. In a time-frequency relation coordinate plane, time piecewise points, namely abscissa piecewise points are uniformly distributed, and 2n+3 time piecewise points can be expressed as follows:

[00007]$t_s={\left[-\frac{T}{2},T_{2.Math..Math.n},\mathrm{.Math.}.Math.,T_{21},0,T_{11},\mathrm{.Math.}.Math.,T_{1.Math..Math.n},\frac{T}{2}\right]}^T,\mathrm{where}$$T_{1.Math..Math.i}=\frac{T}{2.Math..Math.n+2}.Math.i,T_{2.Math..Math.i}=-\frac{T}{2.Math..Math.n+2}.Math.i,$

(i=1, 2, . . . , n) are known quantities. In the time-frequency relation coordinate plane, 2n frequency control points are given and expressed as follows: B.sub.c=[B.sub.2n, . . . , B.sub.21, B.sub.11, . . . , B.sub.1n].sup.T, then corresponding 2n+3 frequency piecewise points can be expressed as follows:

[00008]$B_s={\left[-\frac{B}{2},B_{2.Math..Math.n},\mathrm{.Math.}.Math.,B_{21},0,B_{11},\mathrm{.Math.}.Math.,B_{1.Math..Math.n},\frac{B}{2}\right]}^T.$

When the piecewise points in the time-frequency relation coordinate plane are given, the time-frequency relation of the multi-degree-of-freedom FM signal can be expressed by using the piecewise linear function as follows:

[00009]$\begin{array}{cc}f\left(t\right)=\{\begin{array}{cc}{B}_{2.Math..Math.n}+k_{2.Math..Math.n}\left(t-T_{2.Math..Math.n}\right)& -T/2t<T_{2.Math..Math.n}\\ \mathrm{.Math.}& \mathrm{.Math.}\\ B_{21}+k_{21}\left(t-T_{21}\right)& T_{22}t<T_{21}\\ k_{20}.Math.t& T_{21}t<0\\ k_{10}.Math.t& 0t<T_{11}\\ B_{11}+k_{11}\left(t-T_{11}\right)& T_{11}t<T_{12}\\ \mathrm{.Math.}& \mathrm{.Math.}\\ B_{1.Math..Math.n}+k_{1.Math..Math.n}\left(t-T_{1.Math..Math.n}\right)& T_{1.Math..Math.n}t<T/2\end{array}& \left(1\right)\end{array}$

[0044] In the function (1), k.sub.1, and k.sub.2i represent frequency modulation rates of stages of the piecewise linear function, and the time-frequency relation is illustrated in FIG. 1.

[00010]$\begin{array}{cc}{k}_{10}=\frac{B_{11}}{T_{11}},\mathrm{.Math.}.Math.,\underset{1i<n}{K_{1.Math..Math.i}}=\frac{B_{1.Math.\left(i+1\right)}-B_{1.Math..Math.i}}{T_{1.Math.\left(i+1\right)}-T_{1.Math..Math.i}},\mathrm{.Math.}.Math.,k_{1.Math..Math.n}=\frac{\frac{B}{2}-B_{1.Math..Math.n}}{\frac{T}{2}-T_{1.Math..Math.n}}& \left(2\right)\\ k_{20}=\frac{B_{21}}{T_{21}},\mathrm{.Math.}.Math.,\underset{1i<n}{K_{2.Math..Math.i}}=\frac{B_{2.Math.\left(i+1\right)}-B_{2.Math..Math.i}}{T_{2.Math.\left(i+1\right)}-T_{2.Math..Math.i}},\mathrm{.Math.}.Math.,k_{2.Math..Math.n}=\frac{\frac{B}{2}+B_{2.Math..Math.n}}{\frac{T}{2}+T_{2.Math..Math.n}}.& \left(3\right)\end{array}$

[0045] The time domain function of the multi-degree-of-freedom FM signal can be expressed as follows:

s(t)=A exp(j2.sub.0.sup.if(t)dt) (4)

In the function (4), A is signal amplitude, and f(t) is the time-frequency function.

[0046] At block 102, constraint condition and objective function of multi-degree-of-freedom FM signal optimization are established.

[0047] Specifically, the constraint condition and the objective function of the multi-degree-of-freedom FM signal optimization are established based on the main lobe width (MW) and the peak sidelobe ratio (PSLR) in the performance of the autocorrelation function of the multi-degree-of-freedom FM signal.

[0048] Specifically, ideal performance of the autocorrelation function of the multi-degree-of-freedom FM signal includes: the main lobe as narrow as possible, the peak sidelobe ratio (PSLR) as low as possible and the rapidly reduced envelope of the sidelobe fluctuation. However, the three ideal performances cannot be satisfied at the same time. Generally, the MW of the autocorrelation function is defined to be less than the peak by 3 dB. It can be known from the definitions in the formulas (1) to (4) that, a specific form of the multi-degree-of-freedom FM signal depends on the control point. Once the number of the frequency control points is determined, the time piecewise points t.sub.s are uniformly distributed on the time axis and are known quantities. The multi-degree-of-freedom FM signal can be defined by the vector B.sub.s including 2n frequency control points. A mathematical model for the multi-degree-of-freedom FM signal optimization, which is based on two key performances of impact response function of point target of an SAR image: the PSLR and the MW, is provided. The mathematical model is for the nonlinear constrained optimization. The requirement of the nonlinear constrained optimization is to perform the sidelobe suppression as much as possible without widening the main lobe, which can be expressed as follows:

[00011]$\begin{array}{cc}\underset{B_c}{\min }.Math.f\left(\mathrm{PSLR}\right),\mathrm{such}.Math..Math.\mathrm{that}.Math..Math.c\left(\mathrm{MW}\right)0,-B/2B_{c}B/2& \left(5\right)\end{array}$

[0049] In the formula (5), c(MW) is nonlinear constraint of the MW.

[0050] At block 103, an algorithm model of augmented Lagrangian genetic algorithm is determined based on the constraint condition and the objective function.

[0051] Specifically, the augmented Lagrangian genetic algorithm is a generalized form of genetic algorithm, and is advanced algorithm for solving the complex constrained optimization by combining the genetic algorithm and generalized Lagrange algorithm.

[0052] The augmented Lagrangian genetic algorithm is mathematically expressed as follows:

[00012]$\begin{array}{cc}\left(x,,s,\right)=f\left(x\right)-\underset{i=1}{\overset{m}{.Math.}}.Math.{}_i.Math.s_i.Math.\log \left(s_i-c_i\left(x\right)\right)+\underset{i=m+1}{\overset{\mathrm{mt}}{.Math.}}.Math.{}_i.Math.{\mathrm{ceq}}_i\left(x\right)+\frac{}{2}.Math.\underset{\mathrm{mt}}{\overset{i=\mathrm{mt}}{.Math.}}.Math.{{\mathrm{ceq}}_i\left(x\right)}^2& \left(6\right)\end{array}$

[0053] In the formula (6), .sub.i is Lagrange multiplier and is a non-negative number, s.sub.i is a non-negative number and represents a total offset to ensure a nonzero antilogarithm of the logarithm, is penalty factor, ceq.sub.i(x) and c.sub.i(x) are equality constraint and nonlinear inequality constraint respectively, f(x) is fitness function, in is the number of the nonlinear constraints, and mt is the total number of the constraints.

[0054] It is required to perform the sidelobe suppression as much as possible without widening the main lobe, which can be expressed as follows:

(B.sub.c,,s)=f(B.sub.c)s log(sc)(B.sub.c)) (7)

[0055] In the formula (7), f(B.sub.c)=PSLR(B.sub.c), where PSLR(B.sub.c) is the PSLR of the signal obtained according to the frequency control points B.sub.c; c(B.sub.c)=MW(B.sub.c), where MW(B.sub.c) is the MW of the multi-degree-of-freedom FM signal obtained according to the frequency control points B.sub.c; and a is a constraint value of the main lobe.

[0056] Specifically, the specific solution problem is divided into two parts by using the augmented Lagrangian genetic algorithm. One part is traditional genetic algorithm, and the other part is augmented Lagrangian algorithm. The augmented Lagrangian algorithm is used for solving the constraint problem, in which and s are continuously updated according to the formula (7).

[0057] The optimization problem is simulated as a dynamic optimization process of the natural selection of survival of the fittest by using the genetic algorithm. In search space, a chromosome represents a variable to be determined in the specific solution problem, and generally the genetic algorithm includes operations of selection, crossover and mutation of the chromosome. First, the variable is encoded according to the solution problem, fitness of the chromosomes are calculated according to an objective function value, and the chromosomes are selected according to a certain rule. Second, the selected chromosome pairs are crossed according to the mating probability r.sub.c to generate offspring. Finally, the mutation operation is performed on genes on the chromosome according to the certain mutation probability r.sub.m, and a new individual is generated in the search space of the variable. In a whole iterative updating optimization process, the chromosome individual with high fitness has high probability to be selected for generating the offspring, and the individual with low fitness is replaced by a better offspring.

[0058] In regard to the multi-degree-of-freedom FM signal optimization, in an encoding process of the problem, each chromosome is regarded as the frequency control point vector B.sub.c including 2n frequency components in the embodiment of the disclosure. An initialization signal is selected, so that the frequency control point vector B.sub.c and parameters of the algorithm can be initialized.

[0059] At block 104, characteristic parameters of first iteration of the algorithm model are initialized, where the characteristic parameters of the first iteration at least include the Lagrange multiplier of the first iteration and the offset of the first iteration.

[0060] Specifically, roulette wheel selection (RWS) rule is used to select male parents to perform crossover process in the disclosure. It is assumed that the total number of the selected chromosomes is K. K random numbers are firstly generated for the K chromosomes in the crossover process. When the random numbers corresponding to the chromosomes are less than the crossover probability r.sub.c, it is indicated that the chromosomes are selected for the crossover operation. The crossover operation is performed at one cross point herein, and a position of the cross point is randomly generated. The male parents exchange genes at the cross point to generate a new chromosome. The mutation operation refers to an operation of changing the genes on the chromosome, and the genes to be mutated are randomly selected. In the multi-degree-of-freedom FM signal optimization, the number of genes on each chromosome is 2n, then the total number of the genes can be expressed as follows: L=2Kn. The number M of the mutated genes depends on mutation probability, and can be specifically expressed as follows: M=r.sub.mL. The mutation operation is performed on M genes randomly selected from L genes, and mutation operation can be expressed as follows:

p.sub.k(i)=p.sub.k(i)*(1+rand) (8)

[0061] Algorithm parameters are initialized. The following initialized algorithm parameters are set: the number of chromosomes n, the crossover probability r.sub.c, the mutation probability r.sub.m, the Lagrange multiplier A and the offset s.

[0062] At block 105, an initialization signal of the multi-degree-of-freedom FM signal is acquired, and initialization frequency control points of the initialization signal are determined based on the initialization signal.

[0063] The window function P(f) is selected as a power spectrum of the multi-degree-of-freedom FM signal. Group delay function is calculated as follows:

[00013]$\begin{array}{cc}T\left(f\right)={}_0^f.Math.\frac{P\left(f\right)}{C}.Math.\mathrm{df}.Math..Math.\mathrm{where}& \left(9\right)\\ C=\frac{1}{T}.Math.{}_{-B/2}^{B/2}.Math.P\left(f\right).Math.\mathrm{df}& \left(10\right)\end{array}$

[0064] And the time-frequency function of the multi-degree-of-freedom FM signal can be calculated as follows:

f(t)=T.sup.1(f) (11)

[0065] Generally, in the function (11), an analytical form of an inverse function cannot be obtained. Therefore, the function (11) is solved by using numerical method in the embodiment of the disclosure. After the time-frequency relation is obtained, B.sub.s can be obtained according to the time piecewise points t.sub.s, and then the frequency control points B.sub.c are obtained, and the initialization operation is completed.

[0066] At block 106, fitness of the initialization frequency control points, the Lagrange multiplier of second iteration and the offset of the second iteration are calculated based on the characteristic parameters of the first iteration and the initialization frequency control points by using the algorithm model; and frequency control points of the second iteration are determined from the initialization frequency control points based on the fitness by using a selection strategy.

[0067] Specifically, the selection strategy includes RWS rule and crossover and mutation rule.

[0068] Specifically, the multi-degree-of-freedom FM signal corresponding to the control point represented by each chromosome is determined according to the formulas (1) to (4), and autocorrelation operation is performed on the multi-degree-of-freedom FM signal. MW and PSLR of the output signal are calculated. Fitness of the chromosome is calculated according to the formula (7). According to the formula (7), and offset s of the next optimization are calculated by using the augmented Lagrangian algorithm. Then the chromosomes are selected according to the RWS rule, and crossover operation is performed at one cross point according to the crossover probability. Finally, the mutation operation in the formula (8) is performed according to the mutation probability r.sub.m.

[0069] At block 107, frequency control points of third iteration are determined based on the frequency control points of the second iteration, the Lagrange multiplier of the second iteration and the offset of the second iteration by using the algorithm model and the selection strategy. Based on the algorithm model and the selection strategy, iterations of the Lagrange multiplier and the offset s are continuously performed, until the augmented Lagrangian genetic algorithm converges to obtain an optimized multi-degree-of-freedom FM signal.

[0070] Specifically, the selection strategy includes RWS rule and crossover and mutation rule.

[0071] Specifically, based on the obtained next generation of the chromosomes, namely the frequency control points of the second iteration, the fitness, and offset s of the next optimization are calculated with reference to the operation at block 106 by using the formula (7). And then the roulette selection, crossover and mutation operations are successively performed. That is, iterations of the Lagrange multiplier and the offset s are continuously performed by cycling the operation at block 106, until the augmented Lagrangian genetic algorithm converges. The convergence of the algorithm indicates that multi-degree-of-freedom FM signal design and optimization are completed.

[0072] At block 108, signal transmission is performed by adopting the optimized multi-degree-of-freedom FM signal.

[0073] In conclusion, according to the embodiment of the disclosure, the operations of an optimization process of a multi-degree-of-freedom FM signal are illustrated as follows.

TABLE-US-00001 Optimization process of multi-degree-of-freedom FM signal Block 1 Algorithm parameters are initialized. The following initialized algorithm parameters are set: the number of chromosomes n, the crossover probability r.sub.c, the mutation probability r.sub.m, the Lagrange multiplier and the offset s. Block 2 Based on the principle of stationary phase, an initialization multi-degree-of-freedom FM signal is generated, and initialization frequency control points B.sub.c, namely the initialization chromosome, are obtained. Block 3 The multi-degree-of-freedom FM signal corresponding to the control point represented by each chromosome is determined according to the formulas (1) to (4); the autocorrelation operation is performed on the multi-degree-of-freedom FM signal; MW and PSLR of the output signal are calculated; fitness of the chromosome is calculated according to the formula (7); and and offset s of the next optimization are calculated according to the formula (7) by using the augmented Lagrangian algorithm. Block 4 The chromosomes are selected according to the RWS rule. Block 5 Crossover operation is performed at one cross point according to the crossover probability r.sub.c. Block 6 Mutation operation in the formula (8) is performed according to the mutation probability r.sub.m. Block 7 The operations at block 3 to block 7 are cycled until the algorithm converges.

[0074] The function of the disclosure will be further described in detail in conjunction with specific examples below.

[0075] A multi-degree-of-freedom FM signal is generated according to the principle of stationary phase by using 35 dB of Taylor window function as the power spectrum, and initialization frequency control points of the signal are obtained. The pulse width of the signal is 10 us, the bandwidth is 100 MHz, and the sampling frequency is 120 MHz. The result of the signal optimization is illustrated in Table 1:

TABLE-US-00002 TABLE 1 Waveform PSLR MW Initialization waveform 34.3 dB 1.19 Optimized waveform 40.3 dB 1.17

[0076] The MW is in unit of sampling number. It can be seen from the Table 1 that compared with the initialization signal, the process of the multi-degree-of-freedom FM signal optimization in the disclosure enables the sidelobe of the multi-degree-of-freedom FM signal to be suppressed by 6 dB without widening the main lobe. In addition, it is to be noted that Taylor window is the window function with excellent performance, and the PSLR and MW indexes of the initialization waveform are good enough. According to the method, the sidelobe can further be suppressed by 6 dB without widening the main lobe, thereby the effectiveness of the disclosure can be adequately confirmed.

[0077] FIGS. 3-5 illustrate a time-frequency relation curve, a time domain waveform and output waveform of autocorrelation function of the optimized multi-degree-of-freedom FM signal. FIG. 5 is a time domain waveform diagram of the optimized signal. Compared with the LFM signal with the same time width and an equivalent sidelobe level generated by windowing, the optimized signal has capability of maintaining 1.2 dB of the system SNR, which has a great practical engineering significance.

[0078] Therefore, the multi-degree-of-freedom FM signal can be effectively designed and optimized by using the method in the disclosure.

[0079] The embodiments of the disclosure provide a device for designing and optimizing a multi-degree-of-freedom FM signal. As illustrated in FIG. 6, the device includes a time-frequency function establishing unit 601, an algorithm model determining unit 602, an algorithm parameter initialization unit 603, a frequency control point initialization unit 604 and an iteration unit 605.

[0080] The time-frequency function establishing unit 601 is configured to determine a time domain function of the multi-degree-of-freedom FM signal. Specifically, the time-frequency function of the multi-degree-of-freedom FM signal is established in the form of piecewise linear function, and then the time domain function of the multi-degree-of-freedom FM signal is determined based on the time-frequency function.

[0081] The algorithm model determining unit 602 is configured to establish constraint condition and objective function of multi-degree-of-freedom FM signal optimization. The algorithm model determining unit 602 is further configured to determine an algorithm model of augmented Lagrangian genetic algorithm based on the constraint condition and the objective function. Specifically, the algorithm model determining unit 602 is further configured to establish the constraint condition and the objective function of the multi-degree-of-freedom FM signal optimization, based on the main lobe width (MW) and the peak sidelobe ratio (PSLR) in the performance of the autocorrelation function of the multi-degree-of-freedom FM signal.

[0082] The algorithm parameter initialization unit 603 is configured to initialize characteristic parameters of first iteration of the algorithm model, where the characteristic parameters of the first iteration at least include the Lagrange multiplier of the first iteration and the offset of the first iteration.

[0083] The frequency control point initialization unit 604 is configured to acquire an initialization signal of the multi-degree-of-freedom FM signal, and determine initialization frequency control points of the initialization signal based on the initialization signal.

[0084] The iteration unit 605 is configured to calculate fitness of the initialization frequency control points, the Lagrange multiplier of second iteration and the offset of the second iteration based on the characteristic parameters of the first iteration and the initialization frequency control points by using the algorithm model. The iteration unit 605 is further configured to determine frequency control points of the second iteration from the initialization frequency control points based on the fitness by using a selection strategy. The iteration unit 605 is further configured to determine frequency control points of third iteration based on the frequency control points of the second iteration, the Lagrange multiplier of the second iteration and the offset of the second iteration by using the algorithm model and the selection strategy. The iteration unit 605 is further configured to continuously perform the iterations of the Lagrange multiplier and the offset based on the algorithm model and the selection strategy, until the augmented Lagrangian genetic algorithm converges to obtain an optimized multi-degree-of-freedom FM signal. Signal transmission is performed by adopting the optimized multi-degree-of-freedom frequency-modulation signal.

[0085] Specifically, the selection strategy includes roulette wheel selection (RWS) rule and crossover and mutation rule.

[0086] It should be understood by those skilled in the art that, the realized functions of the units in the device for designing and optimizing the multi-degree-of-freedom FM signal as illustrated in FIG. 6 can be understood with reference to the related description of the preceding method for designing and optimizing the multi-degree-of-freedom FM signal. The functions of the units in the device for designing and optimizing the multi-degree-of-freedom FM signal as illustrated in FIG. 6 may be realized by programs running on a processor, and may also be realized by a specific logical circuit.

[0087] The embodiments of the disclosure further provide a computer-readable storage medium having stored thereon computer-executable instructions, which when being executed by a processor, to execute the operations of the above method.

[0088] The above is merely the embodiments of the disclosure and is not intended to limit the scope of protection of the disclosure.