CN114615122A - 通信信号的频偏确定方法及装置 - Google Patents
通信信号的频偏确定方法及装置 Download PDFInfo
- Publication number
- CN114615122A CN114615122A CN202210114597.7A CN202210114597A CN114615122A CN 114615122 A CN114615122 A CN 114615122A CN 202210114597 A CN202210114597 A CN 202210114597A CN 114615122 A CN114615122 A CN 114615122A
- Authority
- CN
- China
- Prior art keywords
- frequency offset
- sequence
- change rate
- trs
- determining
- 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
- 238000004891 communication Methods 0.000 title claims abstract description 84
- 238000000034 method Methods 0.000 title claims abstract description 84
- 230000008859 change Effects 0.000 claims abstract description 145
- 230000006870 function Effects 0.000 claims description 48
- 238000004364 calculation method Methods 0.000 claims description 39
- 239000011159 matrix material Substances 0.000 claims description 21
- 238000012545 processing Methods 0.000 claims description 21
- 238000004590 computer program Methods 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 description 83
- 230000005540 biological transmission Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 11
- 238000007476 Maximum Likelihood Methods 0.000 description 9
- 125000004122 cyclic group Chemical group 0.000 description 7
- 230000001419 dependent effect Effects 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 230000001413 cellular effect Effects 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 230000004927 fusion Effects 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 239000013598 vector Substances 0.000 description 4
- 238000007792 addition Methods 0.000 description 3
- 238000005562 fading Methods 0.000 description 3
- 201000003042 peeling skin syndrome Diseases 0.000 description 3
- 229920001467 poly(styrenesulfonates) Polymers 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000004929 transmission Raman spectroscopy Methods 0.000 description 3
- 238000009795 derivation Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000010267 cellular communication Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2655—Synchronisation arrangements
- H04L27/2657—Carrier synchronisation
- H04L27/266—Fine or fractional frequency offset determination and synchronisation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本申请提供一种通信信号的频偏确定方法及装置,其中方法包括:基于跟踪参考信号TRS确定通信信号的频偏变化率估计值;基于所述频偏变化率估计值确定所述通信信号的频偏估计值。本申请实施例提供的通信信号的频偏确定方法及装置,通过TRS信号确定频偏变化率估计值,然后根据频偏变化率估计值确定频偏估计值,以便于精确地跟踪通信信号的频偏,提高载波频率精度。
Description
技术领域
本申请涉及通信技术领域,尤其涉及一种通信信号的频偏确定方法及装置。
背景技术
当前正交频分复用(Orthogonal Frequency Division Multiplexing,OFDM)系统中成熟的频谱估计算法大多将通信信号的多普勒频偏视为恒定值,而在蜂窝高速移动场景和星地融合场景这类高动态场景下,多普勒频偏具有快速时变性。为了在高动态场景下得到更精确的频偏估计值,不仅需要估计接收信号的载波频率,还需要估计接收信号的频偏变化率。
目前OFDM系统中对多普勒频偏变化率估计的关注度不高,有二阶差分(SecondOrder Difference,SOD)算法、多普勒频偏联合其一阶变化率的最大似然估计(MaximumLikelihood,ML)算法和基于第五代移动通信(the 5th Generation MobileCommunication,5G)主同步信号(primary synchronization signal,PSS)的一阶频偏变化率估计算法三种算法。
SOD算法基于同一OFDM符号做频偏变化率估计,频偏在一个OFDM符号内的变化非常细微,因此SOD算法的估计精度较低。ML算法的精度依赖于区间长度和步长大小,并且算法复杂度高。基于PSS的频谱变化率估计算法提升了估计精度,但仍然达不到10-8这个量级。
发明内容
针对现有技术存在的上述问题,本申请提供一种通信信号的频偏确定方法及装置。
第一方面,本申请提供一种通信信号的频偏确定方法,包括:
基于跟踪参考信号TRS确定通信信号的频偏变化率估计值;
基于所述频偏变化率估计值确定所述通信信号的频偏估计值。
可选地,所述基于跟踪参考信号TRS确定通信信号的频偏变化率估计值,包括:
对第一接收TRS信号的第一接收序列和所述第一接收TRS信号对应的第一发送序列进行共轭相关处理,获取第一序列;
对第二接收TRS信号的第二接收序列和所述第二接收TRS信号对应的第二发送序列进行共轭相关处理,获取第二序列;
基于所述第一序列和所述第二序列确定第一频偏变化率估计值;
所述第一接收TRS信号和所述第二接收TRS信号位于相邻帧中的相同位置,所述第一接收序列和所述第二接收序列不同。
可选地,所述基于所述第一序列和所述第二序列确定第一频偏变化率估计值,包括:
对所述第一序列和所述第二序列进行共轭相关处理,获得第三序列;
对所述第三序列中的前N/2项和所述第三序列中的后N/2项进行共轭相关处理,确定所述第一频偏变化率估计值,N为正整数。
可选地,所述获得第三序列之后,还包括:
基于所述第三序列的相位信息构造第一损失函数;
通过最小二乘法对所述第一损失函数的回归系数进行求解,确定第二频偏变化率估计值。
可选地,所述基于所述第三序列的相位信息构造第一损失函数满足如下计算公式:
其中,J1表示所述第一损失函数,N表示快速傅里叶变换点数,A表示回归系数矩阵,S表示自变量矩阵,D(n)表示所述第三序列,L表示所述第一接收TRS信号和所述第二接收TRS信号之间的采样点数,α表示一阶频偏变化率,表示所述第一接收TRS信号和所述第二接收TRS信号的相位差。
可选地,所述确定第二频偏变化率估计值之后,还包括:
基于所述第一接收TRS信号和所述第二接收TRS信号的幅度设置权值,更新所述第一损失函数为第二损失函数;
通过最小二乘法对所述第二损失函数的回归系数进行求解,确定第三频偏变化率估计值。
可选地,所述基于所述第一接收TRS信号和所述第二接收TRS信号的幅度设置权值,满足如下计算公式:
其中,wn表示所述权值,y(ls-1,n)表示所述第一接收TRS信号的第一接收序列,y(ls,n)表示所述第二接收TRS信号的第二接收序列,ls-1表示s-1时刻的第l个正交频分复用符号,ls表示s时刻的第l个正交频分复用符号,N表示快速傅里叶变换点数。
可选地,所述第三频偏变化率估计值满足如下计算公式:
第二方面,本申请还提供一种通信信号的频偏确定装置,包括:
第一确定模块,用于基于跟踪参考信号TRS确定通信信号的频偏变化率估计值;
第二确定模块,用于基于所述频偏变化率估计值确定所述通信信号的频偏估计值。
第三方面,本申请还提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如第一方面所述任一种通信信号的频偏确定方法的步骤。
本申请还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现如第一方面所述任一种通信信号的频偏确定方法的步骤。
本申请还提供一种计算机程序产品,包括计算机程序,所述计算机程序被处理器执行时实现如第一方面所述任一种通信信号的频偏确定方法的步骤。
本申请实施例提供的通信信号的频偏确定方法及装置,通过TRS信号确定频偏变化率估计值,然后根据频偏变化率估计值确定频偏估计值,以便于精确地跟踪通信信号的频偏,提高载波频率精度。
附图说明
为了更清楚地说明本申请或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以基于这些附图获得其他的附图。
图1是本申请实施例提供的通信信号的频偏确定方法的流程示意图之一;
图2是本申请实施例提供的TRS在时频资源格上的位置示意图;
图3是本申请实施例提供的不同算法在不同信噪比下的性能对比示意图之一;
图4是本申请实施例提供的不同算法在不同信噪比下的性能对比示意图之二;
图5是本申请实施例提供的基于变化率的频偏确定方法和传统静态频偏确定方法的频偏跟踪性能的对比示意图;
图6是本申请实施例提供的通信信号的频偏确定装置的结构示意图;
图7是本申请实施例提供的电子设备的结构示意图。
具体实施方式
为了更好地对本申请实施例中的技术方案进行描述,下面对相关知识进行介绍。
(1)多普勒频偏
多普勒频偏(Doppler Shift)是指当移动台以恒定的速率沿某一方向移动时,由于传播路程差的原因,会造成相位和频率的变化,通常将这种变化称为多普勒频偏。
(2)蜂窝高速移动场景和星地融合场景下的多普勒频偏
在蜂窝高速移动场景和星地融合场景下,通信信号的多普勒频偏不仅大,而且具有变化率。如表1所示,以低轨卫星场景为例,表1中给出了低轨卫星场景下不同高度、载波频率下最大多普勒频偏和最大多普勒频偏变化率参考值。
表1:低轨卫星场景下的多普勒频偏及其变化率
星地融合通信沿用了地面蜂窝通信中的OFDM系统,对频率偏移非常敏感,小数倍频偏会破坏子载波之间的正交性,造成子载波间的干扰,严重降低系统性能;整数倍频偏不会引起子载波间的干扰,但是会导致接收端恢复的频域数据码元序列循环位移,使得系统的误比特率升高。
当前OFDM系统中成熟的频偏估计算法大多都是将通信信号的频偏视为恒定值,这类算法在静止或低动态的场景下较为适用。而在蜂窝高速移动场景和星地融合场景这类高动态场景中,频偏具有时变性,因此,为了在高动态场景下得到更精确的频偏估计值,不仅需要估计接收信号的载波频率,还需要估计接收信号的频偏变化率,以便于精确地跟踪通信信号的频偏、提高载波频率精度。
目前OFDM系统中对多普勒频偏变化率估计的关注度不高,有SOD算法、ML算法和基于PSS的频偏变化率估计算法三种算法。
(3)归一化频偏变化率
归一化频偏变化率α可以表示为如下公式:
其中,α表示相邻采样点之间的多普勒频偏数值相对于子载波间隔的归一化,Δfdoppler表示每秒内多普勒频率的变化数值,Ts表示采样间隔,Δfsubcarrier表示子载波间隔。
(4)SOD算法
SOD算法通过两次差分,从循环前缀中提取出相位,算法如下:
令y(n)为接收OFDM符号循环前缀中的信号,其对应于OFDM符号后半部分接收信号中的y(n+N),则y(n)和y(n+N)可以表示为如下公式:
(5)ML算法
定义I和I′,分别表示第l个符号的循环前缀和第l个符号后Ncp个数据,I和I′分别表示为如下计算公式:
I={0,1,...,Ncp-1}
I′={N,N+1...,N+Ncp-1}
则OFDM符号中采样点的特征满足如下计算公式:
最大似然函数可以被构造为满足如下计算公式:
由于∏nf(x(n))不受ε0和α影响,因此上式可以被转化为如下计算公式:
其中,f(x(n),x(n+N))为二维的复高斯分布概率密度函数,满足如下计算公式:
其中,ρ为x(n)与x(n+N)的相关表示,满足如下计算公式:
可以得到|ρ|满足如下计算公式:
f(x(n))和f(x(n+N))均为一维高斯分布,分别表示为如下计算公式:
因此,将以上一维高斯分布与二维高斯分布代入Λ(ε0,α),可以得到如下计算公式:
去掉上式中与ε0和α独立的分量,最大似然函数可以简化为如下计算公式:
(6)基于PSS的频偏变化率估计算法
假设在同一个OFDM符号上多普勒频偏视作不变,利用两个相邻且间隔为5ms的PSS信号获取频偏变化率。
令相邻两个PSS接收信号ym-1(n)和ym(n)表示分别表示为如下计算公式:
εm=εm-1+αL
其中,L表示相邻两个PSS的采样点数。将相邻PSS进行共轭相关处理,满足如下计算公式:
将R(n)的前后半段共轭相关得到如下计算公式:
最终可以求得频偏变化率估计值满足如下公式:
(7)三种不同算法的对比
以表1中多普勒频偏变化率为8.16kHz/s的情况为例,假设子载波间隔为15kHz,子载波数目为1024,基于归一化频偏变化率的计算公式确定归一化频偏变化率如下:
归一化频偏变化率为3.5×10-8,即相邻采样点间频率变化数值为子载波间隔的3.5×10-8倍,
①SOD算法是基于同一OFDM符号中循环前缀和数据部分之间的频率变化来设计算法的,假设归一化频偏变化率为3.5×10-8,那么在同一OFDM符号内相邻采样点之间的归一化频率变化为3.5×10-8,频偏在一个OFDM符号内的变化非常细微,因此SOD算法的计算精度较低。在一次仿真中设置归一化频偏变化率量级为10-1,假设归一化频偏变化率为0.1,子载波间隔为2232Hz,子载波数目为4096,则频偏变化率的值约为2.04GHz/s,与表1中示出的频偏变化率不相符。
②ML算法是在一定的区间内,设置步长来搜索最佳的频偏变化率估计值,算法复杂度高,且计算结果与设置的区间的长度、步长的大小相关联。
③基于PSS的频偏变化率估计算法通过在更大的时间尺度上累积频偏变化率对参考信号相位的影响,使其对小范围的频偏变化率更加敏感,提高了频偏变化率的估计精度,但是仍然达不到10-8的量级。
针对现有技术存在的上述问题,本申请实施例提供一种通信信号的频偏确定方法及装置,利用TRS信号来估计频偏变化率,进一步估计频偏,降低算法复杂度的同时提升频偏估计精度。
为使本申请的目的、技术方案和优点更加清楚,下面将结合本申请中的附图,对本申请中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
图1是本申请实施例提供的通信信号的频偏确定方法的流程示意图之一,如图1所示,本申请实施例提供一种通信信号的频偏确定方法,执行主体为用于通信信号的频偏确定的装置,该方法至少包括以下步骤:
步骤101、基于跟踪参考信号确定通信信号的频偏变化率估计值。
具体地,终端在配置的时候会在相应的时频位置去测量信道状态信息参考信号(Channel State Information-Reference Signal,CSI-RS),获取信道状态信息。用于时频跟踪的CSI-RS信号经过特殊配置,被称作跟踪参考信号(Tracing Reference Signal,TRS)。TRS是一种多周期的CSI-RS。
当终端处于空闲状态或初始状态时,使用主同步信号(primary synchronizationsignal,PSS)和辅同步信号(secondary synchronization signal,SSS)实现时频同步,在无线资源控制(Radio Resource Control,RRC)连接建立后,对终端配置TRS,用于跟踪阶段的时频同步。
图2是本申请实施例提供的TRS在时频资源格上的位置示意图,如图2所示,在时间维度,一个时隙(Slot)内包含两个TRS,在一个时隙内两个TRS间隔是4个OFDM符号,频域上间隔4个子载波,TRS的发送周期为2uXpslot,其中,u为子载波间隔配置,u∈{0,1,2,3,4},Xp∈{10,20,40,80}。
本申请实施例中通过TRS信号来确定通信信号的频偏变化率估计值。
步骤102、基于频偏变化率估计值确定通信信号的频偏估计值。
具体地,为了在高动态场景下得到更精确的频偏估计值,不仅需要估计接收信号的载波频率,还需要估计接收信号的频偏变化率,以便于精确地跟踪通信信号的频偏,提高载波频率精度。
第l个OFDM符号第n个子载波上的多普勒频偏估计值ε(l,n)满足如下计算公式:
第l个OFDM符号第n个子载波上的时域接收TRS信号的接收序列y(l,n)的表达式如下:
其中,N表示快速傅里叶变换(Fast Fourier Transform,FFT)的点数,H(k)表示信道频域响应值,X(l,k)表示与接收TRS信号相对应的频域发送信号的发送序列,ε(l,n)表示频偏。相对应地,与接收TRS信号相对应的本地时域发送信号的发送序列x(l,n)的表达式如下:
其中,N表示FFT点数,X(l,k)表示与接收TRS信号相对应的频域发送信号的发送序列。
本申请实施例提供的通信信号的频偏确定方法,通过TRS信号确定频偏变化率估计值,然后根据频偏变化率估计值确定频偏估计值,以便于精确地跟踪通信信号的频偏,提高载波频率精度。
可选地,基于TRS信号确定通信信号的频偏变化率估计值,包括:
对第一接收TRS信号的第一接收序列和第一接收TRS信号对应的第一发送序列进行共轭相关处理,获取第一序列;
对第二接收TRS信号的第二接收序列和第二接收TRS信号对应的第二发送序列进行共轭相关处理,获取第二序列;
基于第一序列和第二序列确定第一频偏变化率估计值;
第一接收TRS信号和第二接收TRS信号位于相邻帧中的相同位置,第一接收序列和第二接收序列不同。
具体地,在5G帧结构中利用相邻帧中相同位置的TRS信号来确定频偏变化率估计值。
可选地,增大时间尺度,在当前帧中联合前一帧中对应位置的TRS信号来确定频偏变化率,例如TRS的发送周期为10ms。
第一接收TRS信号为当前帧中的接收TRS信号,第二接收TRS信号为相邻帧中相同位置的接收TRS信号,相邻帧中对应位置的TRS信号的序列不同,即第一接收序列和第二接收序列不同。
当前帧中第一接收TRS信号的第一接收序列y(ls-1,n)的表达式如下:
其中,N表示FFT点数,H(k)表示信道频域响应值,X(l,k)表示与第一接收TRS信号相对应的频域发送信号的发送序列,εs-1表示s-1时刻第一接收TRS信号的频偏,表示s-1时刻第一接收TRS信号的初始相位。
前一帧中相同位置的第二接收TRS信号的第二接收序列y(ls,n)的表达式如下:
其中,N表示FFT点数,H(k)表示信道频域响应值,X(l,k)表示与第二接收TRS信号相对应的频域发送信号的发送序列,εs表示s时刻第二接收TRS信号的频偏,表示s时刻第二接收TRS信号的初始相位。
第一接收TRS信号和第二接收TRS信号的频偏之间满足如下计算公式:
εs=εs-1+αL
其中,α表示一阶频偏变化率,L表示第一接收TRS信号和第二接收TRS信号之间的采样点数。
对第一接收TRS信号的第一接收序列y(ls-1,n)和第一接收TRS信号对应的第一发送序列x(ls-1,n)进行共轭相关处理,获取第一序列X1,满足如下计算公式:
X1=x*(ls-1,n)y(ls-1,n)
对第二接收TRS信号的第二接收序列y(ls,n)和第二接收TRS信号对应的第二发送序列x(ls,n)进行共轭相关处理,获取第二序列X2,满足如下计算公式:
X2=x*(ls,n)y(ls,n)
将第一接收TRS信号的第一接收序列和第一接收TRS信号对应的第一发送序列做共轭相关得到第一序列,将第二接收TRS信号的第二接收序列和第二接收TRS信号对应的第二发送序列做共轭相关得到第二序列之后,对得到的第一序列和第二序列进行处理,确定第一频偏变化率估计值。
可选地,基于第一序列和第二序列确定第一频偏变化率估计值,包括:
对第一序列和第二序列进行共轭相关处理,获得第三序列;
对第三序列中的前N/2项和第三序列中的后N/2项进行共轭相关处理,确定第一频偏变化率估计值,N为正整数。
具体地,本申请实施例提供一种基于TRS共轭相关的频偏变化率估计算法(TRS共轭相关算法),具体包括:
获取第一序列X1和第二序列X2之后,将第一序列和第二序列做共轭相关,得到第三序列D(n),满足如下计算公式:
其中,C为实数。
对第三序列D(n)的前N/2项和第三序列中的后N/2项进行共轭相关处理,满足如下计算公式:
其中,C′为实数,N为正整数。
本申请实施例提供的通信信号的频偏确定方法,利用5G帧结构中相邻帧中相同位置的TRS信号,通过对第一TRS接收信号的第一接收序列和对应的第一发送序列共轭相关得到第一序列,通过对第二TRS接收信号的第二接收序列和对应的第二发送序列共轭相关得到第二序列之后,将第一序列和第二序列共轭相关得到第三序列,对第三序列的前后半段共轭相关,获得第一频偏变化率估计值,降低了频偏变化率估计过程的计算复杂度,同时提升了计算精度。
可选地,获得第三序列之后,还包括:
基于第三序列的相位信息构造第一损失函数;
通过最小二乘法对第一损失函数的回归系数进行求解,确定第二频偏变化率估计值。
具体地,最小二乘法(Least Squares Method,LS)通过最小化误差的平方和寻找数据的最佳函数匹配。利用第三序列的相位信息构造第一损失函数,然后通过最小二乘法对第一损失函数的回归系数进行优化求解,确定第二频偏变化率估计值。
最小二乘法的原理如下:假设自变量矩阵设自变量矩阵S中包含u个输入的特征向量,第i个特征向量的表达式如下:
s(i)=[1 si,1 si,2 ... si,(p-1)]T,p=1,2,3,…,u
其中,上标T表示转置。相对应地,自变量矩阵S的表达式如下:
特征向量中每个元素对应一个回归系数,系数矩阵A的表达式如下:
A=[a0 a1 a2 ... ap-1]T
第i个输入特征向量对应因变量矩阵Υ中第i个因变量γ(i),因变量矩阵Υ的表达式如下:
Υ=[γ(1) γ(2) γ(3) ... γ(u-1)]T
应因变量矩阵Υ、自变量矩阵S和系数矩阵A之间存在线性关系,表达式如下:
Υ=SA
最小二乘法的最终目的是求解出回归系数矩阵A,使得损失函数取最小值。系统损失函数J的表达式如下:
使得损失函数J取最小值的回归系数矩阵A即为求解得到的回归系数A的估计值。因此损失函数J为关于a的函数,其表达式如下:
J=F(a0,a1,a2,...,ap-1)
本申请实施例提供一种基于最小二乘法的频偏变化率估计算法(LS算法),具体包括:
将最小二乘法应用到频偏变化率的估计过程中,利用第三序列的相位信息构造第一损失函数,然后通过最小二乘法对第一损失函数的回归系数进行优化求解,确定第二频偏变化率估计值。第三序列D(n)的表达式如下:
其中,C为实数,n∈{0,1,...N-1},因此可以利用第三序列的相位构造线性方程如下:
γ=[Dl,0 Dl,1 Dl,2 ... Dl,N-1]T
自变量矩阵S的表达式如下:
回归系数矩阵A的表达式如下:
因此,第一损失函数的表达式如下:
其中,J1表示第一损失函数,N表示FFT点数,A表示回归系数矩阵,S表示自变量矩阵,D(n)表示第三序列,L表示第一接收TRS信号和第二接收TRS信号之间的采样点数,α表示一阶频偏变化率,表示第一接收TRS信号和第二接收TRS信号的相位。
令
解得第二频偏变化率估计值的表达式如下:
本申请实施例提供的通信信号的频偏确定方法,通过TRS共轭相关算法得到第三序列后,利用第三序列的相位构造线性方程,通过最小二乘法求解出第二频偏变化率估计值,进一步提升了频偏确定过程中频偏变化率估计的计算精度。
可选地,确定第二频偏变化率估计值之后,还包括:
基于第一接收TRS信号和第二接收TRS信号的幅度设置权值,更新第一损失函数为第二损失函数;
通过最小二乘法对第二损失函数的回归系数进行求解,确定第三频偏变化率估计值。
具体地,本申请实施例提供一种基于加权最小二乘法的频偏变化率估计算法(WEIGHTED-LS算法),具体包括:
在通过最小二乘法求解出第二频偏变化率估计值之后,考虑到噪声对通信信号的影响,可以根据接收TRS信号的幅度信息设置权值(Weighted),根据第一接收TRS信号和第二接收TRS信号的幅度平方和获得加权最小二乘的权值,进一步提升频偏变化率估计值的估计精度。
权值wn满足如下计算公式:
其中,N表示FFT点数。
基于第一接收TRS信号的第二接收TRS信号的幅度信息设置权值wn,满足如下计算公式:
其中,wn表示权值,y(ls-1,n)表示第一接收TRS信号的第一接收序列,y(ls,n)表示第二接收TRS信号的第二接收序列,ls-1表示s-1时刻的第l个OFDM,ls表示s时刻的第l个正交频分复用符号,N表示FFT点数。
在这种情况下,更新第一损失函数为第二损失函数,第二损失函数满足如下计算公式:
其中,J2表示第二损失函数,N表示FFT点数,L表示第一接收TRS信号和第二接收TRS信号之间的采样点数,α表示一阶频偏变化率,表示第一接收TRS信号和第二接收TRS信号的相位差,D(n)表示第三序列。
本申请实施例提供的通信信号的频偏确定方法,通过相邻接收帧中TRS信号幅度平方和获得加权最小二乘的权值,进一步提升了频偏确定过程中频偏变化率估计的计算精度。
下面对本申请实施例提出的通信信号的频偏确定方法的估计精度进行分析说明,图3是本申请实施例提供的不同算法在不同信噪比下的性能对比示意图之一,如图3所示,本申请实施例对比了在高斯加性白噪声(Additive White Gaussian Noise,AWGN)信道下SOD算法、基于PSS的频偏变化率估计算法(PSS共轭相关算法)、基于TRS共轭相关的频偏变化率估计算法(TRS共轭相关算法)、基于最小二乘法的频偏变化率估计算法(LS算法)、基于加权最小二乘法的频偏变化率估计算法(WEIGHTED-LS算法)的性能,以均方误差(MeanSquared Error,MSE)作为评价指标,横坐标表示信噪比,其单位为dB,纵坐标为MSE。
MSE满足如下计算公式:
表2:仿真参数
从图3中可以看出,在AWGN信道下,随着信噪比的增大,各算法的性能均有所提升。其中,SOD算法的算法性能非常差,精度远达不到1×10-8。PSS共轭相关算法通过在更大的时间尺度上累积多普勒频偏变化率对参考信号相位的影响,但性能表现仍然较差。
进一步地,对TRS共轭相关算法、LS算法和WEIGHTED-LS算法的性能示意图在纵坐标上进行放大。可以看出,本申请实施例提供的三种频偏变化率估计算法在低信噪比下MSE均能达到10-1。当MSE达到10-16时,WEIGHTED-LS算法的性能优于LS算法4dB,LS算法的性能优于TRS共轭相关算法2dB。
图4是本申请实施例提供的不同算法在不同信噪比下的性能对比示意图之二,如图4所示,本申请实施例对比了在衰落信道下SOD算法、基于PSS的频偏变化率估计算法(PSS共轭相关算法)、基于TRS共轭相关的频偏变化率估计算法(TRS共轭相关算法)、基于最小二乘法的频偏变化率估计算法(LS算法)、基于加权最小二乘法的频偏变化率估计算法(WEIGHTED-LS算法)的性能,以MSE作为评价指标,横坐标表示信噪比,其单位为dB,纵坐标为MSE。仿真参数如表2所示。
从图4中可以看出,在衰落信道下,随着信噪比的增大,各算法的性能均有所提升。在MSE达到10-16时,TRS共轭相关算法的信噪比为14dB,相较于AWGN信道下MSE达到10-1时的信噪比高出4dB。
当MSE达到10-16时,WEIGHTED-LS算法的性能优于LS算法2dB,LS算法的性能优于TRS共轭相关算法2dB。
图5是本申请实施例提供的基于变化率的频偏确定方法和传统静态频偏确定方法的频偏跟踪性能的对比示意图,如图5所示,横坐标表示时间(Time),其单位为毫秒(ms),纵坐标为归一化频偏。在基于变化率的频偏确定方法中,用于跟踪频偏变化率取值为衰落信道下WEIGHTED-LS算法在信噪比为14dB时的平均频偏变化率估计值。
从图5中可以看出,在频偏具有时变性的高动态通信场景下,使用传统静态频偏确定方法将导致频偏估计误差随时间的变化而不断增大;使用基于频偏变化率的频偏确定方法能够不断跟踪时变性频偏,在高动态通信场景下更为适用。
下面对本申请实施例提出的通信信号的频偏确定方法的算法复杂度进行分析说明,1次复数乘法相当于4次实数乘法和2次加法,1次复数加法相当于2次实数加法。算法复杂度如表3所示,TRS共轭相关算法的计算复杂度最低,LS算法的复杂度相对有所提高,而WEIGHTED-LS算法中进一步提高了算法的复杂度。可以看出,提升频偏变化率估计值的精度牺牲了算法的复杂度。
表3:算法复杂度
下面对本申请提供的通信信号的频偏确定装置进行描述,下文描述的通信信号的频偏确定装置与上文描述的通信信号的频偏确定方法可相互对应参照。
图6是本申请实施例提供的通信信号的频偏确定装置的结构示意图,如图6所示,本申请实施例提供一种通信信号的频偏确定装置,该装置包括:
第一确定模块601,用于基于跟踪参考信号TRS确定通信信号的频偏变化率估计值;
第二确定模块602,用于基于所述频偏变化率估计值确定所述通信信号的频偏估计值。
可选地,所述第一确定模块还用于:
对第一接收TRS信号的第一接收序列和所述第一接收TRS信号对应的第一发送序列进行共轭相关处理,获取第一序列;
对第二接收TRS信号的第二接收序列和所述第二接收TRS信号对应的第二发送序列进行共轭相关处理,获取第二序列;
基于所述第一序列和所述第二序列确定第一频偏变化率估计值;
所述第一接收TRS信号和所述第二接收TRS信号位于相邻帧中的相同位置,所述第一接收序列和所述第二接收序列不同。
可选地,所述基于所述第一序列和所述第二序列确定第一频偏变化率估计值,包括:
对所述第一序列和所述第二序列进行共轭相关处理,获得第三序列;
对所述第三序列中的前N/2项和所述第三序列中的后N/2项进行共轭相关处理,确定所述第一频偏变化率估计值,N为正整数。
可选地,所述获得第三序列之后,还包括:
基于所述第三序列的相位信息构造第一损失函数;
通过最小二乘法对所述第一损失函数的回归系数进行求解,确定第二频偏变化率估计值。
可选地,所述基于所述第三序列的相位信息构造第一损失函数满足如下计算公式:
其中,J1表示所述第一损失函数,N表示快速傅里叶变换点数,A表示回归系数矩阵,S表示自变量矩阵,D(n)表示所述第三序列,L表示所述第一接收TRS信号和所述第二接收TRS信号之间的采样点数,α表示一阶频偏变化率,表示所述第一接收TRS信号和所述第二接收TRS信号的相位差。
可选地,所述确定第二频偏变化率估计值之后,还包括:
基于所述第一接收TRS信号和所述第二接收TRS信号的幅度设置权值,更新所述第一损失函数为第二损失函数;
通过最小二乘法对所述第二损失函数的回归系数进行求解,确定第三频偏变化率估计值。
可选地,所述基于所述第一接收TRS信号和所述第二接收TRS信号的幅度设置权值,满足如下计算公式:
其中,wn表示所述权值,y(ls-1,n)表示所述第一接收TRS信号的第一接收序列,y(ls,n)表示所述第二接收TRS信号的第二接收序列,ls-1表示s-1时刻的第l个正交频分复用符号,ls表示s时刻的第l个正交频分复用符号,N表示快速傅里叶变换点数。
可选地,所述第三频偏变化率估计值满足如下计算公式:
在此需要说明的是,本申请实施例提供的上述装置,能够实现上述方法实施例所实现的所有方法步骤,且能够达到相同的技术效果,在此不再对本实施例中与方法实施例相同的部分及有益效果进行具体赘述。
图7是本申请实施例提供的电子设备的结构示意图,如图7所示,该电子设备可以包括:处理器(processor)701、通信接口(Communications Interface)702、存储器(memory)703和通信总线704,其中,处理器701,通信接口702,存储器703通过通信总线704完成相互间的通信。处理器701可以调用存储器703中的逻辑指令,以执行通信信号的频偏确定方法,该方法包括:
基于跟踪参考信号TRS确定通信信号的频偏变化率估计值;
基于所述频偏变化率估计值确定所述通信信号的频偏估计值。
此外,上述的存储器703中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
另一方面,本申请还提供一种计算机程序产品,所述计算机程序产品包括计算机程序,计算机程序可存储在非暂态计算机可读存储介质上,所述计算机程序被处理器执行时,计算机能够执行上述各方法所提供的通信信号的频偏确定方法,该方法包括:
基于跟踪参考信号TRS确定通信信号的频偏变化率估计值;
基于所述频偏变化率估计值确定所述通信信号的频偏估计值。
又一方面,本申请还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现以执行上述各方法提供的通信信号的频偏确定方法,该方法包括:
基于跟踪参考信号TRS确定通信信号的频偏变化率估计值;
基于所述频偏变化率估计值确定所述通信信号的频偏估计值。
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以基于实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是:以上实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。
Claims (10)
1.一种通信信号的频偏确定方法,其特征在于,包括:
基于跟踪参考信号TRS确定通信信号的频偏变化率估计值;
基于所述频偏变化率估计值确定所述通信信号的频偏估计值。
2.根据权利要求1所述的通信信号的频偏确定方法,其特征在于,所述基于跟踪参考信号TRS确定通信信号的频偏变化率估计值,包括:
对第一接收TRS信号的第一接收序列和所述第一接收TRS信号对应的第一发送序列进行共轭相关处理,获取第一序列;
对第二接收TRS信号的第二接收序列和所述第二接收TRS信号对应的第二发送序列进行共轭相关处理,获取第二序列;
基于所述第一序列和所述第二序列确定第一频偏变化率估计值;
所述第一接收TRS信号和所述第二接收TRS信号位于相邻帧中的相同位置,所述第一接收序列和所述第二接收序列不同。
3.根据权利要求2所述的通信信号的频偏确定方法,其特征在于,所述基于所述第一序列和所述第二序列确定第一频偏变化率估计值,包括:
对所述第一序列和所述第二序列进行共轭相关处理,获得第三序列;
对所述第三序列中的前N/2项和所述第三序列中的后N/2项进行共轭相关处理,确定所述第一频偏变化率估计值,N为正整数。
4.根据权利要求3所述的通信信号的频偏确定方法,其特征在于,所述获得第三序列之后,还包括:
基于所述第三序列的相位信息构造第一损失函数;
通过最小二乘法对所述第一损失函数的回归系数进行求解,确定第二频偏变化率估计值。
6.根据权利要求4所述的通信信号的频偏确定方法,其特征在于,所述确定第二频偏变化率估计值之后,还包括:
基于所述第一接收TRS信号和所述第二接收TRS信号的幅度设置权值,更新所述第一损失函数为第二损失函数;
通过最小二乘法对所述第二损失函数的回归系数进行求解,确定第三频偏变化率估计值。
9.一种通信信号的频偏确定装置,其特征在于,包括:
第一确定模块,用于基于跟踪参考信号TRS确定通信信号的频偏变化率估计值;
第二确定模块,用于基于所述频偏变化率估计值确定所述通信信号的频偏估计值。
10.一种电子设备,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1至8任一项所述通信信号的频偏确定方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210114597.7A CN114615122B (zh) | 2022-01-30 | 2022-01-30 | 通信信号的频偏确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210114597.7A CN114615122B (zh) | 2022-01-30 | 2022-01-30 | 通信信号的频偏确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114615122A true CN114615122A (zh) | 2022-06-10 |
CN114615122B CN114615122B (zh) | 2023-11-03 |
Family
ID=81859730
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210114597.7A Active CN114615122B (zh) | 2022-01-30 | 2022-01-30 | 通信信号的频偏确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114615122B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115715003A (zh) * | 2023-01-09 | 2023-02-24 | 四川创智联恒科技有限公司 | 一种高速移动通信系统的时偏调整方法及系统 |
US20230079512A1 (en) * | 2021-09-16 | 2023-03-16 | Apple Inc. | Frequency offset delta tracking for nr connected mode discontinuous reception carrier aggregation |
CN116633749A (zh) * | 2023-07-24 | 2023-08-22 | 北京智芯微电子科技有限公司 | 频偏估计方法、频偏补偿方法及其系统、存储介质和芯片 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060133527A1 (en) * | 2004-12-11 | 2006-06-22 | Heejung Yu | Residual frequency, phase, timing offset and signal amplitude variation tracking apparatus and methods for OFDM systems |
US20090154601A1 (en) * | 2007-12-12 | 2009-06-18 | Harris Corporation | Communications device and related method with improved acquisition estimates of frequency offset and phase error |
WO2016019657A1 (zh) * | 2014-08-07 | 2016-02-11 | 中兴通讯股份有限公司 | 一种频偏估计的方法和装置、存储介质 |
US20180323918A1 (en) * | 2017-05-05 | 2018-11-08 | Mediatek Inc. | Tracking Reference Signal And Framework Thereof In Mobile Communications |
CN110708267A (zh) * | 2019-09-30 | 2020-01-17 | 北京中科晶上科技股份有限公司 | 频偏信息估计值确定方法 |
CN111554321A (zh) * | 2020-04-20 | 2020-08-18 | 北京达佳互联信息技术有限公司 | 降噪模型训练方法、装置、电子设备及存储介质 |
CN112398764A (zh) * | 2020-11-05 | 2021-02-23 | 上海擎昆信息科技有限公司 | 一种联合dmrs和ptrs的频偏估计方法及系统 |
CN112422473A (zh) * | 2020-11-09 | 2021-02-26 | 上海微波技术研究所(中国电子科技集团公司第五十研究所) | 联合频偏估计策略的系统及方法 |
-
2022
- 2022-01-30 CN CN202210114597.7A patent/CN114615122B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060133527A1 (en) * | 2004-12-11 | 2006-06-22 | Heejung Yu | Residual frequency, phase, timing offset and signal amplitude variation tracking apparatus and methods for OFDM systems |
US20090154601A1 (en) * | 2007-12-12 | 2009-06-18 | Harris Corporation | Communications device and related method with improved acquisition estimates of frequency offset and phase error |
WO2016019657A1 (zh) * | 2014-08-07 | 2016-02-11 | 中兴通讯股份有限公司 | 一种频偏估计的方法和装置、存储介质 |
US20180323918A1 (en) * | 2017-05-05 | 2018-11-08 | Mediatek Inc. | Tracking Reference Signal And Framework Thereof In Mobile Communications |
CN110708267A (zh) * | 2019-09-30 | 2020-01-17 | 北京中科晶上科技股份有限公司 | 频偏信息估计值确定方法 |
CN111554321A (zh) * | 2020-04-20 | 2020-08-18 | 北京达佳互联信息技术有限公司 | 降噪模型训练方法、装置、电子设备及存储介质 |
CN112398764A (zh) * | 2020-11-05 | 2021-02-23 | 上海擎昆信息科技有限公司 | 一种联合dmrs和ptrs的频偏估计方法及系统 |
CN112422473A (zh) * | 2020-11-09 | 2021-02-26 | 上海微波技术研究所(中国电子科技集团公司第五十研究所) | 联合频偏估计策略的系统及方法 |
Non-Patent Citations (1)
Title |
---|
赵黎;冯新宇;: "基于训练序列的OFDM符号定时与频偏联合估计", 西安工程大学学报, no. 05 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20230079512A1 (en) * | 2021-09-16 | 2023-03-16 | Apple Inc. | Frequency offset delta tracking for nr connected mode discontinuous reception carrier aggregation |
US11792881B2 (en) * | 2021-09-16 | 2023-10-17 | Apple Inc. | Frequency offset delta tracking for NR connected mode discontinuous reception carrier aggregation |
CN115715003A (zh) * | 2023-01-09 | 2023-02-24 | 四川创智联恒科技有限公司 | 一种高速移动通信系统的时偏调整方法及系统 |
CN116633749A (zh) * | 2023-07-24 | 2023-08-22 | 北京智芯微电子科技有限公司 | 频偏估计方法、频偏补偿方法及其系统、存储介质和芯片 |
CN116633749B (zh) * | 2023-07-24 | 2023-12-29 | 北京智芯微电子科技有限公司 | 频偏估计方法、频偏补偿方法及其系统、存储介质和芯片 |
Also Published As
Publication number | Publication date |
---|---|
CN114615122B (zh) | 2023-11-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110290581B (zh) | 一种5g系统中的快速时频同步方法及终端 | |
CN114615122B (zh) | 通信信号的频偏确定方法及装置 | |
CN1802831B (zh) | 用于ofdm信号的自适应相位补偿的方法和装置 | |
US9571305B2 (en) | Channel estimation by time-domain parameter extraction | |
US8824269B2 (en) | Method and system for a reference signal (RS) timing loop for OFDM symbol synchronization and tracking | |
CN112714086B (zh) | 一种频偏估计方法及基站 | |
CN101018219B (zh) | 一种空频信号处理方法 | |
CN101534287A (zh) | 移动通信系统中进行载波频率偏移校正的方法和装置 | |
CN113678416B (zh) | 用于联合信道和相位噪声估计的循环先导序列 | |
US20120213315A1 (en) | Process for estimating the channel in a ofdm communication system, and receiver for doing the same | |
WO2021143644A1 (zh) | 一种载波相位跟踪方法及装置 | |
WO2009095410A2 (en) | Method for estimating a channel, and decoder, terminal, and computer program | |
CN111131119A (zh) | 一种正交频分复用系统高精度定时偏移估计的方法及装置 | |
CN114301745B (zh) | 载波频偏和采样频偏的确定方法及装置 | |
CN111131120B (zh) | 一种基于ml同步的高精度定时偏移估计方法及装置 | |
Noschese et al. | A low-complexity approach for time of arrival estimation in OFDM systems | |
CN108989261B (zh) | 一种通信系统的定时同步方法、装置及相关设备 | |
CN107835141B (zh) | 自相关与互相关结合的多段重复序列ofdm同步算法 | |
RU2658335C1 (ru) | Способ совместной оценки канала связи и мягкой демодуляции для COFDM-сигналов и устройство для его реализации | |
CN114584449B (zh) | 一种基于CAZAC序列的改进M-part定时同步与频偏方法 | |
CN111212014B (zh) | 一种基于跳频规则的nb-iot定时提前量估计方法 | |
Manzoor et al. | A novel noise power and SNR estimation in WiMAX systems | |
CN117118791B (zh) | 频偏估计方法、装置、计算机设备和存储介质 | |
Canet et al. | Low complexity time synchronization algorithm for OFDM systems with repetitive preambles | |
US20230396476A1 (en) | Radio transmission device and radio reception device |
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 |