CN112953863A - 一种适用于超低信噪比多径传输环境的信道参数估计方法 - Google Patents

一种适用于超低信噪比多径传输环境的信道参数估计方法 Download PDF

Info

Publication number
CN112953863A
CN112953863A CN202110177061.5A CN202110177061A CN112953863A CN 112953863 A CN112953863 A CN 112953863A CN 202110177061 A CN202110177061 A CN 202110177061A CN 112953863 A CN112953863 A CN 112953863A
Authority
CN
China
Prior art keywords
multipath
vector
decision
path
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
Application number
CN202110177061.5A
Other languages
English (en)
Other versions
CN112953863B (zh
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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN202110177061.5A priority Critical patent/CN112953863B/zh
Publication of CN112953863A publication Critical patent/CN112953863A/zh
Application granted granted Critical
Publication of CN112953863B publication Critical patent/CN112953863B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L25/03178Arrangements involving sequence estimation techniques
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/309Measuring or estimating channel quality parameters
    • H04B17/336Signal-to-interference ratio [SIR] or carrier-to-interference ratio [CIR]
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/0204Channel estimation of multiple channels
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03878Line equalisers; line build-out devices

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Quality & Reliability (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
  • Noise Elimination (AREA)

Abstract

本发明属于无线通信技术领域,公开了一种适用于超低信噪比多径传输环境的信道参数估计方法,包括:产生两组PN序列PN1与PN2,利用PN1与PN2对训练序列进行组帧;根据训练序列的帧结构特点,定义并计算接收信号多径分量检测的判决向量;利用接收信号的判决值,计算多径检测判决门限;对接收信号进行检测,估计超过判决门限多径的路径时延与路径衰减;估计上述各多径分量的相位偏移。本发明主要解决了在低信噪比多径信道下,如何快速准确的估计出信道多径分量基本参数的问题。本发明对训练序列帧结构进行设计,在保证估计精度的同时,大大降低了计算复杂度,可适用于接收机的精定时同步、信道均衡、多径信号合并等部分。

Description

一种适用于超低信噪比多径传输环境的信道参数估计方法
技术领域
本发明属于无线通信技术领域,尤其涉及一种适用于超低信噪比多径传输环境的信道参数估计方法。
背景技术
目前,业内常用的现有技术是这样的:
信号在传输过程中,由于障碍物的存在,会经过多点反射,从多条不同路径到达接收端,这种多径信号的幅度、相位和到达时间均不相同,它们相互叠加会产生电平衰落和时延扩展。由于多径衰落的存在,直接利用无线信道进行通信,存在严重的码间干扰。而信道参数估计可以从接收数据中将信号经过的信道模型参数估计出来,再通过信道均衡技术对接收信号进行补偿,从而可以消除或减弱多径效应带来的码间干扰问题。因此,快速准确地对信道模型参数进行估计,在多径信道传输中是很有必要的。
根据先验信息的不同,信道参数估计方法可分为基于参考信号的估计、半盲估计及盲估计三类,基于参考信号的估计适用于突发传输系统,通过对接收训练序列与本地训练序列的比较,获得信道传输特征,而盲估计和半盲信道估计算法无需或者需要较短的训练序列,频谱效率高,但是盲估计和半盲估计方法的计算复杂度较高,且可能出现相位模糊、误差传播、收敛慢或陷入局部极小等问题。Rontogiannis提出了一种通过延迟解耦优化方法来估计多径信道冲激响应,Meng Ye通过提取滑动和曲线的突变信息,提出了一种改进的多径信道滑动窗口时延估计方法,AsifIqbal提出了一种基于时域循环相关的信道冲激响应估计技术,Yu Yang提出了一种基于匹配滤波和最小二乘法的联合改进的多径参数估计算法,使用Chirp信号或m序列估计路径数和每个路径的时延,并通过最小二乘法来估计每个路径的衰减,Yong-Hwa Kim引入了时延估计方案和衰减因子估计方案来分析多径参数,Dmitriy Shutin开发了空间交替广义期望最大化方法的稀疏变分贝叶斯扩展,用于估计频率和空间选择性无线信道中相关多径分量的参数。
综上所述,现有技术存在的问题是:
对于低信噪比多径信道,更高的信道参数估计精度,意味着更长的训练序列,或更高的迭代次数,而过长的序列或过多的迭代次数意味着增加大量的计算复杂度,导致其在实际通信系统中的应用成本难以承受,并且估计时间也会大大增加,造成通信系统的整体性能很低。
解决上述技术问题的难度:
(1)在低信噪比环境下,信号往往湮没在噪声当中,如何避免噪声对信道参数估计性能的影响,从而获得准确的多径信道基本参数,是亟需解决的第一个难点。
(2)随着训练序列长度的增加,相关窗口相应变得更长,同一时间需要的计算量大大增加,如何设计训练序列帧结构,使得系统在提高训练序列长度的同时,几乎不增加估计方法的复杂度,是亟需解决的第二个难点。
解决上述技术问题的意义:如果可以在保证系统在低信噪比下信道参数估计精度的同时,降低多径信道参数估计的实现复杂度,则可以快速获得精确的多径信道基本参数,大大降低通信系统的应用成本。
发明内容
针对现有技术存在的问题,本发明提供了一种适用于超低信噪比多径传输环境的信道参数估计方法。
本发明是这样实现的,一种适用于超低信噪比多径传输环境的信道参数估计方法,包括:
产生两组PN序列PN1与PN2,利用PN1与PN2对训练序列进行组帧,可以在估计精度不变的情况下,减小相关窗口的长度,从而极大降低计算复杂度;
根据训练序列的帧结构特点,定义并计算接收信号多径分量检测的判决向量;
利用接收信号的判决值,计算多径检测判决门限;
对接收信号进行检测,将得到的判决向量与判决门限进行比较,估计超过判决门限多径的路径时延与路径衰减;
估计上述各多径分量的相位偏移,得到接收信号中各多径分量的基本参数。
进一步,所述产生两组伪随机PN序列PN1与PN2,利用PN1与PN2对训练序列进行组帧具体包括:产生两组PN序列PN1:a=[a0,a1,...,an,...,aN-1]T与PN2:b=[b0,b1,...,bm,...,bM-1]T,其中PN1的长度为N,PN2的长度为M,要求M≥Kmax,其中
Figure BDA0002941003380000031
表示最大多径时延对应的符号个数,τmax表示信道的最大多径时延,Ts表示符号周期,
Figure BDA0002941003380000032
表示向下取整,序列a与b均表示模值为1的复序列:
Figure BDA0002941003380000033
其中,exp(·)表示指数运算,当b中的每一项均为±1,即θm∈{0,π}时,可以在基本不损失估计性能的情况下,极大降低估计复杂度;
发射机利用两个PN序列a与b进行下述的组帧操作,将a中的每一项an与b相乘,得到组帧输出的第n个向量:
xn=anb=[anb0,anb1,...,anbm,...,anbM-1]T
其中[·]T表示转置运算,xn表示组帧输出的第n个向量,由于a的长度为N,故相乘后可得到N个向量x0,x1,...,xN-1
将这些向量按顺序前后拼接起来,即可得到最终组帧后的训练序列输出x=[x0,x1,...,xnM+m,...,xNM-1]T,其可用下述公式表示:
Figure BDA0002941003380000034
其中,xnM+m=anbm表示第n个组帧向量的第m项。
进一步,所述根据训练序列的帧结构特点,定义并计算接收信号多径分量检测的判决向量具体包括:将组帧后的训练序列x经过低信噪比多径信道,得到接收输入信号r:
Figure BDA0002941003380000041
由于存在多径效应,接收信号向量为多个路径的信号经过不同延时后叠加得到的,造成符号间干扰,因此接收信号向量比发射信号向量更长,增加的长度为Kmax
定义判决向量为J=[J0,J1,...,Jk,...,JM-1]T,其中每一项Jk通过下述计算方法得到,从接收信号中截取长度为NM的向量rk
rk=[rk,rk+1,...,rk+NM-1]T
再将该向量分为N组,每组长度为M,则第k个截取向量的第n组可表示为:
rk,n=[rk+nM,rk+nM+1,...,rk+nM+m,...,rk+(n+1)M-1]T
其中,k=0,1,...,M-1,n=0,1,...,N-1,m=0,1,...,M-1,将第k个截取向量的每一组转置后分别与PN序列b的共轭相乘,然后求计算结果的模值,并将得到的N个结果相加,其可用下述公式表示:
Figure BDA0002941003380000042
其中,Jk表示判决向量中的第k项,|·|表示取模,(·)*表示取共轭,则判决向量J可用下述公式表示:
Figure BDA0002941003380000043
进一步,所述利用接收信号的判决值,计算多径检测判决门限具体包括:设判决门限向量γ=[γ01,...,γk,...,γM-1]T,其长度同样为M,由判决向量J的定义可知,该判决向量的每一项表示接收向量rk,n与PN序列b的互相关模值,而rk,n的帧格式中本身包含序列b,故会在二者对齐时存在一个明显的相关峰值,判决门限向量计算每个判决值周围点的平均值,再扩大一定倍数后,与该点的判决值进行比较,得到该点位置的判决结果,判决门限γk可用下述公式表示:
Figure BDA0002941003380000051
其中,u表示放大倍数,v表示取平均的范围。
步骤四,根据多径信道的定义,接收信号可用下述公式表示:
Figure BDA0002941003380000052
其中,rnM+m表示第n个组帧向量的第m个符号对应的接收信号,L表示多径数量,αl表示第l条路径的衰减,
Figure BDA0002941003380000053
表示第l条路径的相位,Kl表示第l条路径的时延,wnM+m表示噪声,由定义可知,Kl≤Kmax<M,故经过延时后第n个发射信号向量的起始位置,仍然落在该接收信号向量的持续时间范围内,即:
Figure BDA0002941003380000054
其中,|m-Kl|<M,设接收信号在l=l0处存在一多径,相应的路径衰减为
Figure BDA0002941003380000055
路径时延为
Figure BDA0002941003380000056
相位为
Figure BDA0002941003380000057
则判决向量中第
Figure BDA0002941003380000058
项可表示为:
Figure BDA0002941003380000059
其中,|·|表示取模运算,可以看到,由于an不参与后续相关运算,故可以直接提取出来,且|an|=1,因此在不增加复杂度的情况下,将判决向量的值扩大了N倍,在相关运算部分,只有当
Figure BDA0002941003380000061
即l=l0时,对应的多径才具有自相关特性,相关结果进一步将判决值扩大了M倍,而其他多径分量由于不具有相关性,故不会被放大,因此仅会存在l0多径分量一个峰值;
上面的分析可知,每个多径分量在判决向量中存在一个相关峰值,故可以比较判决向量与判决门限向量的大小,来找到超过判决门限的多径分量,则路径时延估计值
Figure BDA0002941003380000062
可用下述公式表示:
Figure BDA0002941003380000063
其中,
Figure BDA0002941003380000064
表示所有超过判决门限的路径时延估计值的集合,相应的路径衰减的估计值
Figure BDA0002941003380000065
可用下述公式表示:
Figure BDA0002941003380000066
其中,
Figure BDA0002941003380000067
表示所有超过判决门限的路径衰减估计值的集合。
进一步,所述估计各多径分量的相位偏移具体包括:对于第l0个路径时延估计值
Figure BDA0002941003380000068
及路径衰减估计值
Figure BDA0002941003380000069
求接收信号与本地PN序列的相关结果,相关结果的幅角即为第l0个路径的相位估计值
Figure BDA00029410033800000610
其可用下述公式表示:
Figure BDA00029410033800000611
其中,arg[·]表示求幅角运算,可以看到,由于取消了模值运算,原本的标量相加变为矢量相加,故模值大的矢量在计算中占更大的权重,而训练序列的自相关性,将第l0条多径分量对应的相关值扩大了M倍,故在计算结果中的权重远大于噪声及其他多径分量,因此可以通过求相关结果的幅角来估计第l0条多径的相位,则路径相位的估计值
Figure BDA0002941003380000071
可用下述公式表示:
Figure BDA0002941003380000072
其中,
Figure BDA0002941003380000073
表示所有超过判决门限的路径相位估计值的集合。
本发明的另一目的在于提供一种实施权利要求所述适用于超低信噪比多径传输环境的信道参数估计方法的低信噪比多径信道下的低复杂度时域信道参数估计系统。
综上所述,本发明的优点及积极效果为:
本发明提供了一种适用于超低信噪比多径传输环境的信道参数估计方法,主要解决了在低信噪比多径信道下,如何快速准确的估计出信道多径分量基本参数的问题。本发明对训练序列帧结构进行设计,在保证估计精度的同时,大大降低了计算复杂度,同时,本发明根据当前信道状态,动态调整判决门限的大小,保证在不同信道环境下的工作性能。
本发明可在极低信噪比多径信道的恶劣环境下,准确估计出接收强度高的若干多径分量的基本参数,避免了多径效应引起的符号间干扰,同时本发明将训练序列分为两组进行组帧,在保证估计精度的同时,降低了相关窗口的长度,极大地降低了系统的复杂度。本发明所述方法可适用于接收机的精定时同步、信道均衡、多径信号合并等部分,在消除多径效应引起的符号间干扰的同时,还可利用不同路径接收信号之间具有一定的独立性,对多径信号进行分集合并处理,从而有效提升系统的整体性能。
附图说明
为了更清楚地说明本申请实施例的技术方案,下面将对本申请实施例中所需要使用的附图做简单的介绍。
图1是本发明实施例提供的信道参数估计方法流程图。
图2是本发明实施例提供的训练序列组帧方式示意图。
图3是本发明实施例提供的接收机信道参数估计及判决流程图。
图4是本发明实施例提供的路径时延估计准确率仿真结果示意图。
图5是本发明实施例提供的路径衰减估计精度仿真结果示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
针对现有技术存在的问题,本发明提供了一种适用于超低信噪比多径传输环境的信道参数估计方法,下面结合附图对本发明的应用原理作详细的描述。
如图1所示,本发明实施例提供的适用于超低信噪比多径传输环境的信道参数估计方法包括以下步骤:
S101:产生两组PN序列PN1与PN2,利用PN1与PN2对训练序列进行组帧;
S102:根据训练序列的帧结构特点,定义并计算接收信号多径分量检测的判决向量;
S103:利用接收信号的判决值,计算多径检测判决门限;
S104:对接收信号进行检测,估计超过判决门限多径的路径时延与路径衰减;
S105:估计上述各多径分量的相位偏移。
本发明实施例提供的多径信道下短突发直接序列扩频系统的时间与频率同步方法具体包括以下步骤:
步骤S101中,产生两组PN序列PN1:a=[a0,a1,...,an,...,aN-1]T与PN2:b=[b0,b1,...,bm,...,bM-1]T,其中PN1的长度为N,PN2的长度为M,要求M≥Kmax,其中
Figure BDA0002941003380000081
表示最大多径时延对应的符号个数,τmax表示信道的最大多径时延,Ts表示符号周期,
Figure BDA0002941003380000082
表示向下取整,序列a与b均表示模值为1的复序列:
Figure BDA0002941003380000083
其中,exp(·)表示指数运算,当b中的每一项均为±1,即θm∈{0,π}时,可以在基本不损失估计性能的情况下,极大降低估计复杂度;
发射机利用两个PN序列a与b进行下述的组帧操作,将a中的每一项an与b相乘,得到组帧输出的第n个向量:
xn=anb=[anb0,anb1,...,anbm,...,anbM-1]T
其中[·]T表示转置运算,xn表示组帧输出的第n个向量,由于a的长度为N,故相乘后可得到N个向量x0,x1,...,xN-1
本发明采用的训练序列组帧方式如图2所示,将这些向量按顺序前后拼接起来,即可得到最终组帧后的训练序列输出x=[x0,x1,...,xnM+m,...,xNM-1]T,其可用下述公式表示:
Figure BDA0002941003380000091
其中,xnM+m=anbm表示第n个组帧向量的第m项。
步骤S102中,本发明所述信道参数估计与判决方法如图3所示,将组帧后的训练序列x经过低信噪比多径信道,得到接收输入信号r:
Figure BDA0002941003380000092
由于存在多径效应,接收信号向量为多个路径的信号经过不同延时后叠加得到的,造成符号间干扰,因此接收信号向量比发射信号向量更长,增加的长度为Kmax
定义判决向量为J=[J0,J1,...,Jk,...,JM-1]T,其中每一项Jk通过下述计算方法得到,从接收信号中截取长度为NM的向量rk
rk=[rk,rk+1,...,rk+NM-1]T
再将该向量分为N组,每组长度为M,则第k个截取向量的第n组可表示为:
rk,n=[rk+nM,rk+nM+1,...,rk+nM+m,...,rk+(n+1)M-1]T
其中,k=0,1,...,M-1,n=0,1,...,N-1,m=0,1,...,M-1,将第k个截取向量的每一组转置后分别与PN序列b的共轭相乘,然后求计算结果的模值,并将得到的N个结果相加,其可用下述公式表示:
Figure BDA0002941003380000101
其中,Jk表示判决向量中的第k项,|·|表示取模,(·)*表示取共轭,则判决向量J可用下述公式表示:
Figure BDA0002941003380000102
步骤S103中,设判决门限向量γ=[γ01,...,γk,...,γM-1]T,其长度同样为M,由判决向量J的定义可知,该判决向量的每一项表示接收向量rk,n与PN序列b的互相关模值,而rk,n的帧格式中本身包含序列b,故会在二者对齐时存在一个明显的相关峰值,判决门限向量计算每个判决值周围点的平均值,再扩大一定倍数后,与该点的判决值进行比较,得到该点位置的判决结果,判决门限γk可用下述公式表示:
Figure BDA0002941003380000103
其中,u表示放大倍数,v表示取平均的范围。
步骤S104中,根据多径信道的定义,接收信号可用下述公式表示:
Figure BDA0002941003380000104
其中,rnM+m表示第n个组帧向量的第m个符号对应的接收信号,L表示多径数量,αl表示第l条路径的衰减,
Figure BDA0002941003380000105
表示第l条路径的相位,Kl表示第l条路径的时延,wnM+m表示噪声,由定义可知,Kl≤Kmax<M,故经过延时后第n个发射信号向量的起始位置,仍然落在该接收信号向量的持续时间范围内,即:
Figure BDA0002941003380000111
其中,|m-Kl|<M,设接收信号在l=l0处存在一多径,相应的路径衰减为
Figure BDA0002941003380000112
路径时延为
Figure BDA0002941003380000113
相位为
Figure BDA0002941003380000114
则判决向量中第
Figure BDA0002941003380000115
项可表示为:
Figure BDA0002941003380000116
其中,|·|表示取模运算,可以看到,由于an不参与后续相关运算,故可以直接提取出来,且|an|=1,因此在不增加复杂度的情况下,将判决向量的值扩大了N倍,在相关运算部分,只有当
Figure BDA0002941003380000117
即l=l0时,对应的多径才具有自相关特性,相关结果进一步将判决值扩大了M倍,而其他多径分量由于不具有相关性,故不会被放大,因此仅会存在l0多径分量一个峰值;
上面的分析可知,每个多径分量在判决向量中存在一个相关峰值,故可以比较判决向量与判决门限向量的大小,来找到超过判决门限的多径分量,则路径时延估计值
Figure BDA0002941003380000118
可用下述公式表示:
Figure BDA0002941003380000119
其中,
Figure BDA00029410033800001110
表示所有超过判决门限的路径时延估计值的集合,相应的路径衰减的估计值
Figure BDA00029410033800001111
可用下述公式表示:
Figure BDA00029410033800001112
其中,
Figure BDA0002941003380000121
表示所有超过判决门限的路径衰减估计值的集合。
步骤S105中,对于第l0个路径时延估计值
Figure BDA0002941003380000122
及路径衰减估计值
Figure BDA0002941003380000123
求接收信号与本地PN序列的相关结果,相关结果的幅角即为第l0个路径的相位估计值
Figure BDA0002941003380000124
其可用下述公式表示:
Figure BDA0002941003380000125
其中,arg[·]表示求幅角运算,可以看到,由于取消了模值运算,原本的标量相加变为矢量相加,故模值大的矢量在计算中占更大的权重,而训练序列的自相关性,将第l0条多径分量对应的相关值扩大了M倍,故在计算结果中的权重远大于噪声及其他多径分量,因此可以通过求相关结果的幅角来估计第l0条多径的相位,则路径相位的估计值
Figure BDA0002941003380000126
可用下述公式表示:
Figure BDA0002941003380000127
其中,
Figure BDA0002941003380000128
表示所有超过判决门限的路径相位估计值的集合。
本发明所述方法的具体应用方案包括但不限于以下两点:接收机的精定时同步、多径信号分集合并。
接收机在进行符号定时同步时,可能会受到噪声等干扰的影响,造成定时位置偏移,从而降低接收性能,而定时位置的偏移可以看做是接收信号中每一个多径分量均增加了相同的时延,这个时延会在本发明所述方法的路径时延估计部分估计出来,再通过信道均衡等方法对时延进行补偿,从而可以消除符号定时位置偏移对系统性能的影响。
在多径信道中,不同时间到达的接收信号路径分量之间具有一定的独立性,如果可以将这些多径分量分离出来,再按一定的方法进行合并,则可以通过时间分集技术,提高接收信号的可靠性,而本发明所述方法可以估计出每一个多径分量的路径时延、衰减与相位,接收机可以根据这些估计参数,将各多径分量分别提取出来,从而获得时间分集增益。
以上两点仅为本发明所述方法的具体应用方案实例而已,不作为本发明应用的限制条件。
下面结合仿真对本发明的应用效果作详细的描述。
1)仿真条件:
本发明的仿真实验使用MATLAB R2016a仿真软件,PN1选择长度为4的ZC复序列,PN2选择长度为255的m序列,仿真信噪比范围为-20dB~0dB。
2)仿真内容及结果分析:
仿真1:在多径衰落信道下,用本发明所述信道参数估计方法估计出的多径信道时延,与MengYe提出的多径信道时延估计方法,其估计准确率对比如图4所示,可以看到,当信噪比仅为-10dB时,本发明所述信道时延估计准确率已经近似达到100%,因此本发明所述方法具有很强的抗干扰能力与信道适应性,而Meng所提出的方法比本发明的方法估计准确率稍差,且每次估计均需要计算MN长度相关窗口的相关值,计算复杂度远高于本发明所述方法。
仿真2:在多径衰落信道下,用本发明所述信道参数估计方法估计出的多径信道衰减,其估计均方误差如图5所示,可以看到,即使处在低信噪比信道环境中,本发明所述信道衰减估计方法仍然具有很高的估计精度。
结合图4和图5可见,本发明的适用于超低信噪比多径传输环境的信道参数估计方法,在非常恶劣的信道环境中,仍然可以保证很高的估计准确率与估计精度,有效避免了多径效应引起的码间干扰,且本发明仅需进行一次卷积运算与若干乘加运算即可实现,计算复杂度极低。综上所述,本发明所述方案对于低信噪比多径信道环境下的信道参数估计方案具有较优的整体性能。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (7)

1.一种适用于超低信噪比多径传输环境的信道参数估计方法,其特征在于,所述适用于超低信噪比多径传输环境的信道参数估计方法包括:
步骤一,产生两组PN序列PN1与PN2,利用PN1与PN2对训练序列进行组帧;
步骤二,根据训练序列的帧结构特点,定义并计算接收信号多径分量检测的判决向量;
步骤三,利用接收信号的判决值,计算多径检测判决门限;
步骤四,对接收信号进行检测,估计超过判决门限多径的路径时延与路径衰减;
步骤五,估计上述各多径分量的相位偏移。
2.如权利要求1所述的适用于超低信噪比多径传输环境的信道参数估计方法,其特征在于,所述步骤一具体包括:产生两组PN序列PN1:a=[a0,a1,...,an,...,aN-1]T与PN2:b=[b0,b1,...,bm,...,bM-1]T,其中PN1的长度为N,PN2的长度为M,要求M≥Kmax,其中
Figure FDA0002941003370000011
表示最大多径时延对应的符号个数,τmax表示信道的最大多径时延,Ts表示符号周期,
Figure FDA0002941003370000012
表示向下取整,序列a与b均表示模值为1的复序列:
an=exp(jθn),n=0,1,...,N-1
bm=exp(jθm),m=0,1,...,M-1;
其中,exp(·)表示指数运算,当b中的每一项均为±1,即θm∈{0,π}时,可以在基本不损失估计性能的情况下,极大降低估计复杂度;
发射机利用两个PN序列a与b进行下述的组帧操作,将a中的每一项an与b相乘,得到组帧输出的第n个向量:
xn=anb=[anb0,anb1,...,anbm,...,anbM-1]T
其中,[·]T表示转置运算,xn表示组帧输出的第n个向量,由于a的长度为N,故相乘后可得到N个向量x0,x1,...,xN-1
将这些向量按顺序前后拼接起来,即可得到最终组帧后的训练序列输出x=[x0,x1,...,xnM+m,...,xNM-1]T,其可用下述公式表示:
Figure FDA0002941003370000021
其中,xnM+m=anbm表示第n个组帧向量的第m项。
3.如权利要求1所述的适用于超低信噪比多径传输环境的信道参数估计方法,其特征在于,所述步骤二具体包括:将组帧后的训练序列x经过低信噪比多径信道,得到接收输入信号r:
Figure FDA0002941003370000022
由于存在多径效应,接收信号向量为多个路径的信号经过不同延时后叠加得到的,造成符号间干扰,因此接收信号向量比发射信号向量更长,增加的长度为Kmax
定义判决向量为J=[J0,J1,...,Jk,...,JM-1]T,其中每一项Jk通过下述计算方法得到,从接收信号中截取长度为NM的向量rk
rk=[rk,rk+1,...,rk+NM-1]T
再将该向量分为N组,每组长度为M,则第k个截取向量的第n组可表示为:
rk,n=[rk+nM,rk+nM+1,...,rk+nM+m,...,rk+(n+1)M-1]T
其中,k=0,1,...,M-1,n=0,1,...,N-1,m=0,1,...,M-1,将第k个截取向量的每一组转置后分别与PN序列b的共轭相乘,然后求计算结果的模值,并将得到的N个结果相加,其可用下述公式表示:
Figure FDA0002941003370000023
其中,Jk表示判决向量中的第k项,|·|表示取模,(·)*表示取共轭,则判决向量J可用下述公式表示:
Figure FDA0002941003370000031
4.如权利要求1所述的适用于超低信噪比多径传输环境的信道参数估计方法,其特征在于,所述步骤三具体包括:设判决门限向量γ=[γ01,...,γk,...,γM-1]T,其长度同样为M,由判决向量J的定义可知,该判决向量的每一项表示接收向量rk,n与PN序列b的互相关模值,而rk,n的帧格式中本身包含序列b,故会在二者对齐时存在一个明显的相关峰值,判决门限向量计算每个判决值周围点的平均值,再扩大一定倍数后,与该点的判决值进行比较,得到该点位置的判决结果,判决门限γk可用下述公式表示:
Figure FDA0002941003370000032
其中,u表示放大倍数,v表示取平均的范围。
5.如权利要求1所述的适用于超低信噪比多径传输环境的信道参数估计方法,其特征在于,所述步骤四具体包括:根据多径信道的定义,接收信号可用下述公式表示:
Figure FDA0002941003370000033
其中,rnM+m表示第n个组帧向量的第m个符号对应的接收信号,L表示多径数量,αl表示第l条路径的衰减,
Figure FDA0002941003370000034
表示第l条路径的相位,Kl表示第l条路径的时延,wnM+m表示噪声,由定义可知,Kl≤Kmax<M,故经过延时后第n个发射信号向量的起始位置,仍然落在该接收信号向量的持续时间范围内,即:
Figure FDA0002941003370000035
其中,|m-Kl|<M,设接收信号在l=l0处存在一多径,相应的路径衰减为
Figure FDA0002941003370000036
路径时延为
Figure FDA0002941003370000037
相位为
Figure FDA0002941003370000038
则判决向量中第
Figure FDA0002941003370000039
项可表示为:
Figure FDA0002941003370000041
其中,|·|表示取模,可以看到,由于an不参与后续相关运算,故可以直接提取出来,且|an|=1,因此在不增加复杂度的情况下,将判决向量的值扩大了N倍,在相关运算部分,只有当
Figure FDA0002941003370000042
即l=l0时,对应的多径才具有自相关特性,相关结果进一步将判决值扩大了M倍,而其他多径分量由于不具有相关性,故不会被放大,因此仅会存在l0多径分量一个峰值;
上面的分析可知,每个多径分量在判决向量中存在一个相关峰值,故可以比较判决向量与判决门限向量的大小,来找到超过判决门限的多径分量,则路径时延估计值
Figure FDA0002941003370000043
可用下述公式表示:
Figure FDA0002941003370000044
其中,
Figure FDA0002941003370000045
表示所有超过判决门限的路径时延估计值的集合,相应的路径衰减的估计值
Figure FDA0002941003370000046
可用下述公式表示:
Figure FDA0002941003370000047
其中,
Figure FDA0002941003370000048
表示所有超过判决门限的路径衰减估计值的集合。
6.如权利要求1所述的适用于超低信噪比多径传输环境的信道参数估计方法,其特征在于,所述步骤五具体包括:对于第l0个路径时延估计值
Figure FDA0002941003370000049
及路径衰减估计值
Figure FDA00029410033700000410
求接收信号与本地PN序列的相关结果,相关结果的幅角即为第l0个路径的相位估计值
Figure FDA00029410033700000411
其可用下述公式表示:
Figure FDA0002941003370000051
其中,arg[·]表示求幅角运算,可以看到,由于取消了模值运算,原本的标量相加变为矢量相加,故模值大的矢量在计算中占更大的权重,而训练序列的自相关性,将第l0条多径分量对应的相关值扩大了M倍,故在计算结果中的权重远大于噪声及其他多径分量,因此可以通过求相关结果的幅角来估计第l0条多径的相位,则路径相位的估计值
Figure FDA0002941003370000052
可用下述公式表示:
Figure FDA0002941003370000053
其中,
Figure FDA0002941003370000054
表示所有超过判决门限的路径相位估计值的集合。
7.一种应用权利要求1~6任意一项所述适用于超低信噪比多径传输环境的信道参数估计方法的无线通信平台。
CN202110177061.5A 2021-02-09 2021-02-09 一种适用于超低信噪比多径传输环境的信道参数估计方法 Active CN112953863B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110177061.5A CN112953863B (zh) 2021-02-09 2021-02-09 一种适用于超低信噪比多径传输环境的信道参数估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110177061.5A CN112953863B (zh) 2021-02-09 2021-02-09 一种适用于超低信噪比多径传输环境的信道参数估计方法

Publications (2)

Publication Number Publication Date
CN112953863A true CN112953863A (zh) 2021-06-11
CN112953863B CN112953863B (zh) 2022-12-02

Family

ID=76244621

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110177061.5A Active CN112953863B (zh) 2021-02-09 2021-02-09 一种适用于超低信噪比多径传输环境的信道参数估计方法

Country Status (1)

Country Link
CN (1) CN112953863B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113890577A (zh) * 2021-09-28 2022-01-04 电子科技大学 基于信号相似度的快速分集方法
CN115767756A (zh) * 2022-11-02 2023-03-07 南京邮电大学 一种低信噪比下重要多径分量的辨识方法、系统及介质
WO2024032607A1 (zh) * 2022-08-10 2024-02-15 维沃移动通信有限公司 帧结构确定方法、装置、通信设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1925470A (zh) * 2005-08-29 2007-03-07 中兴通讯股份有限公司 一种时分双工移动通信系统的多径判决方法
WO2007030972A1 (fr) * 2005-09-15 2007-03-22 Zte Corporation Procede de traitement de bande de base pour ameliorer un rapport signal sur bruit selon un echantillonnage multiple
CN103338166A (zh) * 2013-07-01 2013-10-02 北京大学 一种改进的信道估计方法
US20180191533A1 (en) * 2016-05-06 2018-07-05 Amlogic (Shanghai) Co., Ltd. Channel estimation method
CN109981513A (zh) * 2019-01-24 2019-07-05 西安电子科技大学 室内高速大容量mimo-ofdm系统的时间与频率同步方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1925470A (zh) * 2005-08-29 2007-03-07 中兴通讯股份有限公司 一种时分双工移动通信系统的多径判决方法
WO2007030972A1 (fr) * 2005-09-15 2007-03-22 Zte Corporation Procede de traitement de bande de base pour ameliorer un rapport signal sur bruit selon un echantillonnage multiple
CN103338166A (zh) * 2013-07-01 2013-10-02 北京大学 一种改进的信道估计方法
US20180191533A1 (en) * 2016-05-06 2018-07-05 Amlogic (Shanghai) Co., Ltd. Channel estimation method
CN109981513A (zh) * 2019-01-24 2019-07-05 西安电子科技大学 室内高速大容量mimo-ofdm系统的时间与频率同步方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
QIUHONG LIU等: "《Joint Blind Estimation of PN Codes and Channels for Long-Code DSSS Signals in Multiple Paths at Low SNR》", 《IEEE》 *
吴钊等: "基于宽带突发单载波频域均衡传输的时域精细信道估计方法", 《电子与信息学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113890577A (zh) * 2021-09-28 2022-01-04 电子科技大学 基于信号相似度的快速分集方法
WO2024032607A1 (zh) * 2022-08-10 2024-02-15 维沃移动通信有限公司 帧结构确定方法、装置、通信设备及存储介质
CN115767756A (zh) * 2022-11-02 2023-03-07 南京邮电大学 一种低信噪比下重要多径分量的辨识方法、系统及介质

Also Published As

Publication number Publication date
CN112953863B (zh) 2022-12-02

Similar Documents

Publication Publication Date Title
CN112953863B (zh) 一种适用于超低信噪比多径传输环境的信道参数估计方法
KR100757696B1 (ko) 개선된 신호 포착과 프로세싱을 지닌 코드 분할 다중 접속 시스템 및 오퍼레이션 방법
KR20070007029A (ko) 부호 특정 결합을 이용한 ds-cdma 간섭 억제를 위한방법 및 장치
KR20060047917A (ko) 잡음 전력 추정 장치, 잡음 전력 추정 방법 및 신호 검출장치
US6614836B1 (en) Biased-corrected rake receiver for direct sequence spread spectrum waveform
JP5038410B2 (ja) 汎用rake受信機における干渉推定のための方法と装置
JPWO2007088953A1 (ja) 等化装置及び等化方法
US8290105B2 (en) Signal reception device and method of signal reception timing detection
CN109450489A (zh) 一种扩频水声通信的导频序列干扰抵消方法
CN111953380B (zh) 基于范数拟合的非周期长码直扩信号时延估计方法及系统
Ahmed et al. Reduced-rank adaptive multiuser detection in hybrid direct-sequence time-hopping ultrawide bandwidth systems
JP2003218826A (ja) 直交周波数分割多重信号の受信方式及び受信機
CN113259291A (zh) 利用水声连续信号动态多普勒跟踪实现的相位补偿方法
CN109150776B (zh) 一种短码直扩信号伪码序列盲估计方法
CN107465472B (zh) 一种基于路径合成的多径时延估计方法
CN116418364A (zh) 一种基于扩频序列信号的低虚警率检测方法
KR20060066255A (ko) 상관 계수의 대칭성을 이용한 상관 장치 및 그 방법
Manzoor et al. Novel SNR estimation technique in Wireless OFDM systems
CN111327547B (zh) 一种具有阈值去噪处理的时域信道估计装置
EP2076969A1 (en) A mehtod for simplifying the calculatinos for pre-whitening in a g-rake receiver,
JP2005151396A (ja) 受信装置および受信制御方法
US20190393910A1 (en) Timing estimation device and timing estimation method
CN112688889B (zh) 一种无人机测控系统的信道估计方法和装置
CN110768687B (zh) 一种mc-cdma接收机多径、异址干扰联合抑制方法及系统
CN110233811B (zh) 一种改进的ufmc系统干扰消除的联合方案

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