CN101552630A - 一种基于航空通信信道下的波束成形方法 - Google Patents
一种基于航空通信信道下的波束成形方法 Download PDFInfo
- Publication number
- CN101552630A CN101552630A CN 200910084480 CN200910084480A CN101552630A CN 101552630 A CN101552630 A CN 101552630A CN 200910084480 CN200910084480 CN 200910084480 CN 200910084480 A CN200910084480 A CN 200910084480A CN 101552630 A CN101552630 A CN 101552630A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- msup
- mtr
- 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
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000004891 communication Methods 0.000 title claims abstract description 19
- 239000013598 vector Substances 0.000 claims abstract description 36
- 239000011159 matrix material Substances 0.000 claims abstract description 23
- 238000005070 sampling Methods 0.000 claims abstract description 12
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000005562 fading Methods 0.000 abstract description 11
- 238000005516 engineering process Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 241000209094 Oryza Species 0.000 description 2
- 235000007164 Oryza sativa Nutrition 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 235000009566 rice Nutrition 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001268 conjugating effect Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Landscapes
- Radio Transmission System (AREA)
Abstract
本发明涉及一种基于航空通信信道下的波束成形方法,属于信号处理和通信领域。该方法首先估计基带信号的自相关矩阵的逆矩阵,并计算出协方差向量。然后,将协方差向量与自相关矩阵的逆矩阵相乘,得到最小均方误差的权值向量;同时,求出贝叶斯权值向量。最后,将最小均方误差的权值向量与贝叶斯权值向量相加,即可得出空域滤波器所需要的权值。该方法能够在多径衰落以及强干扰环境中,充分利用后验概率信息,在采样同步偏差存在的情况下形成最佳空域波形,并输出接近最佳的信扰比。
Description
技术领域
本发明涉及一种数字波束成形方法,特别涉及一种基于航空通信信道的数字波束成形方法,属于信号处理和通信领域。
背景技术
在无线通信中,由于收信机和发信机之间的多条传播路径造成的多径效应是一个普遍存在的问题。通常,航空信道被认为是由一个视距传播(LOS)成分和一个多径散射成分组成,称为莱斯(Rice)衰落信道。其中多径散射成分是由传播环境中物体对电磁波的反射和散射产生的。在航空通信中,飞行器不同的状态会导致不同的信道衰落形式,比如飞行器在飞行途中、跑道滑行以及停场状态导致的信道衰落是不同的。一般来说,在停场状态下,由于停机坪、机库的影响,视距成分是不存在的,这导致了瑞利(Rayleigh)衰落,是最严重的衰落形式。而在其他状态下,接收信号通常包含LOS成分和散射成分,这会导致莱斯衰落。同时在航空通信中还存在不同用户对期望信号的干扰,在战争环境下,还存在来自敌方的同频强干扰。因此应用于这种环境下的通信系统要同时具有抗多径衰落以及抗干扰的性能。
Rake接收方法是目前广泛采用的一种抗衰落方法,它采用分集合并技术,对多个信号分别进行处理、合并从而获得更大的益处。但是由于航空信道中可能存在较强的同频干扰,采用Rake接收并不能很好地消除干扰。采用扩频,跳频技术可以使通信系统具有很强的抗干扰能力,但是该技术本身并不能消除多径衰落的影响。当期望信号存在多个非相关多径时,波束成形技术能够保留其中的一路信号,而把零陷对准其它信号,这样,阵列能够减小由非相关多径带来的干扰,也可以消除不同到达角度上(DOA)的干扰信号。
传统的波束成形的实现方法是由多个辐射元排阵构成天线阵列,用波束形成网络向阵列单元激励所需的振幅和相位,以形成不同形状的成形波束。它的优点是可对波束数目和形状进行灵活控制,并可控制波束作快速扫描;但结构较复杂,造价高。已知的另一种做法是利用数字芯片对采样和下变频得到的数字基带信号作一系列的复数矩阵运算,这种做法所需的计算量较大,且处理过程带来一定的误差,使得处理结果难以达到空域滤波器的误差指标。同时,应用于通信系统中的波束成形系统不同于传统雷达中所采用波束天线阵列,由于雷达中回波的存在,可以方便估计期望信号的DOA,虽然在通信系统中可采用已有的多重信号分类(MUSIC)方法和旋转不变信号参数估计(ESPRIT)方法估计通信系统中期望信号的DOA,但是实现复杂度高,且误差较大。已知的最小均方误差(MMSE)波束成形方法可以避免DOA估计,降低实现的复杂度,但是传统的MMSE方法对接收端采样同步偏差较为敏感,该偏差会导致阵列输出信扰噪声比(SINR)下降,该方法在采样同步偏差为10的情况下的输出SINR在图4中给出。
发明内容
本发明的目的在于提供一种基于航空通信信道下的波束成形方法,该方法能够有效解决航空通信信道中的多径衰落以及干扰等问题,降低传统方法的实现成本和复杂度,并且能够克服传统数字波束成形方法的DOA估计不准以及采样同步偏差问题。
本发明所采用的技术方案如下:
一种基于航空通信信道下的波束成形方法,一种基于航空通信信道下的波束成形方法,其实现步骤包括自相关矩阵处理、协方差向量处理、权值计算输出,其流程图如图1所示,包括以下步骤:
步骤一:利用长度为N的训练序列xn,其中xn为M×1维的列向量,M为天线阵列的单元个数,估计基带信号的自相关矩阵和协方差向量
估计公式为:
步骤三:求取贝叶斯(Bayesian)权值向量w2,方法如下:
假设实际信号的到来方向u服从先验密度函数q(u)的随机变量,这里考虑q(u)是由L个离散点u1,u2,…,uL所组成的先验参数空间,10≤L≤30。在干扰存在的情况下,
其中,
其中,L为先验参数空间的维数。
③根据②中的a求Bayesian权值向量w2为:
步骤四:将w1与贝叶斯w2相加,即可得出空域滤波器所需要的权值w:
w=w1+w2 (7)
至此,就完成了波束的成形。
有益效果
本发明方法对比现有技术,能够在采样同步偏差存在的情况下形成期望的空间响应波形,输出的SINR接近于最佳情况,降低了传统方法的实现成本和复杂度
附图说明
图1为本发明的波束成形方法的流程图;
图2为理论上最佳权值所产生的波束特性曲线示意图;
图3为实现的波束成形器所产生的波束特性曲线示意图;
图4为输出的信扰比随训练序列变化特性曲线示意图。
具体实施方式
下面结合附图及实施例对本发明的优选实施方式做进一步详细说明。
①计算出各个采样时刻向量的积,即Rn=xnxn H。
该过程仅用到复数乘法运算,由于Rn满足Rn H=Rn,所以在计算时可以只计算该矩阵的主对角线以及主对角线以上的矩阵元素,这样能够节省M(M-1)/2个复数乘法器以及复数加法器。
②将①中计算的Rn由
计算自相关矩阵该过程用到复数求和运算,由于N为2的整数幂,上式的除法运算可通过将求和结果右移log2 N位来实现,该矩阵中的其他元素可以通过对上三角矩阵中的相应元素求共轭得到,即rij=rji *
(8)
其中hm (l)为行向量nl中的元素。
步骤四:由(7)计算最终的权值w。
至此,就完成了波束的成形。
实施例
设定接收端从天线经多路接收机下变频的中频信号,接收端的天线阵列单元数量为4,中频频率为70MHz,在AD中经过64MHz的带通采样得到6MHz的信号,送入FPGA中进行处理。FPGA系统主时钟为32MHz,经下变频到基带输出I,Q两路数据到波束成形模块,I,Q两路均为16bit量化的数据信息。现利用本发明方法对I,Q两路为16bit量化的数据信息进行波束成形处理。
利用长度为64的训练序列,计算出协方差向量
w1=[-0.0054-j0.0016 0.0043+j0.0004 0.004-j0.0038 0.0077-j0.0019]T
步骤三:求取贝叶斯权值向量w2,方法如下:
a(u0)=[1 0.9511+j0.3090 0.8090+j0.5878 0.5878+j0.8090]T
a(u1)=[1 0.9686+j0.2487 0.8763+j0.4818 0.7290+j0.6845]T
a(u2)=[1 0.9823+j0.1874 0.9298+j0.3681 0.8443+j0.5358]T
a(u3)=[1 0.9921+j0.1253 0.9686+j0.2487 0.9298+j0.3681]T
a(u4)=[1 0.9980+j0.0628 0.9921+j0.1253 0.9823+j0.1874]T
a(u5)=[1 1 1 1]T
a(u6)=[1 0.9980-j0.0628 0.9921-j0.1253 0.9823-j0.1874]T
a(u7)=[1 0.9921-j0.1253 0.9686-j0.2487 0.9298-j0.3681]T
a(u8)=[1 0.9823-j0.1874 0.9298-j0.3681 0.8443-j0.5358]T
a(u9)=[1 0.9686-j0.2487 0.8763-j0.4818 0.7290-j0.6845]T
a(u10)=[1 0.9511-j0.3090 0.8090-j0.5878 0.5878-j0.8090]T
根据求出a:
则Bayesian权值向量w2为:
最后,将MMSE权值向量w1与Bayesian权值向量w2相加,即可得出w:
w=w1+w2
=[0.3497-j0.0128 0.7214+j0.0268 0.7127-j0.0282
0.3452+j0.0134]T
至此,就完成了波束的成形。
将计算出的权值w输入matlab中计算波束的方向特性图,图2显示了理论上的最佳权值所产生的波束特性曲线,而图3示出了所设计的波束成形器所计算的权值所产生的波束特性曲线。图4示出了所设计的波束成形器所输出信扰比随训练序列长度变化的曲线,并对比传统的MMSE方法在采样同步误差为10个采样点的情况下的输出SINR,可见在本实施例中,本发明的方法可以带来5dB以上的好处。
Claims (1)
1、一种基于航空通信信道下的波束成形方法,其特征在于包括以下步骤:
估计公式为:
步骤三:求取贝叶斯权值向量w2,方法如下:
假设实际信号的到来方向u服从先验密度函数q(u)的随机变量,这里考虑q(u)是由L个离散点u1,u2,…,uL所组成的先验参数空间,10≤L≤30;在干扰存在的情况下,
其中,
(·)H表示矩阵或向量的共轭转置运算,表示天线之间的相位偏移;
其中,L为先验参数空间的维数;
③根据②中的a求贝叶斯(Bayesian)权值向量w2为:
步骤四:将最小均方误差的权值向量w1与贝叶斯权值向量w2相加,即可得出空域滤波器所需要的权值w:
w=w1+w2
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910084480 CN101552630B (zh) | 2009-05-19 | 2009-05-19 | 一种基于航空通信信道下的波束成形方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910084480 CN101552630B (zh) | 2009-05-19 | 2009-05-19 | 一种基于航空通信信道下的波束成形方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101552630A true CN101552630A (zh) | 2009-10-07 |
CN101552630B CN101552630B (zh) | 2013-01-16 |
Family
ID=41156648
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200910084480 Expired - Fee Related CN101552630B (zh) | 2009-05-19 | 2009-05-19 | 一种基于航空通信信道下的波束成形方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101552630B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102497644A (zh) * | 2011-12-22 | 2012-06-13 | 电子科技大学 | 一种低复杂度正交迭代波束成形的方法 |
CN109039405A (zh) * | 2018-07-04 | 2018-12-18 | 西安电子科技大学 | 一种基于干扰子空间旋转的跳空通信方法 |
CN111726157A (zh) * | 2020-06-17 | 2020-09-29 | 中山大学 | 面向航空通信的自适应波束赋形方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100399721C (zh) * | 2005-01-11 | 2008-07-02 | 中国人民解放军理工大学 | 基于发送辅助选择用户反馈的正交预波束成形传输方法 |
-
2009
- 2009-05-19 CN CN 200910084480 patent/CN101552630B/zh not_active Expired - Fee Related
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102497644A (zh) * | 2011-12-22 | 2012-06-13 | 电子科技大学 | 一种低复杂度正交迭代波束成形的方法 |
CN102497644B (zh) * | 2011-12-22 | 2014-08-13 | 电子科技大学 | 一种低复杂度正交迭代波束成形的方法 |
CN109039405A (zh) * | 2018-07-04 | 2018-12-18 | 西安电子科技大学 | 一种基于干扰子空间旋转的跳空通信方法 |
CN109039405B (zh) * | 2018-07-04 | 2021-03-23 | 西安电子科技大学 | 一种基于干扰子空间旋转的跳空通信方法 |
CN111726157A (zh) * | 2020-06-17 | 2020-09-29 | 中山大学 | 面向航空通信的自适应波束赋形方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101552630B (zh) | 2013-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101369014B (zh) | 应用于多输入多输出雷达的双边约束自适应波束形成方法 | |
JP6278961B2 (ja) | レーダ装置 | |
CN110850445B (zh) | 一种基于空时采样协方差求逆的脉冲干扰抑制方法 | |
CN101881822B (zh) | 一种针对共享谱雷达同频干扰的抑制方法 | |
KR101555996B1 (ko) | 저피탐 표적 탐지를 위한 주파수 다이버시티를 이용한 공분산 행렬 합성 장치 및 방법 | |
CN106646394A (zh) | 一种fda‑mimo雷达中基于循环优化的发射波束矩阵设计方法 | |
JP5073809B2 (ja) | チャネル情報予測システム及びチャネル情報予測方法 | |
CN105319545B (zh) | 提高stap检测性能的mimo-ofdm雷达波形设计方法 | |
Fascista et al. | A Pseudo Maximum likelihood approach to position estimation in dynamic multipath environments | |
CN103760519A (zh) | 高分辨率doa估计方法及系统 | |
CN114095318B (zh) | 智能超表面辅助的混合构型毫米波通信系统信道估计方法 | |
JP4481336B2 (ja) | チャネル情報予測システム及びチャネル情報予測方法 | |
CN101572574B (zh) | 基于最小二乘—最小均方的智能天线自适应干扰抑制方法 | |
CN103777197A (zh) | 单基地mimo雷达中降维传播算子的方位估计方法 | |
CN110361760A (zh) | 基于子空间跟踪的gnss接收机多波束指向抗干扰方法 | |
CN111880198B (zh) | 基于交替极化敏感阵列的空时极化抗干扰方法 | |
CN104345300A (zh) | 杂波空时谱线性补偿的机载非正侧视阵雷达stap方法 | |
CN110824414A (zh) | 一种到达角估计的装置及方法 | |
Kumari et al. | Low resolution sampling for joint millimeter-wave MIMO communication-radar | |
CN105487054A (zh) | 提高基于mimo-ofdm雷达stap最差检测性能的稳健波形设计方法 | |
Qian et al. | Joint design for co-existence of MIMO radar and MIMO communication system | |
CN114114188A (zh) | 一种低旁瓣的fda雷达通信一体化波形设计方法 | |
CN101552630A (zh) | 一种基于航空通信信道下的波束成形方法 | |
CN103701515A (zh) | 一种数字多波束形成方法 | |
CN104346532A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130116 Termination date: 20140519 |