CN107070818A - 基于毫米波mimo系统的多普勒频偏估计方法及装置 - Google Patents

基于毫米波mimo系统的多普勒频偏估计方法及装置 Download PDF

Info

Publication number
CN107070818A
CN107070818A CN201610887151.2A CN201610887151A CN107070818A CN 107070818 A CN107070818 A CN 107070818A CN 201610887151 A CN201610887151 A CN 201610887151A CN 107070818 A CN107070818 A CN 107070818A
Authority
CN
China
Prior art keywords
mrow
msub
msup
theta
transmitting
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
CN201610887151.2A
Other languages
English (en)
Other versions
CN107070818B (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.)
Beijing University of Posts and Telecommunications
Original Assignee
Beijing University of Posts and Telecommunications
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 Beijing University of Posts and Telecommunications filed Critical Beijing University of Posts and Telecommunications
Priority to CN201610887151.2A priority Critical patent/CN107070818B/zh
Priority to PCT/CN2016/106880 priority patent/WO2018068365A1/zh
Publication of CN107070818A publication Critical patent/CN107070818A/zh
Application granted granted Critical
Publication of CN107070818B publication Critical patent/CN107070818B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202—Channel estimation
    • H04L25/022—Channel estimation of frequency response
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413—MIMO systems
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202—Channel estimation
    • H04L25/0222—Estimation of channel variability, e.g. coherence bandwidth, coherence time, fading frequency
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Radio Transmission System (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明实施例提供了一种基于毫米波MIMO系统的多普勒频偏估计方法及装置，方法包括：根据毫米波MIMO系统中的接收端接收到的信号，构建针对多普勒频偏的最大似然估计模型；根据所述最大似然估计模型，确定所述多普勒频偏的求解表达式；根据所述求解表达式对所述多普勒频偏进行估计。相比于现有技术，应用本发明实施例，可以对毫米波MIMO系统的多普勒频偏进行估计。

Description

基于毫米波MIMO系统的多普勒频偏估计方法及装置
技术领域
本发明涉及无线通信技术领域，特别是涉及基于毫米波MIMO系统的多普勒频偏估计方法及装置。
背景技术
无线通信系统所涉及的信道通常为多径时变衰落信道，其接收信号的幅度和相位会随时间发生变化。衰落信道变化的快慢取决于信道多普勒频偏，多普勒频偏越大，信道变化越快。这就需要实时估计出多普勒频偏，根据多普勒频偏来动态调整系统参数，以获得最优的接收性能。其中，多普勒频偏是指由接收端或移动端的移动造成的频率的变化。多普勒频偏的估计在系统参数的选择、优化和自适应方法都有广泛的应用。
目前，多普勒频偏的估计方法主要有基于信道自相关特性的估计、基于电平通过率的估计、基于开关分集的估计等等。但是每种估计方法都仅应用于各自对应的应用场景，应用范围有限。比如，基于60GHz CS-OFDM MIMO系统的实时频率同步以及相位偏移自动追踪的频偏估计方法，只适用于本地振荡器所针对的场景即无直射径的场景。针对平坦衰落的迭代频偏估计方法，只针对平坦衰落信道的场景。
由于每种估计方法都仅适用于其对应的应用场景，因此上述的估计方法均不能被应用于毫米波MIMO系统，并且现有技术中也没有针对于毫米波MIMO系统的多普勒频偏估计方法，因此，如何对毫米波MIMO系统的多普勒频偏进行估计是亟待解决的问题。
发明内容
本发明实施例的目的在于提供基于毫米波MIMO系统的多普勒频偏估计方法及装置，以对毫米波MIMO系统的多普勒频偏进行估计。
为达到上述目的，本发明实施例还提供了一种基于毫米波MIMO系统的多普勒频偏估计方法，包括：
根据毫米波MIMO系统中的接收端接收到的信号，构建针对多普勒频偏的最大似然估计模型；
根据所述最大似然估计模型，确定所述多普勒频偏的求解表达式；
根据所述求解表达式对所述多普勒频偏进行估计。
可选地，所述根据毫米波MIMO系统中的接收端接收到的信号，构建针对多普勒频偏的最大似然估计模型，包括：
根据信噪比SINR，确定波束成形权重向量；
根据毫米波MIMO系统中的接收端接收到的信号y(n)以及所述波束成形权重向量，构建针对fd、θp、βp的代价函数，其中，fd为移动接收端相对于发射端产生的多普勒频偏，θp为移动接收端相对于发射端的位置角度，βp为移动接收端的衰落损耗；
根据所述代价函数以及所述波束成形权重向量，构建针对多普勒频偏的最大似然估计模型。
可选地，
其中，w(θ)为波束成形权重向量，θ为接收端相对于发射端的位置角度，H为信道矩阵，x(n)为发射端的发射信号的矩阵，aR(θp)为移动接收端的方向向量， 为aT(θp)的转置矩阵，aT(θp)为移动发射端的方向向量，λ为发射信号的波长，nT为发射端的发射天线数，nR为接收端的接收天线数，dT为发射端天线各元素之间的距离，dR为接收端天线各元素之间的距离，Rin为干扰噪声相关矩阵，σ为噪音正态分布函数，为接收天线的单位矩阵，θi为第i个静止接收端相对于发射端的位置角度，βi为第i个静止接收端的衰落损耗，L表示静止接收端个数，为aT(θi)的转置矩阵，aT(θi)为第i个静止发射端的方向向量，aR(θi)为第i个静止接收端的方向向量，
所确定的波束成形权重向量
其中，aR(θ)为接收端的方向向量， 为Rin的逆矩阵；
其中，n表示发射端发射信号的符号，N表示发射端发射信号的总符号数，v(n)为接收端的噪声的矩阵；
构建的针对fd、θp、βp的代价函数
其中，β为接收端的衰落损耗，为aT(θ)的转置矩阵，aT(θ)为发射端的方向向量，
构建的针对多普勒频偏的最大似然估计模型
其中，为aT(θ)的共轭转置矩阵。
可选地，根据所述最大似然估计模型，确定的所述多普勒频偏的求解表达式为：
其中，
为达到上述目的，本发明实施例还提供了一种基于毫米波MIMO系统的多普勒频偏估计装置，包括：
构建模块，用于根据毫米波MIMO系统中的接收端接收到的信号，构建针对多普勒频偏的最大似然估计模型；
确定模块，用于根据所述最大似然估计模型，确定所述多普勒频偏的求解表达式；
估计模块，用于根据所述求解表达式对所述多普勒频偏进行估计。
可选地，所述构建模块包括：
第一确定子单元，用于根据信噪比SINR，确定波束成形权重向量；
第一构建子单元，用于根据毫米波MIMO系统中的接收端接收到的信号y(n)以及所述波束成形权重向量，构建针对fd、θp、βp的代价函数，其中，fd为移动接收端相对于发射端产生的多普勒频偏，θp为移动接收端相对于发射端的位置角度，βp为移动接收端的衰落损耗；
第二构建子单元，用于根据所述代价函数以及所述波束成形权重向量，构建针对多普勒频偏的最大似然估计模型。
可选地，
其中，w(θ)为波束成形权重向量，θ为接收端相对于发射端的位置角度，H为信道矩阵，x(n)为发射端的发射信号的矩阵，aR(θp)为移动接收端的方向向量， 为aT(θp)的转置矩阵，aT(θp)为移动发射端的方向向量，λ为发射信号的波长，nT为发射端的发射天线数，nR为接收端的接收天线数，dT为发射端天线各元素之间的距离，dR为接收端天线各元素之间的距离，Rin为干扰噪声相关矩阵，σ为噪音正态分布函数，为接收天线的单位矩阵，θi为第i个静止接收端相对于发射端的位置角度，βi为第i个静止接收端的衰落损耗，L表示静止接收端个数，为aT(θi)的转置矩阵，aT(θi)为第i个静止发射端的方向向量，aR(θi)为第i个静止接收端的方向向量，
所确定的波束成形权重向量
其中，aR(θ)为接收端的方向向量， 为Rin的逆矩阵；
其中，n表示发射端发射信号的符号，N表示发射端发射信号的总符号数，v(n)为接收端的噪声的矩阵；
构建的针对fd、θp、βp的代价函数
其中，β为接收端的衰落损耗，为aT(θ)的转置矩阵，aT(θ)为发射端的方向向量，
构建的针对多普勒频偏的最大似然估计模型
其中，为aT(θ)的共轭转置矩阵。
可选地，所述确定模块确定出的求解表达式为：
其中，
应用本发明实施例提供的基于毫米波MIMO系统的多普勒频偏估计方法及装置，可以实现对毫米波MIMO系统的多普勒频偏的估计。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为毫米波MIMO系统波束成形结构示意图；
图2为本发明实施例所提供的基于毫米波MIMO系统的多普勒频偏估计方法的流程示意图；
图3为本发明实施例所提供的基于毫米波MIMO系统的多普勒频偏估计方法中构建针对多普勒频偏的最大似然估计模型过程的流程示意图；
图4为本发明实施例所提供的基于毫米波MIMO系统的多普勒频偏估计装置的结构示意图；
图5为本发明实施例所提供的基于毫米波MIMO系统的多普勒频偏估计装置中构建模块的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
为了对毫米波MIMO系统的多普勒频偏进行估计，本发明实施例提供了一种基于毫米波MIMO系统的多普勒频偏估计方法及装置，以下进行详细说明。
需要说明的是，毫米波MIMO系统波束成型结构可以如图1所示，假设该毫米波MIMO系统波所占用频段为28GHZ，发射端则选择了nT根发射天线和nT个发射频增益，接收端则采用了nR根接收天线和nR个接收频增益。
毫米波MIMO系统的具体工作流程如下：发射端的各路发射符号经过基带信号处理之后进行射频链路，进而各路发射符号之间进行发射预编码，编码完成后成为发射信号并被送至天线阵列进行发射，发射信号经过毫米波MIMO信道之后被接收天线阵列接收到，进行与发射端顺序相反的操作之后就得到了接收信号。
在实际应用中，发射端的每个均匀线性阵列可以都包含8个水平元素，接收端的每个均匀线性阵列可以都包含4个水平元素，发射端的天线孔径的波束半宽度可以为大约水平10°，垂直20°，接收端的天线孔径的波束半宽度可以为大约水平25°，垂直60°。对于射频增益，发射端可以为21dBi，接收端的射频增益可以为8dBi。
如图2所示，本发明实施例提供的一种基于毫米波MIMO系统的多普勒频偏估计方法，可以包括以下步骤：
S210，根据毫米波MIMO系统中的接收端接收到的信号，构建针对多普勒频偏的最大似然估计模型；
具体地，为了实现对毫米波MIMO系统中的多普勒频偏的估计，可以根据毫米波MIMO系统中的接收端接收到的信号，构建针对多普勒频偏的最大似然估计模型。
S220，根据所述最大似然估计模型，确定所述多普勒频偏的求解表达式；
具体地，为了降低算法复杂度，可以对所述最大似然估计模型进行化简，进而确定出多普勒频偏的求解表达式。例如，可以将最大似然估计模型简化成二维快速傅里叶变换(Fast Fourier Transform，FFT)求解形式的求解表达式。
S230，根据所述求解表达式对所述多普勒频偏进行估计；
具体地，确定出多普勒频偏的求解表达式后，可以采用二分法、梯度下降法等方法，对求解表达式进行求解，获得多普勒频偏的估计值，进而实现对多普勒频偏的估计。
进一步地，如图3所示，根据毫米波MIMO系统中的接收端接收到的信号，构建针对多普勒频偏的最大似然估计模型的过程，可以包括以下步骤：
S310，根据信噪比SINR，确定波束成形权重向量；
具体地，信干噪比的表达式为：
其中，w(θ)为波束成形权重向量，θ为接收端相对于发射端的位置角度，H为信道矩阵，x(n)为发射端的发射信号的矩阵，aR(θp)为移动接收端的方向向量， 为aT(θp)的转置矩阵，aT(θp)为移动发射端的方向向量，λ为发射信号的波长，nT为发射端的发射天线数，nR为接收端的接收天线数，dT为发射端天线各元素之间的距离，dR为接收端天线各元素之间的距离，Rin为干扰噪声相关矩阵，σ为噪音正态分布函数，为接收天线的单位矩阵，θi为第i个静止接收端相对于发射端的位置角度，βi为第i个静止接收端的衰落损耗，L表示静止接收端个数，为aT(θi)的转置矩阵，aT(θi)为第i个静止发射端的方向向量，aR(θi)为第i个静止接收端的方向向量，
为构建多普勒频偏的最大似然估计模型，需要最大化信噪比，如果要最大化信噪比，则要满足并且波束成形权重向量还满足wHaR(θ)＝1，进而可以推导出波束成形权重向量的表达式。
因此，根据上述条件，所确定的波束成形权重向量其中，aR(θ)为接收端的方向向量， 为Rin的逆矩阵。
S320，根据毫米波MIMO系统中的接收端接收到的信号y(n)以及所述波束成形权重向量，构建针对fd、θp、βp的代价函数，其中，fd为移动接收端相对于发射端产生的多普勒频偏，θp为移动接收端相对于发射端的位置角度，βp为移动接收端的衰落损耗；
具体地，接收端的接收信号的表达式为：
其中，n表示发射端发射信号的符号，N表示发射端发射信号的总符号数，v(n)为接收端的噪声的矩阵；
然后，根据毫米波MIMO系统中的接收端接收到的信号y(n)以及所述波束成形权重向量，构建针对fd、θp、βp的代价函数(fd,θp,βp)；
具体地，最终所确定的针对fd、θp、βp的代价函数(fd,θp,βp)的表达式为：
其中，β为接收端的衰落损耗，为aT(θ)的转置矩阵，aT(θ)为发射端的方向向量，
S330，根据所述代价函数以及所述波束成形权重向量，构建针对多普勒频偏的最大似然估计模型；
具体地，首先，利用β最小化代价函数(fd,θp,βp)得到第一等式；
具体地，所得到的第一等式为：
其中，为aT(θ)的共轭转置矩阵；
然后，将所述第一等式代入到代价函数(fd,θp,βp)中后进行最小化处理，得到第一代价函数(fd,θ)1；
具体地，所得到第一代价函数(fd,θ)1的表达式为：
然后，将所确定的波束成形权重向量w(θ)代入到第一代价函数(fd,θ)1中，得到第二代价函数(fd,θ)2；
具体地，所得到第二代价函数(fd,θ)2的表达式为：
最后，对第二代价函数(fd,θ)2进行化简，得到针对多普勒频偏的最大似然估计模型(fd,θ)；
即，所构建的针对多普勒频偏的最大似然估计模型(fd,θ)为：
进一步地，根据所构建的针对多普勒频偏的最大似然估计模型(fd,θ)，可以确定出多普勒频偏的求解表达式。
具体地，令则(fd,θ)中的数学公式转换为第一表达式G(n)1为：
将aT(θp)和aR(θR)的表达式代入第一表达式G(n)1中，则第一表达式G(n)1可以转换为第二表达式G(n)2为：
其中，rk(n)为的结果中不包括指数部分的部分。
令则第二表达式G(n)2可以转换为第三表达式G(n)3为：
将第三表达式G(n)3展开，可以得到第四表达式G(n)4：
对第四表达式G(n)4进行整理，可以得到G(n)的一般表达式如下：
如果考虑MIMO发射端及接收端两个相邻天线的距离为发射信号的二分之一波长，即以及因此
利用对第四表达式G(n)4进行简化，得到第五表达式G(n)5为：
因此，当时，G(n)的一般表达式为：
对G(n)的一般表达式求其期望值可以得到以下表达式：
令则E{G(n)}可以化简为如下形式：
最后，将化简后的E{G(n)}的表达式代入到针对多普勒频偏的最大似然估计模型(fd,θ)中并进行化简，确定出多普勒频偏的求解表达式
所确定出的多普勒频偏的求解表达式为：
可以理解，对于参数估计问题，克拉美罗界(Cramer Rao Low Bound，CRLB)为任何无偏估计量的方差确定了一个下限。即由于不可能求得方差小于下限的无偏估计量，因此为比较无偏估计量的性能提供了一个标准。而且当无偏估计量达不到CRLB时也可以渐进达到这个下界。
为衡量多普勒频偏fd、移动接收端相对于发射端的位置角度θp以及移动接收端的衰落损耗βp的估计，现推导出这三个参数的克拉美劳下界。
首先，用信号项s表示接收信号y(n)中出现的表达式将v(n)简化成v，因此接收信号y(n)的表达式将会变成y＝s+n。
取时间段为n到n+N-1的接收信号y组成样本向量，可以表示为：y＝[yT(n),yT(n+1),…,yT(n+N-1)]。
将移动接收端的衰落损耗βp看成实部和虚部形式，分别用和表示，则将所有待估计参数用向量的形式表示为：
假设噪声样本v(n)到v(n+N-1)之间是互不相关的，同时忽略干扰带来的影响，参照费希尔信息矩阵(Fisher Information Matrix，FIM)推广的Slepian-Bangs公式，可以得出关于估计向量Ω的费希尔信息矩阵，如下：
其中：
是nR×4的矩阵，是4×nR的矩阵，
其中，□R和□T的表达式为：
应用本实施例，实现了对毫米波MIMO系统的多普勒频偏的估计。
如图4所示，本发明实施例还提供了一种基于毫米波MIMO系统的多普勒频偏估计装置，包括：
构建模块410，用于根据毫米波MIMO系统中的接收端接收到的信号，构建针对多普勒频偏的最大似然估计模型；
确定模块420，用于根据所述最大似然估计模型，确定所述多普勒频偏的求解表达式；
估计模块430，用于根据所述求解表达式对所述多普勒频偏进行估计。
进一步地，如图5所示，所述构建模块410可以包括：
第一确定子单元510，用于根据信噪比SINR，确定波束成形权重向量；
第一构建子单元520，用于根据毫米波MIMO系统中的接收端接收到的信号y(n)以及所述波束成形权重向量，构建针对fd、θp、βp的代价函数，其中，fd为移动接收端相对于发射端产生的多普勒频偏，θp为移动接收端相对于发射端的位置角度，βp为移动接收端的衰落损耗；
第二构建子单元530，用于根据所述代价函数以及所述波束成形权重向量，构建针对多普勒频偏的最大似然估计模型。
进一步地，第一确定子单元510根据的噪比可以为
其中，w(θ)为波束成形权重向量，θ为接收端相对于发射端的位置角度，H为信道矩阵，x(n)为发射端的发射信号的矩阵，aR(θp)为移动接收端的方向向量， 为aT(θp)的转置矩阵，aT(θp)为移动发射端的方向向量，λ为发射信号的波长，nT为发射端的发射天线数，nR为接收端的接收天线数，dT为发射端天线各元素之间的距离，dR为接收端天线各元素之间的距离，Rin为干扰噪声相关矩阵，σ为噪音正态分布函数，为接收天线的单位矩阵，θi为第i个静止接收端相对于发射端的位置角度，βi为第i个静止接收端的衰落损耗，L表示静止接收端个数，为aT(θi)的转置矩阵，aT(θi)为第i个静止发射端的方向向量，aR(θi)为第i个静止接收端的方向向量，
第一确定子单元510所确定的波束成形权重向量
其中，aR(θ)为接收端的方向向量， 为Rin的逆矩阵。
第一构建子单元520根据的毫米波MIMO系统中的接收端接收到的信号y(n)的表达式可以为：
其中，n表示发射端发射信号的符号，N表示发射端发射信号的总符号数，v(n)为接收端的噪声的矩阵；
第一构建子单元520构建的针对fd、θp、βp的代价函数为：
其中，β为接收端的衰落损耗，为aT(θ)的转置矩阵，aT(θ)为发射端的方向向量，
第二构建子单元530构建的针对多普勒频偏的最大似然估计模型为：其中，为aT(θ)的共轭转置矩阵。
进一步地，确定模块420确定出的求解表达式可以为：
其中，
应用本实施例，实现了对毫米波MIMO系统的多普勒频偏的估计。
对于装置实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
本说明书中的各个实施例均采用相关的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
以上所述仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等，均包含在本发明的保护范围内。

Claims (8)

1.一种基于毫米波MIMO系统的多普勒频偏估计方法，其特征在于，包括：
根据毫米波MIMO系统中的接收端接收到的信号，构建针对多普勒频偏的最大似然估计模型；
根据所述最大似然估计模型，确定所述多普勒频偏的求解表达式；
根据所述求解表达式对所述多普勒频偏进行估计。
2.根据权利要求1所述的方法，其特征在于，所述根据毫米波MIMO系统中的接收端接收到的信号，构建针对多普勒频偏的最大似然估计模型，包括：
根据信噪比SINR，确定波束成形权重向量；
根据毫米波MIMO系统中的接收端接收到的信号y(n)以及所述波束成形权重向量，构建针对fd、θp、βp的代价函数，其中，fd为移动接收端相对于发射端产生的多普勒频偏，θp为移动接收端相对于发射端的位置角度，βp为移动接收端的衰落损耗；
根据所述代价函数以及所述波束成形权重向量，构建针对多普勒频偏的最大似然估计模型。
3.根据权利要求2所述的方法，其特征在于，
<mrow> <mi>S</mi> <mi>I</mi> <mi>N</mi> <mi>R</mi> <mo>=</mo> <mfrac> <mrow> <mo>|</mo> <msub> <mi>&amp;beta;</mi> <mi>p</mi> </msub> <msup> <mo>|</mo> <mn>2</mn> </msup> <mi>E</mi> <mo>{</mo> <mo>|</mo> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mi>d</mi> </msub> <mi>n</mi> </mrow> </msup> <msup> <mi>w</mi> <mi>H</mi> </msup> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <msub> <mi>a</mi> <mi>R</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <msubsup> <mi>a</mi> <mi>T</mi> <mi>T</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>}</mo> </mrow> <mrow> <msup> <mi>w</mi> <mi>H</mi> </msup> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msub> <mi>w</mi> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>;</mo> </mrow>
其中，w(θ)为波束成形权重向量，θ为接收端相对于发射端的位置角度，H为信道矩阵，x(n)为发射端的发射信号的矩阵，aR(θp)为移动接收端的方向向量， 为aT(θp)的转置矩阵，aT(θp)为移动发射端的方向向量，λ为发射信号的波长，nT为发射端的发射天线数，nR为接收端的接收天线数，dT为发射端天线各元素之间的距离，dR为接收端天线各元素之间的距离，Rin为干扰噪声相关矩阵，σ为噪音正态分布函数，为接收天线的单位矩阵，θi为第i个静止接收端相对于发射端的位置角度，βi为第i个静止接收端的衰落损耗，L表示静止接收端个数，为aT(θi)的转置矩阵，aT(θi)为第i个静止发射端的方向向量，aR(θi)为第i个静止接收端的方向向量，
所确定的波束成形权重向量
其中，aR(θ)为接收端的方向向量， 为Rin的逆矩阵；
<mrow> <mi>y</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>&amp;beta;</mi> <mi>p</mi> </msub> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mi>d</mi> </msub> <mi>n</mi> </mrow> </msup> <msub> <mi>a</mi> <mi>R</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <msubsup> <mi>a</mi> <mi>T</mi> <mi>T</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>L</mi> </munderover> <msub> <mi>&amp;beta;</mi> <mi>i</mi> </msub> <msub> <mi>a</mi> <mi>R</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msubsup> <mi>a</mi> <mi>T</mi> <mi>T</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>v</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>,</mo> <mi>n</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>N</mi> <mo>;</mo> </mrow>
其中，n表示发射端发射信号的符号，N表示发射端发射信号的总符号数，v(n)为接收端的噪声的矩阵；
构建的针对fd、θp、βp的代价函数
其中，β为接收端的衰落损耗，为aT(θ)的转置矩阵，aT(θ)为发射端的方向向量，
构建的针对多普勒频偏的最大似然估计模型
其中，为aT(θ)的共轭转置矩阵。
4.根据权利要求3所述的方法，其特征在于，根据所述最大似然估计模型，确定的所述多普勒频偏的求解表达式为：
<mrow> <msub> <mover> <mi>f</mi> <mo>~</mo> </mover> <mi>d</mi> </msub> <mo>=</mo> <munder> <mi>argmax</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>c</mi> </msub> <mo>,</mo> <msub> <mi>f</mi> <mi>d</mi> </msub> <mo>)</mo> </mrow> </munder> <mo>|</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <msub> <mi>n</mi> <mi>R</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>+</mo> <mi>&amp;gamma;</mi> <mrow> <mo>(</mo> <msub> <mi>n</mi> <mi>T</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </munderover> <mi>g</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>,</mo> <mi>m</mi> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mi>c</mi> </msub> <mi>m</mi> </mrow> </msup> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mi>d</mi> </msub> <mi>n</mi> </mrow> </msup> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow>
其中，
5.一种基于毫米波MIMO系统的多普勒频偏估计装置，其特征在于，包括：
构建模块，用于根据毫米波MIMO系统中的接收端接收到的信号，构建针对多普勒频偏的最大似然估计模型；
确定模块，用于根据所述最大似然估计模型，确定所述多普勒频偏的求解表达式；
估计模块，用于根据所述求解表达式对所述多普勒频偏进行估计。
6.根据权利要求5所述的装置，其特征在于，所述构建模块包括：
第一确定子单元，用于根据信噪比SINR，确定波束成形权重向量；
第一构建子单元，用于根据毫米波MIMO系统中的接收端接收到的信号y(n)以及所述波束成形权重向量，构建针对fd、θp、βp的代价函数，其中，fd为移动接收端相对于发射端产生的多普勒频偏，θp为移动接收端相对于发射端的位置角度，βp为移动接收端的衰落损耗；
第二构建子单元，用于根据所述代价函数以及所述波束成形权重向量，构建针对多普勒频偏的最大似然估计模型。
7.根据权利要求6所述的装置，其特征在于，
<mrow> <mi>S</mi> <mi>I</mi> <mi>N</mi> <mi>R</mi> <mo>=</mo> <mfrac> <mrow> <mo>|</mo> <msub> <mi>&amp;beta;</mi> <mi>p</mi> </msub> <msup> <mo>|</mo> <mn>2</mn> </msup> <mi>E</mi> <mo>{</mo> <mo>|</mo> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mi>d</mi> </msub> <mi>n</mi> </mrow> </msup> <msup> <mi>w</mi> <mi>H</mi> </msup> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <msub> <mi>a</mi> <mi>R</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <msubsup> <mi>a</mi> <mi>T</mi> <mi>T</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>}</mo> </mrow> <mrow> <msup> <mi>w</mi> <mi>H</mi> </msup> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msub> <mi>w</mi> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>;</mo> </mrow>
其中，w(θ)为波束成形权重向量，θ为接收端相对于发射端的位置角度，H为信道矩阵，x(n)为发射端的发射信号的矩阵，aR(θp)为移动接收端的方向向量， 为aT(θp)的转置矩阵，aT(θp)为移动发射端的方向向量，λ为发射信号的波长，nT为发射端的发射天线数，nR为接收端的接收天线数，dT为发射端天线各元素之间的距离，dR为接收端天线各元素之间的距离，Rin为干扰噪声相关矩阵，σ为噪音正态分布函数，为接收天线的单位矩阵，θi为第i个静止接收端相对于发射端的位置角度，βi为第i个静止接收端的衰落损耗，L表示静止接收端个数，为aT(θi)的转置矩阵，aT(θi)为第i个静止发射端的方向向量，aR(θi)为第i个静止接收端的方向向量，
所确定的波束成形权重向量
其中，aR(θ)为接收端的方向向量， 为Rin的逆矩阵；
<mrow> <mi>y</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>&amp;beta;</mi> <mi>p</mi> </msub> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mi>d</mi> </msub> <mi>n</mi> </mrow> </msup> <msub> <mi>a</mi> <mi>R</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <msubsup> <mi>a</mi> <mi>T</mi> <mi>T</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>L</mi> </munderover> <msub> <mi>&amp;beta;</mi> <mi>i</mi> </msub> <msub> <mi>a</mi> <mi>R</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msubsup> <mi>a</mi> <mi>T</mi> <mi>T</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>v</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>,</mo> <mi>n</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>N</mi> <mo>;</mo> </mrow>
其中，n表示发射端发射信号的符号，N表示发射端发射信号的总符号数，v(n)为接收端的噪声的矩阵；
构建的针对fd、θp、βp的代价函数
其中，β为接收端的衰落损耗，为aT(θ)的转置矩阵，aT(θ)为发射端的方向向量，
构建的针对多普勒频偏的最大似然估计模型
其中，为aT(θ)的共轭转置矩阵。
8.根据权利要求7所述的装置，其特征在于，所述确定模块确定出的求解表达式为：
<mrow> <msub> <mover> <mi>f</mi> <mo>~</mo> </mover> <mi>d</mi> </msub> <mo>=</mo> <munder> <mi>argmax</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>c</mi> </msub> <mo>,</mo> <msub> <mi>f</mi> <mi>d</mi> </msub> <mo>)</mo> </mrow> </munder> <mo>|</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <msub> <mi>n</mi> <mi>R</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>+</mo> <mi>&amp;gamma;</mi> <mrow> <mo>(</mo> <msub> <mi>n</mi> <mi>T</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </munderover> <mi>g</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>,</mo> <mi>m</mi> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mi>c</mi> </msub> <mi>m</mi> </mrow> </msup> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mi>d</mi> </msub> <mi>n</mi> </mrow> </msup> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow>
其中，
CN201610887151.2A 2016-10-11 2016-10-11 基于毫米波mimo系统的多普勒频偏估计方法及装置 Active CN107070818B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201610887151.2A CN107070818B (zh) 2016-10-11 2016-10-11 基于毫米波mimo系统的多普勒频偏估计方法及装置
PCT/CN2016/106880 WO2018068365A1 (zh) 2016-10-11 2016-11-23 一种基于毫米波mimo系统的多普勒频偏估计方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610887151.2A CN107070818B (zh) 2016-10-11 2016-10-11 基于毫米波mimo系统的多普勒频偏估计方法及装置

Publications (2)

Publication Number Publication Date
CN107070818A true CN107070818A (zh) 2017-08-18
CN107070818B CN107070818B (zh) 2019-05-03

Family

ID=59616929

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610887151.2A Active CN107070818B (zh) 2016-10-11 2016-10-11 基于毫米波mimo系统的多普勒频偏估计方法及装置

Country Status (2)

Country Link
CN (1) CN107070818B (zh)
WO (1) WO2018068365A1 (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108828504A (zh) * 2018-03-30 2018-11-16 西安电子科技大学 基于部分相关波形的mimo雷达目标方向快速估计方法
CN109699069A (zh) * 2018-12-20 2019-04-30 桂林电子科技大学 一种基于三维波束形成的联合多普勒频偏和载波频偏的估计方法
CN113472396A (zh) * 2020-03-30 2021-10-01 北京大学 一种基于多普勒效应的上行空分多址系统及实现方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101873152A (zh) * 2010-07-12 2010-10-27 西安电子科技大学 基于最佳带宽滤波预处理的多普勒频移估计方法
US20160020936A1 (en) * 2014-07-17 2016-01-21 Conversant Intellectual Property Management Inc. System and method for frequency synchronization of doppler-shifted subcarriers

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101873152A (zh) * 2010-07-12 2010-10-27 西安电子科技大学 基于最佳带宽滤波预处理的多普勒频移估计方法
US20160020936A1 (en) * 2014-07-17 2016-01-21 Conversant Intellectual Property Management Inc. System and method for frequency synchronization of doppler-shifted subcarriers

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
乔静: "《快时变MIMO-OFDM系统频偏和信道联合估计技术设计》", 30 September 2014 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108828504A (zh) * 2018-03-30 2018-11-16 西安电子科技大学 基于部分相关波形的mimo雷达目标方向快速估计方法
CN109699069A (zh) * 2018-12-20 2019-04-30 桂林电子科技大学 一种基于三维波束形成的联合多普勒频偏和载波频偏的估计方法
CN109699069B (zh) * 2018-12-20 2021-04-02 桂林电子科技大学 一种基于三维波束形成的联合多普勒频偏和载波频偏的估计方法
CN113472396A (zh) * 2020-03-30 2021-10-01 北京大学 一种基于多普勒效应的上行空分多址系统及实现方法

Also Published As

Publication number Publication date
WO2018068365A1 (zh) 2018-04-19
CN107070818B (zh) 2019-05-03

Similar Documents

Publication Publication Date Title
EP2037594B1 (en) Method and apparatus for communicating in a wireless personal area network using adaptive beamforming
US10794984B2 (en) System, method and computer-readable medium for estimating direction of arrival of a signal incident on at least one antenna array
EP3348001B1 (en) Systems and methods for beam selection for hybrid beamforming
US7831232B2 (en) Multiple input multiple output communication apparatus
US20140334564A1 (en) Method and system for providing low-complexity hybrid precoding in wireless communication systems
CN106559367A (zh) 基于低秩张量分解的mimo‑ofdm系统毫米波信道估计方法
US20180337717A1 (en) Communication device and methods thereof
EP3522391B1 (en) Selecting beams based on channel measurements
US10320461B2 (en) System, method and computer-readable medium for estimating direction of arrival of a signal incident on at least one antenna array
CN107070818B (zh) 基于毫米波mimo系统的多普勒频偏估计方法及装置
CN103905105A (zh) 一种双流波束赋形方法和装置
CN111446999A (zh) 基于多臂强盗的位置辅助波束对准方法及其系统
CN104168047A (zh) 利用压缩感知的单端时域波束搜索方法
CN104168046B (zh) 利用压缩感知的单端频域波束搜索方法
KR102494846B1 (ko) 무향칼만필터 기반의 빔 추적 시스템 및 방법
CN104218984B (zh) 利用压缩感知的双端频域波束搜索方法
Nishimoto et al. Performance evaluation of cross-polarized antenna selection over 2 GHz measurement-based channel models
CN111418163B (zh) 无线装置以及无线通信控制方法
Thallapalli et al. Channel Estimation under Dual Wideband Effects for THz Massive MIMO Systems
US20240235897A1 (en) Communication method based on quantum channel estimation and base station
CN114362799B (zh) 波束赋形的方法及毫米波通信系统、电子设备、存储介质
US20230025961A1 (en) Beam selection apparatus and method in wireless communication system
CN110912600B (zh) 一种通信方法、装置、设备及存储介质
US20230388158A1 (en) Ai-augmented channel estimation
Inagaki et al. DOA estimation of desired signals by cyclic ESPRIT based on noise subspace and its performance improvement

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