CN104270328B - A kind of signal to noise ratio real-time estimation method - Google Patents

A kind of signal to noise ratio real-time estimation method Download PDF

Info

Publication number
CN104270328B
CN104270328B CN201410590639.XA CN201410590639A CN104270328B CN 104270328 B CN104270328 B CN 104270328B CN 201410590639 A CN201410590639 A CN 201410590639A CN 104270328 B CN104270328 B CN 104270328B
Authority
CN
China
Prior art keywords
signal
data
mean power
noise ratio
passband
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.)
Expired - Fee Related
Application number
CN201410590639.XA
Other languages
Chinese (zh)
Other versions
CN104270328A (en
Inventor
王甲峰
尹显东
文豪
李蕾
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CHENGDU JIUXIN ELECTRONIC TECHNOLOGY CO LTD
Institute of Electronic Engineering of CAEP
Original Assignee
CHENGDU JIUXIN ELECTRONIC TECHNOLOGY CO LTD
Institute of Electronic Engineering of CAEP
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by CHENGDU JIUXIN ELECTRONIC TECHNOLOGY CO LTD, Institute of Electronic Engineering of CAEP filed Critical CHENGDU JIUXIN ELECTRONIC TECHNOLOGY CO LTD
Priority to CN201410590639.XA priority Critical patent/CN104270328B/en
Publication of CN104270328A publication Critical patent/CN104270328A/en
Application granted granted Critical
Publication of CN104270328B publication Critical patent/CN104270328B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The present invention relates to a kind of real-time estimation method, receiving terminal receives passband data signal and is, output result is converted to by quadrature frequency conversion, according to some specific input dataMean powerAnd output dataMean power, obtain receivingkThe estimate of signal to noise ratio during individual dataFor:;The present invention receives signal total mean power and useful signal mean power and the relation of noise average power, solves useful signal mean power and noise average power, and then estimate signal to noise ratio using before and after quadrature frequency conversion;The present invention is not only suitable for single carrier digital modulation signals, is also applied for multi-carrier digital modulated signal;Mass data need not be cached, the features such as amount of calculation is small by having compared with other methods of estimation, with complexity is low, wide adaptability, efficiency high the features such as, it is particularly suitable for hardware realization, it is adaptable to the technical field such as cooperative communication, the analysis of non-co-operation signal, electromagnetic environment monitor.

Description

A kind of signal to noise ratio real-time estimation method
Technical field
The present invention relates to the method for estimation of signal to noise ratio, a kind of specifically signal to noise ratio of digital modulation signals side of estimation in real time Method.
Background technology
Signal to noise ratio is an important parameter for weighing channel quality, and many applications are required for signal to noise ratio as priori, Such as the adaptive coding and modulating in cooperative communication, bit error rate estimation, Turbo decodings;And analyzed for non-co-operation signal, letter Make an uproar than being the important references parameter for selecting Modulation Identification method, and be the important evidence for assessing data after demodulating confidence level;Separately Outside, the normalized premise of signal power in the accurate estimation of signal to noise ratio or demodulating process, only realizes the normalization of signal, The correct judgement of symbol can be realized.
At present, signal-noise ratio estimation method mainly has the method for estimation based on symbol square, segmentation symbol square signal-to-noise ratio (SNR) estimation, base Method of estimation, maximum Likelihood in Subspace Decomposition, and auto-relativity function method etc..These methods are carrying out noise It is required for caching substantial amounts of data during compared estimate, and estimated accuracy requires that more high required data buffer storage space is bigger, this Undoubtedly add hard-wired complexity;In addition, the above method is only suitable for estimating single carrier number in addition to digital signal processing The signal to noise ratio of word modulated signal, and digital signal processing computation complexity is very big, hardware is realized difficult.
The content of the invention
The present invention proposes a kind of digital modulation signals signal to noise ratio real-time estimation method, directly utilizes natural birth in demodulating process Raw signal, estimates signal to noise ratio in real time, only needs to cache five data in estimation procedure, so as to overcome above method needs A large amount of data cached shortcomings, and method proposed by the invention there is no particular requirement that to digital modulation signals, both be applicable Multi-carrier digital modulated signal is also applied in single carrier digital modulation signals.
Technical scheme is as follows:
A kind of signal to noise ratio real-time estimation method, it is characterised in that real-time estimating step is as follows:
(1) what receiving terminal was received is that the passband data signal crossed by noise pollution is x (n), wherein,
N=i+1, i+2 ..., i+L;The L is data length, the i.e. length of window for signal-to-noise ratio (SNR) estimation, and i is window Sequence number;When the passband data signal received is k-th of input data x (k), the k ∈ n;
(2) so x (k) mean power is:
Wherein Px(k-1) it is the mean power of x (k-1) when receiving -1 data of kth;
(3) after for received passband data signal x (n), carry out quadrature frequency conversion and passband signal be converted into baseband signal, Quadrature frequency conversion output result y (n);After influence of the LPF to baseband signal is ignored, the y (n) can be with equivalently represented For:Y (n)=b (n)+v (n) * h (n), wherein b (n) are complex baseband signal, and v (n) is white complex gaussian noise, under h (n) is orthogonal The shock response of low pass filter in conversion process, wherein h (n) data length are K;Therefore, when receiving k-th of input number After x (k), obtained output data is y (k);
(4) so y (k) mean power is:
(5) according to x (k) mean power Px(k) with y (k) mean power Py(k) believe when obtaining receiving k-th of data Make an uproar than estimateFor:
Wherein
(6) k=k+1 is taken, (1)-(5) is repeated, until k=L, takesEstimation procedure terminates.
Further, the quadrature frequency conversion described in step (3) includes two steps of shift frequency and LPF, LPF Effect be to filter out the high-frequency signal produced during shift frequency.
From above-mentioned estimation flow, the present invention only needs additional buffered c, k, Px(k)、Py(k) andFive mediants According to, and because quadrature frequency conversion is the general process of signal demodulation, so complexity increase is several as caused by algorithm for estimating It can ignore.
Beneficial effects of the present invention are as follows:
The present invention is put down using signal total mean power before and after quadrature frequency conversion, is received with useful signal mean power and noise The relation of equal power, solves useful signal mean power and noise average power, and then estimate signal to noise ratio;The present invention was both fitted For single carrier digital modulation signals, multi-carrier digital modulated signal is also applied for;Have compared with other methods of estimation and be not required to Cache mass data, the features such as amount of calculation is small, with complexity is low, wide adaptability, efficiency high the features such as, be particularly suitable for hardware Realize, it is adaptable to the technical field such as cooperative communication, the analysis of non-co-operation signal, electromagnetic environment monitor.
Brief description of the drawings
Fig. 1 is estimation steps flow chart of the invention;
Fig. 2 is that the signal to noise ratio in the specific embodiment of the invention estimates dynamic process figure in real time.
Embodiment
A kind of signal to noise ratio real-time estimation method disclosed by the invention, as shown in figure 1, its real-time estimating step is as follows:
(1) what receiving terminal was received is that the passband data signal crossed by noise pollution is x (n), wherein,
N=i+1, i+2 ..., i+L;The L is data length, the i.e. length of window for signal-to-noise ratio (SNR) estimation, and i is window Sequence number;When the passband data signal received is k-th of input data x (k), the k ∈ n;
(2) so x (k) mean power is:
Wherein Px(k-1) it is the mean power of x (k-1) when receiving -1 data of kth;
(3) after for received passband data signal x (n), carry out quadrature frequency conversion and passband signal be converted into baseband signal, Quadrature frequency conversion output result y (n);After influence of the LPF to baseband signal is ignored, the y (n) can be with equivalently represented For:Y (n)=b (n)+v (n) * h (n), wherein b (n) are complex baseband signal, and v (n) is white complex gaussian noise, under h (n) is orthogonal The shock response of low pass filter in conversion process, wherein h (n) data length are K;Therefore, when receiving k-th of input number After x (k), obtained output data is y (k);
(4) so y (k) mean power is:
(5) according to x (k) mean power Px(k) with y (k) mean power Py(k) believe when obtaining receiving k-th of data Make an uproar than estimateFor:
Wherein
(6) k=k+1 is taken, (1)-(5) is repeated, until k=L, takesEstimation procedure terminates.
Specifically above-mentioned signal-noise ratio estimation method is verified by taking a kind of 16QAM signals as an example.
Simulation parameter is:Character rate 1/16sps;Carrier frequency 1/4Hz;Sample rate 1Hz;Forming filter is root liter Cosine filter, form factor is 0.35;Data length is 16000 (1000 symbols);Signal to noise ratio is 10dB;Low pass filter H (n) is FIR filter, and exponent number is 16, is 3/8Hz by frequency, specific coefficient is:0.0000,0.0048,0.0080- 0.0089, -0.0429, -0.0290,0.0973,0.2834,0.3745,0.2834,0.0973, -0.0290, -0.0429, - 0.0089 0.0080,0.0048,0.0000.
The real-time estimating step of signal to noise ratio is as follows:
(1) P is takenx(0)=Py(0)=0, L=16000, utilizeWherein
C=0.3255 is calculated, and takes k=1
(2) data x (k) is received;
(3) calculate
(4) k-th of output data y (k) after quadrature frequency conversion is received;
(5) calculate
(6) calculate
(7) k=k+1 is made, (1)-(6) are obtained until k=L
As shown in Figure 2, it is shown that the dynamic process that signal to noise ratio is estimated in real time.

Claims (2)

1. a kind of signal to noise ratio real-time estimation method, it is characterised in that real-time estimating step is as follows:
(1) what receiving terminal was received is that the passband data signal crossed by noise pollution is x (n), wherein,
N=i+1, i+2 ..., i+L;The L is data length, the i.e. length of window for signal-to-noise ratio (SNR) estimation, and i is window sequence Number;When the passband data signal received is k-th of input data x (k), the k ∈ n;
(2) so x (k) mean power is:
P x ( k ) = 1 k Σ n = 1 k | x ( n ) | 2 = k - 1 k [ 1 k - 1 Σ n = 1 k - 1 | x ( n ) | 2 + | x ( k ) | 2 k - 1 ] = 1 k [ ( k - 1 ) P x ( k - 1 ) + | x ( k ) | 2 ] ,
Wherein Px(k-1) it is the mean power of x (k-1) when receiving -1 data of kth;
(3) after for received passband data signal x (n), carry out quadrature frequency conversion and passband signal is converted into baseband signal, it is orthogonal Down coversion output result y (n);After influence of the LPF to baseband signal is ignored, the y (n) is equivalently represented to be:Y (n)= B (n)+v (n) * h (n), wherein b (n) are complex baseband signal, and v (n) is white complex gaussian noise, and h (n) is during quadrature frequency conversion The data length of the shock response of low pass filter, wherein h (n) is K;Therefore, after k-th of input data x (k) is received, Obtained output data is y (k);
(4) so y (k) mean power is:
(5) according to x (k) mean power Px(k) with y (k) mean power Py(k) signal to noise ratio when obtaining receiving k-th of data EstimateFor:
ρ ^ ( k ) = 10 lg P y ( k ) - cP x ( k ) P x ( k ) - P y ( k ) ,
Wherein
(6) k=k+1 is taken, (1)-(5) is repeated, until k=L, takesEstimation procedure terminates.
2. a kind of signal to noise ratio real-time estimation method according to claim 1, it is characterised in that:Described in step (3) just Down coversion is handed over to include two steps of shift frequency and LPF, the effect of LPF is to filter out the high frequency letter produced during shift frequency Number.
CN201410590639.XA 2014-10-29 2014-10-29 A kind of signal to noise ratio real-time estimation method Expired - Fee Related CN104270328B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410590639.XA CN104270328B (en) 2014-10-29 2014-10-29 A kind of signal to noise ratio real-time estimation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410590639.XA CN104270328B (en) 2014-10-29 2014-10-29 A kind of signal to noise ratio real-time estimation method

Publications (2)

Publication Number Publication Date
CN104270328A CN104270328A (en) 2015-01-07
CN104270328B true CN104270328B (en) 2017-07-25

Family

ID=52161820

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410590639.XA Expired - Fee Related CN104270328B (en) 2014-10-29 2014-10-29 A kind of signal to noise ratio real-time estimation method

Country Status (1)

Country Link
CN (1) CN104270328B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106330362A (en) * 2016-08-25 2017-01-11 中国电子科技集团公司第十研究所 Data assisted signal to noise ratio estimation method
CN110876187A (en) * 2018-09-03 2020-03-10 普天信息技术有限公司 Uplink power control method and device
CN117556246B (en) * 2024-01-09 2024-03-19 电信科学技术第五研究所有限公司 Method for separating single wave signal from carrier mixed signal

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7190741B1 (en) * 2002-10-21 2007-03-13 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Real-time signal-to-noise ratio (SNR) estimation for BPSK and QPSK modulation using the active communications channel
CN101977169A (en) * 2010-11-09 2011-02-16 西安电子科技大学 Time domain parameter blind evaluation method of OFDM (Orthogonal Frequency Division Multiplexing) signals
CN102137053A (en) * 2011-05-06 2011-07-27 中国工程物理研究院电子工程研究所 Method for estimating signal to noise ratio of BPSK (Binary Phase Shift Keying) signal
CN102307166A (en) * 2011-08-31 2012-01-04 成都久鑫电子科技有限公司 SNR (signal to noise ratio) estimation method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7190741B1 (en) * 2002-10-21 2007-03-13 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Real-time signal-to-noise ratio (SNR) estimation for BPSK and QPSK modulation using the active communications channel
CN101977169A (en) * 2010-11-09 2011-02-16 西安电子科技大学 Time domain parameter blind evaluation method of OFDM (Orthogonal Frequency Division Multiplexing) signals
CN102137053A (en) * 2011-05-06 2011-07-27 中国工程物理研究院电子工程研究所 Method for estimating signal to noise ratio of BPSK (Binary Phase Shift Keying) signal
CN102307166A (en) * 2011-08-31 2012-01-04 成都久鑫电子科技有限公司 SNR (signal to noise ratio) estimation method

Also Published As

Publication number Publication date
CN104270328A (en) 2015-01-07

Similar Documents

Publication Publication Date Title
CN102724147B (en) A kind of channel estimation methods of underwater sound OFDM
CN103200139B (en) A kind of ofdm signal bandwidth blind estimation
CN112733811B (en) Method for identifying underwater sound signal modulation modes based on improved dense neural network
CN101406019A (en) Method and apparatus for calculating likelihood metric of a received signal in a digital communication system
CN110138459A (en) Sparse underwater sound orthogonal frequency division multiplexing channel estimation methods and device based on base tracking denoising
CN111884685A (en) Synchronous demodulation method and device for digital communication signal
CN104270328B (en) A kind of signal to noise ratio real-time estimation method
CN108270702A (en) Turbo iteration equalizing detection methods based on MCMC
CN111245580B (en) Signal-to-noise ratio calculation system and method based on hardware logic circuit
CN102307166B (en) SNR (signal to noise ratio) estimation method
CN107947899B (en) Multi-user signal-to-noise ratio blind estimation method in single-carrier interleaved frequency division multiple access system
CN109167744B (en) Phase noise joint estimation method
CN116506273A (en) Novel MPSK modulation signal identification and classification method
CN115378776A (en) MFSK modulation identification method based on cyclic spectrum parameters
CN115664905A (en) Wi-Fi equipment identification system and method based on multi-domain physical layer fingerprint characteristics
CN102315835A (en) Method for estimating roll-off coefficient of forming filter
CN115913849A (en) Electromagnetic signal identification method based on one-dimensional complex value residual error network
Feng et al. An accurate modulation recognition method of QPSK signal
CN104753842B (en) The signal modulation mode recognition methods differentiated based on peak
CN107248964B (en) Method for estimating code rate of underlay frequency spectrum sharing time-frequency overlapping signals
CN111917674A (en) Modulation identification method based on deep learning
CN109639618B (en) Low-complexity soft output demodulation method suitable for high-order quadrature amplitude modulation signal
CN111935044B (en) PSK (phase Shift keying) and QAM (Quadrature amplitude modulation) signal modulation identification method based on EMD (empirical mode decomposition)
CN108449120B (en) Low-complexity QAM-MIMO detection method based on differential measurement
CN114615113B (en) Multi-symbol BEM channel estimation method for reducing noise of base coefficient by utilizing preamble sequence

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170725

Termination date: 20181029

CF01 Termination of patent right due to non-payment of annual fee