CN109194422A - A kind of SNR estimation method based on subspace - Google Patents
A kind of SNR estimation method based on subspace Download PDFInfo
- Publication number
- CN109194422A CN109194422A CN201811024519.8A CN201811024519A CN109194422A CN 109194422 A CN109194422 A CN 109194422A CN 201811024519 A CN201811024519 A CN 201811024519A CN 109194422 A CN109194422 A CN 109194422A
- Authority
- CN
- China
- Prior art keywords
- subspace
- signal
- covariance matrix
- estimation method
- dimension
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000011159 matrix material Substances 0.000 claims abstract description 38
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 12
- 238000004891 communication Methods 0.000 claims abstract description 7
- 239000013598 vector Substances 0.000 claims description 11
- 238000007476 Maximum Likelihood Methods 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 abstract description 9
- 238000012545 processing Methods 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000013459 approach Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/309—Measuring or estimating channel quality parameters
- H04B17/336—Signal-to-interference ratio [SIR] or carrier-to-interference ratio [CIR]
Landscapes
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Optical Communication System (AREA)
Abstract
The SNR estimation method based on subspace that the invention discloses a kind of, comprising the following steps: step 1: the estimated value for receiving signal and solving covariance matrix obtained by communication equipment；Step 2: solving the dimension of signal subspace；Step 3: Eigenvalues Decomposition being carried out to the estimated value for the covariance matrix that step 1 obtains, and according to signal subspace dimension descending arrayed feature value；Step 4: noise power and signal power being solved according to the characteristic value that step 3 obtains, to obtain the estimated value of SNR.A kind of the advantages of SNR estimation method based on subspace provided by the invention, is: the application use estimation method calculation amount it is smaller, processing speed faster, estimation performance be better than presently, there are other algorithms.
Description
Technical Field
The invention relates to the technical field of signal processing, in particular to an SNR estimation method based on subspace.
Background
Mars are the planets in the solar system that are most detected by humans. Due to the fact that the distance between the mars and the earth is far, signal transmission time is prolonged, and signal attenuation is large. Meanwhile, the relative distance changes rapidly, and the load and power consumption of the detector are also very limited. These real situations all cause great trouble to the ground fire communication. A short communication arc results in less effective communication time and therefore a higher effective data rate. Moreover, real-time changes of the relative distance and the relative attitude between the devices can also cause real-time changes of the signal-to-noise ratio (SNR) of the signals received by the receiving end. This requires real-time adjustment of the code rate to increase it as much as possible. SNR estimation uses parameter estimation theory to calculate the noise-to-power ratio of a signal and provides the required channel quality information through modulation mode switching, rate adaptation, power control and channel allocation methods. Furthermore, many parameter estimation algorithms require SNR as a precondition to optimize performance.
The SNR estimation method includes a time domain method and a frequency domain method. The time domain method is classified into a Data Aided (DA) method and a non-data aided (NDA) method. The DA method has higher estimation accuracy than the NDA method, but it requires insertion of a periodic preamble sequence, which reduces transmission efficiency. In the time domain method, the DA-based SNR estimation method includes a Minimum Mean Square Error (MMSE), a maximum likelihood estimation (ML), a discrete symbol matrix estimation (SSME), a signal-to-noise separation of high-order cumulants, and the like. The NDA-based SNR estimation method includes a second-order fourth-order matrix estimation method (M2M4), a Signal-to-variance Ratio estimation (SVR), a Squared Signal-to-noise variance Ratio estimation (SNV), a Data Fitting estimation (DF), and the like. The classical frequency domain method is based on the flat characteristic of a White Noise power spectrum, is suitable for the SNR estimation of an Additive White Gaussian Noise (AWGN) channel, and is not suitable for the color Noise environment.
Disclosure of Invention
The technical problem to be solved by the invention is to provide a subspace-based SNR estimation method with better performance.
The invention solves the technical problems through the following technical scheme:
a subspace-based SNR estimation method, comprising the steps of:
step 1: solving an estimated value of a covariance matrix through a received signal obtained by communication equipment;
step 2: solving the dimension of the signal subspace;
and step 3: carrying out eigenvalue decomposition on the estimated value of the covariance matrix obtained in the step 1, and arranging the eigenvalues in a descending order according to the signal subspace dimension;
and 4, step 4: and (4) solving the noise power and the signal power according to the characteristic value obtained in the step (3) so as to obtain an estimated value of the SNR.
Preferably, the method for solving the estimated value of the covariance matrix in step 1 is:
taking L receiving signals from t time forward to form a receiving signal vector r (t) at t time, and then
r(t)＝[rt,rt-1,…,rt-L+1]T
The covariance matrix can be obtained as
Rrr＝E[r(t)rH(t)]
Assuming K invariant observation vectors, the estimated value of the covariance matrix is
Where k is {0,1, …, L-1 }.
Preferably, the dimension p of the signal subspace is estimated in step 2 using a minimum description length criterion MDL,
k and the subspace dimension p satisfy the following relationship
p＝argkmin MDL(k)。
Preferably, the eigenvalue decomposition method of step 3 is:
assuming that the signal power and the noise power are uncorrelated, the covariance matrix estimate can be expressed as:
Rrr＝Rxx+Rnn
according to the eigenvalue decomposition theory, RrrCan be decomposed into
Rrr＝AΣAH
Wherein A is composed of orthogonal feature vectors; diagonal matrix Σ ═ diag (b)i) From the characteristic values b of the matrixiIs formed therein
b1≥b2≥…bL
The autocorrelation matrix of white noise satisfies
Wherein,is the variance of the noise power;
the eigenvalue of the covariance matrix obtained by descending order is
Wherein,representing the signal power of the ith eigenvector, p representing the dimension of the signal subspace, the first p eigenvectors constituting the signal subspace, and the last L-p eigenvectors constituting the noise subspace.
Using the estimated value of the covariance matrix obtained in step 1The above decomposition is carried out to obtain the characteristic value RrrIs estimated from the maximum likelihood.
Preferably, the SNR estimate is
ρ＝10lg(Ps/Pn)
Wherein,
Pnas noise power, PsIs the signal power.
The subspace-based SNR estimation method provided by the invention has the advantages that: the estimation method adopted by the application has the advantages of small calculation amount, higher processing speed and better estimation performance than other existing algorithms.
Drawings
FIG. 1 is a diagram illustrating the performance of a subspace-based SNR estimation method provided by an embodiment of the present invention for different SNR cases;
FIG. 2 is a diagram illustrating SNR estimation performance of a subspace-based SNR estimation method under different snapshot conditions according to an embodiment of the present invention;
fig. 3 is a comparison of estimated performance of different algorithms.
Detailed Description
In order that the objects, aspects and advantages of the present invention will become more apparent, the present invention will be further described in detail with reference to the accompanying drawings in which specific embodiments are shown.
In a complex AWGN channel, the received signal may be represented as
rk＝xk+nk
Wherein nk is complex white Gaussian noise with zero mean value under sequential sampling, and its variance isThe noise column vector consisting of L times of sampling data before t time can be expressed as
n(t)＝[nt,nt-1,…,nt-L+1]T
Wherein (C)TRepresenting transpose, the second moment of the noise vector can be represented as
Where E [ ] represents the mathematical expectation, the superscript H represents the conjugate transpose, and I is an L identity matrix.
xkAt a carrier frequency fcThe lower band pass signal is represented as
Where θ represents the initial phase, fsRepresenting the sampling frequency.
Is a complex equivalent baseband signal that, for a Multiple Phase Shift Keying (MPSK) signal,
wherein A ∈ R+Representing the amplitude of the signal.
In performing a Multiple Quadrature Amplitude Modulation (MQAM) signal,
wherein A islAnd e C represents the amplitude corresponding to the l-th phase.
Let PsAnd PnRepresenting signal power and noise power, respectively, the SNR can be expressed as
ρ＝10lg(Ps/Pn)
Based on the above analysis, the emphasis of SNR estimation is to find the signal power PsSum noise power PnThe subspace-based SNR method provided by the invention comprises the following steps:
step 1: solving an estimated value of a covariance matrix through a received signal obtained by communication equipment;
taking L receiving signals from t time forward to form a receiving signal vector r (t) at t time, and then
r(t)＝[rt,rt-1,…,rt-L+1]T
The covariance matrix can be obtained as
Rrr＝E[r(t)rH(t)]
In practice, a covariance matrix R is obtainedrrIs very difficult and therefore we can only estimate with a finite length of the received signal.
Assuming K invariant observation vectors, the estimated value of the covariance matrix is
Where k is {0,1, …, L-1 }.
Step 2: solving the dimension of the signal subspace;
the minimum description length criterion MDL is heard in the preferred embodiment to estimate the number of dimensions p of the signal subspace,
k and the subspace dimension p satisfy the following relationship
p＝argkmin MDL(k)。
And step 3: carrying out eigenvalue decomposition on the estimated value of the covariance matrix obtained in the step 1, and arranging the eigenvalues in a descending order according to the signal subspace dimension;
assuming that the signal power and the noise power are uncorrelated, the covariance matrix estimate can be expressed as:
Rrr＝Rxx+Rnn
according to the eigenvalue decomposition theory, RrrCan be decomposed into
Rrr＝AΣAH
Wherein A is composed of orthogonal feature vectors; diagonal matrix Σ ═ diag (b)i) From the characteristic values b of the matrixiIs formed therein
b1≥b2≥…bL
The autocorrelation matrix R of white noisennSatisfy the requirement of
Wherein,is the variance of the noise power;
the eigenvalue of the covariance matrix obtained by descending order is
Wherein,representing the signal power of the ith eigenvector, P representing the dimension of the signal subspace, the first P eigenvectors constituting the signal subspace, the last L-P eigenvectors constituting the noise subspace, the noise power P of whichnFrom L to pWhen the dimension p of the signal subspace is constant, the SNR estimation can be passedAnd PnTo obtain an estimate of (c).
Using the estimated value of the covariance matrix obtained in step 1By performing the decomposition, the characteristic value R can be obtainedrrIs estimated from the maximum likelihood.
And 4, step 4: and (4) solving the noise power and the signal power according to the characteristic value obtained in the step (3) so as to obtain an estimated value of the SNR.
It can be known from the previous analysis
ρ＝10lg(Ps/Pn)
Wherein,
Pnas noise power, PsIs the signal power.
The calculated quantity of the above algorithm is O (L)2N+L3) A sub-multiplication, where N is the number of symbols.
The performance of the algorithm provided by the embodiment of the invention is analyzed by using MATLAB simulation, and the result is as follows:
FIG. 1 is an estimation performance of a subspace-based SNR estimation method under different SNR cases; it can be seen from fig. 1 that the subspace approach has better SNR estimation performance. As the SNR increases, the estimate of the SNR can approach the true value.
FIG. 2 is a diagram of SNR estimation performance for different snapshot scenarios for a subspace-based SNR estimation method; it can be seen that the estimation of SNR gets better and better as the number of snapshots increases.
Fig. 3 is a comparison of the estimated performance of the different algorithms (N-300); comparison was made with DF, M2M4, SVR and SNV algorithms. It can be seen from the figure that the estimation performance of the subspace approach is better than other algorithms.
The above-mentioned embodiments, objects, technical solutions and advantages of the present invention are further described in detail, it should be understood that the above-mentioned embodiments are only examples of the present invention, and are not intended to limit the present invention, and any modifications, equivalent substitutions, improvements, etc. made by those skilled in the art without departing from the spirit and principles of the present invention should fall within the protection scope defined by the claims of the present invention.
Claims (5)
1. A subspace-based SNR estimation method, characterized in that: the method comprises the following steps:
step 1: solving an estimated value of a covariance matrix through a received signal obtained by communication equipment;
step 2: solving the dimension of the signal subspace;
and step 3: carrying out eigenvalue decomposition on the estimated value of the covariance matrix obtained in the step 1, and arranging the eigenvalues in a descending order according to the signal subspace dimension;
and 4, step 4: and (4) solving the noise power and the signal power according to the characteristic value obtained in the step (3), thereby obtaining an estimated value of the SNR.
2. The subspace-based SNR estimation method according to claim 1, wherein: the method for solving the estimation value of the covariance matrix in the step 1 comprises the following steps:
taking L receiving signals from t time forward to form a receiving signal vector r (t) at t time, and then
r(t)＝[rt,rt-1,…,rt-L+1]T
The covariance matrix can be obtained as
Rrr＝E[r(t)rH(t)]
Assuming K invariant observation vectors, the estimated value of the covariance matrix is
Where k is {0,1, …, L-1 }.
3. A subspace-based SNR estimation method according to claim 2, characterized in that: the dimension p of the signal subspace is estimated in step 2 using the minimum description length criterion MDL,
k and the subspace dimension p satisfy the following relationship
p＝argkmin MDL(k)。
4. A subspace-based SNR estimation method according to claim 3, characterized in that: the eigenvalue decomposition method in step 3 comprises the following steps:
assuming that the signal power and the noise power are uncorrelated, the covariance matrix estimate can be expressed as:
Rrr＝Rxx+Rnn
decomposition according to eigenvaluesTheory, RrrCan be decomposed into
Rrr＝AΣAH
Wherein A is composed of orthogonal feature vectors; diagonal matrix Σ ═ diag (b)i) From the eigenvalues b of the matrixiIs formed therein
b1≥b2≥…bL
The autocorrelation matrix of white noise satisfies
Wherein,is the variance of the noise power;
the eigenvalue of the covariance matrix obtained by descending order is
Wherein,representing the signal power of the ith eigenvector, p representing the dimension of the signal subspace, the first p eigenvectors constituting the signal subspace, and the last L-p eigenvectors constituting the noise subspace.
Using the estimated value of the covariance matrix obtained in step 1Performing the decomposition to obtain a characteristic value RrrIs estimated from the maximum likelihood.
5. The subspace-based SNR estimation method according to claim 4, wherein: SNR estimation value of
ρ＝10lg(Ps/Pn)
Wherein,
Pnas noise power, PsIs the signal power.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811024519.8A CN109194422B (en) | 2018-09-04 | 2018-09-04 | SNR estimation method based on subspace |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811024519.8A CN109194422B (en) | 2018-09-04 | 2018-09-04 | SNR estimation method based on subspace |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109194422A true CN109194422A (en) | 2019-01-11 |
| CN109194422B CN109194422B (en) | 2021-06-22 |
Family
ID=64912133
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811024519.8A Active CN109194422B (en) | 2018-09-04 | 2018-09-04 | SNR estimation method based on subspace |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109194422B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111294126A (en) * | 2020-03-13 | 2020-06-16 | 南京航空航天大学 | Subspace tracking method for SNR estimation |
| CN113364539A (en) * | 2021-08-09 | 2021-09-07 | 成都华日通讯技术股份有限公司 | Blind estimation method for signal-to-noise ratio of digital signal in frequency spectrum monitoring equipment |
| WO2023246915A1 (en) * | 2022-06-21 | 2023-12-28 | 锐捷网络股份有限公司 | Method and apparatus for determining noise power, and device and storage medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102938254A (en) * | 2012-10-24 | 2013-02-20 | 中国科学技术大学 | Voice signal enhancement system and method |
| CN103854660A (en) * | 2014-02-24 | 2014-06-11 | 中国电子科技集团公司第二十八研究所 | Four-microphone voice enhancement method based on independent component analysis |
| US9721582B1 (en) * | 2016-02-03 | 2017-08-01 | Google Inc. | Globally optimized least-squares post-filtering for speech enhancement |
| CN107316648A (en) * | 2017-07-24 | 2017-11-03 | 厦门理工学院 | A kind of sound enhancement method based on coloured noise |
-
2018
- 2018-09-04 CN CN201811024519.8A patent/CN109194422B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102938254A (en) * | 2012-10-24 | 2013-02-20 | 中国科学技术大学 | Voice signal enhancement system and method |
| CN103854660A (en) * | 2014-02-24 | 2014-06-11 | 中国电子科技集团公司第二十八研究所 | Four-microphone voice enhancement method based on independent component analysis |
| US9721582B1 (en) * | 2016-02-03 | 2017-08-01 | Google Inc. | Globally optimized least-squares post-filtering for speech enhancement |
| CN107316648A (en) * | 2017-07-24 | 2017-11-03 | 厦门理工学院 | A kind of sound enhancement method based on coloured noise |
Non-Patent Citations (1)
| Title |
|---|
| 牛铜: "基于子空间的语音增强算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111294126A (en) * | 2020-03-13 | 2020-06-16 | 南京航空航天大学 | Subspace tracking method for SNR estimation |
| CN111294126B (en) * | 2020-03-13 | 2021-08-03 | 南京航空航天大学 | Subspace tracking method for SNR estimation |
| CN113364539A (en) * | 2021-08-09 | 2021-09-07 | 成都华日通讯技术股份有限公司 | Blind estimation method for signal-to-noise ratio of digital signal in frequency spectrum monitoring equipment |
| CN113364539B (en) * | 2021-08-09 | 2021-11-16 | 成都华日通讯技术股份有限公司 | Blind estimation method for signal-to-noise ratio of digital signal in frequency spectrum monitoring equipment |
| WO2023246915A1 (en) * | 2022-06-21 | 2023-12-28 | 锐捷网络股份有限公司 | Method and apparatus for determining noise power, and device and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109194422B (en) | 2021-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109194422B (en) | SNR estimation method based on subspace | |
| US7079988B2 (en) | Method for the higher-order blind identification of mixtures of sources | |
| Mühlhaus et al. | Automatic modulation classification for MIMO systems using fourth-order cumulants | |
| CN1983910B (en) | Method for detecting signal in multi-antenna digital wireless telecommunication system | |
| KR100663525B1 (en) | Interference power measurement apparatus and method required space-time beam forming | |
| Grossi et al. | Adaptive detection and localization exploiting the IEEE 802.11 ad standard | |
| Li et al. | Joint estimation of channel parameters for MIMO communication systems | |
| US20090060008A1 (en) | System and method for blind source separation of signals using noise estimator | |
| CN106716866B (en) | Ping-pong beamforming | |
| US20060176971A1 (en) | Multi input multi output wireless communication reception method and apparatus | |
| CN106941463A (en) | A kind of single-bit quantification mimo system channel estimation methods and system | |
| Rodríguez-Fernández et al. | A compressive sensing-maximum likelihood approach for off-grid wideband channel estimation at mmWave | |
| CN114726686B (en) | Uniform area array millimeter wave large-scale MIMO channel estimation method | |
| CN109547376A (en) | A kind of Modulation Identification method of the gauss hybrid models based on modulating characteristic | |
| EP2293482B1 (en) | Signal detecting method and device for multi-code multi-transmitting multi-receiving system | |
| Bahloul et al. | An efficient likelihood-based modulation classification algorithm for multiple-input multiple-output systems | |
| Pourkhaatoun et al. | High-resolution independent component analysis-based time-of-arrival estimation for line-of-sight multipath environments | |
| US10498558B1 (en) | Symbol detection in shared wireless channel | |
| CN108566227B (en) | Multi-user detection method | |
| Bhatti et al. | Spectrum sensing using principal component analysis | |
| Rice et al. | Phase 1 final report: Preamble assisted equalization for aeronautical telemetry (PAQ) | |
| CN111901058A (en) | Multi-antenna auxiliary broadband spectrum sensing method based on sub-nyquist sampling | |
| Ikhlef et al. | Constant modulus algorithms using hyperbolic Givens rotations | |
| Causey et al. | Blind multiuser detection using linear prediction | |
| Salam et al. | A unified practical approach to modulation classification in cognitive radio using likelihood-based techniques |
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 |