CN110837075A - 一种低复杂度的极化参数估计跟踪装置及方法 - Google Patents

一种低复杂度的极化参数估计跟踪装置及方法 Download PDF

Info

Publication number
CN110837075A
CN110837075A CN201911109436.3A CN201911109436A CN110837075A CN 110837075 A CN110837075 A CN 110837075A CN 201911109436 A CN201911109436 A CN 201911109436A CN 110837075 A CN110837075 A CN 110837075A
Authority
CN
China
Prior art keywords
polarization
axis
submodule
sub
parameter estimation
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.)
Pending
Application number
CN201911109436.3A
Other languages
English (en)
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.)
Hangzhou Dianzi University
Hangzhou Electronic Science and Technology University
Original Assignee
Hangzhou Electronic Science and Technology 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 Hangzhou Electronic Science and Technology University filed Critical Hangzhou Electronic Science and Technology University
Priority to CN201911109436.3A priority Critical patent/CN110837075A/zh
Publication of CN110837075A publication Critical patent/CN110837075A/zh
Pending legal-status Critical Current

Links

Images

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/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)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明公开了一种低复杂度的极化参数估计跟踪装置及方法，装置包括天线信号处理模块和极化参数估计模块，所述天线信号处理模块包括若干子阵信号处理模块，将混合双极化天线阵列接收到的模拟信号波束成型并转换为数字信号；所述极化参数估计模块用于将天线信号处理模块输出进行信号处理，估计来波极化参数。本发明采用平面混合双极化阵列来对信号进行接收，同时充分利用来波到达角自适应跟踪与估计以提高波束成型接收信噪比，从而提高极化参数估计的精度和收敛速度，而且在计算过程中仅仅涉及普通的四则运算，计算量小便于快速实现。

Description

一种低复杂度的极化参数估计跟踪装置及方法
技术领域
本发明属于信号处理技术领域，涉及一种低复杂度的极化参数估计跟踪装 置及方法。
背景技术
混合天线阵列作为一种支持未来毫米波高速通信的新型自适应天线结构， 由于其在平衡硬件成本和性能、阵列增益和收发机尺寸等方面的潜在优势，受 到了阵列信号处理领域的广泛关注。它由多个模拟子阵列组成，每个子阵列都 包含一个数字处理链路。一方面，使用模拟子阵可以有效地产生天线增益，并 缓解射频链路的空间限制问题。另一方面，多个数字链路可以提供复用能力、 高灵活性、适应性和性能优化。
传统的单极化天线只能接收到入射波的单方向电场分量，不能捕捉到完整 的来波极化状态。特别是当天线的极化方向与入射波的极化方向垂直时，接收 信号的强度将完全丧失。因此，双极化接收天线是最佳接收的必要条件。为了 获得最佳接收，来波的极化状态参数需要进行估计来对接收信号极化分级合并。 混合双极化天线阵列的集成将充分利用两者的优点来提高系统性能。
基于双极化阵列的信号极化参数估计是双极化阵列信号处理的一个重要研 究内容，在过去三十年，收到广泛的关注，并提出了许多有效的极化参数估计 方法。传统的最大似然估计方法由于需要相当大的计算量，在实际应用中并不 能被采用。经典的MUSIC和ESPRIT方法可以获得比较高的参数估计精度，但 是他们要求协方差矩阵和奇异值分解，其带来的计算量与天线数量的立方成比 例增长，故引入了高代价的系统硬件实现成本和复杂度。因此，这些算法仅适 用于小数量的数字天线阵列，并不能真正应用于未来大规模的天线阵列。
发明内容
为解决上述问题，本发明的目的在于提供一种低复杂度的极化参数估计跟 踪装置及方法。采用平面混合双极化阵列来对信号进行接收，同时充分利用来 波到达角自适应跟踪与估计以提高波束成型接收信噪比，从而提高极化参数估 计的精度和收敛速度，而且在计算过程中仅仅涉及普通的四则运算，计算量小 便于快速实现。
为实现上述目的，本发明的技术方案为：
一种低复杂度的极化参数估计跟踪装置，包括天线信号处理模块和极化参 数估计模块，其中，
所述天线信号处理模块包括若干子阵信号处理模块，将混合双极化天线阵 列接收到的模拟信号波束成型并转换为数字信号；
所述极化参数估计模块用于将天线信号处理模块输出进行信号处理，估计 来波极化参数，包括互相关计算子模块、子阵相位控制子模块、数字波束成型 子模块、极化参数估计子模块和最大比合并子模块，所述互相关计算子模块对 所有空间相邻的子阵信号处理模块的输出信号在两个正交极化方向分别进行互 相关运算，进而估计到达角度信息；所述子阵相位控制子模块将互相关计算子 模块输出的结果反馈控制子阵信号处理模块；所述数字波束成型子模块将互相 关计算子模块输出的结果对多个子阵信号处理模块输出的信号分别进行加权波 束成型；所述极化参数估计子模块对数字波束成型子模块输出信号和互相关计 算子模块输出的结果进行计算，输出来波极化参数估计；所述最大比合并子模块计算最大比合并系数，对数字波束成型子模块输出信号进行最大比合并。
优选地，所述子阵信号处理模块包括依次连接的子阵天线子模块、子阵相 移器子模块和子阵信号转换子模块，所述子阵天线子模块包括若干个由沿x轴 和y轴极化方向的阵子组成的双极化阵子，用于接收来波模拟信号；所述子阵 相移器子模块包括若干个模拟相移器，用于调节双极化阵子的接收信号相位， 由子阵相位控制子模块反馈设置；所述子阵信号转换子模块包括模拟信号合并 器，用于合并模拟相移器的输出，和将合并信号转换到数字基带。
基于上述目的，本发明还提供了一种低复杂度的极化参数估计跟踪方法， 其特征在于，包括以下步骤：
S10，对混合双极化阵列沿着x轴和y轴极化方向上的信号分别进行相邻 子阵输出信号的相关运算；
S20，合并后获得到达角信息估计；
S30，利用到达角信息估计对数字波束成型的加权系数进行设置；
S40，利用沿着x轴和y轴极化方向上的数字波束成型输出信号估计极化 参数；
S50，利用到达角信息估计和极化参数估计对数字波束成型输出信号进行 加权合并，判断迭代是否达到要求；
S60，否，则利用到达角信息估计对所有相移器进行更新；
S70，是，则根据系统性能要求，输出极化参数估计结果。
优选地，所述S10为在互相关计算子模块中，计算所有相邻的子阵信号处 理模块输出的互相关，对两个极化方向的信号求和，利用下式(1)和(2)分 别对沿x轴和u轴方向子阵互相关结果进行滤波，得到第i个时刻滤波后的互相关
Figure BDA0002272285600000031
和
Figure BDA0002272285600000032
Figure BDA0002272285600000033
Figure BDA0002272285600000034
其中，0＜1为滤波系数；c轴为x轴或y轴；Mx和My分别为沿x轴和 y轴方向子阵的数量；mx＝0,1,……,Mx-1、my＝0,1,……,My-1； 为第myMx+mx个子阵在c轴极化方向第i个时刻的输出采样信号， ()*为表示共轭。
优选地，所述S20为在互相关计算子模块中，估算到达角信息
Figure BDA0002272285600000036
和其中，arg{}表示取幅角，
Figure BDA0002272285600000038
和
Figure BDA0002272285600000039
包含来 波的天顶角和方位角信息，
Figure BDA0002272285600000041
和
Figure BDA0002272285600000042
分别为沿x轴和y轴方向子阵的间距，假设
Figure BDA0002272285600000043
由
Figure BDA0002272285600000044
和
Figure BDA0002272285600000045
得到
Figure BDA0002272285600000046
和
Figure BDA0002272285600000047
θ为来波的天顶角，φ为来波的方位角。
优选地，所述S30为在数字波束成型子模块中，根据下式(3)对第myMx+ mx个子阵信号处理模块沿着c轴极化方向上的输出信号进行加权，加权值为
Figure BDA0002272285600000048
产生数字波束成型信号sc[i],由式(4)得到，
Figure BDA0002272285600000049
Figure BDA00022722856000000410
优选地，所述S40为在极化参数估计子模块中，由下式(5)和(6)对数 字波束成型子模块的输出互相关和总功率分别进行滤波，得到第i个时刻滤波后 的互相关和总功率为，
Figure BDA00022722856000000411
Figure BDA00022722856000000412
其中，0＜β＜1为滤波系数，
Figure BDA00022722856000000413
为子阵输出的加性零均值高斯白噪声功率， 由式(7)分别得到第i个时刻的互相关与功率比值，
Figure BDA00022722856000000414
其中，为
Figure BDA00022722856000000416
的实部，Vc (i)为
Figure BDA00022722856000000417
的虚部，j为虚数单位， 依据x轴和y轴极化方向接收信号的功率大小，x轴对应式(8)，y轴对应式 (9)，进行第i个时刻的极化状态参数估计，
Figure BDA00022722856000000418
这里，
记
Figure BDA0002272285600000053
Figure BDA0002272285600000054
这里，
Figure BDA0002272285600000055
Figure BDA0002272285600000056
记
其中，极化参数信息包括极化辅助角γ和极化相位差η。
优选地，所述S50为在最大比合并子模块中，由下式(10)计算第i个时刻 c轴极化方向的最大比合并系数
Figure BDA0002272285600000058
Figure BDA0002272285600000059
记
Figure BDA00022722856000000511
输出最大比合并信号
Figure BDA00022722856000000512
在极化参数估计子模块中，判断 迭代是否达到要求。
优选地，所述S60为迭代未达到要求，在子阵相位控制子模块中，由下式 (11)对第nyNx+nx个相移器第i+1个时刻的相移值进行设置和样点更新，并 返回S10，
Figure BDA0002272285600000061
其中，Nx和Ny分别为每个子阵中沿x轴和y轴方向阵子的数量，nx＝ 0,1,……,Nx-1,ny＝0,1,……,Ny-1。
优选地，所述S70为迭代达到要求，跳出循环，输出S40中极化状态参数 估计γ(I)和η(I)，其中，I为迭代要求的样点数量。
与现有技术相比，本发明的有益效果如下：利用混合双极化阵列的相邻子 阵输出信号进行相关运算，和对数字波束成型输出信号进行普通的四则运算来 估计极化参数，使得估计的复杂度大幅度降低。另外，由于对相移器的自适应 调整，使得波束逐渐调整到来波方向，从而提高了极化参数估计精度，收敛速 度和合并信号信噪比。同时，在极化状态时变情况下，本发明可以自适应地跟 踪极化状态。
附图说明
图1为本发明实施例的低复杂度的极化参数估计跟踪装置的结构示意图；
图2为本发明一实施例的低复杂度的极化参数估计跟踪装置的混合双极化 天线阵列示意图；
图3为本发明实施例的低复杂度的极化参数估计跟踪方法的步骤流程图；
图4为本发明实施例的低复杂度的极化参数估计跟踪的均方误差MSE性能 对迭代次数的关系曲线图；
图5为本发明实施例的低复杂度的极化参数估计跟踪的合并信号输出功率 对迭代次数的关系曲线图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实 施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅 仅用以解释本发明，并不用于限定本发明。
相反，本发明涵盖任何由权利要求定义的在本发明的精髓和范围上做的替 代、修改、等效方法以及方案。进一步，为了使公众对本发明有更好的了解， 在下文对本发明的细节描述中，详尽描述了一些特定的细节部分。对本领域技 术人员来说没有这些细节部分的描述也可以完全理解本发明。
装置实施例1
参见图1，所示为本发明实施例的低复杂度的极化参数估计跟踪装置的结 构示意图一种低复杂度的极化参数估计跟踪装置，包括天线信号处理模块11， 用于将图2所示的混合双极化天线阵列接收到的模拟信号波束成型并转换为数 字信号；
极化参数估计模块12，用于将天线信号处理模块11输出进行信号处理， 估计来波极化参数。
天线信号处理模块11包括多个子阵信号处理模块111，子阵信号处理模块 111进一步包括多个由沿x轴和y轴极化方向的阵子组成的双极化阵子，用于接 收来波模拟信号；子阵相移器子模块1112包括多个用于调节阵子接收信号相位 的模拟相移器，由子阵相位控制子模块122反馈进行设置；子阵信号转换子模 块1113包括合并相移器输出的模拟信号合并器，并将合并信号转换到数字基带， 这里包括将射频或中频信号下变频到基带。
极化参数估计模块12包括互相关计算子模块121，用于对所有空间相邻的 子阵信号处理模块111的输出信号在两个正交极化方向分别进行互相关运算， 进而估计到达角度信息；子阵相位控制子模块122，用于将互相关计算子模块 121输出的结果反馈控制子阵相移器子模块1112中的模拟相移器；数字波束成 型子模块123，用于将互相关计算子模块121输出的结果对多个子阵信号处理 模块111输出的信号(沿x轴和y轴极化方向)分别进行加权波束成型；极化参 数估计子模块124，用于对数字波束成型子模块123输出信号和互相关计算子 模块121输出的结果进行计算，输出来波极化参数估计；最大比合并子模块125， 用于计算最大比合并系数，对数字波束成型子模块123输出信号进行最大比合 并。
参见图3，为本发明的方法实施例，一种低复杂度的极化参数估计跟踪方 法，包括以下步骤：
S10，对混合双极化阵列沿着x轴和y轴极化方向上的信号分别进行相邻 子阵输出信号的相关运算；
S20，合并后获得到达角信息估计；
S30，利用到达角信息估计对数字波束成型的加权系数进行设置；
S40，利用沿着x轴和y轴极化方向上的数字波束成型输出信号估计极化 参数；
S50，利用到达角信息估计和极化参数估计对数字波束成型输出信号进行 加权合并，判断迭代是否达到要求；
S60，否，则利用到达角信息估计对所有相移器进行更新；
S70，是，则根据系统性能要求，输出极化参数估计结果。
上述方法基于以下公式或定义，
假设采用MxNx×MyNy个双极化天线构成的均匀平面方阵接收一个来波信 号
Figure BDA0002272285600000082
其中Mx和My分别表示沿x轴和y轴方向子阵的数量；Nx和Ny分别表示 每个子阵中沿x轴和y轴方向阵子的数量。令第m＝myMx+mx个子阵对应的模拟 波束成型在c(x或y)轴极化方向输出信号表示为
Figure BDA0002272285600000081
其中，mx＝0,1,……,Mx-1；my＝0,1,……,My-1；Ec表示c轴极化方向的 天线响应；Ps(θ,φ)表示子阵的方向图；λ和(θ，φ)分别表示来波信号的波 长，天顶角和方位角；和
Figure BDA0002272285600000092
分别表示沿x轴和y轴方向子阵的间距；
Figure BDA0002272285600000093
表 示对应子阵输出的加性零均值高斯白噪声，其功率表示为
Figure BDA0002272285600000094
需要估计的极化参数信息包括极化辅助角γ和极化相位差η包含在Ec中，他 们的关系可以表示为，
Ex＝sinγcosθcosφejη-cosγsinφ；
Ey＝sinγcosθsinφejη+cosγcosφ；
首先对模拟波束成型输出信号进行采样得到
Figure RE-GDA0002363624780000095
然后通过数字加权获得数字波束成型信号
Figure RE-GDA0002363624780000097
最后对sc[i] 进行加权合并得到
Figure RE-GDA0002363624780000098
为c轴极化方向的合并系数。
S10为在互相关计算子模块中，计算所有相邻的子阵信号处理模块输出的 互相关，对两个极化方向的信号求和，利用下式(1)和(2)分别对沿x轴和y轴 方向子阵互相关结果进行滤波，得到第i个时刻滤波后的互相关
Figure BDA0002272285600000099
和
Figure BDA00022722856000000910
Figure BDA00022722856000000911
其中，0＜μ＜1为滤波系数；c轴为x轴或y轴；Mx和My分别为沿x轴和 y轴方向子阵的数量；mx＝0,1,……,Mx-1、my＝0,1,……,My-1；
Figure BDA00022722856000000913
为第myMx+mx个子阵在c轴极化方向第i个时刻的输出采样信号， ()*为表示共轭。
S20为在互相关计算子模块中，估算到达角信息
Figure BDA00022722856000000914
和
Figure BDA00022722856000000915
其中，arg{}表示取幅角，
Figure BDA00022722856000000916
和包含来波的天顶角和方位 角信息，
Figure BDA0002272285600000101
和分别为沿x轴和y轴方向子阵的间距，假设
Figure BDA0002272285600000103
由
Figure BDA0002272285600000104
和
Figure BDA0002272285600000105
得到和
Figure BDA0002272285600000107
θ为来波的天顶角，φ 为来波的方位角。
S30为在数字波束成型子模块中，根据下式(3)对第myMx+mx个子阵信 号处理模块沿着c轴极化方向上的输出信号进行加权，加权值为
Figure BDA0002272285600000108
产 生数字波束成型信号sc[i],由式(4)得到，
Figure BDA0002272285600000109
Figure BDA00022722856000001010
S40为在极化参数估计子模块中，由下式(5)和(6)对数字波束成型子 模块的输出互相关和总功率分别进行滤波，得到第i个时刻滤波后的互相关和总 功率为，
Figure BDA00022722856000001011
Figure BDA00022722856000001012
其中，0＜β＜1为滤波系数，
Figure BDA00022722856000001013
为子阵输出的加性零均值高斯白噪声功率， 由式(7)分别得到第i个时刻的互相关与功率比值，
Figure BDA00022722856000001014
其中，为的实部，jVc (i)为
Figure BDA00022722856000001017
的虚部，j为虚数单位， 依据x轴和y轴极化方向接收信号的功率大小，x轴对应式(8)，y轴对应式 (9)，进行第i个时刻的极化状态参数估计，
Figure BDA00022722856000001018
这里，
Figure BDA0002272285600000111
Figure BDA0002272285600000112
记
Figure BDA0002272285600000113
Figure BDA0002272285600000114
这里，
Figure BDA0002272285600000115
Figure BDA0002272285600000116
记
Figure BDA0002272285600000117
其中，极化参数信息包括极化辅助角γ和极化相位差η。
S50为在最大比合并子模块中，由下式(10)计算第i个时刻c轴极化方向 的最大比合并系数
Figure BDA0002272285600000118
Figure BDA0002272285600000119
记
输出最大比合并信号在极化参数估计子模块中，判断 迭代是否达到要求。
S60为迭代未达到要求，在子阵相位控制子模块中,由下式(11)对第 nyNx+nx个相移器第i+1个时刻的相移值进行设置和样点更新，并返回S10，
Figure BDA0002272285600000121
其中，Nx和Ny分别为每个子阵中沿x轴和y轴方向阵子的数量，nx＝ 0,1,……,Nx-1,ny＝0,1,……,Ny-1。
S70为迭代达到要求，跳出循环，输出S40中极化状态参数估计γ(I)和η(I)， 其中，I为迭代要求的样点数量。
为评价本发明所提出的极化参数的估计方案，对其均方误差(MSE)性能 和最大比合并信号输出信号噪声功率比(SNR)进行了计算机仿真，还将结果 与传统的混合单极化阵列方案进行比较。为了进行比较，采用Mx＝My＝ 2,Nx＝Ny＝4,d＝λ/2。
图4为本发明提供的方案极化参数估计MSE性能对迭代次数的计算机仿 真结果。为了获得统计性能，仿真假设到达角度和极化参数在设定范围内均匀 分布，并执行10000次的独立仿真。从图4可以看出，波形B11为利用式(8) 的极化相位差η，B12为利用式(9)的极化相位差η，波形B21为利用式(8)的极 化辅助角γ，B22为利用式(9)的极化辅助角γ，本发明的方案随着迭代次数的 增加，估计的MSE性能得到明显改善。
图5为本发明提供的方案与现有技术的最大比合并信号输出SNR对迭代次 数的计算机仿真结果。这里假设来波的极化参数γ＝π/4和η＝-π/4，到达角 信息
Figure BDA0002272285600000122
和uy＝0，极化波在沿x轴和y轴极化方向的阵子上接收功率比 为0.5。从图5可以看出，波形B31为本发明的最大比合并信号输出SNR，波 形B32为现有技术中只有y轴极化的合并信号输出SNR，形B33为现有技术中 只有x轴极化的合并信号输出SNR，本发明的方案由于采用了两个极化方向信 号的最大比合并，与现有技术相比获得了更大的接收SNR。另外，本发明的方 案随着迭代次数的增加，估计的角度信息和极化参数精确度提高，因而合并后 的信号SNR性能得到明显改善。
上述所提发明适用于高速数据通信，例如卫星通信和雷达感知。所公开的 装置非常适合毫米波至太赫兹波极化参数的自适应跟踪估计。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发 明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明 的保护范围之内。

Claims (10)

1.一种低复杂度的极化参数估计跟踪装置，其特征在于，包括天线信号处理模块和极化参数估计模块，其中，
所述天线信号处理模块包括若干子阵信号处理模块，将混合双极化天线阵列接收到的模拟信号波束成型并转换为数字信号；
所述极化参数估计模块用于将天线信号处理模块输出进行信号处理，估计来波极化参数，包括互相关计算子模块、子阵相位控制子模块、数字波束成型子模块、极化参数估计子模块和最大比合并子模块，所述互相关计算子模块对所有空间相邻的子阵信号处理模块的输出信号在两个正交极化方向分别进行互相关运算，进而估计到达角度信息；所述子阵相位控制子模块将互相关计算子模块输出的结果反馈控制子阵信号处理模块；所述数字波束成型子模块将互相关计算子模块输出的结果对多个子阵信号处理模块输出的信号分别进行加权波束成型；所述极化参数估计子模块对数字波束成型子模块输出信号和互相关计算子模块输出的结果进行计算，输出来波极化参数估计；所述最大比合并子模块计算最大比合并系数，对数字波束成型子模块输出信号进行最大比合并。
2.根据权利要求1所述的低复杂度的极化参数估计跟踪装置，其特征在于，所述子阵信号处理模块包括依次连接的子阵天线子模块、子阵相移器子模块和子阵信号转换子模块，所述子阵天线子模块包括若干个由沿x轴和y轴极化方向的阵子组成的双极化阵子，用于接收来波模拟信号；所述子阵相移器子模块包括若干个模拟相移器，用于调节双极化阵子的接收信号相位，由子阵相位控制子模块反馈设置；所述子阵信号转换子模块包括模拟信号合并器，用于合并模拟相移器的输出，和将合并信号转换到数字基带。
3.一种采用权利要求1或2所述装置的低复杂度的极化参数估计跟踪方法，其特征在于，包括以下步骤：
S10，对混合双极化阵列沿着x轴和y轴极化方向上的信号分别进行相邻子阵输出信号的相关运算；
S20，合并后获得到达角信息估计；
S30，利用到达角信息估计对数字波束成型的加权系数进行设置；
S40，利用沿着x轴和y轴极化方向上的数字波束成型输出信号估计极化参数；
S50，利用到达角信息估计和极化参数估计对数字波束成型输出信号进行加权合并，判断迭代是否达到要求；
S60，否，则利用到达角信息估计对所有相移器进行更新；
S70，是，则根据系统性能要求，输出极化参数估计结果。
4.根据权利要求3所述的方法，其特征在于，所述S10为在互相关计算子模块中，计算所有相邻的子阵信号处理模块输出的互相关，对两个极化方向的信号求和，利用下式(1)和(2)分别对沿x轴和y轴方向子阵互相关结果进行滤波，得到第i个时刻滤波后的互相关
Figure FDA0002272285590000021
和
Figure FDA0002272285590000022
Figure FDA0002272285590000023
Figure FDA0002272285590000024
其中，0＜μ＜1为滤波系数；c轴为x轴或y轴；Mx和My分别为沿x轴和y轴方向子阵的数量；mx＝0,1,……,Mx-1、my＝0,1,……,My-1；
Figure FDA0002272285590000025
为第myMx+mx个子阵在c轴极化方向第i个时刻的输出采样信号，()*为表示共轭。
5.根据权利要求3所述的方法，其特征在于，所述S20为在互相关计算子模块中，估算到达角信息和其中，arg{}表示取幅角，
Figure FDA0002272285590000028
和包含来波的天顶角和方位角信息，
Figure FDA00022722855900000210
和
Figure FDA00022722855900000211
分别为沿x轴和y轴方向子阵的间距，假设
Figure FDA00022722855900000212
由和
Figure FDA00022722855900000214
得到
Figure FDA0002272285590000031
和
Figure FDA0002272285590000032
θ为来波的天顶角，φ为来波的方位角。
6.根据权利要求5所述的方法，其特征在于，所述S30为在数字波束成型子模块中，根据下式(3)对第myMx+mx个子阵信号处理模块沿着c轴极化方向上的输出信号进行加权，加权值为
Figure FDA0002272285590000033
产生数字波束成型信号sc[i],由式(4)得到，
Figure FDA0002272285590000034
Figure FDA0002272285590000035
。
7.根据权利要求6所述的方法，其特征在于，所述S40为在极化参数估计子模块中，由下式(5)和(6)对数字波束成型子模块的输出互相关和总功率分别进行滤波，得到第i个时刻滤波后的互相关和总功率为，
Figure FDA0002272285590000036
Figure FDA0002272285590000037
其中，0＜β＜1为滤波系数，
Figure FDA0002272285590000038
为子阵输出的加性零均值高斯白噪声功率，由式(7)分别得到第i个时刻的互相关与功率比值，
Figure FDA0002272285590000039
其中，
Figure FDA00022722855900000310
为
Figure FDA00022722855900000311
的实部，Vc (i)为
Figure FDA00022722855900000312
的虚部，j为虚数单位，依据x轴和y轴极化方向接收信号的功率大小，x轴对应式(8)，y轴对应式(9)，进行第i个时刻的极化状态参数估计，
Figure FDA00022722855900000313
这里，
Figure FDA0002272285590000042
记
Figure FDA0002272285590000043
Figure FDA0002272285590000044
这里，
Figure FDA0002272285590000045
Figure FDA0002272285590000046
记
Figure FDA0002272285590000047
其中，极化参数信息包括极化辅助角γ和极化相位差η。
8.根据权利要求7所述的方法，其特征在于，所述S50为在最大比合并子模块中，由下式(10)计算第i个时刻c轴极化方向的最大比合并系数
Figure FDA0002272285590000048
Figure FDA0002272285590000049
记
Figure FDA00022722855900000410
Figure FDA00022722855900000411
输出最大比合并信号
Figure FDA0002272285590000051
在极化参数估计子模块中，判断迭代是否达到要求。
9.根据权利要求8所述的方法，其特征在于，所述S60为迭代未达到要求，在子阵相位控制子模块中,由下式(11)对第nyNx+nx个相移器第i+1个时刻的相移值进行设置和样点更新，并返回S10，
Figure FDA0002272285590000052
其中，Nx和Ny分别为每个子阵中沿x轴和y轴方向阵子的数量，nx＝0,1,……,Nx-1,ny＝0,1,……,Ny-1。
10.根据权利要求8所述的方法，其特征在于，所述S70为迭代达到要求，跳出循环，输出S40中极化状态参数估计γ(I)和η(I)，其中，I为迭代要求的样点数量。
CN201911109436.3A 2019-11-13 2019-11-13 一种低复杂度的极化参数估计跟踪装置及方法 Pending CN110837075A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911109436.3A CN110837075A (zh) 2019-11-13 2019-11-13 一种低复杂度的极化参数估计跟踪装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911109436.3A CN110837075A (zh) 2019-11-13 2019-11-13 一种低复杂度的极化参数估计跟踪装置及方法

Publications (1)

Publication Number Publication Date
CN110837075A true CN110837075A (zh) 2020-02-25

Family

ID=69574854

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911109436.3A Pending CN110837075A (zh) 2019-11-13 2019-11-13 一种低复杂度的极化参数估计跟踪装置及方法

Country Status (1)

Country Link
CN (1) CN110837075A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112731278A (zh) * 2020-12-28 2021-04-30 杭州电子科技大学 一种部分极化信号的角度与极化参数欠定联合估计方法
CN112887040A (zh) * 2021-01-18 2021-06-01 杭州电子科技大学 一种极化参数估计装置及方法
SE2100129A1 (en) * 2021-08-31 2023-03-01 Saab Ab Multi-channel active array system and Method for obtaining positional information of an object

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015194042A1 (ja) * 2014-06-20 2015-12-23 日産自動車株式会社 電波画像生成装置及び電波画像生成方法
WO2017008307A1 (en) * 2015-07-16 2017-01-19 Nec Corporation Method and apparatus for performing beamforming
CN106848546A (zh) * 2017-01-24 2017-06-13 哈尔滨工业大学（威海） 一种宽带双极化天线阵列装置及高分辨测向方法
CN109738853A (zh) * 2018-11-21 2019-05-10 南京航空航天大学 一种电磁矢量互质阵基于旋转不变性的角度与极化估计方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015194042A1 (ja) * 2014-06-20 2015-12-23 日産自動車株式会社 電波画像生成装置及び電波画像生成方法
WO2017008307A1 (en) * 2015-07-16 2017-01-19 Nec Corporation Method and apparatus for performing beamforming
CN106848546A (zh) * 2017-01-24 2017-06-13 哈尔滨工业大学（威海） 一种宽带双极化天线阵列装置及高分辨测向方法
CN109738853A (zh) * 2018-11-21 2019-05-10 南京航空航天大学 一种电磁矢量互质阵基于旋转不变性的角度与极化估计方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
HANG LI等: "Adaptive AoA and Polarization Estimation for Receiving Polarized mmWave Signals", 《IEEE WIRELESS COMMUNICATIONS LETTERS》 *
HANG LI等: "Enhanced AoA Estimation Using Localized Hybrid Dual-Polarized Arrays", 《 2019 IEEE 90TH VEHICULAR TECHNOLOGY CONFERENCE》 *
HANG LI等: "Joint Adaptive AoA and Polarization Estimation Using Hybrid Dual-Polarized Antenna Arrays", 《IEEE ACCESS 》 *
HANG LI等: "Low-Complexity Multiuser Receiver for Massive Hybrid Array mmwave Communications", 《IEEE TRANSACTIONS ON COMMUNICATIONS》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112731278A (zh) * 2020-12-28 2021-04-30 杭州电子科技大学 一种部分极化信号的角度与极化参数欠定联合估计方法
CN112731278B (zh) * 2020-12-28 2023-11-03 杭州电子科技大学 一种部分极化信号的角度与极化参数欠定联合估计方法
CN112887040A (zh) * 2021-01-18 2021-06-01 杭州电子科技大学 一种极化参数估计装置及方法
CN112887040B (zh) * 2021-01-18 2022-10-25 杭州电子科技大学 一种极化参数估计装置及方法
SE2100129A1 (en) * 2021-08-31 2023-03-01 Saab Ab Multi-channel active array system and Method for obtaining positional information of an object
WO2023033693A1 (en) * 2021-08-31 2023-03-09 Saab Ab Multi-channel active array system and method for obtaining positional information of an object
SE545126C2 (en) * 2021-08-31 2023-04-11 Saab Ab Multi-channel active array system and Method for obtaining positional information of an object

Similar Documents

Publication Publication Date Title
CN108462521B (zh) 自适应阵列天线的抗干扰实现方法
CN108051782B (zh) 基于子阵划分的大规模相控阵差波束形成系统
CN110837075A (zh) 一种低复杂度的极化参数估计跟踪装置及方法
CN109787671B (zh) 一种混合波束成形装置及方法
CN112737647A (zh) 一种基于混合连接的毫米波大规模mimo预编码方法
CN110708103B (zh) 一种无需预延迟的宽带波束形成方法
CN110824414A (zh) 一种到达角估计的装置及方法
CN112636773A (zh) 基于数字频域补偿的宽带时域波束成形方法
CN113162670A (zh) 数字多波束校正与合成方法
CN111859644A (zh) 一种圆形相控阵列天线的接收波束形成与方位扫描方法
CN111983556A (zh) 一种到达角估计的装置及方法
CN108667489B (zh) 多波束波形发射方法及系统
CN113572503B (zh) 一种低复杂度的基于gp的改进混合波束成形方法
CN113067616B (zh) 基于时间调制数字超表面的dbf相控阵系统
CN111162878B (zh) 一种基于子阵降维带约束的多域联合抗干扰方法
CN112887040B (zh) 一种极化参数估计装置及方法
Okkonen Uniform linear adaptive antenna array beamforming implementation with a wireless open-access research platform
CN106842147B (zh) 一种解决栅瓣干扰问题的数字波束形成方法
CN111898087B (zh) 阵列天线子向量循环约束优化波束形成系统及方法
CN112526510A (zh) 基于方向图分集的单通道角度超分辨方法及系统
Zhang et al. A low-complexity Laguerre wideband beamformer
CN115833894B (zh) 基于子阵的数字模拟合成自适应抗干扰方法
Sun et al. Shared aperture multibeam forming of time-modulated linear array
RU2255396C2 (ru) Способ энергетической оптимизации моноимпульсных антенных решеток с совместным формированием лучей
Yu et al. Adaptive beamforming technique for large-scale arrays with various subarray selections

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20200225

RJ01 Rejection of invention patent application after publication