CN105891771B - 一种提高估计精度的基于连续分布的角度估计方法与设备 - Google Patents
一种提高估计精度的基于连续分布的角度估计方法与设备 Download PDFInfo
- Publication number
- CN105891771B CN105891771B CN201610210465.9A CN201610210465A CN105891771B CN 105891771 B CN105891771 B CN 105891771B CN 201610210465 A CN201610210465 A CN 201610210465A CN 105891771 B CN105891771 B CN 105891771B
- Authority
- CN
- China
- Prior art keywords
- angle
- mrow
- msub
- phi
- mtd
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000005070 sampling Methods 0.000 claims abstract description 108
- 238000013461 design Methods 0.000 claims abstract description 33
- 239000011159 matrix material Substances 0.000 claims description 60
- 239000013598 vector Substances 0.000 claims description 58
- 230000005540 biological transmission Effects 0.000 claims description 20
- 238000012549 training Methods 0.000 claims description 18
- 238000007781 pre-processing Methods 0.000 claims description 16
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 239000000654 additive Substances 0.000 claims description 4
- 230000000996 additive effect Effects 0.000 claims description 4
- 238000005562 fading Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 2
- 241000022852 Letis Species 0.000 claims 1
- 238000004891 communication Methods 0.000 abstract description 14
- 238000005516 engineering process Methods 0.000 abstract description 10
- 238000012545 processing Methods 0.000 abstract description 7
- 238000010586 diagram Methods 0.000 description 8
- 238000004088 simulation Methods 0.000 description 7
- 230000001427 coherent effect Effects 0.000 description 3
- 230000006835 compression Effects 0.000 description 3
- 238000007906 compression Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000001413 cellular effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000012772 sequence design Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/04—Details
- G01S3/12—Means for determining sense of direction, e.g. by combining signals from directional antenna or goniometer search coil with those from non-directional antenna
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radio Transmission System (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明提供了一种提高估计精度的基于连续分布的角度估计方法与设备，属于宽带无线通信技术和多天线技术领域。所述设备包括发端码本设计模块、收端码本设计模块和信道估计模块，信道估计模块包括预处理模块和角度迭代估计器；预处理模块在角度域将连续分布的角度信息分解为整数倍采样角度与小数倍采样角度，角度迭代估计器借鉴Turbo译码原理实现角度估计。角度估计方法包括发送端码本设计、接收端码本设计和信道估计。本发明在天线数目较少的条件下即可实现高精度角度估计，极大地降低了装置部署天线的成本，降低了角度估计的算法复杂度，减少了接收端的处理时间，适用于大规模天线阵列下的单径毫米波通信系统，具有很好的推广应用前景。
Description
技术领域
本发明涉及一种提高角度估计精度的检测方法，尤其是涉及毫米波天线阵列出发角/到达角估计的方法，属于宽带无线通信技术和多天线技术领域。
背景技术
为了满足下一代移动通信系统(5G)对高速数据传输率和大容量的需求，需要寻找更大传输带宽的频谱资源和研究高频谱效率的传输技术。毫米波频段的可用传输带宽可超过1GHz，是提供更大传输带宽的潜在频谱资源，因而毫米波通信技术有着广阔的发展前景。毫米波段的通信频段高，信号路径损耗十分严重，相比于微波通信，毫米波通信的信号衰减高达20-40dB左右，链路损耗问题十分突出。由于毫米波的波长较短，毫米波通信系统可以在收发两端同时部署大规模天线阵列，利用波束赋形、相干接收技术获得高的阵列增益来补偿严重的路径损耗。但波束赋形与相干接收都需要准确的信道状态信息，尤其是发送角和到达角信息。同时，由于毫米波频段的射频链路成本高昂，为节省成本，采用大规模天线阵列的通信系统中使用的射频链路数通常远小于天线阵元数目。在射频链路数目远少于天线阵元的硬件约束条件下，传统的空间超分辨率角度估计方法(例如MUSIC/ESPRINT等算法)不能使用。
其次，传统的角度估计方法在角度空间离散化时，通常假定角度分布在离散的采样点上，而实际角度分布是连续的。因此，为了满足角度分布在离散的采样点上的假设，在给定载波频率及天线阵元间隔条件下，只能通过增加阵元数来实现。然而，增加天线阵元数不仅增加了硬件成本，同时增大风阻，不利用实际应用。当角度分布不在离散的采样点上时，传统的角度估计算法性能显著下降。
对于高频段的毫米波蜂窝系统，为了克服高频段信号衰减严重的缺点，可以通过在收发两端配置大规模天线阵列并采用自适应波束赋形与相干接收技术获得大的阵列增益来提升链路的传输性能，而超分辨的角度估计技术是核心。现有技术主要有：
[1]中国专利申请：一种基于空时数据的高分辨目标方位的估计方法，公开号为104392114A，公开日期为2015年3月4日。该方案通过将天线输出的数据在时间域、空间域和延迟域形成相关矩阵，对相关矩阵进行去噪处理，利用循环特征分解方法计算信号子空间，从信号和噪声子空间的关系求解噪声子空间，得到噪声子空间的投影矩阵，基于投影矩阵反解目标方位和俯仰信息。
[2]中国专利申请：高精度阵列天线接收系统角度估计的装置及其方法，公开号为102394686A，公开日期为2012年3月28日。该方案通过选择梯度下降法或牛顿法，对各个阵元设置自适应权值，采用MUSIC算法来估计接收信号的到达角度，将与该角度对应的权值参数作为自适应迭代的初值。计算系统的输出信号功率对权值参数进行自适应迭代；搜索最优的权值向量获得相应的角度值。
[3]中国专利申请：一种可扩展的用于均匀圆阵二维到达角的快速估计算法，公开号为104931923A，公开日期为2015年9月23日。该方案提出一种基于均匀圆阵的基于网格的迭代旋转不变技术估计型号参数的算法来获得二维到达角的超分辨估计。采用空间网格划分，利用循环补偿和迭代使用经典ESPRIT算法来进行估计。
[4]自适应压缩感知(ACS)算法的角度估计技术，基于二分法逐次减半角度空间，使用空间匹配滤波器搜索多径所在的角度范围，即收发两端多次交互发送训练序列，通过多次迭代，实现角度估计。
现有的四种天线阵列角度方法在毫米波信道角度估计应用中都存在局限性。[1][2][3]中所设计的方法虽然属于超分辨的角度估计方法，但都是基于MUSIC/ESPRIT类算法，而MUSIC/ESPRIT算法需要在数字域使用各天线的直接输出。当存在射频链路约束时，各天线的输出经过接收合并后才能在数字域被信道估计模块进行处理，因而MUSIC/ESPRIT类算法无法使用。[4]中给出的基于自适应压缩感知算法的角度估计技术的估计精度受限于物理天线数，且完成估计所需迭代次数较多，相应的估计时延和训练序列的开销较大。因此，如何在射频链路数有限的情况下，设计超分辨率的角度估计方法，在不增加天线阵元数目基础上，实现任意角度分辨率和通用低复杂的精确角度估计，是目前学术界和产业界都非常关注的热点。
发明内容
本发明为了克服传统角度估计方法对于角度仅仅分布在离散的采样点上的假设的限制，从角度连续分布的实际出发，提供了一种提高估计精度的基于连续分布的角度估计方法及设备，针对大规模均匀线性天线阵列下的毫米波通信系统，在不增加天线阵元数目的基础上，实现了高精度的角度估计。
本发明的基于连续分布的角度估计设备，包括发端码本设计模块、收端码本设计模块和信道估计模块。信道估计模块包括预处理模块和角度迭代估计器。
发端码本设计模块在连续的R个时隙采用相同的波束赋形矩阵发送相同的训练序列，其中R为正整数，训练序列x为全1向量，波束赋形矩阵WB中某一个列向量为第i0行元素为1、其他行元素为0的单位列向量，WB的其他列向量均为0向量。
收端码本设计模块设计在连续的R个时隙采用不同的合并矩阵接收数据，其中，合并矩阵由随机单位矢量构成，任意两个中的单位矢量均不同，r＝1,2...R。
预处理模块对接收信号进行预处理，在角度域将连续分布的角度信息分解为整数倍采样角度与小数倍采样角度。
角度迭代估计器通过在整数倍采样角度与小数倍采样角度两部分之间相互迭代，进行角度估计。
本发明的基于连续分布的角度估计方法，实现步骤如下：
步骤1，发送端码本设计，具体为：发送端在连续的R个时隙采用相同的波束赋形矩阵发送相同的训练序列；训练序列x为全1向量，R为正整数；波束赋形矩阵WB中某一个列向量为第i0行元素为1、其他行元素为0的单位列向量，WB的其他列向量均为0向量。
步骤2，接收端码本设计，具体为：接收端在连续的R个时隙采用不同的合并矩阵接收数据，由随机单位矢量构成，任意两个中的单位矢量均不同，r＝1,2...R。
步骤3，进行信道估计，具体为：(3.1)接收端对接收信号进行预处理，在角度域将连续分布的角度信息分解为整数倍采样角度与小数倍采样角度；(3.2)通过在整数倍采样角度与小数倍采样角度两部分之间相互迭代，进行角度估计。
所述的(3.1)中，设天线到达角和发送角分别表示为θM和θB，然后转换得到角度转换值φM和φB为：
对区间进行离散化，分为N等分，N为天线数目；
则将角度转换值φ表示为其中，φk为整数倍采样角度，是与φ相邻最近的第k个离散角度采样点，k＝1,2,...,N；Δ为小数倍采样角度，是φ与离散点φk之间的偏差；当φ取φB时，N为基站天线数目NB，当φ取φM时，N为移动站天线数目NM。
本发明的优点与积极效果在于：
(1)克服传统角度估计方法对于角度仅仅分布在离散的采样点上的假设的限制，通过在接收端进行预处理，在角度域将连续分布的角度信息分解为与其最邻近的离散角度采样点估计和其与离散点的偏差估计两部分，即整数倍采样角度信息和小数倍采样角度信息。这种处理方式可以不受天线阵列尺寸的限制，在天线数目较少的条件下，即可实现高精度角度估计，可以超出整数倍均匀角度估计的分辨率下限，极大地降低了装置部署天线的成本。
(2)采用传统的收发端训练方式，包括训练序列、发端波束赋形矩阵、収端合并矩阵的设计，即可以使接收端提取的感知矩阵能满足约束等距性条件；
(3)采用常规的压缩重构算法即可完成角度估计，极大地降低了角度估计的算法复杂度，减少了接收端的处理时间；
(4)借鉴Turbo译码原理，通过整数倍采样角度与小数倍采样角度两部分之间相互迭代，来实现可靠的连续的角度估计。
本发明适用于大规模天线阵列下的单径毫米波通信系统，克服了现有角度估计方法资源开销大，角度分辨率限制的问题，故具有很好的推广应用前景。
附图说明
图1是本发明的基于连续分布的角度估计设备的实现框图；
图2是本发明的角度估计方法中接收端对接受信号预处理的示意图；
图3是本发明的角度迭代估计器的实现框图；
图4是本发明实施例中不同阵元数下，角度估计误差和均匀整数倍角度估计误差下限仿真图；
图5是本发明实施例中不同阵元数与其相邻最近的整数倍离散点估计错误概率仿真对比图。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明提供了一种提高估计精度的基于连续分布的角度估计方法及其设备。针对大规模均匀线性天线阵列下的毫米波通信系统，充分挖掘毫米波信道的稀疏结构特性，基于压缩感知技术，设计了训练序列辅助的角度估计方法。通过在接收端对接收信号进行预处理，在角度域将连续分布的角度信息分解为与其最接近的离散角度采样点估计和其与离散角度采样点的偏差估计两部分，即整数倍采样角度信息与小数倍采样角度信息。收发两端设计恰当的训练序列，使接收端接收的训练信号模型满足约束等距性条件，采用基于随机测量的压缩重构方法能够估计出整数倍采样角度；当整数倍采样角度已知条件下，通过最小二乘(LS：least square)算法即可估计出小数倍采样角度。借鉴Turbo译码原理，设计了在整数倍采样角度与小数倍采样角度两部分之间相互迭代，来实现连续的角度估计。因此，本发明在不增加天线数目的基础上，便可以实现高精度的角度估计。
毫米波信道中，直射径能量远大于非直射径的能量和，基于最大化信道容量准则的波束赋形如只给直射径分配功率而忽略非直射径，信道吞吐量仅有轻微的损失，因此对于配置了大规模天线阵列的毫米波通信系统，可以仅仅考虑直射径的角度估计而忽略其余非直射径。
设毫米波通信系统的基站(BS：Base Station)配置NB个天线和KB个射频链路，移动站(MS：Moblie Station)有NM个天线阵元和KM个射频链路。当BS与MS之间的距离远大于天线阵列尺寸时，天线阵列上各阵元的接收信号衰落幅值近似一致，仅存在相位差。信道矩阵可表示为其中，θM和θB分别为天线到达角(AOA)和发送角(AOD)，β是衰落系数，aM(·)为移动站天线阵列的导向矢量，aB(·)为基站天线阵列的导向矢量，记将天线到达角和发送角进行数值转换，得到的数值φM和φB分别为：
则：
如图1所示，本发明的基于连续分布的角度估计设备包含三大模块：发端码本设计模块，收端码本设计模块，信道估计模块。信道估计模块包括预处理模块和角度迭代估计器。发端码本设计模块设计在连续的时隙内采用相同的波束赋形矩阵发送相同的训练序列。收端码本设计模块设计在连续的时隙内采用不同的合并矩阵接收数据。预处理模块对接收信号进行预处理，在角度域将连续分布的角度信息分解为整数倍采样角度与小数倍采样角度。角度迭代估计器借鉴Turbo译码原理，通过在整数倍采样角度与小数倍采样角度两部分之间相互迭代，进行角度估计。
对应地，本发明的基于连续分布的角度估计方法涉及三个部分：发送端码本设计；接收端码本设计；信道估计。发送端的码本设计包括训练序列设计、发端赋形矩阵设计；接收端的码本设计包含収端合并矩阵设计。设计原则是使角度估计的信号模型满足基于随机测量的压缩重构算法可重构条件。下面本发明举例说明三个模块/部分的实现。
第一步，发送端码本设计，即发端码本设计模块的实现。
(1.1)发送端在R个时隙内发送相同的全一向量R为正整数。
(1.2)发送端设计波束赋形矩阵WB，
设WB的第n0个列向量为 为第i0个元素为1、其他元素为0的NB维的单位列向量，n0∈{1,2,...,KB}，i0取值集合为集合PB从集合RB中随机选取；WB的其余列向量均为0向量。
(1.3)发送端在连续的R个时隙采用相同的波束赋形矩阵WB发送相同的训练序列x。
第二步，接收端码本设计，即收端码本设计模块的实现。
(2.1)接收端在连续的R个时隙采用不同的合并矩阵接收数据，合并矩阵由KM个随机单位矢量构成，即的第m个列向量其中，m＝1,2,...,KM，为第个元素为1，其他元素为零的NM维单位列向量，取值集合为任意两个合并矩阵中的单位矢量均不相同，从Rr中随机选取，其中R1＝{1,2,...,NM}，R2表示从R1选了KM个数以后剩下的集合，Rr表示从R1选了(r-1)KM个数以后剩下的集合。
(2.2)发送信号通过BS端波束赋形处理和MS端接收合并处理后接收信号为：
其中，为接收信号，为加性高斯白噪声，为R个时隙收到的合并在一块的总合并矩阵
第三步，进行信道估计，即信道估计模块的实现。信道估计分为两步：接收端的预处理，角度迭代估计。
(3.1)接收端的预处理，即信道估计模块中预处理模块的实现。
接收端对接收信号进行预处理，在角度域将连续分布的角度信息分解为整数倍采样角度信息与小数倍采样角度信息。具体实现如下。
(3.1.1)对连续角度φ所属的区间进行离散化，分为N等分，N为天线数目，φ取φB或φM，对应N取值为NB或NM。将φ分为与其相邻最近的第k(k＝1,2,...,N)个离散角度采样点φk和其与离散点之间的偏差Δ两部分，如图2所示，即整数倍采样角度与小数倍采样角度：
(3.1.2)信道矩阵可表示为：
其中，为φM相邻最近的第k(k＝1,2,...,NM)个整数倍采样角度，为φB相邻最近的第q(q＝1,2,...,NB)个整数倍采样角度。与φM的偏差记做ΔM，与φB的偏差记做ΔB，和分别表示包含到达角或发送角小数倍采样角度的一个对角矩阵，称为偏差矩阵，具体形式为：
本发明实施例在说明时，偏差矩阵与偏差的字符表示不同之处在于，前者用粗体字符表示，后者用非粗体字符，例如偏差矩阵ΔM与偏差ΔM。
(3.1.3)根据均匀整数倍角度估计的信道稀疏处理方式，选取离散傅里叶变换为一组正交基，将天线阵列间的毫米波连续信道矩阵H中离散部分进行稀疏表征，相应形式为：
式中，和为离散傅里叶变换矩阵，为整数倍采样角度的虚拟信道矩阵。
(3.1.4)经过预处理后接收导频序列的形式为：
式中，为仅包含到达角度小数倍采样角度的对角矩阵，是稀疏度为1的仅包含到达角的整数倍采样角度的虚拟信道向量。
(3.1.5)MS端作为接收端重复(3.1.1)～(3.1.4)过程，估计MS端的AoA信息，即BS端的AoD信息。
(3.2)接收端的角度迭代估计，即接收端的角度迭代估计器的实现。
借鉴turbo译码原理，设计了角度估计迭代器，通过在整数倍采样角度与小数倍采样角度两部分之间相互迭代，来实现连续分布的角度估计，迭代过程如图3所示，具体描述如下。
(3.2.1)初始化估计器的参数，包括到达角/发送角的小数倍采样角度，初始偏差Δ＝0；
(3.2.2)整数倍采样角度的估计：在小数倍采样角度已知的条件下，由接收信号公式通过压缩
y＝WΔFg+z (6)
其中，W是压缩感知重构算法的采样矩阵，对应公式(5)中的Δ·F是压缩感知重构算法的变换矩阵，对应公式(5)中的由偏差Δ可得到偏差矩阵。
采取传统的压缩感知重构算法，例如正交匹配追踪算法、Turbo估计算法等，即可估计整数倍采样角度的虚拟信道稀疏向量g。
(3.2.3)小数倍采样角度的估计：在整数倍采样角度的虚拟信道稀疏向量g已知的条件下，由对角矩阵的点乘特性，由LS算法即可估计出小数倍采样角度矩阵Δ。
(3.2.4)检查迭代终止条件，若不满足，将Δ作为初始值继续第二部迭代过程。若满足进入(3.2.5)。
(3.2.5)终止迭代，输出小数倍采样角度矩阵的估计值和整数倍采样角度的虚拟信道稀疏向量
(3.2.6)根据整数倍采样角度的虚拟信道向量非零元素所在的行号k，得到天线接收端到达角所对应的与其相邻最近的整数倍采样点：
(3.2.7)根据偏差矩阵通过对角元素的线性拟合，可以得到天线的发送角和到达角所对应的小数倍采样角度的估计值
(3.2.8)根据所求得的与其相邻最近的整数倍采样点和其与整数倍采样点的角度偏差，得到估计的连续角度如下：
利用上述等式，接收端估计出天线发送角和到达角
下面结合图4和图5，说明本发明方法进行多次仿真实施试验的结果。仿真实验结果图是在加性白高斯噪声信道状况下，随机生成100000次的链路级仿真实施试验图，为了简便，其中收发端天线数目，且射频链路数设定相等，即NB＝NM＝N，R＝5。
图4中曲线为本发明方法实施例中在不同阵元数下角度估计误差曲线，横坐标为发射信噪比(SNR)，纵坐标为平均角度估计误差。仿真实验表明，在发射信噪比从0dB增加到20dB的过程中，本发明实施例中角度估计方法的角度估计误差逐渐降低。当N＝512时，本发明的角度估计方法在SNR＝3dB左右开始超过基于均匀整数倍角度估计方法的角度估计误差下限值，此时基于均匀整数倍角度估计方法的角度估计误差下限值如图中lower boundof 512 base in discrete angle所示的线段。当N＝1024时，本发明的角度估计方法在SNR＝7dB左右开始超过基于均匀整数倍角度估计方法的角度估计误差下限值，基于均匀整数倍角度估计方法的角度估计误差下限值如图中lower bound of 1024 base in discreteangle所示的线段。当SNR＝20dB时，天数数目为512，本发明的与误差下限相差0.003406；天线数目1024，本发明的与误差下限相差0.000728。随着阵元数N增加，角度分辨率增大，角度估计误差也越小。
图5是本发明的角度估计方法中与其相邻最近的均匀整数倍离散点估计错误概率仿真图，横坐标为发射信噪比(SNR)，纵坐标为均匀整数倍离散点估计错误概率。在发射信噪比从0dB增加到20dB的过程中，随着阵元数的增加，对应于压缩感知中的稀疏率越低，其重构性能也越好，与其相邻最近的均匀整数倍离散点估计错误概率越低。因此，通过实验证明本发明能在不增加天线阵元数目的基础上，实现更高精度的角度估计的目的。
Claims (6)
1.一种提高估计精度的基于连续分布的角度估计方法，其特征在于，实现步骤如下：
步骤1，发送端码本设计，具体为：发送端在连续的R个时隙采用相同的波束赋形矩阵发送相同的训练序列；
训练序列x为全1向量，R为正整数；波束赋形矩阵WB中某一个列向量为第i0行元素为1、其他行元素为0的单位列向量，WB的其他列向量均为0向量；
步骤2，接收端码本设计，具体为：接收端在连续的R个时隙采用不同的合并矩阵接收数据，由随机单位矢量构成，任意两个中的单位矢量均不同，r＝1,2...R；
步骤3，进行信道估计，具体为：(3.1)接收端对接收信号进行预处理，在角度域将连续分布的角度信息分解为整数倍采样角度与小数倍采样角度；(3.2)通过在整数倍采样角度与小数倍采样角度两部分之间相互迭代进行角度估计；
所述的(3.1)中，设天线到达角和发送角分别表示为θM和θB，然后转换得到角度转换值φM和φB为：
<mrow>
<msub>
<mi>&phi;</mi>
<mi>M</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>sin&theta;</mi>
<mi>M</mi>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>,</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>&phi;</mi>
<mi>B</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>sin&theta;</mi>
<mi>B</mi>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>,</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>)</mo>
<mo>;</mo>
</mrow>
对区间进行离散化，分为N等分，N为天线数目；
则将角度转换值φ表示为其中，φk为整数倍采样角度，是与φ相邻最近的第k个离散角度采样点，k＝1,2,...,N；Δ为小数倍采样角度，是φ与离散点φk之间的偏差；当φ取φB时，N为基站天线数目NB，当φ取φM时，N为移动站天线数目NM。
2.根据权利要求1所述的一种提高估计精度的基于连续分布的角度估计方法，其特征在于，所述的(3.1)中，设为φM相邻最近的第k(k＝1,2,...,NM)个整数倍采样角度，为φB相邻最近的第q(q＝1,2,...,NB)个整数倍采样角度，与φM的偏差记做ΔM，与φB的偏差记做ΔB，记包含到达角或发送角小数倍采样角度的一个对角矩阵为偏差矩阵，表示为ΔM和ΔB如下：
<mrow>
<msub>
<mi>&Delta;</mi>
<mi>M</mi>
</msub>
<mo>=</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;&Delta;</mi>
<mi>M</mi>
</msub>
</mrow>
</msup>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mi>B</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>&Delta;</mi>
<mi>M</mi>
</msub>
</mrow>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>&Delta;</mi>
<mi>B</mi>
</msub>
<mo>=</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;&Delta;</mi>
<mi>B</mi>
</msub>
</mrow>
</msup>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mi>B</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>&Delta;</mi>
<mi>B</mi>
</msub>
</mrow>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
信道矩阵其中，β是衰落系数，aM(·)为移动站天线阵列的导向矢量，aB(·)为基站天线阵列的导向矢量；
选取离散傅里叶变换为一组正交基，将信道矩阵H中的离散部分进行稀疏表征，则H表示为：其中，FB和FM为离散傅里叶变换矩阵，G为整数倍采样角度的虚拟信道矩阵；
经过预处理后接收信号y表示为：其中，g是稀疏度为1的仅包含到达角的整数倍采样角度的虚拟信道向量；为R个时隙收到的合并得到矩阵，z为加性高斯白噪声。
3.根据权利要求1或2所述的一种提高估计精度的基于连续分布的角度估计方法，其特征在于，所述的(3.2)中，实现过程包括：
在小数倍采样角度已知的条件下对整数倍采样角度的估计时，采取压缩感知重构算法估计整数倍采样角度的虚拟信道稀疏向量g；
在整数倍采样角度的虚拟信道稀疏向量g已知的条件下，对小数倍采样角度的估计时，采用最小二乘算法估计出小数倍采样角度矩阵Δ；
当满足迭代终止条件时，输出小数倍采样角度矩阵的估计值和整数倍采样角度的虚拟信道稀疏向量
根据非零元素所在的行号k，得到天线接收端到达角所对应的与其相邻最近的整数倍采样点φk，
根据通过对角元素的线性拟合，得到天线的发送角和到达角所对应的小数倍采样角度的估计值
最后根据下式得到估计的连续角度
<mrow>
<mfrac>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mover>
<mi>&theta;</mi>
<mo>~</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
<mn>2</mn>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</mfrac>
<mo>+</mo>
<mover>
<mi>&Delta;</mi>
<mo>~</mo>
</mover>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>.</mo>
</mrow>
4.一种提高估计精度的基于连续分布的角度估计设备，其特征在于，包括发端码本设计模块、收端码本设计模块和信道估计模块；信道估计模块包括预处理模块和角度迭代估计器；
发端码本设计模块在连续的R个时隙采用相同的波束赋形矩阵发送相同的训练序列，其中R为正整数，训练序列x为全1向量，波束赋形矩阵WB中某一个列向量为第i0行元素为1、其他行元素为0的单位列向量，WB的其他列向量均为0向量；
收端码本设计模块设计在连续的R个时隙采用不同的合并矩阵接收数据，其中，合并矩阵由随机单位矢量构成，任意两个中的单位矢量均不同，r＝1,2...R；
预处理模块对接收信号进行预处理，在角度域将连续分布的角度信息分解为整数倍采样角度与小数倍采样角度；
角度迭代估计器通过在整数倍采样角度与小数倍采样角度两部分之间相互迭代，进行角度估计。
5.根据权利要求4所述的一种提高估计精度的基于连续分布的角度估计设备，其特征在于，所述的预处理模块将角度转换值φ分为整数倍采样角度φk和小数倍采样角度Δ，具体是：对区间进行离散化，分为N等分，N为天线数目；设φk是与φ相邻最近的第k个离散角度采样点，k＝1,2,...,N；Δ是φ与离散点φk之间的偏差；
角度转换值φ表示为：
设天线到达角和发送角分别表示为θM和θB，对应的角度转换值φM和φB分别为：
<mrow>
<msub>
<mi>&phi;</mi>
<mi>M</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>sin&theta;</mi>
<mi>M</mi>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>,</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>&phi;</mi>
<mi>B</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>sin&theta;</mi>
<mi>B</mi>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>,</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>)</mo>
<mo>;</mo>
</mrow>
当φ取φB时，N为基站天线数目NB，当φ取φM时，N为移动站天线数目NM；
设为φM相邻最近的第k(k＝1,2,...,NM)个整数倍采样角度，为φB相邻最近的第q(q＝1,2,...,NB)个整数倍采样角度，与φM的偏差记做ΔM，与φB的偏差记做ΔB，记包含到达角或发送角小数倍采样角度的一个对角矩阵为偏差矩阵，表示为ΔM和ΔB如下：
<mrow>
<msub>
<mi>&Delta;</mi>
<mi>M</mi>
</msub>
<mo>=</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;&Delta;</mi>
<mi>M</mi>
</msub>
</mrow>
</msup>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mi>B</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>&Delta;</mi>
<mi>M</mi>
</msub>
</mrow>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>&Delta;</mi>
<mi>B</mi>
</msub>
<mo>=</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;&Delta;</mi>
<mi>B</mi>
</msub>
</mrow>
</msup>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mi>B</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>&Delta;</mi>
<mi>B</mi>
</msub>
</mrow>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
信道矩阵其中，β是衰落系数，aM(·)为移动站天线阵列的导向矢量，aB(·)为基站天线阵列的导向矢量；
选取离散傅里叶变换为一组正交基，将信道矩阵H中的离散部分进行稀疏表征，则H表示为：其中，FB和FM为离散傅里叶变换矩阵，G为整数倍采样角度的虚拟信道矩阵；
经过预处理后接收信号y表示为：其中，g是稀疏度为1的仅包含到达角的整数倍采样角度的虚拟信道向量；为R个时隙收到的合并得到矩阵，z为加性高斯白噪声。
6.根据权利要求4或5所述的一种提高估计精度的基于连续分布的角度估计设备，其特征在于，所述的角度迭代估计器，在小数倍采样角度已知的条件下对整数倍采样角度的估计时，采取压缩感知重构算法估计整数倍采样角度的虚拟信道稀疏向量g；在整数倍采样角度的虚拟信道稀疏向量g已知的条件下，对小数倍采样角度的估计时，采用最小二乘算法估计出小数倍采样角度矩阵Δ；
设满足迭代终止条件时，输出的小数倍采样角度矩阵的估计值为输出的整数倍采样角度的虚拟信道稀疏向量为根据非零元素所在的行号k，得到天线接收端到达角所对应的与其相邻最近的整数倍采样点φk；根据通过对角元素的线性拟合，得到估计值最后估计出连续角度为：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610210465.9A CN105891771B (zh) | 2016-04-06 | 2016-04-06 | 一种提高估计精度的基于连续分布的角度估计方法与设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610210465.9A CN105891771B (zh) | 2016-04-06 | 2016-04-06 | 一种提高估计精度的基于连续分布的角度估计方法与设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105891771A CN105891771A (zh) | 2016-08-24 |
| CN105891771B true CN105891771B (zh) | 2018-01-02 |
Family
ID=57012994
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610210465.9A Active CN105891771B (zh) | 2016-04-06 | 2016-04-06 | 一种提高估计精度的基于连续分布的角度估计方法与设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105891771B (zh) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10178697B2 (en) * | 2017-02-22 | 2019-01-08 | Mitsubishi Electric Research Laboratories | Systems and methods for separating collided packets in shared wireless channels |
| CN107426119B (zh) * | 2017-05-31 | 2020-02-07 | 东南大学 | 毫米波通信的十字形信道估计方法、装置及系统 |
| CN108449098B (zh) * | 2018-03-22 | 2020-08-21 | 京东方科技集团股份有限公司 | 信号接收方法及装置 |
| CN108494449B (zh) * | 2018-03-23 | 2021-03-30 | 中国人民解放军陆军工程大学 | 大规模mimo高速移动场景下角度域信道追踪方法 |
| CN108512787B (zh) * | 2018-04-13 | 2021-09-24 | 电子科技大学 | 大规模mimo系统的超精细信道估计方法 |
| CN108933745B (zh) * | 2018-07-16 | 2020-07-10 | 北京理工大学 | 一种基于超分辨率角度和时延估计的宽带信道估计方法 |
| CN109738854B (zh) * | 2018-12-14 | 2020-07-10 | 北京邮电大学 | 一种天线阵列来波方向的到达角估计方法 |
| CN110187304B (zh) * | 2019-05-21 | 2021-05-04 | 泰凌微电子(上海)股份有限公司 | 一种信号到达角估计方法及装置 |
| TWI711835B (zh) * | 2019-08-12 | 2020-12-01 | 群邁通訊股份有限公司 | Aod估算方法及裝置 |
| CN112422458B (zh) * | 2019-08-23 | 2022-03-25 | 中兴通讯股份有限公司 | 信道估计方法、装置和计算机存储介质 |
| CN111541015B (zh) * | 2020-04-07 | 2023-01-24 | 南京市德赛西威汽车电子有限公司 | 一种改善天线角度分辨率的方法及天线 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1535048A (zh) * | 2003-03-31 | 2004-10-06 | 大唐移动通信设备有限公司 | 一种估计固定波束空间到达方向的方法 |
| CN1917396A (zh) * | 2005-08-18 | 2007-02-21 | 中兴通讯股份有限公司 | 一种波达角估计和波束赋形的方法 |
| CN104698430A (zh) * | 2015-03-27 | 2015-06-10 | 北京邮电大学 | 一种用于提高精度的基于虚拟天线阵列的角度估计方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100749451B1 (ko) * | 2005-12-02 | 2007-08-14 | 한국전자통신연구원 | Ｏｆｄｍ 기지국 시스템에서의 스마트 안테나 빔 형성 방법및 장치 |
-
2016
- 2016-04-06 CN CN201610210465.9A patent/CN105891771B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1535048A (zh) * | 2003-03-31 | 2004-10-06 | 大唐移动通信设备有限公司 | 一种估计固定波束空间到达方向的方法 |
| CN1917396A (zh) * | 2005-08-18 | 2007-02-21 | 中兴通讯股份有限公司 | 一种波达角估计和波束赋形的方法 |
| CN104698430A (zh) * | 2015-03-27 | 2015-06-10 | 北京邮电大学 | 一种用于提高精度的基于虚拟天线阵列的角度估计方法 |
Non-Patent Citations (4)
| Title |
|---|
| 基于扩展码滤波法的联合角度一时延估计;王维新 等;《信号处理》;20051231;第21卷(第6期);第668-671，652页 * |
| 多经信道环境下LTE及LTE-A系统中的多普勒频偏估计;刘航;《中国优秀硕士学位论文全文数据库 信息科技辑》;20150815(第8期);第I136-405页 * |
| 大规模MIMO系统的信道与DOA估计研究;胡安中;《中国博士学位论文全文数据库 信息科技辑》;20150415(第4期);第I136-14页 * |
| 毫米波通信中波束赋形性能仿真平台的设计与实现;朱晨;《中国优秀硕士学位论文全文数据库 信息科技辑》;20150815(第8期);第I136-616页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105891771A (zh) | 2016-08-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105891771B (zh) | 一种提高估计精度的基于连续分布的角度估计方法与设备 | |
| CN107315160B (zh) | 基于内插虚拟阵列信号原子范数最小化的互质阵列波达方向估计方法 | |
| CN112910578B (zh) | 一种针对毫米波3d mimo信道的路径参数提取方法 | |
| Huang et al. | Frequency-domain AoA estimation and beamforming with wideband hybrid arrays | |
| JP2020509683A (ja) | シンボルを復号する方法及びシンボルを受信し、復号する受信機 | |
| CN102175989B (zh) | 一种非相干分布式信号二维波达角的测定方法 | |
| CN110022274B (zh) | 一种毫米波mimo-ofdm系统的联合信道与载波频率偏移估计方法 | |
| CN109738854A (zh) | 一种天线阵列来波方向的到达角估计方法 | |
| CN108683619B (zh) | 一种低复杂度的大规模mimo信道参数估计方法 | |
| CN114095318B (zh) | 智能超表面辅助的混合构型毫米波通信系统信道估计方法 | |
| CN109861933B (zh) | 一种基于music算法和预编码的毫米波mimo信道估计方法 | |
| CN110824414A (zh) | 一种到达角估计的装置及方法 | |
| CN111983556A (zh) | 一种到达角估计的装置及方法 | |
| CN113189592A (zh) | 考虑幅相互耦误差的车载毫米波mimo雷达测角方法 | |
| Lin et al. | 3D wideband mmWave localization for 5G massive MIMO systems | |
| CN114726686B (zh) | 一种均匀面阵毫米波大规模mimo信道估计方法 | |
| CN109946663B (zh) | 一种线性复杂度的Massive MIMO目标空间方位估计方法和装置 | |
| Zhou et al. | Individual Channel Estimation for RIS-aided Communication Systems-A General Framework | |
| CN110719127A (zh) | 一种具有恒模约束的毫米波mimo系统波束成形方法 | |
| CN114520699B (zh) | 可重构智能表面辅助毫米波通信的低复杂度信道估计方法 | |
| CN110636018A (zh) | 一种网格补偿大规模mimo信道估计方法 | |
| CN114966523A (zh) | 大规模mimo测向中基于幂迭代算法的doa估计方法 | |
| CN110213186A (zh) | 基于混合结构的esprit的到达角与信噪比估计方法 | |
| Adrian-Ionut et al. | A speed convergence Least Squares Constant Modulus Algorithm for smart antenna beamforming | |
| CN115766350A (zh) | 大规模mimo系统中一种同时信道估计与定位方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |