CN112887040B - 一种极化参数估计装置及方法 - Google Patents
一种极化参数估计装置及方法 Download PDFInfo
- Publication number
- CN112887040B CN112887040B CN202110061089.2A CN202110061089A CN112887040B CN 112887040 B CN112887040 B CN 112887040B CN 202110061089 A CN202110061089 A CN 202110061089A CN 112887040 B CN112887040 B CN 112887040B
- Authority
- CN
- China
- Prior art keywords
- parameter estimation
- polarization parameter
- submodule
- polarization
- sub
- 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
- 230000010287 polarization Effects 0.000 title claims abstract description 147
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000012545 processing Methods 0.000 claims abstract description 54
- 238000003491 array Methods 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 18
- 239000000654 additive Substances 0.000 claims description 9
- 230000000996 additive effect Effects 0.000 claims description 9
- 238000004088 simulation Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000010363 phase shift Effects 0.000 claims description 3
- 238000004891 communication Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/20—Monitoring; Testing of receivers
- H04B17/29—Performance testing
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种极化参数估计的装置及方法,该装置至少包括:天线信号处理模块:用于将阵列天线接收到的模拟信号波束成型,并转换为数字信号;极化参数估计模块:用于将天线信号处理模块输出进行信号处理估计来波的极化参数。相对于传统的互相关能量比极化参数估计方法只利用x轴或者y轴极化方向的接收信号能量进行极化参数估计,特别地,当所采用极化方向的接收信号能量低时,极化参数估计的精度会大大下降。本发明采用平面混合双极化阵列来对信号进行接收,联合利用来波到达角估计信息对来波信号进行极化参数估计,从而获得最大比合并信号,同时联合利用两个正交的极化方向接收信号能量进行极化参数估计,从而提高极化参数估计的精度。
Description
技术领域
本发明涉及阵列信号处理领域,尤其涉及一种高精度的极化参数估计装置及方法。
背景技术
全数字天线阵列可以应用在无线通信***中获得最优性能,但是其在未来大规模天线***中引入的高成本及高实现复杂度使其难以实现。混合天线阵列作为一种支持未来毫米波高速通信的新型天线结构,在平衡硬件成本和性能、阵列增益和收发机尺寸等方面具有优势,因此受到了阵列信号处理领域的广泛关注。混合阵列由多个模拟子阵列组成,其中每个天线接收信号的相位由相移器进行调整,每个子阵列都连接一个数字处理链路。一方面,使用模拟子阵可以有效地产生天线增益,并减小射频链路的硬件成本。另一方面,多个数字链路可以提供灵活地复用能力并对***性能进行联合优化。
传统天线并不考虑极化特性,因此只能接收到入射波的单方向电场分量,无法捕捉到完整的来波极化信息。尤其当天线的极化方向与入射波的极化方向垂直时,接收信号的强度为零。因此,双极化接收天线是最佳接收的必要条件。为了获得最佳接收,来波的极化信息需要进行估计来对接收信号进行极化最大比合并。混合双极化天线阵列的集成可以联合两者的优点来对***性能进行大幅提升。
基于双极化阵列的信号极化参数估计在过去三十年里一直获得信号处理领域的广泛关注,期间许多经典及有效的估计方法被提出并深入研究。性能最优的最大似然估计方法要求相当大的计算量,因此在实际应用中无法被采用。经典的MUSIC和ESPRIT方法可以获得比较高的参数估计精度,但是协方差矩阵和奇异值分解引入的计算量与天线数量的立方成正比,使得这些算法仅适用于小数量的天线阵列,无法真正应用于未来大规模的天线阵列。近来被提出的基于混合双极化阵列的极化参数估计方法仅利用x轴或者y轴极化方向的接收信号能量进行极化参数估计,使得另一维的接收信号没有有效利用,导致估计的精度损失。
故,针对现有技术存在的技术缺陷,实有必要提出一种解决方案克服现有技术存在的技术问题。
发明内容
有鉴于此,本发明提供一种极化参数估计的装置及方法,相对于传统的互相关能量比极化参数估计方法只利用x轴或者y轴极化方向的接收信号能量进行极化参数估计,传统方法当所采用极化方向的接收信号能量低时,极化参数估计的精度会大大下降。本发明采用平面混合双极化阵列来对信号进行接收,联合利用来波到达角估计信息对来波信号进行极化参数估计,从而获得最大比合并信号。本方法联合利用两个正交的极化方向接收信号能量进行极化参数估计,从而提高极化参数估计的精度。
为达到上述目的,本发明的技术方案具体是这样实现的:
一种极化参数估计的装置,该装置包括:
天线信号处理模块:用于将阵列天线接收到的模拟信号波束成型,并转换为数字信号;
极化参数估计模块:用于将天线信号处理模块输出进行信号处理估计来波的极化参数;
其中,天线信号处理模块包括多个子阵信号处理模块,所述子阵信号处理模块包括:
子阵天线子模块:包含多个由沿x轴和y轴极化方向的阵子组成的双极化阵子,用于接收来波模拟信号;
子阵相移器子模块:包含多个用于调节阵子接收信号相位的模拟相移器,其根据来波到达角信息进行设置;
子阵信号转换子模块:包含用于合并相移器输出的模拟信号合并器,并将合并信号转换到数字基带,包括将射频或中频信号下变频到基带;
所述极化参数估计模块包括:
波达角估计子模块:用于对子阵信号处理模块在两个正交极化方向的输出信号进行处理,进而估计到达角度信息;根据到达角估计信息反馈对子阵相移器子模块中的模拟相移器进行设置;
极化参数估计子模块:用于对子阵信号处理模块在两个正交极化方向的输出信号进行处理,获得两组可选的极化参数估计值,再根据功率计算子模块反馈确定期望的极化参数估计值并输出;
最大比合并子模块:用于通过波达角估计子模块和极化参数估计子模块输出计算最大比合并系数,对子阵信号处理模块输出信号进行最大比合并;根据功率计算子模块反馈确定期望的最大比合并系数,同时输出期望的最大比合并信号;
功率计算子模块:用于计算最大比合并子模块输出的两路信号功率,通过比较功率大小确定期望的信号,同时反馈给极化参数估计子模块和最大比合并子模块。
本发明还公开了一种极化参数估计的方法,该方法包括以下步骤:
对混合双极化阵列沿着x轴和y轴极化方向上的接收信号进行处理,获得到达角信息估计;
利用估计的到达角信息对相移器进行设置,得到每个子阵沿着x轴和y轴极化方向上的波束成型输出信号;
利用波束成型输出信号估计极化参数,获得两组可选估计;
利用到达角信息估计和可选极化参数估计对波束成型输出信号进行加权合并,从两组可选估计中确定期望极化参数估计和最大比输出信号。
作为进一步的改进方案,该方法由极化参数估计模块执行,其进一步包括以下步骤:
步骤S4:在极化参数估计子模块中,利用等式(5)进行极化参数估计,得到两组可选估计
步骤S5:在最大比合并子模块中,利用等式(6)计算最大比合并系数
其中:
步骤S6:在功率计算子模块中,利用等式(7)对最大比合并子模块123的输出信号sk[i](k=1,2)的功率进行滤波,得到第i个时刻滤波后的平均功率:
其中,0<β<1是滤波系数;根据仿真选择最优的滤波系数;利用等式(8)确定期望索引,
这里I1表示迭代样点的要求数量;
本发明的技术效果如下:
传统的互相关能量比极化参数估计方法只利用x轴或者y轴极化方向的接收信号能量进行极化参数估计,传统方法当所采用极化方向的接收信号能量低时,极化参数估计的精度会大大下降。相对于现有技术,本发明正是由于联合利用沿着x轴和y轴极化方向上的波束成型输出信号进行处理来估计极化参数,使得估计的精度大幅度提高。
附图说明
图1为本发明实施例的极化参数估计装置的结构示意图;
图2为本发明实施例的极化参数估计装置的混合双极化天线阵列示意图;
图3为本发明实施例的极化参数估计方法的步骤流程图;
图4为本发明实施例的极化参数估计均方误差MSE性能对接收信噪比的关系曲线图;
具体实施方式
为使本发明的目的,技术方案及优点更加清楚明白,以下将参照附图并举实施例,对本发明作进一步说明。
参见图1,所示为本发明提供的极化参数估计装置第一实施例的流程与结构示意图,包括:
天线信号处理模块11:用于将混合双极化天线阵列接收到的模拟信号波束成型并转换为数字信号,参见图2,所示为本发明实施例的极化参数估计装置的混合双极化天线阵列示意图;
极化参数估计模块12:用于将天线信号处理模块11输出信号进行处理,估计来波极化参数。
进一步天线信号处理模块11包括:
多个子阵信号处理模块111。
进一步子阵信号处理模块111包括:
子阵天线子模块1111:包含多个由沿x轴和y轴极化方向的阵子组成的双极化阵子,用于接收来波模拟信号;
子阵相移器子模块1112:包含多个用于调节阵子接收信号相位的模拟相移器,根据来波到达角信息进行设置;
子阵信号转换子模块1113:包含用于合并相移器输出的模拟信号合并器,并将合并信号转换到数字基带,这里包括将射频或中频信号下变频到基带;
进一步极化参数估计模块12包括:
波达角估计子模块121:用于对子阵信号处理模块111在两个正交极化方向的输出信号进行处理,进而估计到达角度信息。根据到达角估计信息反馈对子阵相移器子模块1112中的模拟相移器进行设置;
极化参数估计子模块122:用于对子阵信号处理模块111在两个正交极化方向的输出信号进行处理,获得两组可选的极化参数估计值,再根据功率计算子模块124反馈确定期望的极化参数估计值并输出;
最大比合并子模块123:用于通过波达角估计子模块121和极化参数估计子模块122输出计算最大比合并系数,对子阵信号处理模块111输出信号进行最大比合并。根据功率计算子模块124反馈确定期望的最大比合并系数,同时输出期望的最大比合并信号;
功率计算子模块124:用于计算最大比合并子模块123输出的两路信号功率,通过比较功率大小确定期望的信号,同时反馈给极化参数估计子模块122和最大比合并子模块123;
本发明提供的第二实施例是一种极化参数估计方法31,由极化参数估计模块12执行。
为了便于理解本实施例,下面给出一些预备知识。
1)假设采用MxNx×MyNy个双极化天线构成的均匀平面方阵接收一个来波信号其中Mx和My分别表示沿x轴和y轴方向子阵的数量;Nx和Ny分别表示每个子阵中沿x轴和y轴方向阵子的数量。令第m=myMx+mx(m=0,1,......,MxMy-1;mx=0,1,......,Mx-1;my=0,1,......,My-1)个子阵对应的模拟波束成型在c(x或y)轴极化方向输出信号表示为其中Ec表示c轴极化方向的天线响应;Ps(θ,φ)表示子阵的方向图;λ和(θ,φ)分别表示来波信号的波长,天顶角和方位角;和分别表示沿x轴和y轴方向子阵的间距;表示对应子阵输出的加性零均值高斯白噪声,其功率表示为
需要估计的极化参数信息1)极化辅助角γ和2)极化相位差η包含在Ec中,他们的关系可以表示为
Ex=sinγcosθcosφejη-cosγsinφ
Ey=sinγcosθsinφejη+cosγcosφ
本发明还公开了一种极化参数估计方法,其步骤流程如图3所示,包括:
步骤S4:在极化参数估计子模块122中,利用等式(5)进行极化参数估计,得到两组可选估计
步骤S5:在最大比合并子模块123中,利用等式(6)计算最大比合并系数
这里
步骤S6:在功率计算子模块124中,利用等式(7)对最大比合并子模块123的输出信号sk[i](k=1,2)的功率进行滤波,得到第i个时刻滤波后的平均功率
这里0<β<1是滤波系数。实际应用中,可以根据仿真选择最优的滤波系数。利用等式(8)确定期望索引,
这里I1表示迭代样点的要求数量。
需要指出的是,利用混合双极化阵列对到达角信息估计的方法采用现有技术方法,本发明重点在于提出极化参数估计装置及方法,因此具体的到达角信息估计的方法并未列出。
为评价本发明所提出的极化参数的估计方案,对极化参数估计均方误差(MSE)性能进行了计算机仿真,同时还将结果与只利用x轴或者y轴极化方向的接收信号能量进行极化参数估计的性能进行比较。在仿真中,我们采用Mx=4,My=1,Nx=8,Ny=1,d=λ2的混合双极化阵列接收信号。
图4为本发明提供的方案极化参数估计MSE性能对阵子的平均接收信噪比的计算机仿真结果。为了获得统计性能,仿真假设到达角度和极化参数在设定范围内均匀分布,并执行50000次的独立仿真。假设步骤32中精确地获得到达角估计。一次估计采用的样点数量I=I1=8,功率滤波系数β=0.25。从图4可以看出,与只利用x轴或者y轴极化方向的接收信号能量进行极化参数估计的方法比较,本发明方案提出的极化参数估计方法的MSE性能得到明显改善。
上述所提发明及所公开的装置适用于采用双极化天线收发的数据通信***,例如卫星通信。
前述仅描述了本发明的一些实施例,本领域的技术人员可以在不偏离本发明的范围和精神的情况下对其进行修改和改变,所述实施例是说明性的而不是限制性的。如果这些修改和改变属于本发明权利要求及其等同技术范围之内,则本发明也意图包含这些修改和改变在内。
Claims (2)
1.一种极化参数估计的装置,其特征在于,至少包括:
天线信号处理模块:用于将阵列天线接收到的模拟信号波束成型,并转换为数字信号;
极化参数估计模块:用于将天线信号处理模块输出进行信号处理估计来波的极化参数;
其中,天线信号处理模块包括多个子阵信号处理模块,所述子阵信号处理模块包括:
子阵天线子模块:包含多个由沿x轴和y轴极化方向的阵子组成的双极化阵子,用于接收来波模拟信号;
子阵相移器子模块:包含多个用于调节阵子接收信号相位的模拟相移器,其根据来波到达角信息进行设置;
子阵信号转换子模块:包含用于合并相移器输出的模拟信号合并器,并将合并信号转换到数字基带,包括将射频或中频信号下变频到基带;
所述极化参数估计模块包括:
波达角估计子模块:用于对子阵信号处理模块在两个正交极化方向的输出信号进行处理,进而估计到达角度信息;根据到达角估计信息反馈对子阵相移器子模块中的模拟相移器进行设置;
极化参数估计子模块:用于对子阵信号处理模块在两个正交极化方向的输出信号进行处理,获得两组可选的极化参数估计值,再根据功率计算子模块反馈确定期望的极化参数估计值并输出;
最大比合并子模块:用于通过波达角估计子模块和极化参数估计子模块输出计算最大比合并系数,对子阵信号处理模块输出信号进行最大比合并;根据功率计算子模块反馈确定期望的最大比合并系数,同时输出期望的最大比合并信号;
功率计算子模块:用于计算最大比合并子模块输出的两路信号功率,通过比较功率大小确定期望的信号,同时反馈给极化参数估计子模块和最大比合并子模块;
极化参数估计模块执行以下步骤:
步骤S4:在极化参数估计子模块中,利用等式(5)进行极化参数估计,得到两组可选估计
步骤S5:在最大比合并子模块中,利用等式(6)计算最大比合并系数
其中:
步骤S6:在功率计算子模块中,利用等式(7)对最大比合并子模块的输出信号sk[i](k=1,2)的功率进行滤波,得到第i个时刻滤波后的平均功率:
其中,0<β<1是滤波系数;根据仿真选择最优的滤波系数;利用等式(8)确定期望索引,
这里I1表示迭代样点的要求数量;
2.一种极化参数估计的方法,其特征在于,该方法包括以下步骤:
对混合双极化阵列沿着x轴和y轴极化方向上的接收信号进行处理,获得到达角信息估计;
利用估计的到达角信息对相移器进行设置,得到每个子阵沿着x轴和y轴极化方向上的波束成型输出信号;
利用波束成型输出信号估计极化参数,获得两组可选估计;
利用到达角信息估计和可选极化参数估计对波束成型输出信号进行加权合并,从两组可选估计中确定期望极化参数估计和最大比输出信号;
该方法由极化参数估计模块执行,其进一步包括以下步骤:
步骤S4:在极化参数估计子模块中,利用等式(5)进行极化参数估计,得到两组可选估计
步骤S5:在最大比合并子模块中,利用等式(6)计算最大比合并系数
其中:
步骤S6:在功率计算子模块中,利用等式(7)对最大比合并子模块的输出信号sk[i](k=1,2)的功率进行滤波,得到第i个时刻滤波后的平均功率:
其中,0<β<1是滤波系数;根据仿真选择最优的滤波系数;利用等式(8)确定期望索引,
这里I1表示迭代样点的要求数量;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110061089.2A CN112887040B (zh) | 2021-01-18 | 2021-01-18 | 一种极化参数估计装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110061089.2A CN112887040B (zh) | 2021-01-18 | 2021-01-18 | 一种极化参数估计装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112887040A CN112887040A (zh) | 2021-06-01 |
CN112887040B true CN112887040B (zh) | 2022-10-25 |
Family
ID=76048806
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110061089.2A Active CN112887040B (zh) | 2021-01-18 | 2021-01-18 | 一种极化参数估计装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112887040B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110824414A (zh) * | 2019-10-14 | 2020-02-21 | 杭州电子科技大学 | 一种到达角估计的装置及方法 |
CN110837075A (zh) * | 2019-11-13 | 2020-02-25 | 杭州电子科技大学 | 一种低复杂度的极化参数估计跟踪装置及方法 |
CN111983556A (zh) * | 2020-08-19 | 2020-11-24 | 杭州电子科技大学 | 一种到达角估计的装置及方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8754810B2 (en) * | 2009-02-02 | 2014-06-17 | Commonwealth Scientific And Industrial Research Organisation | Hybrid adaptive antenna array |
CN108337031A (zh) * | 2017-12-22 | 2018-07-27 | 南京理工大学 | 大规模mimo中基于根值music的混合到达角估计 |
-
2021
- 2021-01-18 CN CN202110061089.2A patent/CN112887040B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110824414A (zh) * | 2019-10-14 | 2020-02-21 | 杭州电子科技大学 | 一种到达角估计的装置及方法 |
CN110837075A (zh) * | 2019-11-13 | 2020-02-25 | 杭州电子科技大学 | 一种低复杂度的极化参数估计跟踪装置及方法 |
CN111983556A (zh) * | 2020-08-19 | 2020-11-24 | 杭州电子科技大学 | 一种到达角估计的装置及方法 |
Non-Patent Citations (1)
Title |
---|
Adaptive AoA and Polarization Estimation for Receiving Polarized mmWave Signals;Hang Li;《IEEE》;20190612;第2-3节 * |
Also Published As
Publication number | Publication date |
---|---|
CN112887040A (zh) | 2021-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108462521B (zh) | 自适应阵列天线的抗干扰实现方法 | |
CN101369014B (zh) | 应用于多输入多输出雷达的双边约束自适应波束形成方法 | |
CN105891771B (zh) | 一种提高估计精度的基于连续分布的角度估计方法与设备 | |
CN109787671B (zh) | 一种混合波束成形装置及方法 | |
CN104991236B (zh) | 一种单基地mimo雷达非圆信号相干源波达方向估计方法 | |
CN110515033B (zh) | 一种基于Toeplitz矩阵恢复的欠通道测向***及方法 | |
CN112737647A (zh) | 一种基于混合连接的毫米波大规模mimo预编码方法 | |
WO2002052677A1 (en) | Base station, base station module and method for direction of arrival estimation | |
CN110824414A (zh) | 一种到达角估计的装置及方法 | |
CN107276649B (zh) | 低精度模数转换与混合预编码结合的无反馈波束训练方法 | |
CN110837075A (zh) | 一种低复杂度的极化参数估计跟踪装置及方法 | |
CN109861933B (zh) | 一种基于music算法和预编码的毫米波mimo信道估计方法 | |
CN111983556B (zh) | 一种到达角估计的装置及方法 | |
CN103037494B (zh) | 多波束天线功率控制方法 | |
CN111859644A (zh) | 一种圆形相控阵列天线的接收波束形成与方位扫描方法 | |
CN104868946B (zh) | 自适应加权的子阵级混合mimo‑相控阵***的干扰抑制方法 | |
CN113067616B (zh) | 基于时间调制数字超表面的dbf相控阵*** | |
CN113644941A (zh) | 一种基于大规模mimo接收阵列结构的模糊相位快速消除方法 | |
Li et al. | Adaptive AoA and polarization estimation for receiving polarized mmWave signals | |
CN112887040B (zh) | 一种极化参数估计装置及方法 | |
Okkonen | Uniform linear adaptive antenna array beamforming implementation with a wireless open-access research platform | |
CN113917388A (zh) | 一种使用混合adc的doa估计和定位方法 | |
CN109738857B (zh) | 基于交叉偶极子阵列的非圆信号定位快速估计方法及装置 | |
CN115189722B (zh) | 适于大规模混合模数mimo的最大化扇区接收功率的doa估计方法 | |
CN112526510A (zh) | 基于方向图分集的单通道角度超分辨方法及*** |
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 |