CN113221059B - Fast conjugate gradient direction finding algorithm without constructing covariance matrix - Google Patents
Fast conjugate gradient direction finding algorithm without constructing covariance matrix Download PDFInfo
- Publication number
- CN113221059B CN113221059B CN202010721858.2A CN202010721858A CN113221059B CN 113221059 B CN113221059 B CN 113221059B CN 202010721858 A CN202010721858 A CN 202010721858A CN 113221059 B CN113221059 B CN 113221059B
- Authority
- CN
- China
- Prior art keywords
- vector
- signal
- equation
- initial value
- subspace
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
- G06F2218/10—Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention belongs to the technical field of array signal processing, and particularly relates to a fast conjugate gradient direction finding algorithm which can remarkably reduce the calculation complexity without constructing a covariance matrix under the condition of not losing the precision.
Description
The technical field is as follows:
the invention belongs to the technical field of array signal processing, and particularly relates to a fast conjugate gradient direction finding algorithm which can remarkably reduce the computational complexity without constructing a covariance matrix under the condition of not losing the precision.
Background art:
estimation of the direction of arrival of a signal is an important research topic often encountered in applications such as radar, sonar, wireless communication, and passive positioning. The subspace algorithm represented by the multiple signal classification and rotation invariant subspace is provided, the leap of the traditional spatial spectrum estimation to super-resolution angle measurement is realized, but the huge calculation amount of the MUSIC algorithm and the lower estimation precision of the ESPRIT algorithm hinder the engineering progress of the super-resolution algorithm. To address this problem, the emergence of a Krylov subspace-based Conjugate Gradient (CG) algorithm was motivated.
The CG algorithm is one of the most classical Direction of Arrival (DOA) estimation methods based on Krylov subspace. Different from the traditional characteristic subspace algorithm, the CG algorithm utilizes the Krylov subspace to replace the signal subspace, and the effective algorithm of the rapid DOA estimation is realized by gradient iteration. Since the CG algorithm is at the beginning of iterationBoth the initial value setting and the iteration process need to calculate an Array Covariance Matrix (ACM), and the calculation complexity of the ACM is about O (M) 2 N), if the calculation of ACM is avoided during the gradient iteration, the computational complexity will be greatly reduced, especially at large fast beat numbers. More importantly, most existing CG algorithms maintain good estimation performance while reducing computational complexity compared to the mucci algorithm. Therefore, a Krylov subspace-based CG technique should be chosen over the conventional approach.
The algorithm enriches the theoretical connotation of CG, but needs to estimate a signal subspace by means of iteration in practical engineering application. As is well known, the setting of the initial value of the iteration affects the efficiency and accuracy of the algorithm, and improper handling may result in the iteration failing to converge or converging to a wrong solution.
The invention content is as follows:
aiming at the problem of high complexity of a CG algorithm caused by the conventional calculation of the ACM, the invention provides a fast conjugate gradient direction finding algorithm which reduces the calculation complexity of the CG algorithm by replacing the ACM with a received data matrix through improving a wienerhoff equation and realizes the purpose of remarkably reducing the calculation complexity without constructing a covariance matrix under the condition of not losing the precision.
The invention can be achieved by the following measures:
a fast conjugate gradient direction finding algorithm without constructing a covariance matrix is characterized in that after a radiation source signal is received, a wiener Hough equation is improved by using a least square method to obtain a least square wiener Hough equation, the least square wiener Hough equation is a normal equation, a reference signal only related to received data is constructed by using the normal equation, an iteration initial value is set, conjugate gradient iteration is carried out, a signal subspace is estimated according to an information source number, and a signal wave arrival direction is obtained.
The method for improving the wiener Hough equation by using the least square method specifically comprises the following steps of:
step 1-1: wienerhoff equation RW = r xd Simplifying the wienerhoff equation, R and R xd Respectively expressed as:wherein X is a matrix of received data signals, and has X = [ X (0), X (1), \ 8230;, X (M-1)]D is reference signal, d = [ d (0), d (1), \8230;, d (M-1)];
The above formula is substituted into the wienerhoff equation
After simplification, X is obtained H W=d H ;
Step 1-2: defining a least square wienerhoff equation as r + X H W=d H Wherein r is a least square residual vector, and the least square wiener Hough equation is a normal equation;
step 1-3: to avoid computing the array covariance matrix R and to obtain an improved transition vector from the least squares wiener Hoff equation of v = X H p;
Step 1-4: defining the intermediate variable normal equation residual as q = Xr, where q is related only to the array received data and the reference signal.
The construction of the reference signal related to the received data only specifically comprises the following steps:
step 2-1: defining the reference signal asWherein a (theta) is a search steering vector, and theta e [ -90 DEG, 90 DEG]D is a reference signal;
step 2-2: the data matrix of the received signal of the k array element is expressed asWhen theta = theta i When i belongs to {1, \8230;, K }, the initial value q of the residual error of the normal equation cg,0 Is shown as
It can be known that the normal equation residual error initial value q cg,0 For linear combination of signal steering vectors, then q cg,0 The K-dimensional Krylov subspace is a subset of the signal subspace, denoted as
span{Q cg,K }≡span{A(θ)},
Wherein
Q cg,K ={q cg,0 ,q cg,1 ,…,q cg,K-1 };
Step 2-3: when theta is not equal to theta i I belongs to {1, \8230;, K }, and the initial value q of the residual error of the normal equation cg,0 Expressed as:
at this time, the normal equation residual vector is also a linear combination of the guide vectors, is combined by K signal vectors and a search guide vector, namely, a K + 1-dimensional Krylov subspace spanned by the normal equation residual vector at this time is also a subset of a subspace spanned by the signal guide vectors, and is expressed as
span{Q cg,K+1 }≡span{A(θ),a(θ)},
Wherein
Q cg,K+1 ={q cg,0 ,q cg,1 ,…,q cg,K-1 ,q cg,K }。
The method specifically comprises the following steps of setting an iteration initial value and carrying out conjugate gradient iteration:
step 3-1: in the setting of the initial value of the conjugate gradient, the initial weight vector is: w is a 0 =0, search vector initial value: p is a radical of 1 =-r 1 = -d, where d is a reference signal;
step 3-2: the least square residual vector initial value in the initial value setting is as follows: r is 0 =d H -X H W, search vector initial value: p is a radical of 1 =q 0 ,γ 0 =||q 0 || 2 The initial value of the normal equation residual vector is:q 0 =Xr 0 ;
step 3-3: when the iteration number is k =1,2, \8230, M, an intermediate transition variable is obtained: v. of k =X H p k And determining the step size: alpha is alpha k =γ i-1 /||v k || 2 To obtain a weight vector: w is a k =w k-1 +α k p k Using the least squares residual: r is k =r k-1 -α k v k To obtain the normal equation residual: q. q.s k =Xr k Determining a conjugate direction vector:
the method for estimating the signal subspace according to the information source number and obtaining the signal direction of arrival specifically comprises the following steps:
step 4-1: given the number of sources K, only K iterations are needed to obtain the Krylov subspace equivalent to the signal subspace: w is a group of cg,K ={w cg,0 ,w cg,1 ,…,w cg,K-1 The spectrum peak search function of the MUSIC algorithm is as follows:where a (theta) is search guide vector, theta belongs to (-90 DEG, 90 DEG), and U N Is a noise subspace;
step 4-2: if the source number is known a priori, the iteration number is selected to be K =1,2, \8230;, K, the new algorithm spectral peak search function is:wherein I M Is an identity matrix, W cg,K Is a Krylov subspace equivalent to the signal subspace.
The invention receives the radiation source signal and expresses as follows: assuming that L mutually independent array elements form a Uniform Linear Array (ULA) at equal intervals d, K far-field narrow-band signals in space are considered to be incident to the array. Wherein, assuming that K is known a priori, d satisfies d ≦ λ 2 to avoid phase ambiguity, λ is the wavelength of the narrowband signal, and the array receiving radiation source signal is:where a (θ) is an array flow pattern matrix of dimension L × K, S (t) is an incident signal vector of dimension K × 1, N (t) is an additive white gaussian noise vector of dimension L × 1, α (θ) is a column vector of a (θ), which can be expressed as:
the method realizes that the problem of high complexity caused by setting an iteration initial value and calculating the ACM in the iteration process is avoided by improving the wiener Hough equation and replacing the ACM with the received data matrix, and provides technical support for the engineering realization of the direction of arrival.
Description of the drawings:
FIG. 1 is a flow chart of the present invention.
FIG. 2 is a comparison graph of the spectral peak search of the present invention and the MUSIC algorithm, CG algorithm, where M =10, SNR =10dB, N =100, K =2, θ = 1 =20°,θ 2 =40°。
FIG. 3 is a graph of RMSE as a function of input SNR, where M =10, the signal-to-noise ratio is varied from-20 dB to 0dB, N =200, K =2, and θ, in 5dB steps 1 =10°,θ 2 =30°。
FIG. 4 is a graph comparing the computational efficiency of the present invention and the MUSIC algorithm, CG algorithm, RV-CG algorithm, where M is from 10 to 70, the step size is 10,SNR =10dB, N =200, K =2, θ = 1 =10°,θ 2 =30°。
The specific implementation mode is as follows:
aiming at the problem that initial values and an iteration process both need to calculate ACM to cause high algorithm complexity, the invention provides a method which does not need to construct a covariance matrix DOA estimator, replaces the ACM by a received data matrix through improving a wiener Hough equation to reduce the calculation complexity of a CG algorithm, and theoretical analysis and test results show that the method can remarkably reduce the calculation complexity without losing precision and can effectively avoid wrong solutions in the iteration process, thereby providing theoretical reference for actual engineering of CG.
As shown in the attached figure 1, the invention relates to a new conjugate gradient method based on a Krylov subspace without constructing a covariance matrix, which comprises the following specific steps:
a first step of receiving a radiation source signal using an antenna array, the first step comprising the steps of: (1) Assuming that L mutually independent array elements form a Uniform Linear Array (ULA) at equal intervals d, K far-field narrow-band signals in space are considered to be incident to the array. Wherein, assuming that K is known a priori, d is satisfied to be ≦ λ/2 to avoid phase ambiguity, λ being the wavelength of the narrowband signal. The array receives the radiation source signals as follows:
where a (θ) is an array flow pattern matrix of dimension L × K, S (t) is an incident signal vector of dimension K × 1, N (t) is an additive white gaussian noise vector of dimension L × 1, α (θ) is a column vector of a (θ), which can be expressed as:
and a second step of improving the wienerhoff equation by using a least square method, wherein the second step comprises the following steps:
(1) Wienerhoff equation RW = r xd The classical conjugate gradient solution process is:
first, setting an initial weight vector as w 0 =0, initial value of search vector is p 1 =-r 1 = d, when k =1,2 \ 8230, M, the step size is determined as:the transition vector is: v = Rp k And updating the weight vector: w is a k =w k-1 +α k p k And obtaining a gradient vector: r is k+1 =r k +α k Rd k Determining a conjugate direction vector:searching for a direction vector p k Is R-conjugated, i.e. for all i ≠ j, there is
Memory matrix P k =[p 1 ,p 2 ,…,p k ]Is an M x k dimensional matrix composed of search direction vectors, then P k Is column full rank. The conjugate gradient algorithm obtained by the analysis can be converged to a real solution in the k +1 gradient direction and P step at most by M steps of iteration k Spanned k-dimensional space orthogonality:
after k iterations, the negative gradient direction vector r 1 ,r 2 ,…,r k Are orthogonal to each other. Search direction p 1 ,p 2 ,…,p k Spanned space, spanned space in gradient direction, matrix R and initial gradient vector R 1 The Krylov subspace formed is equal: thereby obtaining a signal subspace or a noise subspace, and thus performing DOA estimation.
(2) Reduction of the wienerhoff equation, R and R xd Can be respectively expressed as:
wherein X is a matrix of received data signals, and has X = [ X (0), X (1), \ 8230;, X (M-1)]D is reference signal, and has d = [ d (0), d (1), \8230;, d (M-1)]。
Therefore, the above formula is introduced into the wienerhoff equationAfter simplification, X is obtained H W=d H 。
(3) At the moment, X in the wienerhoff equation does not meet the structure of the Hermite conjugate matrix, so in order to meet the requirement on conjugate gradient iteration in the Krylov subspace, the least square wienerhoff equation is defined to be r + X H W=d H Wherein r is a least square residual vector, and the least square wiener Hough equation is a normal equation.
(4) In the conventional conjugate gradient algorithm, the transition vector is defined as v = Rp,
wherein p is a search vector, to avoid calculating the array covariance matrix R, and to obtain an improved transition vector according to the least squares wiener hough equation, v = X H p。
(5) Defining the residual of the normal equation of the intermediate variable as q = Xr, wherein q is only related to the array received data and the reference signal.
A third step of constructing a reference signal related to only the received data, the third step comprising the steps of:
(1) Defining the reference signal asWhere a (theta) is a search steering vector, theta ∈ [ -90 DEG, 90 DEG]D is a reference signal;
(2) The data matrix of the received signal of the kth array element is expressed asWhen theta = theta i When i belongs to {1, \8230;, K }, the initial value q of the residual error of the normal equation cg,0 Is shown as
From the above formula, it can be seen that the residual error initial value q of the normal equation cg,0 For linear combination of signal steering vectors, then q cg,0 The K-dimensional Krylov subspace is a signal subspaceIs represented as
span{Q cg,K Is ≡ span { A (θ) }, where Q cg,K ={q cg,0 ,q cg,1 ,…,q cg,K-1 }。
(3) When theta is not equal to theta i I ∈ {1, \8230;, K }, the initial value q of the residual error of the normal equation cg,0 Expressed as:
the above equation shows that the normal equation residual vector is also a linear combination of the guide vectors, and is a combination of K signal vectors and the search guide vector, that is, the K + 1-dimensional Krylov subspace spanned by the normal equation residual vector is also a subset of the subspace spanned by the signal guide vectors, and is expressed as
span{Q cg,K+1 Is { identical to ] span { A (theta), a (theta) }, wherein Q cg,K+1 ={q cg,0 ,q cg,1 ,…,q cg,K-1 ,q cg,K }。
A fourth step of setting an iteration initial value and performing conjugate gradient iteration, the fourth step including the steps of:
(1) The initial weight vector in the classical conjugate gradient initial value setting is:
w 0 =0, search vector initial value: p is a radical of formula 1 =-r 1 = -d, where d is a reference signal.
(2) The initial value of the least square residual vector in the initial value setting of the new algorithm is as follows: r is a radical of hydrogen 0 =d H -X H W, search vector initial value: p is a radical of formula 1 =q 0 ,γ 0 =||q 0 || 2 The initial value of the normal equation residual vector is: q. q.s 0 =Xr 0 。
(3) When the number of iterations is k =1,2, \8230, M, an intermediate transition variable is obtained: v. of k =X H p k And determining the step size: alpha is alpha k =γ i-1 /||v k || 2 To obtain a weight vector: w is a k =w k-1 +α k p k Using the least squares residual: r is k =r k-1 -α k v k To obtain the normal equation residual: q. q.s k =Xr k Determining the conjugate direction vector:
and a fifth step of estimating a signal subspace according to the number of the information sources and obtaining DOA, wherein the fifth step comprises the following steps:
(1) Given that the number of the sources is K, only K iterations are needed to obtain a Krylov subspace equivalent to the signal subspace: w cg,K ={w cg,0 ,w cg,1 ,…,w cg,K-1 }。
(2) The spectrum peak search function of the MUSIC algorithm is as follows:where a (theta) is search guide vector, theta belongs to (-90 DEG, 90 DEG), and U N Is the noise subspace.
(3) If the source number is known a priori, the iteration times are selected to be K =1,2, \8230, K, and the new algorithm spectral peak search function is:wherein I M Is an identity matrix, W cg,K Is a Krylov subspace equivalent to the signal subspace.
Example (b):
the performance of the invention is illustrated by example simulations as follows:
simulation conditions are as follows: assuming a ULA array type with 10 array elements and an array element spacing of d = λ 2, the direction of two incident signals is θ 1 =10 ° and θ 2 =30 °. To further evaluate the performance of the present invention, the number of monte carlo experiments was set to 500, and Root Mean Square Error (RMSE) was used as an evaluation index.
Simulation content and results:
simulation 1, the spectral peak search graphs of the present invention, the MUSIC algorithm and the CG algorithm are compared, and the results are shown in fig. 2.
Setting the initial reference signal to d in fig. 2 0 = Xa (θ)/| | Xa (θ) |, where θ ∈ (-90 °,90 °). The graph shows that the space spectrum peak of the least square CG algorithm is sharper than the space spectrum peaks of the classic CG algorithm and the MUSIC algorithm, the estimation performance is good, and the calculation complexity of the WCC-CG algorithm is lower than that of the classic CG algorithm, so that the WCC-CG algorithm can be better applied to practice.
Simulation 2, comparing the RMSE of the present invention with different algorithms with the input (Signal-to-noise ratio, SNR), shows the results in fig. 3.
It can be seen from fig. 3 that the estimation accuracy of the RV-CG algorithm and the estimation accuracy of the WCC-CG algorithm are both almost close to the estimation accuracy of the MUSIC algorithm, while the RV-CG algorithm and the WCC-CG algorithm are real-valued operations and do not need to calculate the covariance matrix of the received data and the eigenvalue decomposition, so that the computation amount of the RV-CG algorithm and the WCC-CG algorithm is relatively small under the condition that the accuracy is almost the same, and is between-20 dB and-7 dB, the accuracy of the WCC-CG algorithm is better than that of the other two algorithms, and the performance of the WCC-CG algorithm is poorer between-5 dB and 0dB, but is within an acceptable range. Through comparison of the three algorithms, the RV-CG algorithm and the WCC-CG algorithm have certain advantages, conditions required by actual engineering are further met, and the new algorithm is suitable for the conditions of low signal-to-noise ratio and small fast beat number.
And 3, simulating by using different array elements to compare the calculation efficiency of the method with different algorithms. By running MATLAB code in the same environment, the simulation results are given by the PC side with AMD Ryzen 5 3500U with Radon Vega Mobile Gfx 2.10GHz and 8.00 GBRAM. The equivalent evaluation of the calculation efficiency from the viewpoint of CPU time is performed, and the result is shown in fig. 4.
In fig. 4, the calculation time of the algorithm for searching spectral peaks is removed, and only the calculation of the reference signal and the process of obtaining the signal subspace through gradient iteration are recorded. It can be seen from the figure that the computation speed of the CG algorithm and the improved CG algorithm is faster than that of the MUSIC algorithm, and the greater the number of array elements is, the more obvious the advantages of the algorithm are. The WCC-CG algorithm has the fastest calculation speed, and the RV-CG algorithm and the CG algorithm are arranged in the second place.
The invention realizes that the reception data matrix is used for replacing ACM by improving the wienerhoff equation, avoids the problem of high complexity caused by setting an iteration initial value and calculating the ACM in the iteration process, and provides technical support for the engineering realization of the direction of arrival.
Claims (4)
1. A fast conjugate gradient direction-finding algorithm without constructing a covariance matrix is characterized in that after a radiation source signal is received, a wienerhoff equation is improved by using a least square method to obtain a least square wienerhoff equation, the least square wienerhoff equation is a normal equation, a reference signal only related to received data is constructed by using the normal equation, an iteration initial value is set, conjugate gradient iteration is carried out, a signal subspace is estimated according to the number of information sources, and the direction of arrival of the signal is obtained; the method for improving the wienerhoff equation by using the least square method specifically comprises the following steps:
step 1-1: wienerhoff equation RW = r xd Simplifying the wienerhoff equation, R and R xd Respectively expressed as:wherein X is a matrix of received data signals, and has X = [ X (0), X (1), \ 8230;, X (M-1)]D is reference signal, and has d = [ d (0), d (1), \8230;, d (M-1)];
The above formula is introduced into the wienerhoff equation, one
After simplification, X is obtained H W=d H ;
Step 1-2: defining a least square wienerhoff equation as r + X H W=d H Wherein r is a least square residual vector, and the least square wienerhoff equation is a normal equation;
step 1-3: to avoid computing the array covariance matrix R and to obtain an improved transition vector from the least squares wiener Hoff equation of v = X H p;
Step 1-4: defining the residual error of the intermediate variable normal equation as q = Xr, wherein q is only related to array received data and a reference signal;
the construction of the reference signal relating only to the received data comprises in particular the steps of:
step 2-1: defining the reference signal asWherein a (theta) is a search steering vector, and theta e [ -90 DEG, 90 DEG]D is a reference signal;
step 2-2: the data matrix of the received signal of the kth array element is expressed asWhen theta = theta i When i belongs to {1, \8230;, K }, the initial value q of the residual error of the normal equation cg,0 Is shown as
It can be known that the normal equation residual error initial value q cg,0 For linear combination of signal steering vectors, then q cg,0 The K-dimensional Krylov subspace is a subset of the signal subspace, denoted as
span{Q cg,K }≡span{A(θ)},
Wherein
Q cg,K ={q cg,0 ,q cg,1 ,…,q cg,K-1 };
Step 2-3: when theta is not equal to theta i I belongs to {1, \8230;, K }, and the initial value q of the residual error of the normal equation cg,0 Expressed as:
at this time, the normal equation residual vector is also a linear combination of the guide vectors, is combined by K signal vectors and a search guide vector, namely, a K + 1-dimensional Krylov subspace spanned by the normal equation residual vector at this time is also a subset of a subspace spanned by the signal guide vectors, and is expressed as
span{Q cg,K+1 }≡span{A(θ),a(θ)},
Wherein
Q cg,K+1 ={q cg,0 ,q cg,1 ,…,q cg,K-1 ,q cg,K }。
2. The fast conjugate gradient direction-finding algorithm without constructing covariance matrix as claimed in claim 1, wherein the setting of iteration initial value and conjugate gradient iteration specifically comprises the following steps:
step 3-1: in the setting of the initial value of the conjugate gradient, the initial weight vector is: w is a 0 =0, search vector initial value: p is a radical of 1 =-r 1 = -d, where d is a reference signal;
step 3-2: the least square residual vector initial value in the initial value setting is as follows: r is a radical of hydrogen 0 =d H -X H W, search vector initial value: p is a radical of 1 =q 0 ,γ 0 =||q 0 || 2 The initial value of the normal equation residual vector is: q. q of 0 =Xr 0 ;
Step 3-3: when the number of iterations is k =1,2, \8230, M, an intermediate transition variable is obtained: v. of k =X H p k And determining the step size: alpha (alpha) ("alpha") k =γ i-1 /||v k || 2 And obtaining a weight vector: w is a k =w k-1 +α k p k Using the least squares residual: r is a radical of hydrogen k =r k-1 -α k v k Obtaining a normal equation residual: q. q.s k =Xr k Determining the conjugate direction vector:
3. the fast conjugate gradient direction-finding algorithm without constructing covariance matrix as claimed in claim 1, wherein said estimating signal subspace according to source number and obtaining signal direction of arrival specifically comprises the following steps:
step 4-1: given that the number of the sources is K, only K iterations are needed to obtain a Krylov subspace equivalent to the signal subspace: w is a group of cg,K ={w cg,0 ,w cg,1 ,…,w cg,K-1 The spectrum peak search function of the MUSIC algorithm is as follows:where a (theta) is search guide vector, theta belongs to (-90 DEG, 90 DEG), and U N Is a noise subspace;
step 4-2: if the source number is known a priori, the iteration times are selected to be K =1,2, \8230, K, and the new algorithm spectral peak search function is:
wherein I M Is an identity matrix, W cg,K Is a Krylov subspace equivalent to the signal subspace.
4. A fast conjugate gradient direction finding algorithm without constructing covariance matrix as claimed in claim 1 wherein the received radiation source signal is represented as follows: assuming that L mutually independent array elements form a Uniform Linear Array (ULA) by equal spacing d, considering that K far-field narrow-band signals exist in a space and enter the array, wherein K is known a priori, d satisfies d ≦ λ 2 to avoid phase ambiguity, λ is the wavelength of the narrow-band signals, and the array receiving radiation source signals are as follows:where a (θ) is an array flow pattern matrix of dimension L × K, S (t) is an incident signal vector of dimension K × 1, N (t) is an additive white gaussian noise vector of dimension L × 1, α (θ) is a column vector of a (θ), and is expressed as:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010721858.2A CN113221059B (en) | 2020-07-24 | 2020-07-24 | Fast conjugate gradient direction finding algorithm without constructing covariance matrix |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010721858.2A CN113221059B (en) | 2020-07-24 | 2020-07-24 | Fast conjugate gradient direction finding algorithm without constructing covariance matrix |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113221059A CN113221059A (en) | 2021-08-06 |
CN113221059B true CN113221059B (en) | 2023-01-17 |
Family
ID=77085794
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010721858.2A Active CN113221059B (en) | 2020-07-24 | 2020-07-24 | Fast conjugate gradient direction finding algorithm without constructing covariance matrix |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113221059B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102353947A (en) * | 2011-07-08 | 2012-02-15 | 哈尔滨工程大学 | Method for estimating target echo signal subspaces of passive radars based on CSA-MWF (correlation subtraction algorithm-multistage wiener filter) |
CN105699950A (en) * | 2016-04-22 | 2016-06-22 | 西安电子科技大学 | Radar clutter suppression method based on self-adaptive iteration forward and background smooth conjugate gradient |
CN106443594A (en) * | 2016-08-30 | 2017-02-22 | 西安电子科技大学 | Radar antenna array steady beam forming method based on sparse constraint |
CN108981721A (en) * | 2018-07-24 | 2018-12-11 | 浙江大学 | A kind of static infrared earth sensor target angle that determining appearance for micro-nano satellite determines method |
-
2020
- 2020-07-24 CN CN202010721858.2A patent/CN113221059B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102353947A (en) * | 2011-07-08 | 2012-02-15 | 哈尔滨工程大学 | Method for estimating target echo signal subspaces of passive radars based on CSA-MWF (correlation subtraction algorithm-multistage wiener filter) |
CN105699950A (en) * | 2016-04-22 | 2016-06-22 | 西安电子科技大学 | Radar clutter suppression method based on self-adaptive iteration forward and background smooth conjugate gradient |
CN106443594A (en) * | 2016-08-30 | 2017-02-22 | 西安电子科技大学 | Radar antenna array steady beam forming method based on sparse constraint |
CN108981721A (en) * | 2018-07-24 | 2018-12-11 | 浙江大学 | A kind of static infrared earth sensor target angle that determining appearance for micro-nano satellite determines method |
Also Published As
Publication number | Publication date |
---|---|
CN113221059A (en) | 2021-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109061554B (en) | Target arrival angle estimation method based on dynamic update of spatial discrete grid | |
CN108375751B (en) | Multi-source direction-of-arrival estimation method | |
CN111123192B (en) | Two-dimensional DOA positioning method based on circular array and virtual extension | |
CN103902826B (en) | Method for tracking multiple moving targets under impact noise environment | |
CN109375154B (en) | Coherent signal parameter estimation method based on uniform circular array in impact noise environment | |
CN110940949B (en) | Method for estimating DOA of reciprocal array based on quantum penguin search mechanism in strong impact noise environment | |
CN110361691B (en) | Implementation method of coherent source DOA estimation FPGA based on non-uniform array | |
CN112881972B (en) | Direction-of-arrival estimation method based on neural network under array model error | |
CN110244272B (en) | Direction-of-arrival estimation method based on rank-denoising model | |
CN109239646B (en) | Two-dimensional dynamic direction finding method for continuous quantum water evaporation in impact noise environment | |
CN107238812B (en) | Robust dynamic direction finding method based on minimum gap array | |
CN113835063A (en) | Unmanned aerial vehicle array amplitude and phase error and signal DOA joint estimation method | |
CN115236584A (en) | Meter-wave radar low elevation angle estimation method based on deep learning | |
CN113759303B (en) | Gridless angle of arrival estimation method based on particle swarm optimization | |
CN109212466B (en) | Quantum dragonfly evolution mechanism-based broadband direction finding method | |
CN112731280B (en) | ESPRIT-DOA estimation method in inter-mass array mixed noise environment | |
CN109188346B (en) | Single snapshot DOA estimation method for large-scale uniform cylindrical array | |
CN105572629A (en) | Two-dimensional direction measuring method of low operation complexity and applicable to any array structure | |
CN113221059B (en) | Fast conjugate gradient direction finding algorithm without constructing covariance matrix | |
CN113093111B (en) | Uniform circular array two-dimensional coherent signal demodulation method and system based on compressed sensing and genetic algorithm | |
CN116582158A (en) | Large-scale MIMO square matrix information source number and direction of arrival joint estimation method | |
CN113203997B (en) | FPGA-based radar super-resolution direction finding method, system and application | |
CN109683128B (en) | Single-snapshot direction finding method under impact noise environment | |
CN108051773A (en) | EPUMA methods based on lid formula disk criterion estimation number of source | |
Liu et al. | Fast Root-MUSIC algorithm based on Nystrom method and spectral factorization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |