CN101692666A - Method for filtering time domain channel response noises - Google Patents

Method for filtering time domain channel response noises Download PDF

Info

Publication number
CN101692666A
CN101692666A CN200910190542A CN200910190542A CN101692666A CN 101692666 A CN101692666 A CN 101692666A CN 200910190542 A CN200910190542 A CN 200910190542A CN 200910190542 A CN200910190542 A CN 200910190542A CN 101692666 A CN101692666 A CN 101692666A
Authority
CN
China
Prior art keywords
time domain
domain channel
channel response
response
noises
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
Application number
CN200910190542A
Other languages
Chinese (zh)
Other versions
CN101692666B (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.)
LEAGUER MICROELECTRONICS CO., LTD.
Wuxi Jingxin Microelectronics Co., Ltd.
Original Assignee
SHENZHEN LIHE MICROELECTRONICS CO Ltd
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 SHENZHEN LIHE MICROELECTRONICS CO Ltd filed Critical SHENZHEN LIHE MICROELECTRONICS CO Ltd
Priority to CN2009101905429A priority Critical patent/CN101692666B/en
Publication of CN101692666A publication Critical patent/CN101692666A/en
Application granted granted Critical
Publication of CN101692666B publication Critical patent/CN101692666B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Noise Elimination (AREA)

Abstract

The invention discloses a method for filtering time domain channel response noises. Based on time domain synchronous orthogonal frequency division multiplexing system, and on the basis of eliminating inter-symbol interferences and inter-carrier interferences between a pseudorandom sequence and a data frame body part, the method adopts a minimum mean square error algorithm and simultaneously combines the dynamic evaluation of the noise margin of a time domain channel response to accurately differentiate an effective path and an invalid path in a time domain channel impulse response, can filter most of channel noises and accurately estimate channel characteristics under various channel environments, and greatly improves the performance of a system receiver.

Description

A kind of method of filtering time domain channel response noises
Technical field
The present invention relates to a kind of method of filtering time domain channel response noises.
Background technology
In time-domain synchronization OFDM (TDS-OFDM) system, the signal frame that transmits at transmitting terminal can be divided into two parts of non-overlapping copies in time, i.e. pseudorandom PN sequence { c I, k} K=0 M-1And data division { s I, k} K=0 N-1, as depicted in figs. 1 and 2.Yet, because the existence of multipath effect, at receiving terminal received signal (r I, k} K=0 M+N-1But the two parts by aliasing are constituted: { y I, k} K=0 M+L-1{ x I, k} K=0 N+L-1, { y wherein I, k} K=0 M+L-1The linear convolution result of expression PN sequence and channel impulse response, promptly
Figure G2009101905429D0000011
{ x I, k} K=0 N+L-1Expression data division { s I, k} K=0 N-1With the linear convolution result of channel impulse response, promptly
Figure G2009101905429D0000012
As shown in Figure 3 and Figure 4.
The aliasing of received signal under the situation of not protecting prefix, will produce intersymbol interference ICI and intersymbol interference ISI owing to the multidiameter delay expansion causes.Therefore need take certain measure that interference is eliminated, realize received signal (x I, k} K=0 N+L-1{ y I, k} K=0 M+L-1Separation.
According to Fig. 3 and Fig. 4, can obtain the data { r that the i frame receives I, k} K=0M+N+L-2Can be expressed as:
r i,k=u i,k+n i,k,0≤k≤M+N+L-2 (1)
Wherein,
u i , k = x i - 1 , k + N + y i , k 0 ≤ k ≤ L - 2 y i , k L - 1 ≤ k ≤ M x i , k - M + y i , k M + 1 ≤ k ≤ M + L - 1 x i , k - M M + L ≤ k ≤ N + M x i , k - M + y i + 1 , k - N - M M + N ≤ k ≤ N + M + L - 2 - - - ( 2 )
Therefore, in order to recover the initial data that transmitting terminal sends accurately, channel estimating is carried out channel compensation accurately.In order to obtain time domain channel impulse response accurately, adopt following handling process usually:
1), at first the PN sequence of i frame signal is separated with the Frame body portion, obtains the PN sequence { c of i signal frame I, k} K=0 M-1Channel impulse response with the i frame
Figure G2009101905429D0000021
The linear convolution result, promptly eliminate the i frame PN received signal after intersymbol interference and the intersymbol interference
Figure G2009101905429D0000022
As shown in Figure 5.
2), with i frame PN received signal zero padding to 2048 point after the above-mentioned elimination interference
Figure G2009101905429D0000023
3), receiver synchronously after, obtain i signal frame PN sequence { c I, k} K=0 M-1, the point of zero padding to 2048 simultaneously c ' I, k} K=0 2047
4), according to separation signal
Figure G2009101905429D0000024
With the PN sequence of i frame c ' I, k} K=0 2047, utilize the DFT conversion, obtain initial time domain channel impulse response, promptly
h ^ i , k = IDFT { DFT ( { y ^ i , k &prime; } k = 0 2047 ) DFT ( { c i , k &prime; } k = 0 2047 ) } , 0 &le; k < 2047 - - - ( 3 )
Do not consider The noise according to the time domain channel impulse response that formula (3) calculates, under the bigger situation of noise ratio, can have a strong impact on the performance of system receiver.In order to improve the receptivity of system, must take the method for effective filtering noise, to obtain relatively accurate time domain channel estimated value.
Summary of the invention
The objective of the invention is to propose a kind of method of filtering time domain channel response noises, technical problem to be solved is: the channel response noises of the energy filtering overwhelming majority.
Technical problem of the present invention is solved by following technical scheme:
A kind of method of filtering time domain channel response noises is applied to may further comprise the steps under the time-domain synchronous orthogonal frequency-division multiplexing system framework,
(1) receives current frame signal, and pseudo random sequence is separated with the Frame body portion, eliminate corresponding intersymbol interference and intersymbol interference;
(2), and do Fourier transform with pseudo random sequence part zero padding to 2048 point that receives;
(3) receiver synchronously after, pseudo random sequence zero padding to 2048 point that the present frame that obtains is original, and do fast fourier transform;
(4) on the basis of the Fourier transform results that step (2) and step (3) obtain, utilize least-mean-square error algorithm, obtain initial time domain channel impulse response;
(5) according to initial time domain channel impulse response, the noise power of assessment time domain channel response, and adaptive calculating optimum time domain channel thresholding;
(6), distinguish effective diameter and invalid footpath in the time domain channel response, filtering overwhelming majority noise according to best time domain channel thresholding.
In the optimized technical scheme,
Comprise that also step (7) all is changed to zero with the invalid footpath that obtains in the step (6), further eliminates the noise in the time domain channel response.
The beneficial effect that the present invention is compared with the prior art is:
The present invention in conjunction with the self-adaptive estimation of time domain channel response noises, assesses the multipath length of time domain channel response by the frequency domain least-mean-square error algorithm accurately, can filtering overwhelming majority time domain channel response noises.Further, with invalid multipath zero setting, the noise in the channel response is assessed in further elimination.The method that the present invention proposes is handled simple, realize easily, and it is few to take resource, can improve the receptivity of receiver system greatly.
Description of drawings
Fig. 1 is the transmission signal frame schematic diagram in the background technology;
Fig. 2 is the Time Domain Decomposition schematic diagram of the transmission signal frame in the background technology;
Fig. 3 is the Time Domain Decomposition schematic diagram of the received signal frame in the background technology;
Fig. 4 is the received signal frame schematic diagram in the background technology;
Fig. 5 is the i frame PN received signal after the elimination in the background technology is disturbed;
Fig. 6 is a filtering time domain channel response noises flow chart in the embodiment of the invention;
Also in conjunction with the accompanying drawings the present invention is described in further details below by concrete execution mode.
Embodiment
Among the embodiment, the method for filtering time domain channel response noises comprises the steps:
The first step: receive current frame signal, and pseudo random sequence PN is partly separated with the OFDM frame, eliminate corresponding intersymbol interference and intersymbol interference;
Second step:, and do the FFT conversion with pseudo random sequence part zero padding to 2048 point that receives;
The 3rd step: obtain the original pseudo random sequence of present frame, same zero padding to 2048 point, and do the FFT conversion;
The 4th step: according to the time-domain and frequency-domain transformational relation of linear convolution relation, and utilize MMSE frequency domain filtering algorithm, effectively eliminate partial noise, obtain initial time domain channel impulse response;
The 5th step: according to initial time domain channel impulse response, the noise power of assessment time domain channel response, and adaptive calculating optimum time domain channel thresholding;
The 6th step: accurately judge the multipath length of time domain channel response, accurately distinguish effective diameter and invalid footpath in the time domain channel response simultaneously, remove the various noises in the time domain channel estimated value;
The 7th step: the invalid footpath zero setting in directly time domain channel being estimated, eliminate the most of noise in the time domain channel response, obtain final accurate relatively time domain channel response estimated value.
Wherein, why adopting 2048 points in second and third step, is through the preferred result of inventor.If adopt 1024 points, then count and do not reach effect very little; If adopt more the multiple spot number then cause operand to increase sharply.
Utilize least mean-square error MMSE algorithm in the 4th step, formula (3) is optimized, with the influence of filtering partial noise, promptly
h ^ i , k = IDFT { DFT ( { y ^ i , k &prime; } k = 0 2047 ) DFT ( { c i , k &prime; } k = 0 2047 + 1 SNR ) } , 0 &le; k < 2047 - - - ( 4 )
Wherein, SNR is the signal to noise ratio of received signal.
Although the time domain channel response that formula (4) calculates has been eliminated partial noise, but still need do further noise filtering and handle.Generally speaking, the relation of the time domain channel response that obtains of Shi Ji time domain channel response and formula (4) assessment can be expressed as:
h ^ i , k = h i , k + w ~ k , 0 &le; k &le; L - 1 w ~ k , L &le; k &le; 2047 - - - ( 5 )
In above-mentioned, L represents the length of true multipath channel,
Figure G2009101905429D0000043
Expression and signal are fully independently, the time domain average is zero, variance is σ 2Additive white Gaussian noise.
Since the residing wireless channel of terrestrial DTV be a broadband, at a high speed, the channel of (tens kilometers) on a large scale.Under abominable overland propagation environment, its channel circumstance has very strong multipath interference characteristic, radiofrequency signal is through various reflected by objects, diffraction, refraction, through influencing each other behind the different paths arrival receivers, thereby form the various declines of signal spectrum, cause the length variations of channel multi-path very big.
Therefore, need assess accurately,, therefore utilize adaptive dynamic multi-path thresholding to adapt with it because noise situations is dynamic change to the length of time domain channel response.Because the length of multipath must be less than the length of pseudorandom PN sequence, that is:
L≤N cp=M (6)
So, utilize formula (5) and formula (6) that the average power of noise is assessed accurately:
Figure G2009101905429D0000051
By assessing the noise situations of above-mentioned time domain channel response, adopt certain decision threshold, accurately assess effective multipath length on the one hand, on the other hand, to be considered as noise less than the smaller value of this thresholding, directly zero setting, thus effectively eliminate noise in the time domain channel impulse response.By this decision threshold, obtain more accurate time domain channel estimated value:
Figure G2009101905429D0000053
Therefore, how to obtain best time domain channel threshold value h ThrBecome key.For the domain channel response mean square error minimum that assessment is obtained, promptly
&epsiv; min = min { E [ &Sigma; n = 0 N | H ~ n - H n | 2 ] }
Figure G2009101905429D0000056
Wherein, the signal to noise ratio snr and the β factor are
SNR = E [ | X k | 2 ] &sigma; n 2 - - - ( 10 )
β=E[|X k| 2]E[|X k| -2] (11)
The β factor only depends on the constellation modulation system of employing.
Therefore, best threshold value can be set is
h thr optimal = 2 &beta; SNR = 2 E [ | X k | - 2 ] &sigma; n 2 - - - ( 12 )
At last, utilize formula (8) and formula (12), noise by dynamic evaluation time domain channel impulse response, according to state of signal-to-noise, the noise of adaptive filtering channel response, thereby obtain more accurate time domain channel response, to eliminate the distortion effects of multipath effect, the well various declines of compensate for channel to initial data.
The present invention is based on the TDS-OFDM system, in the intersymbol interference (ISI) and inter-carrier interference (ICI) basis eliminated between PN sequence and the OFDM data division, adopt the MMSE filtering algorithm, while is in conjunction with the dynamic evaluation of the noise gate of time domain channel response, effective diameter and invalid footpath in the time domain channel impulse response have been distinguished accurately, filtering most interchannel noises, estimate the characteristic of channel under the various channel circumstances accurately, greatly improved the performance of system receiver.
Above content be in conjunction with concrete preferred implementation to further describing that the present invention did, can not assert that concrete enforcement of the present invention is confined to these explanations.For the general technical staff of the technical field of the invention, without departing from the inventive concept of the premise, can also make some simple deduction or replace, all should be considered as belonging to protection scope of the present invention.

Claims (2)

1. the method for a filtering time domain channel response noises is applied to it is characterized in that under time-domain synchronization OFDM (TDS-OFDM) system framework: may further comprise the steps,
(1) receives current frame signal, and pseudo random sequence (PN) is separated with the Frame body portion, eliminate corresponding intersymbol interference and intersymbol interference;
(2), and do Fourier transform with pseudo random sequence part zero padding to 2048 point that receives;
(3) receiver synchronously after, the pseudorandom that the present frame that obtains is original (PN) sequence zero padding to 2048 point, and do fast fourier transform;
(4) on the basis of the Fourier transform results that step (2) and step (3) obtain, utilize least-mean-square error algorithm, obtain initial time domain channel impulse response;
(5) according to initial time domain channel impulse response, the noise power of assessment time domain channel response, and adaptive calculating optimum time domain channel thresholding;
(6), distinguish effective diameter and invalid footpath in the time domain channel response, filtering overwhelming majority noise according to best time domain channel thresholding.
2. the method for filtering time domain channel response noises according to claim 1 is characterized in that: comprise that also (7) all are changed to zero with the invalid footpath that obtains in the step (6), further eliminate the noise in the time domain channel response.
CN2009101905429A 2009-10-13 2009-10-13 Method for filtering time domain channel response noises Active CN101692666B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009101905429A CN101692666B (en) 2009-10-13 2009-10-13 Method for filtering time domain channel response noises

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009101905429A CN101692666B (en) 2009-10-13 2009-10-13 Method for filtering time domain channel response noises

Publications (2)

Publication Number Publication Date
CN101692666A true CN101692666A (en) 2010-04-07
CN101692666B CN101692666B (en) 2012-05-02

Family

ID=42081325

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009101905429A Active CN101692666B (en) 2009-10-13 2009-10-13 Method for filtering time domain channel response noises

Country Status (1)

Country Link
CN (1) CN101692666B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102664850A (en) * 2012-04-13 2012-09-12 豪威科技(上海)有限公司 Multi-carrier mode low-complexity channel noise-reducing method and device thereof in wireless local area network (WLAN)
CN102710564A (en) * 2012-06-14 2012-10-03 深圳数字电视国家工程实验室股份有限公司 Channel time domain impulse response filter method and device
CN103457885A (en) * 2012-05-31 2013-12-18 华为技术有限公司 Channel estimation method, channel estimation device and receiver
CN106027431A (en) * 2016-05-06 2016-10-12 晶晨半导体(上海)有限公司 Channel estimation method
WO2018000914A1 (en) * 2016-07-01 2018-01-04 晶晨半导体(上海)股份有限公司 Single-carrier channel estimation method
CN110460551A (en) * 2019-02-24 2019-11-15 苏州工业园区新国大研究院 A kind of channel carrier phase recovery method based on least mean-square error
CN112434415A (en) * 2020-11-19 2021-03-02 中国电子科技集团公司第二十九研究所 Method for implementing heterogeneous radio frequency front end model for microwave photonic array system

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102664850A (en) * 2012-04-13 2012-09-12 豪威科技(上海)有限公司 Multi-carrier mode low-complexity channel noise-reducing method and device thereof in wireless local area network (WLAN)
CN102664850B (en) * 2012-04-13 2014-10-22 豪威科技(上海)有限公司 Multi-carrier mode low-complexity channel noise-reducing method and device thereof in wireless local area network (WLAN)
CN103457885A (en) * 2012-05-31 2013-12-18 华为技术有限公司 Channel estimation method, channel estimation device and receiver
CN103457885B (en) * 2012-05-31 2016-09-14 华为技术有限公司 Channel estimation methods, channel estimating apparatus and receiver
CN102710564A (en) * 2012-06-14 2012-10-03 深圳数字电视国家工程实验室股份有限公司 Channel time domain impulse response filter method and device
CN102710564B (en) * 2012-06-14 2015-03-18 深圳数字电视国家工程实验室股份有限公司 Channel time domain impulse response filter method and device
CN106027431A (en) * 2016-05-06 2016-10-12 晶晨半导体(上海)有限公司 Channel estimation method
WO2017190589A1 (en) * 2016-05-06 2017-11-09 晶晨半导体(上海)股份有限公司 Channel estimation method
CN106027431B (en) * 2016-05-06 2019-05-03 晶晨半导体(上海)股份有限公司 A kind of channel estimation methods
WO2018000914A1 (en) * 2016-07-01 2018-01-04 晶晨半导体(上海)股份有限公司 Single-carrier channel estimation method
CN110460551A (en) * 2019-02-24 2019-11-15 苏州工业园区新国大研究院 A kind of channel carrier phase recovery method based on least mean-square error
CN112434415A (en) * 2020-11-19 2021-03-02 中国电子科技集团公司第二十九研究所 Method for implementing heterogeneous radio frequency front end model for microwave photonic array system

Also Published As

Publication number Publication date
CN101692666B (en) 2012-05-02

Similar Documents

Publication Publication Date Title
CN101692666B (en) Method for filtering time domain channel response noises
CN103095639B (en) Orthogonal frequency division multiplexing (OFDM) underwater acoustic communication parallel iterative inter-carrier interference (ICI) elimination method
AU2004315369C1 (en) Pilot transmission and channel estimation for an OFDM system with excess delay spread
US7609773B2 (en) Method of determining the location of the FFT window and the delay spread for the platinum broadcast channel estimator
EP2387840B1 (en) Interpolated channel estimation for mobile ofdm systems
CN101507219B (en) Method and system for time error estimation for data symbols
KR100967058B1 (en) Method for Estimate Channel in Radio Communication and device thereof
US20070297522A1 (en) Method for Signal Processing and a Signal Processor in an Ofdm System
TWI422193B (en) Channel estimation apparatus and method thereof
US7773683B2 (en) Method and apparatus for ICI cancellation in communication systems
CN107332797B (en) Channel estimation method in power line OFDM communication system
KR101241824B1 (en) A receiver of communication system for orthogonal frequency division multiplexing and Method for mitigate a phase noise in thereof
CN101795246A (en) Method and device for estimating channel
CN101155157A (en) Method and apparatus for processing channel evaluation result based on transformed domain and its receiver
CN107483373A (en) A kind of the LMMSE channel estimation methods and device of the weighting of anti-multipath iteration
US20180191533A1 (en) Channel estimation method
JP2003218826A (en) Method for receiving orthogonal frequency division multiplexed signal and receiver using the same
CN109412987B (en) Channel tracking method of OFDM system
KR101390317B1 (en) Apparatus and method for compensation of channel impulse response estimation error in orthogonal frequency division multiplexing systems
TWI466508B (en) Receiver and signal receiving method thereof applied in ofdm system for removing isi
WO2009082849A1 (en) Method and apparatus for channel estimation
Manzoor et al. Novel SNR estimation technique in Wireless OFDM systems
KR101089502B1 (en) ICI Cancellation Method using Time-Varying Channel Information in OFDM Systems
Mohanraj et al. OMP based Iterative Channel Estimation Algorithm for Modified GFDM System
CN104301262B (en) A kind of single-carrier frequency domain balance ultra-wide broadband system channel estimation methods

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: SHENZHEN LIHE MICROELECTRONIC CO., LTD.

Effective date: 20131012

Owner name: WUXI JINGXIN MICROELECTRONICS CO., LTD.

Free format text: FORMER OWNER: SHENZHEN LIHE MICROELECTRONIC CO., LTD.

Effective date: 20131012

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 518052 SHENZHEN, GUANGDONG PROVINCE TO: 214000 WUXI, JIANGSU PROVINCE

TR01 Transfer of patent right

Effective date of registration: 20131012

Address after: 214000 Jiangsu province Wuxi city Huishan District Huishan Xingzhi Road Economic Development Zone No. 35 WITHUB Park C District 56 floor six

Patentee after: Wuxi Jingxin Microelectronics Co., Ltd.

Patentee after: Shenzhen Lihe Microelectronics Co., Ltd.

Address before: 518052 Research Institute of Tsinghua University, Nanshan District Science Park, Guangdong, Shenzhen province C527

Patentee before: Shenzhen Lihe Microelectronics Co., Ltd.

C56 Change in the name or address of the patentee
CP01 Change in the name or title of a patent holder

Address after: 214000 Jiangsu province Wuxi city Huishan District Huishan Xingzhi Road Economic Development Zone No. 35 WITHUB Park C District 56 floor six

Patentee after: Wuxi Jingxin Microelectronics Co., Ltd.

Patentee after: LEAGUER MICROELECTRONICS CO., LTD.

Address before: 214000 Jiangsu province Wuxi city Huishan District Huishan Xingzhi Road Economic Development Zone No. 35 WITHUB Park C District 56 floor six

Patentee before: Wuxi Jingxin Microelectronics Co., Ltd.

Patentee before: Shenzhen Lihe Microelectronics Co., Ltd.