CN104052555A - 一种ofdm系统下无线信道多径参数估计的方法 - Google Patents
一种ofdm系统下无线信道多径参数估计的方法 Download PDFInfo
- Publication number
- CN104052555A CN104052555A CN201410276031.XA CN201410276031A CN104052555A CN 104052555 A CN104052555 A CN 104052555A CN 201410276031 A CN201410276031 A CN 201410276031A CN 104052555 A CN104052555 A CN 104052555A
- Authority
- CN
- China
- Prior art keywords
- footpath
- pdp
- reference signal
- timing error
- wireless channel
- 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
Landscapes
- Mobile Radio Communication Systems (AREA)
Abstract
本发明公开一种OFDM系统下无线信道多径参数估计的方法,方法包括以下步骤:1、对接收信号做符号定时粗同步,提取接收到的参考信号;2、利用接收端参考信号和本地参考信号得到无线信道功率时延谱,确定多径数量,利用功率时延谱与同步误差的关系,补偿同步误差,并估计第一径强度;3、对上一步补偿后的信道冲击响应做处理,消去第一径,利用功率时延谱与相对时延的关系,估计并补偿后一径相对消去径的相对时延,得到后一径强度;4、重复3步骤直到最后一径相对前一径相对时延以及最后一径强度被估计出;5、最后补偿了同步误差,得到了各径强度以及各径相对第一径时延。本发明实现了准确估计了OFDM系统下无线信道多径参数,抗噪性能优良,复杂度低。
Description
技术领域
本发明涉及无线通信技术领域,特别是涉及一种OFDM系统下无线信道多径参数估计的方法。
背景技术
信道测量是无线通信的基础之一,基于信道测量获得的信道模型是链路仿真、系统仿真的关键依据,是通信系统研发必不可少的工具。目前,信道测量方法主要分为时域测量和频域测量。时域测量主要是周期脉冲测量和伪随机序列相关测量,频域测量一个是线性调频法,一个是多载波法,同时发送不同权值不同相位的不同频点的正弦波,而其优选激励信号就是正交频分复用OFDM信号。
另外在宽带通信系统中,多径问题影响整个系统的性能。信道测量中,准确估计信道的多径强度和多径时延,既对功率时延谱PDP的估计,意义重大,它能够还原最真实的信道,对建立相应的信道模型至关重要。
发明内容
针对以上现有技术的不足,本发明提供一种OFDM系统下无线信道多径参数估计的方法,以提高了对无线信道多径参数估计的准确性,复杂度低抗噪能力强。
一种OFDM系统下无线信道多径参数估计的方法,所述方法包括以下步骤:
1)对接收信号做符号定时粗同步,再提取接收端的参考信号;
2)根据接收端参考信号和本地参考信号得到无线信道PDP,确定多径数量,利用PDP与定时误差的关系,估计补偿定时误差,并准确估计第一径强度及补偿后的CFR;
3)对所述补偿后的信道频率响应CFR进行IFFT得到补偿后的信道冲击响应CIR,并对其消去主径,利用PDP与相对时延的关系,估计并补偿后一径相对消去径的相对时延,估计得到后一径强度;
4)重复3步骤直到最后一径相对前一径相对时延以及最后一径强度被估计出;
5)最后补偿了定时误差,得到了各径强度以及各径相对第一径时延。
所述步骤1对接收信号做符号定时粗同步,再提取接收到的参考信号包括首先对接收数据作符号定时粗同步,将得到的同步头往CP推若干点,然后按照无线帧结构确定包含参考信号信息的数据,去掉CP,提取包含参考信号信息的数据。
所述步骤2利用接收端参考信号和本地参考信号得到无线信道PDP包括将接收端参考信号Rk,t除以本地参考信号Sk,t,得到信道频率响应CFR,即Hk,t,对CFR做M点IDFT得到CIR,即hn,t,n=0,1,2…M-1;对若干时间上的CIR取模平方再平均得到PDP,即
所述估计第一径强度为:确定PDP上第一径到原点距离Ni,根据需要的精度穷举Nf,穷举次数为exhausttimes,穷举范围为[-0.5,0.5],穷举定时误差Ni+Nf,n,n=0,1,2...exhausttimes-1,分别按照定时误差与信号相位关系补偿到CFR上,即重新计算得到exhausttimes个PDPn',判断哪个穷举误差Ni+Nf,n对应的PDP原点值最大,既将第一径搬移至原点并抑制它的泄漏,Ni+N'f下,PDP原点值最大,则定时误差为Δτ0=Ni+N'f,此时PDP原点值为第一径强度P0,并得到补偿后的CFR,
所述步骤3估计并补偿后一径相对消去径的相对时延及估计后一径强度具体为:对补偿后的CFR通过IFFT得到补偿过后的CIR,利用线性插值消去补偿过后的CIR的第一径即,h″0,t=(h′1,t+h'M-1,t)/2;但这里不限于用线性插值的方法,可以选用最小二乘拟合或函数逼近或更好的方法代替。消除了第一径后的CIR,h″n,t可视作新的CIR,它存在定时误差,定时误差等于当前第一径与被消去径的相对时延,用步骤2中的方法估计出定时误差以及新的第一径强度,以及补偿后的CFR,这时得到的定时误差其实是该径相对前一径的相对时延Δτ1。
本发明的有益效果如下:
本发明利用相对时延差与PDP的关系以及在CIR函数上逐次消除各径影响的方法,准确的估计了各径强度以及相对时延,具有良好的抗噪性能,且复杂度低。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明。
图1示出本发明方法流程图。
图2示出LTE主同步信号滑动相关结果。
图3示出未补偿定时误差时的PDP。
图4示出不同穷举误差下补偿后PDP原点强度示意图。
图5示出补偿定时误差后的PDP。
图6示出消去第一径后的PDP。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
如图1所示,本发明提供一种OFDM系统下无线信道多径参数估计的方法,包括以下步骤:
A、对接收信号做符号定时粗同步,再提取接收到的参考信号;
B、利用接收端参考信号和本地参考信号得到无线信道PDP,确定多径数量,再利用PDP与定时误差的关系,估计补偿定时误差,并准确估计第一径强度;
C、对上一步补偿后的CIR做处理,消去主径,利用PDP与相对时延的关系,估计并补偿后一径相对消去径的相对时延,估计得到后一径强度;
D、重复C步骤直到最后一径相对前一径相对时延以及最后一径强度被估计出;
E、最后补偿了定时误差,得到了各径强度以及各径相对第一径时延。
在步骤A中,首先对接收数据作符号定时粗同步,将得到的同步头往CP推若干点,然后按照无线帧结构确定包含参考信号信息的数据,去掉CP,提取包含参考信号信息的数据。对于FFT长度为N的OFDM系统,对包含参考信号的数据做N点FFT,再根据实际系统载波情况提取得到接收端参考信号Rk,t,k=0,1,2...M-1,M为参考信号占用载波数目,t为时间。
在步骤B中,首先将接收端参考信号Rk,t除以本地参考信号Sk,t,得到信道频率响应(CFR),即Hk,t,对CFR做M点IDFT得到CIR,即hn,t,n=0,1,2…M-1;对若干时间上的CIR取模平方再平均得到PDP,即Snapshot为统计的CIR个数,要求取得信道冲击响应在相干时间内。根据PDP确定多径数目,tap。确定PDP上第一径到原点距离Ni,根据需要的精度穷举Nf,穷举次数为exhausttimes,穷举范围为[-0.5,0.5],穷举定时误差Ni+Nf,n,n=0,1,2...exhausttimes-1,分别按照定时误差与信号相位关系补偿到CFR上,即重新计算得到exhausttimes个PDPn'。判断哪个穷举误差Ni+Nf,n对应的PDP原点值最大。Ni+N'f下,PDP原点值最大,则定时误差为Δτ0=Ni+N'f,此时PDP原点值为第一径强度P0,并得到补偿后的CFR,
在步骤C中,先将B步骤中得到的补偿后的CFR通过IFFT得到补偿过后的CIR,再利用线性插值消去的补偿过后的CIR第一径即,h″0,t=(h′1,t+h'M-1,t)/2;但这里不限于用线性插值的方法,可以选用最小二乘拟合或函数逼近或更好的方法代替。将消除了第一径的CIR,h″n,t视作新的CIR,它存在定时误差,定时误差等于当前第一径与前一径的相对时延,用步骤B中的方法估计出定时误差以及新的第一径强度,以及补偿后的CFR,这时得到的定时误差其实是该径相对前一径的相对时延Δτ1。
在步骤D中,重复C步骤直到得到最后一径相对前一径相对时延Δτ0,Δτ1,Δτ2...Δτtap-1,最后一径强度P0,P1,P2...Ptap-1。
在步骤E中,有为各径对第一径相对时延,P0,P1,P2...Ptap-1为各径强度;T=Δτ0为定时误差,此步骤是重复D步骤,将所有的相对时延和强度估计出。
下面以FDD-LTE为例,其中带宽20M,利用信道衰落模块仿真信道,信道设置为三径瑞利信道,参数如下:
表1实验设定瑞利多径信道参数
步骤A:首先将接收基带信号与本地PSS信号做滑动相关,图2为相关结果,找到同步头head=72549,再将同步头向cp内移动cplength/4的长度,cplength=144为CP长度,head=head-cplength/4=72549-144/4=72513。根据LTE协议进行帧同步,得到无线帧帧头,然后按照无线帧结构去除CP,然后将数据进行2048点FFT转到频域,得到时频网格图,按照协议提取参考信号Rk,t,k=0,1,2...199。
在步骤B中,首先将接收端参考信号Rk,t除以本地参考信号Sk,t,得到信道频率响应(CFR),即Hk,t,对CFR做M点IDFT得到CIR,即hn,t,n=0,1,2…M-1;对若干时间上的CIR取模平方再平均得到PDP,即snapshot=600。
图3是存在定时误差时得到的信道PDP,根据多径门限-30dB,多径数目tap=3。确定PDP上第一径到原点距离为21,精度穷举0.1,穷举次数11,穷举范围为[-0.5,0.5],穷举定时误差Ni+Nf,n=21+0.1*n-0.5,n=0,1,2...10,分别按照定时误差与信号相位关系补偿到CFR上,即重新计算得到11个PDPn',结果如图4。
由图4可见Ni+Nf,n=21+0.1*6-0.5=21.1对应的PDP原点值最大,则定时误差为Δτ0=Ni+N'f=21.1,此时PDP原点值为第一径强度P0=8.8854e+07,并得到补偿后的CFR,
在步骤C中,先将B步骤中得到的补偿后的CFR通过IFFT得到补偿过后的CIR,图5为补偿了定时偏差后的PDP.再利用线性插值消去的第一径既,h″0,t=(h′1,t+h'M-1,t)/2;但这里不限于用线性插值的方法,可以选用最小二乘拟合或函数逼近或更好的方法代替。将消除了第一径的CIR,h″n,t视作新的CIR,它存在定时误差,定时误差等于当前第一径与被消去径的相对时延,图6为消除第一径后的PDP。用步骤B中的方法估计出定时误差以及新的第一径强度,以及补偿后的CFR,这时得到的定时误差其实是该径相对前一径的相对时延Δτ1=25.4,P1=8.787e+06。
在步骤D中,重复C步骤第三径相对第二径相对时延Δτ2=74.1以及第三径强度P2=7.516e+05。
在步骤E中,有τ0=0,τ1=25.4,τ2=99.5为各径对第一径相对时延,P0=8.8854e+07,P1=8.787e+06,P2=7.516e+05为各径强度;T=21.1为定时误差。对功率做归一化以及对数化,得参数表如下:
表2估计得瑞利多径信道参数
可见该方法准确估计了信道的多径参数。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (5)
1.一种OFDM系统下无线信道多径参数估计的方法,其特征在于,该方法包括以下步骤:
1)对接收信号做符号定时粗同步,再提取接收端的参考信号;
2)根据接收端参考信号和本地参考信号得到无线信道PDP,确定多径数量,利用PDP与定时误差的关系,估计补偿定时误差,并估计第一径强度及补偿后的CFR;
3)对所述补偿后的CFR进行IFFT得到补偿后的CIR,并对其消去主径,利用PDP与相对时延的关系,估计并补偿后一径相对消去径的相对时延,估计得到后一径强度;
4)重复3步骤直到最后一径相对前一径相对时延以及最后一径强度被估计出;
5)最后补偿定时误差,得到了各径强度以及各径相对第一径时延。
2.根据权利要求1所述的一种OFDM系统下无线信道多径参数估计的方法,其特征在于,所述步骤1对接收信号做符号定时粗同步,再提取接收到的参考信号包括首先对接收数据作符号定时粗同步,将得到的同步头往CP推若干点,然后按照无线帧结构确定包含参考信号信息的数据,去掉CP,提取包含参考信号信息的数据。
3.根据权利要求1所述的一种OFDM系统下无线信道多径参数估计的方法,其特征在于,所述步骤2利用接收端参考信号和本地参考信号得到无线信道PDP包括将接收端参考信号Rk,t除以本地参考信号Sk,t,得到信道频率响应CFR,即Hk,t,对CFR做M点IDFT得到CIR,即hn,t,n=0,1,2…M-1;对若干时间上的CIR取模平方再平均得到PDP,即
4.根据权利要求1所述的一种OFDM系统下无线信道多径参数估计的方法,其特征在于,所述估计第一径强度为:确定PDP上第一径到原点距离Ni,根据需要的精度穷举Nf,穷举次数为exhausttimes,穷举范围为[-0.5,0.5],穷举定时误差Ni+Nf,n,n=0,1,2...exhausttimes-1,分别按照定时误差与信号相位关系补偿到CFR上,即重新计算得到exhausttimes个PDPn',判断哪个穷举误差Ni+Nf,n对应的PDP原点值最大,既将第一径搬移至原点并抑制它的泄漏,Ni+N'f下,PDP原点值最大,则定时误差为Δτ0=Ni+N'f,此时PDP原点值为第一径强度P0,并得到补偿后的CFR,
5.根据权利要求1所述的一种OFDM系统下无线信道多径参数估计的方法,其特征在于,所述步骤3估计并补偿后一径相对消去径的相对时延及估计后一径强度具体为:对补偿后的CFR通过IFFT得到补偿过后的CIR,利用线性插值消去补偿过后的CIR的第一径即,h″0,t=(h′1,t+h'M-1,t)/2;消除了第一径后的CIR,h″n,t视作新的CIR,它存在定时误差,定时误差等于当前第一径与被消去径的相对时延,用步骤2中的方法估计出定时误差以及新的第一径强度,以及补偿后的CFR,这时得到的定时误差其实是该径相对前一径的相对时延Δτ1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410276031.XA CN104052555B (zh) | 2014-06-19 | 2014-06-19 | 一种ofdm系统下无线信道多径参数估计的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410276031.XA CN104052555B (zh) | 2014-06-19 | 2014-06-19 | 一种ofdm系统下无线信道多径参数估计的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104052555A true CN104052555A (zh) | 2014-09-17 |
CN104052555B CN104052555B (zh) | 2016-04-27 |
Family
ID=51504951
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410276031.XA Expired - Fee Related CN104052555B (zh) | 2014-06-19 | 2014-06-19 | 一种ofdm系统下无线信道多径参数估计的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104052555B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106375041A (zh) * | 2016-08-25 | 2017-02-01 | 合肥工业大学 | 用于无线通信信道特征参数提取的便携式装置 |
CN109309537A (zh) * | 2018-09-30 | 2019-02-05 | 电子信息系统复杂电磁环境效应国家重点实验室 | 一种用于无线信道参数测量的信号波形处理方法 |
CN109818885A (zh) * | 2017-11-22 | 2019-05-28 | 电信科学技术研究院 | 一种信道估计方法、装置、设备及计算机可读存储介质 |
CN111595339A (zh) * | 2020-03-30 | 2020-08-28 | 杭州微萤科技有限公司 | 一种提高uwb在室内测距精度的方法及测距系统 |
CN112134823A (zh) * | 2020-09-23 | 2020-12-25 | 中国电子科技集团公司第二十八研究所 | 基于cir峰值的ofdm信号符号同步误差估计方法 |
US11165462B2 (en) * | 2018-11-07 | 2021-11-02 | Samsung Electronics Co., Ltd. | Motion assisted leakage removal for radar applications |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1529440A (zh) * | 2003-09-28 | 2004-09-15 | 中兴通讯股份有限公司 | 码分多址系统多径搜索插值方法 |
CN1592451A (zh) * | 2003-08-28 | 2005-03-09 | 华为技术有限公司 | 一种估计到达时间附加时延误差的方法 |
WO2008000007A1 (de) * | 2006-06-27 | 2008-01-03 | Nxp B.V. | Elektroakustischer wandler |
EP2399372B1 (en) * | 2009-02-17 | 2013-08-07 | Telefonaktiebolaget LM Ericsson (publ) | OFDM receiver having a plurality of FFTS according to G-Rake structure |
-
2014
- 2014-06-19 CN CN201410276031.XA patent/CN104052555B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1592451A (zh) * | 2003-08-28 | 2005-03-09 | 华为技术有限公司 | 一种估计到达时间附加时延误差的方法 |
CN1529440A (zh) * | 2003-09-28 | 2004-09-15 | 中兴通讯股份有限公司 | 码分多址系统多径搜索插值方法 |
WO2008000007A1 (de) * | 2006-06-27 | 2008-01-03 | Nxp B.V. | Elektroakustischer wandler |
EP2399372B1 (en) * | 2009-02-17 | 2013-08-07 | Telefonaktiebolaget LM Ericsson (publ) | OFDM receiver having a plurality of FFTS according to G-Rake structure |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106375041A (zh) * | 2016-08-25 | 2017-02-01 | 合肥工业大学 | 用于无线通信信道特征参数提取的便携式装置 |
CN106375041B (zh) * | 2016-08-25 | 2018-10-23 | 合肥工业大学 | 用于无线通信信道特征参数提取的便携式装置 |
CN109818885A (zh) * | 2017-11-22 | 2019-05-28 | 电信科学技术研究院 | 一种信道估计方法、装置、设备及计算机可读存储介质 |
CN109818885B (zh) * | 2017-11-22 | 2020-07-31 | 电信科学技术研究院 | 一种信道估计方法、装置、设备及计算机可读存储介质 |
CN109309537A (zh) * | 2018-09-30 | 2019-02-05 | 电子信息系统复杂电磁环境效应国家重点实验室 | 一种用于无线信道参数测量的信号波形处理方法 |
CN109309537B (zh) * | 2018-09-30 | 2021-06-01 | 电子信息系统复杂电磁环境效应国家重点实验室 | 一种用于无线信道参数测量的信号波形处理方法 |
US11165462B2 (en) * | 2018-11-07 | 2021-11-02 | Samsung Electronics Co., Ltd. | Motion assisted leakage removal for radar applications |
CN111595339A (zh) * | 2020-03-30 | 2020-08-28 | 杭州微萤科技有限公司 | 一种提高uwb在室内测距精度的方法及测距系统 |
CN111595339B (zh) * | 2020-03-30 | 2022-05-24 | 杭州微萤科技有限公司 | 一种提高uwb在室内测距精度的方法及测距系统 |
CN112134823A (zh) * | 2020-09-23 | 2020-12-25 | 中国电子科技集团公司第二十八研究所 | 基于cir峰值的ofdm信号符号同步误差估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104052555B (zh) | 2016-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104052555B (zh) | 一种ofdm系统下无线信道多径参数估计的方法 | |
CN101425999B (zh) | 正交频分复用接收机的载频偏差同步的方法及装置 | |
CN103346829B (zh) | 兼容lte模式卫星通信初始随机接入两步时延测量方法 | |
CN101938435B (zh) | 一种正交频分复用系统的时偏估计装置及方法 | |
CN104717174B (zh) | 一种复杂多径信道下的ofdm抗干扰同步方法 | |
CN102970271B (zh) | 一种联合估计载波频偏的频率同步方法 | |
CN105007150A (zh) | 低信噪比sc-fde系统同步方法及同步装置 | |
CN103580718A (zh) | 一种低信噪比下的快速时频同步方法 | |
CN101465833A (zh) | 一种正交频分复用信号定时同步方法和装置 | |
CN102868660B (zh) | 一种基于宽带无线通信系统的下行传输同步方法 | |
CN101312447A (zh) | 接收机的整数频偏估计和精细同步方法与装置 | |
CN102387108A (zh) | 一种物理随机接入信道信号的传输方法及装置 | |
CN102571676B (zh) | 正交频分复用系统中帧同步和频偏精确估计的方法 | |
CN102035770B (zh) | 一种利用相关进行的信道估计方法 | |
CN105141562A (zh) | 通信系统及其同步方法 | |
CN102143574B (zh) | 一种适合于IEEE802.16m前导字结构的定时同步方法 | |
CN104254100A (zh) | 一种上行定时提前量的测量方法 | |
CN104242983A (zh) | 一种提高频率同步精度方法和装置 | |
CN105187352A (zh) | 一种基于ofdm前导的整数频偏估计方法 | |
CN113940041B (zh) | 用于改善电缆调制解调器时钟同步的高精度时间戳检测 | |
CN101312378B (zh) | 接收机延时扩展参数的估计方法及延时扩展参数估计单元 | |
CN106789791B (zh) | 基于共轭对称训练序列的移动通信系统载波频偏估计方法 | |
CN104836770A (zh) | 一种基于相关平均与加窗的定时估计方法 | |
CN103188067A (zh) | 一种扩频系统的码片时钟频率偏差误差估计及校正的方法 | |
CN101060366B (zh) | 一种下行同步信号的定位方法及装置 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160427 Termination date: 20170619 |
|
CF01 | Termination of patent right due to non-payment of annual fee |