Method and device for resisting spectrum dispersion interference
阅读说明：本技术 一种抗频谱弥散干扰方法及装置 (Method and device for resisting spectrum dispersion interference ) 是由 尹园威 赵杨 刘利民 韩壮志 史林 于 2021-09-15 设计创作，主要内容包括：本发明公开了一种抗频谱弥散干扰方法及装置,所述方法包括：构建LFM回波信号的第一稀疏基矩阵,构建频谱弥散干扰信号的第二稀疏基矩阵,根据第一稀疏基矩阵和第二稀疏基矩阵构建联合稀疏基矩阵,根据联合稀疏基矩阵对LFM信号进行重构；基于第一干扰参数估计值对干扰参数进行优化,得到第二干扰参数估计值,基于第二干扰参数估计值重新计算第二稀疏基矩阵,将第二干扰参数估计值设置为下次干扰参数优化的第一干扰参数估计值；判断干扰参数进行优化的终止条件,满足终止条件,输出更新后的分数阶域频率分布估计值和LFM信号重构结果,不满足终止条件,重复执行优化参数,本发明通过对干扰参数的迭代优化,提高了频谱弥散干扰下LFM信号重构的精度。(The invention discloses a method and a device for resisting spectrum dispersion interference, wherein the method comprises the following steps: constructing a first sparse basis matrix of the LFM echo signals, constructing a second sparse basis matrix of the spectrum dispersion interference signals, constructing a combined sparse basis matrix according to the first sparse basis matrix and the second sparse basis matrix, and reconstructing the LFM signals according to the combined sparse basis matrix; optimizing the interference parameters based on the first interference parameter estimation value to obtain a second interference parameter estimation value, recalculating a second sparse basis matrix based on the second interference parameter estimation value, and setting the second interference parameter estimation value as the first interference parameter estimation value for next interference parameter optimization; judging the termination condition for optimizing the interference parameters, meeting the termination condition, outputting the updated fractional order domain frequency distribution estimation value and the LFM signal reconstruction result, and repeatedly executing the optimization parameters when the termination condition is not met.)
1. A method for resisting spectral dispersion interference, comprising:
s1, constructing a first sparse basis matrix of an LFM echo signal, constructing a second sparse basis matrix of a spectrum dispersion interference signal according to a first interference parameter estimation value, constructing a combined sparse basis matrix according to the first sparse basis matrix and the second sparse basis matrix, and reconstructing the LFM signal through a reconstruction algorithm according to the combined sparse basis matrix;
s2, optimizing the interference parameters based on the first interference parameter estimation value to obtain a second interference parameter estimation value, recalculating a second sparse basis matrix based on the second interference parameter estimation value, and setting the second interference parameter estimation value as the first interference parameter estimation value of next interference parameter optimization;
and S3, judging a termination condition for optimizing the interference parameters, outputting the updated fractional order domain frequency distribution estimation value and the LFM signal reconstruction result if the termination condition is met, and repeatedly executing the step S2 if the termination condition is not met.
2. The method of claim 1, wherein the first sparse basis matrix is constructed by a fractional Fourier transform method, the first interference parameter estimation value is provided by an STFT algorithm, and the joint sparse basis matrix is constructed by equation 1;
Ψ＝[ΨFr-S,ΨFr-J]formula 1;
therein, ΨFr-SIs a first sparse basis matrix, ΨFr-JIs a second sparse basis matrix.
3. The method according to claim 1, wherein the optimizing the interference parameter based on the first interference parameter estimation value to obtain a second interference parameter estimation value, and recalculating the second sparse basis matrix based on the second interference parameter estimation value specifically includes:
carrying out Taylor expansion through the fractional order frequency of the second sparse basis matrix to obtain a first derivative of the Taylor expansion, and obtaining an approximate representation of the frequency spectrum dispersion interference signal of each slice shown in formula 2 through the first derivative of the Taylor expansion;
sJ＝(ΨF-J+ΔΨF-J) θ + n equation 2;
therein, ΨF-JDictionary matrix, Δ Ψ, for a certain interference slice of the second sparse basis matrixF-JIs an M multiplied by N dimensional matrix, M and N respectively represent the row number and the column number of the sparse dictionary of the current interference slice, theta is a sparse projection vector, and N belongs to [ -N, N]；
Recalculating the second sparse basis matrix by converting equation 2 to the constraint solving problem shown in equation 3;
wherein e ═ δ m/σm，n/σn]TA joint error matrix, Ψ, representing the second sparse basis matrix and the noiseFr-JIs a second sparse basis matrix, ΨΛ Fr-JDenotes ΨFr-JLet column Λ, θΛTaking the Λ element, A, for the sparse projection vector θθ＝[σmΦm,σnIM]。
4. The method according to claim 1, wherein the determining the termination condition for the interference parameter optimization specifically includes:
and taking the change quantity of the projection position of the second sparse basis matrix in the interference parameter optimization process as a termination condition for judging the interference parameter optimization.
5. An anti-spectrum dispersion interference device is characterized by specifically comprising:
the reconstruction module is used for constructing a first sparse basis matrix of the LFM echo signal, constructing a second sparse basis matrix of the spectrum dispersion interference signal according to the first interference parameter estimation value, constructing a combined sparse basis matrix according to the first sparse basis matrix and the second sparse basis matrix, and reconstructing the LFM signal through a reconstruction algorithm according to the combined sparse basis matrix;
the interference parameter optimization module is used for optimizing the interference parameters based on the first interference parameter estimation value to obtain a second interference parameter estimation value, recalculating a second sparse basis matrix based on the second interference parameter estimation value, and setting the second interference parameter estimation value as the first interference parameter estimation value of the next interference parameter optimization;
and the judging module is used for judging the termination condition for optimizing the interference parameters, meeting the termination condition, outputting the updated fractional order domain frequency distribution estimation value and the LFM signal reconstruction result, and repeatedly executing the interference parameter optimizing module if the termination condition is not met.
6. The apparatus according to claim 5, wherein the reconstruction module specifically comprises:
the first sparse basis matrix is constructed by a fractional Fourier transform method;
the first interference parameter estimation value is provided by an STFT algorithm;
the joint sparse basis matrix is constructed through a formula 4;
Ψ＝[ΨFr-S,ΨFr-J]formula 4;
therein, ΨFr-SIs a first sparse basis matrix, ΨFr-JIs a second sparse basis matrix.
7. The apparatus of claim 5, wherein the interference parameter optimization module specifically comprises:
carrying out Taylor expansion through the fractional order frequency of the second sparse basis matrix to obtain a first derivative of the Taylor expansion, and obtaining an approximate representation of the frequency spectrum dispersion interference signal of each slice shown in formula 5 through the first derivative of the Taylor expansion;
sJ＝(ΨF-J+ΔΨF-J) θ + n equation 5;
therein, ΨF-JDictionary matrix, Δ Ψ, for a certain interference slice of the second sparse basis matrixF-JIs an M multiplied by N dimensional matrix, M and N respectively represent the row number and the column number of the sparse dictionary of the current interference slice, theta is a sparse projection vector, and N belongs to [ -N, N]；
Recalculating the second sparse basis matrix by converting equation 5 to the constraint solving problem shown in equation 5;
wherein e ═ δ m/σm，n/σn]TA joint error matrix, Ψ, representing the second sparse basis matrix and the noiseFr-JIs a second sparse basis matrix, ΨΛ Fr-JDenotes ΨFr-JLet column Λ, θΛTaking the Λ element, A, for the sparse projection vector θθ＝[σmΦm,σnIM]。
8. The apparatus according to claim 5, wherein the determining module specifically includes:
and taking the change quantity of the projection position of the second sparse basis matrix in the interference parameter optimization process as a termination condition for judging the interference parameter optimization.
9. An apparatus for resisting spectral dispersion interference, comprising: memory, processor and computer program stored on the memory and executable on the processor, which computer program, when executed by the processor, carries out the steps of the method of spectral dispersion interference rejection according to any one of claims 1 to 4.
10. A computer-readable storage medium, on which an information transfer implementing program is stored, which, when being executed by a processor, implements the steps of the method according to any one of claims 1 to 4.
Technical Field
The invention relates to the technical field of radar anti-interference, in particular to a method and a device for resisting spectrum dispersion interference.
Background
The spectrum dispersion interference is a novel interference means which is specially designed for LFM signals and is commonly used at present and has the characteristics of deception and interference suppression, and the basic mode and the characteristics are as follows: the reconnaissance receiver carries out parameter measurement on the intercepted LFM signal, then slices a complete radar signal into a plurality of segments, and modulates the frequency modulation rate of each segment of slice signal, so that each segment of signal can be expanded to cover the whole bandwidth of the LFM signal, the frequency spectrum dispersion interference signal can also obtain radar signal processing gain, and the suppression and deception effects on specific time frequency points can be realized by changing the interference power and the frequency modulation rate and the frequency spectrum dispersion interference.
The STFT algorithm is widely applied in the traditional radar signal processing process, a rapid algorithm is arranged on the software design, the hardware realization is mature, the average power of an interference signal is usually higher than or even far higher than that of a radar LFM signal, this results in an extremely low signal-to-interference ratio of the echo signal at the front end of the signal processing, the characteristics of the interfering signal dominate the echo signal, so that a coarse estimation of the interfering signal parameters using STFT is both theoretically and practically feasible, however, the fast pilot value of the frequency modulation rate obtained through the STFT usually has a certain deviation from the actual parameters of the interference signal, the magnitude of the deviation directly affects the matching degree of the constructed interference sparse basis matrix and the interference signal, and if the deviation is too large, the interference sparse basis matrix is mismatched with the real interference signal, which may cause that the CS reconstruction algorithm cannot obtain an accurate sparse projection value.
Disclosure of Invention
The invention aims to provide a method and a device for resisting spectrum dispersion interference, and aims to solve the problems in the prior art.
The invention provides a method for resisting spectrum dispersion interference, which comprises the following steps:
s1, constructing a first sparse basis matrix of an LFM echo signal, constructing a second sparse basis matrix of a spectrum dispersion interference signal according to a first interference parameter estimation value, constructing a combined sparse basis matrix according to the first sparse basis matrix and the second sparse basis matrix, and reconstructing the LFM signal through a reconstruction algorithm according to the combined sparse basis matrix;
s2, optimizing the interference parameters based on the first interference parameter estimation value to obtain a second interference parameter estimation value, recalculating a second sparse basis matrix based on the second interference parameter estimation value, and setting the second interference parameter estimation value as the first interference parameter estimation value of next interference parameter optimization;
and S3, judging a termination condition for optimizing the interference parameters, outputting the updated fractional order domain frequency distribution estimation value and the LFM signal reconstruction result if the termination condition is met, and repeatedly executing the step S2 if the termination condition is not met.
Further, a first sparse basis matrix is constructed by a fractional Fourier transform method, a first interference parameter estimation value is provided by an STFT algorithm, and a combined sparse basis matrix is constructed by a formula 1;
Ψ＝[ΨFr-S,ΨFr-J]formula 1;
therein, ΨFr-SIs a first sparse basis matrix, ΨFr-JIs a second sparse basis matrix.
Further, optimizing the interference parameter based on the first interference parameter estimation value to obtain a second interference parameter estimation value, and recalculating a second sparse basis matrix based on the second interference parameter estimation value, which specifically includes:
carrying out Taylor expansion through the fractional order frequency of the second sparse basis matrix to obtain a first derivative of the Taylor expansion, and obtaining an approximate representation of the interference signal of each slice shown in the formula 2 through the first derivative of the Taylor expansion;
sJ＝(ΨF-J+ΔΨF-J) θ + n equation 2;
therein, ΨF-JDictionary matrix, Δ Ψ, for a certain interference slice of the second sparse basis matrixF-JIs an M multiplied by N dimensional matrix, M and N respectively represent the row number and the column number of the sparse dictionary of the current interference slice, theta is a sparse projection vector, and N belongs to [ -N, N]；
Recalculating the second sparse basis matrix by converting equation 2 to the constraint solving problem shown in equation 3;
wherein e ═ δ m/σm，n/σn]TA joint error matrix, Ψ, representing the second sparse basis matrix and the noiseFr-JIs a second sparse basis matrix, ΨΛ Fr-JDenotes ΨFr-JLet column Λ, θΛTaking the Λ element, A, for the sparse projection vector θθ＝[σmΦm,σnIM]。
Further, the determining of the termination condition for optimizing the interference parameter specifically includes:
and taking the change quantity of the projection position of the second sparse basis matrix in the interference parameter optimization process as a termination condition for judging the interference parameter optimization.
The invention provides an anti-spectrum dispersion interference device, which comprises:
the reconstruction module is used for constructing a first sparse basis matrix of the LFM echo signal, constructing a second sparse basis matrix of the spectrum dispersion interference signal according to the first interference parameter estimation value, constructing a combined sparse basis matrix according to the first sparse basis matrix and the second sparse basis matrix, and reconstructing the LFM signal through a reconstruction algorithm according to the combined sparse basis matrix;
the interference parameter optimization module is used for optimizing the interference parameters based on the first interference parameter estimation value to obtain a second interference parameter estimation value, recalculating a second sparse basis matrix based on the second interference parameter estimation value, and setting the second interference parameter estimation value as the first interference parameter estimation value of the next interference parameter optimization;
and the judging module is used for judging the termination condition for optimizing the interference parameters, meeting the termination condition, outputting the updated fractional order domain frequency distribution estimation value and the LFM signal reconstruction result, and repeatedly executing the interference parameter optimizing module if the termination condition is not met.
Further, the reconstruction module specifically includes:
the first sparse basis matrix is constructed by a fractional Fourier transform method;
the first interference parameter estimation value is provided by an STFT algorithm;
the joint sparse basis matrix is constructed through a formula 4;
Ψ＝[ΨFr-S,ΨFr-J]formula 4;
therein, ΨFr-SIs a first sparse basis matrix, ΨFr-JIs a second sparse basis matrix.
Further, the interference parameter optimization module specifically includes:
carrying out Taylor expansion through the fractional order frequency of the second sparse basis matrix to obtain a first derivative of the Taylor expansion, and obtaining an approximate representation of the interference signal of each slice shown in formula 5 through the first derivative of the Taylor expansion;
sJ＝(ΨF-J+ΔΨF-J) θ + n equation 5;
therein, ΨF-JDictionary matrix, Δ Ψ, for a certain interference slice of the second sparse basis matrixF-JIs an M multiplied by N dimensional matrix, M and N respectively represent the row number and the column number of the sparse dictionary of the current interference slice, theta is a sparse projection vector, and N belongs to [ -N, N]；
Recalculating the second sparse basis matrix by converting equation 5 to the constraint solving problem shown in equation 5;
wherein e ═ δ m/σm，n/σn]TA joint error matrix, Ψ, representing the second sparse basis matrix and the noiseFr-JIs a second sparse basis matrix, ΨΛ Fr-JDenotes ΨFr-JLet column Λ, θΛTaking the Λ element, A, for the sparse projection vector θθ＝[σmΦm,σnIM]。
Further, the judging module specifically includes:
and taking the change quantity of the projection position of the second sparse basis matrix in the interference parameter optimization process as a termination condition for judging the interference parameter optimization.
The embodiment of the invention also provides a device for resisting the spectrum dispersion interference, which comprises: a memory, a processor and a computer program stored on the memory and executable on the processor, the computer program, when executed by the processor, implementing the steps of the above-described method of spectral dispersion interference rejection.
The embodiment of the invention also provides a computer-readable storage medium, wherein an implementation program for information transmission is stored on the computer-readable storage medium, and when the implementation program is executed by a processor, the steps of the method for resisting the spectrum dispersion interference are implemented.
By adopting the embodiment of the invention, the interference parameter estimation value is provided through the STFT algorithm, the combined sparse basis matrix of the LFM signal and the spectrum dispersion interference signal is constructed, the reconstruction of the LFM signal is realized, and the reconstruction precision of the LFM signal is improved through the iterative optimization of the interference parameter.
The foregoing description is only an overview of the technical solutions of the present invention, and the embodiments of the present invention are described below in order to make the technical means of the present invention more clearly understood and to make the above and other objects, features, and advantages of the present invention more clearly understandable.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a flow chart of a method for resisting spectrum dispersion interference according to an embodiment of the present invention;
FIG. 2 is a diagram of a first embodiment of the device for resisting spectrum dispersion interference;
fig. 3 is a structural diagram of an anti-spectrum-dispersion interference device according to a second embodiment of the device of the present invention.
Detailed Description
The technical solutions of the present invention will be described clearly and completely with reference to the following embodiments, and it should be understood that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In the description of the present invention, it is to be understood that the terms "center", "longitudinal", "lateral", "length", "width", "thickness", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", "clockwise", "counterclockwise", and the like, indicate orientations and positional relationships based on those shown in the drawings, and are used only for convenience of description and simplicity of description, and do not indicate or imply that the device or element being referred to must have a particular orientation, be constructed and operated in a particular orientation, and thus, should not be considered as limiting the present invention.
Furthermore, the terms "first", "second" and "first" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, features defined as "first", "second", may explicitly or implicitly include one or more of the described features. In the description of the present invention, "a plurality" means two or more unless specifically defined otherwise. Furthermore, the terms "mounted," "connected," and "connected" are to be construed broadly and may, for example, be fixedly connected, detachably connected, or integrally connected; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
Method embodiment
According to an embodiment of the present invention, a method for resisting spectrum dispersion interference is provided, fig. 1 is a flowchart of the method for resisting spectrum dispersion interference according to the embodiment of the present invention, and as shown in fig. 1, the method for resisting spectrum dispersion interference according to the embodiment of the present invention specifically includes:
step S101, a first sparse basis matrix of the LFM echo signal is constructed, a second sparse basis matrix of the spectrum dispersion interference signal is constructed according to the first interference parameter estimation value, a combined sparse basis matrix is constructed according to the first sparse basis matrix and the second sparse basis matrix, the LFM signal is reconstructed through a reconstruction algorithm according to the combined sparse basis matrix, and the step S101 specifically comprises the following steps:
a first sparse basis matrix is constructed by a Pei type discrete fractional Fourier transform method, interference parameters are estimated by an STFT algorithm to obtain the number of slices and the frequency modulation rate, and the number of slices and the frequency modulation rate are substituted into a formula 7 to complete construction of a combined sparse basis matrix.
Ψ＝[ΨFr-S,ΨFr-J]Equation 7;
therein, ΨFr-SIs a first sparse basis matrix, ΨFr-JIs a second sparse basis matrix.
On the basis, signal reconstruction is carried out by using a CS theory, the LFM echo signals are projected through a measurement matrix phi, and an L multiplied by 1 dimensional observation value vector shown in a formula 8 is obtained through measurement;
y ═ Φ s ═ Φ (Ψ θ + n) ═ a θ + z equation 8;
and a phi psi is a L multiplied by 2N-dimensional sensing matrix, N is noise, z is a noise measurement vector, and theta is a sparse projection value of the LFM signal and the spectrum dispersion interference signal. At this time, when signal reconstruction is performed by using the compressive sensing theory, the solution of the LFM echo signal s is converted into the optimal l of the solution constraint of the observation vector y shown in formula 91A norm problem;
where δ is a constant related to noise.
Step S102, optimizing the interference parameters based on the first interference parameter estimation value to obtain a second interference parameter estimation value, recalculating a second sparse basis matrix based on the second interference parameter estimation value, and setting the second interference parameter estimation value as the first interference parameter estimation value of the next interference parameter optimization, wherein the step S102 specifically comprises:
carrying out Taylor expansion through the fractional order frequency of the second sparse basis matrix to obtain a first derivative of the Taylor expansion, and obtaining an approximate representation of the spectrum dispersion interference signal of each slice shown in the formula 10 through the first derivative of the Taylor expansion;
sJ＝(ΨF-J+ΔΨF-J) θ + n equation 10;
therein, ΨF-JDictionary matrix, Δ Ψ, for a certain interference slice of the second sparse basis matrixF-JIs an M multiplied by N dimensional matrix, M and N respectively represent the row number and the column number of the sparse dictionary of the current interference slice, theta is a sparse projection vector, and N belongs to [ -N, N]；
Equation 10 reduces to equation 11;
setting a sparse projection vector theta to contain K solutions, and setting lambda as a column sequence number set of a solution space;
converting formula 11 into a constraint solving problem shown in formula 12;
wherein e ═ δ m/σm，n/σn]TA joint error matrix, Ψ, representing the second sparse basis matrix and the noiseFr-JIs a second sparse basis matrix, ΨΛ Fr-JDenotes ΨFr-JLet column Λ, θΛTaking the Λ element, A, for the sparse projection vector θθ＝[σmΦm,σnIM]。
Obtaining the values shown in formula 13 and formula 14 by formula 12Anda least squares solution of;
wherein, (g)+Representing a pseudo inverse, i.e.
Solving formula 13 requires traversing all solution spaces, and the solving difficulty is high, so psi can be established based on STFTF-JPerforming primary initial matching on the interference signal, and reserving a main peak component in the projection value theta to obtain lambda andas an initial estimation value, further to obtainAnd within ΛThe solution is performed, in which case the least squares solution of equation 13 can be represented by equation 15;
in the formula (I), the compound is shown in the specification,
extracting joint error matricesThe first k element values in (1) are then converted to δ m, which is further substituted into Δ ΨF-JTo psiF-JAnd updating to obtain the updated second sparse basis matrix.
Step S103, judging a termination condition for optimizing the interference parameter, meeting the termination condition, outputting the updated fractional order domain frequency distribution estimation value and the LFM signal reconstruction result, and repeatedly executing step S102, wherein the step S103 specifically comprises the following steps:
the error vector obtained by derivation according to the formula 15 already includes the revised value δ m of the fractional order frequency, and the reconstructed mean square error of the LFM signal has a correlation with the modulation frequency of the interference sparse basis matrix, that is, in the process of approaching to the real interference signal projection position, when δ m is not changed any more or converges to a certain range, the reconstructed error of the LFM signal will also change synchronously, so that the change amount of the second sparse basis matrix projection position in the iterative calculation process is used as the determination condition for determining whether the evaluation algorithm is terminated.
And providing an interference parameter estimation value through an STFT algorithm, constructing a combined sparse basis matrix of the LFM signal and the spectrum dispersion interference signal, performing iterative optimization on the interference parameter, and judging a termination condition of the iterative optimization of the interference parameter through a change amount of a projection position of a second sparse basis matrix to improve the reconstruction precision of the LFM signal.
Apparatus embodiment one
According to an embodiment of the present invention, there is provided an anti-spectrum dispersion interference device, fig. 2 is a structural diagram of the anti-spectrum dispersion interference device according to the embodiment of the present invention, and as shown in fig. 3, the anti-spectrum dispersion interference device according to the embodiment of the present invention specifically includes:
the reconstruction module 20 is configured to construct a first sparse basis matrix of the LFM echo signal, construct a second sparse basis matrix of the spectrum dispersion interference signal according to the first interference parameter estimation value, construct a joint sparse basis matrix according to the first sparse basis matrix and the second sparse basis matrix, and reconstruct the LFM signal according to the joint sparse basis matrix through a reconstruction algorithm, where the reconstruction module 20 is specifically configured to:
a first sparse basis matrix is constructed by a Pei type discrete fractional Fourier transform method, interference parameters are estimated by an STFT algorithm to obtain the number of slices and the frequency modulation rate, and the number of slices and the frequency modulation rate are substituted into a formula 7 to complete construction of a combined sparse basis matrix.
Ψ＝[ΨFr-S,ΨFr-J]Equation 7;
therein, ΨFr-SIs a first sparse basis matrix, ΨFr-JIs a second sparse basis matrix.
On the basis, signal reconstruction is carried out by using a CS theory, the LFM echo signals are projected through a measurement matrix phi, and an L multiplied by 1 dimensional observation value vector shown in a formula 8 is obtained through measurement;
y ═ Φ s ═ Φ (Ψ θ + n) ═ a θ + z equation 8;
and a phi psi is a L multiplied by 2N-dimensional sensing matrix, N is noise, z is a noise measurement vector, and theta is a sparse projection value of the LFM signal and the spectrum dispersion interference signal. At this time, when signal reconstruction is performed by using the compressive sensing theory, the solution of the LFM echo signal s is converted into the optimal l of the solution constraint of the observation vector y shown in formula 91A norm problem;
where δ is a constant related to noise.
The interference parameter optimization module 22 is configured to optimize the interference parameter based on the first interference parameter estimation value to obtain a second interference parameter estimation value, recalculate the second sparse basis matrix based on the second interference parameter estimation value, and set the second interference parameter estimation value as the first interference parameter estimation value for the next interference parameter optimization, where the interference parameter optimization module 22 is specifically configured to:
carrying out Taylor expansion through the fractional order frequency of the second sparse basis matrix to obtain a first derivative of the Taylor expansion, and obtaining an approximate representation of the spectrum dispersion interference signal of each slice shown in the formula 10 through the first derivative of the Taylor expansion;
sJ＝(ΨF-J+ΔΨF-J) θ + n equation 10;
therein, ΨF-JDictionary matrix, Δ Ψ, for a certain interference slice of the second sparse basis matrixF-JIs an M multiplied by N dimensional matrix, M and N respectively represent the row number and the column number of the sparse dictionary of the current interference slice, theta is a sparse projection vector, and N belongs to [ -N, N]；
Equation 10 reduces to equation 11;
setting a sparse projection vector theta to contain K solutions, and setting lambda as a column sequence number set of a solution space;
converting formula 11 into a constraint solving problem shown in formula 12;
wherein e ═ δ m/σm，n/σn]TA joint error matrix, Ψ, representing the second sparse basis matrix and the noiseFr-JIs a second sparse basis matrix, ΨΛ Fr-JDenotes ΨFr-JLet column Λ, θΛTaking the Λ element, A, for the sparse projection vector θθ＝[σmΦm,σnIM]。
Obtaining the values shown in formula 13 and formula 14 by formula 12Anda least squares solution of;
wherein, (g)+Representing a pseudo inverse, i.e.
Solving formula 13 requires traversing all solution spaces, and the solving difficulty is high, so psi can be established based on STFTF-JPerforming primary initial matching on the interference signal, and reserving a main peak component in the projection value theta to obtain lambda andas an initial estimation value, further to obtainAnd within ΛThe solution is performed, in which case the least squares solution of equation 13 can be represented by equation 15;
in the formula (I), the compound is shown in the specification,
extracting joint error matricesThe first k element values in (1) are then converted to δ m, which is further substituted into Δ ΨF-JTo psiF-JAnd updating to obtain the updated second sparse basis matrix.
The judging module 24 is configured to judge a termination condition for optimizing the interference parameter, meet the termination condition, output the updated fractional order domain frequency distribution estimation value and the LFM signal reconstruction result, and repeatedly execute the interference parameter optimizing module if the termination condition is not met, where the judging module 24 is specifically configured to:
and taking the change of the projection position of the second sparse basis matrix in the iterative calculation process as a judgment condition for judging whether the evaluation algorithm is terminated.
Device embodiment II
An embodiment of the present invention provides an apparatus for resisting spectrum dispersion interference, as shown in fig. 3, including: a memory 30, a processor 32 and a computer program stored on the memory 30 and executable on the processor 32, which computer program, when executed by the processor 32, performs the steps as described in the method embodiments.
Device embodiment III
An embodiment of the present invention provides a computer-readable storage medium, on which an implementation program for information transmission is stored, and when the program is executed by a processor 32, the steps as described in the method embodiment are implemented.
The computer-readable storage medium of this embodiment includes, but is not limited to: ROM, RAM, magnetic or optical disks, and the like.
The foregoing description has been directed to specific embodiments of this disclosure. Other embodiments are within the scope of the following claims. In some cases, the actions or steps recited in the claims may be performed in a different order than in the embodiments and still achieve desirable results. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In some embodiments, multitasking and parallel processing may also be possible or may be advantageous.
In the 30 s of the 20 th century, improvements in a technology could clearly be distinguished between improvements in hardware (e.g., improvements in circuit structures such as diodes, transistors, switches, etc.) and improvements in software (improvements in process flow). However, as technology advances, many of today's process flow improvements have been seen as direct improvements in hardware circuit architecture. Designers almost always obtain the corresponding hardware circuit structure by programming an improved method flow into the hardware circuit. Thus, it cannot be said that an improvement in the process flow cannot be realized by hardware physical modules. For example, a Programmable Logic Device (PLD), such as a Field Programmable Gate Array (FPGA), is an integrated circuit whose Logic functions are determined by programming the Device by a user. A digital system is "integrated" on a PLD by the designer's own programming without requiring the chip manufacturer to design and fabricate application-specific integrated circuit chips. Furthermore, nowadays, instead of manually making an Integrated Circuit chip, such Programming is often implemented by "logic compiler" software, which is similar to a software compiler used in program development and writing, but the original code before compiling is also written by a specific Programming Language, which is called Hardware Description Language (HDL), and HDL is not only one but many, such as abel (advanced Boolean Expression Language), ahdl (alternate Hardware Description Language), traffic, pl (core universal Programming Language), HDCal (jhdware Description Language), lang, Lola, HDL, laspam, hardward Description Language (vhr Description Language), vhal (Hardware Description Language), and vhigh-Language, which are currently used in most common. It will also be apparent to those skilled in the art that hardware circuitry that implements the logical method flows can be readily obtained by merely slightly programming the method flows into an integrated circuit using the hardware description languages described above.
The controller may be implemented in any suitable manner, for example, the controller may take the form of, for example, a microprocessor or processor and a computer-readable medium storing computer-readable program code (e.g., software or firmware) executable by the (micro) processor, logic gates, switches, an Application Specific Integrated Circuit (ASIC), a programmable logic controller, and an embedded microcontroller, examples of which include, but are not limited to, the following microcontrollers: ARC 625D, Atmel AT91SAM, Microchip PIC18F26K20, and Silicone Labs C8051F320, the memory controller may also be implemented as part of the control logic for the memory. Those skilled in the art will also appreciate that, in addition to implementing the controller as pure computer readable program code, the same functionality can be implemented by logically programming method steps such that the controller is in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers and the like. Such a controller may thus be considered a hardware component, and the means included therein for performing the various functions may also be considered as a structure within the hardware component. Or even means for performing the functions may be regarded as being both a software module for performing the method and a structure within a hardware component.
The systems, devices, modules or units illustrated in the above embodiments may be implemented by a computer chip or an entity, or by a product with certain functions. One typical implementation device is a computer. In particular, the computer may be, for example, a personal computer, a laptop computer, a cellular telephone, a camera phone, a smartphone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.
For convenience of description, the above devices are described as being divided into various units by function, and are described separately. Of course, the functions of the units may be implemented in the same software and/or hardware or in multiple software and/or hardware when implementing the embodiments of the present description.
One skilled in the art will recognize that one or more embodiments of the present description may be provided as a method, system, or computer program product. Accordingly, one or more embodiments of the present description may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the description may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The description has been presented with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the description. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include forms of volatile memory in a computer readable medium, Random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). Memory is an example of a computer-readable medium.
Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), Read Only Memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
One or more embodiments of the present description may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. One or more embodiments of the specification may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment.
The above description is only an example of this document and is not intended to limit this document. Various modifications and changes may occur to those skilled in the art from this document. Any modifications, equivalents, improvements, etc. which come within the spirit and principle of the disclosure are intended to be included within the scope of the claims of this document.