[发明专利]一种抗阵列系统误差的稳健波束形成方法在审

专利信息
申请号: 201510664765.X 申请日: 2015-10-15
公开(公告)号: CN105306123A 公开(公告)日: 2016-02-03
发明(设计)人: 李文兴;毛晓军;李思;毛云龙;李迎松 申请(专利权)人: 哈尔滨工程大学
主分类号: H04B7/06 分类号: H04B7/06
代理公司: 暂无信息 代理人: 暂无信息
地址: 150001 黑龙江省哈尔滨市南岗区*** 国省代码: 黑龙江;23
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明公开了一种抗阵列系统误差的稳健波束形成方法。首先对阵列天线接收的采样数据协方差矩阵特征分解,根据期望信号导向矢量与期望信号对应特征向量的相关性,估计出期望信号对应的特征向量,进而去掉协方差矩阵中期望信号对应的特征向量,将剩下的特征向量张成的子空间,与Capon谱估计的干扰加噪声协方差矩阵加权组合,作为本发明波束形成方法中的干扰加噪声的协方差矩阵;再利用二次约束二次优化算法得到校正的期望信号导向矢量。本发明能有效避免阵列系统误差引起的性能降低,极大的提高了波束形成器的稳健性,与目前的多数稳健波束形成算法相比,能够更加接近理论输出信干噪比,有更好的输出性能。
搜索关键词: 一种 阵列 系统误差 稳健 波束 形成 方法
【主权项】:
一种抗阵列系统误差的稳健波束形成方法,其特征在于:包括以下步骤,步骤一:构建阵列天线接收数据的采样协方差矩阵:<mrow><mover><mi>R</mi><mo>^</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>K</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><mi>X</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msup><mi>X</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>;</mo></mrow>步骤二:对协方差矩阵进行特征分解:<mrow><mover><mi>R</mi><mo>^</mo></mover><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>&mu;</mi><mi>i</mi></msub><msub><mover><mi>&upsi;</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><msubsup><mover><mi>&upsi;</mi><mo>&OverBar;</mo></mover><mi>i</mi><mi>H</mi></msubsup><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>N</mi></mrow>其中,μi为对协方差矩阵进行特征分解所得到的第i个特征值,为μi所对应的特征向量;步骤三:将估计的期望信号导向矢量向每个特征向量投影,得到:<mrow><mi>p</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mo>|</mo><msubsup><mover><mi>&upsi;</mi><mo>&OverBar;</mo></mover><mi>i</mi><mi>H</mi></msubsup><msub><mover><mi>a</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><msup><mo>|</mo><mn>2</mn></msup><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>N</mi><mo>,</mo></mrow>将p(i)从大到小依次排列,设p(i)的最大值max(p(i))所对应的特征向量为对应的特征值为去除,剩余的特征向量为其分别对应的特征值为得到不包含期望信号分量的协方差矩阵:步骤四:利用干扰区域的Capon空间谱重构干扰加噪声协方差矩阵:<mrow><msub><mover><mi>R</mi><mo>~</mo></mover><mrow><mi>i</mi><mo>+</mo><mi>n</mi></mrow></msub><mo>=</mo><msub><mo>&Integral;</mo><mover><mi>&Theta;</mi><mo>&OverBar;</mo></mover></msub><mfrac><mrow><msup><mi>a</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>a</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow><mrow><msup><mi>a</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><msup><mover><mi>R</mi><mo>^</mo></mover><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>a</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>其中,a(θ)表示θ方向的导向矢量;步骤五:得到干扰加噪声协方差矩阵的估计值:<mrow><msub><mi>R</mi><mrow><mi>i</mi><mo>+</mo><mi>n</mi></mrow></msub><mo>=</mo><mi>&alpha;</mi><msub><mover><mi>R</mi><mo>~</mo></mover><mrow><mi>i</mi><mo>+</mo><mi>n</mi></mrow></msub><mo>+</mo><mi>&beta;</mi><msub><mover><mi>R</mi><mo>&OverBar;</mo></mover><mrow><mi>i</mi><mo>+</mo><mi>n</mi></mrow></msub></mrow>其中α和β为加权因子,<mrow><mi>&alpha;</mi><mo>=</mo><mover><mi>&mu;</mi><mo>^</mo></mover><mo>/</mo><mrow><mo>(</mo><mover><mi>&mu;</mi><mo>^</mo></mover><mo>+</mo><msub><mover><mi>&mu;</mi><mo>^</mo></mover><mn>1</mn></msub><mo>+</mo><msub><mover><mi>&mu;</mi><mo>^</mo></mover><mn>2</mn></msub><mo>+</mo><mo>...</mo><mo>+</mo><msub><mover><mi>&mu;</mi><mo>^</mo></mover><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo><mi>&beta;</mi><mo>=</mo><mo>(</mo><mrow><msub><mover><mi>&mu;</mi><mo>^</mo></mover><mn>1</mn></msub><mo>+</mo><msub><mover><mi>&mu;</mi><mo>^</mo></mover><mn>2</mn></msub><mo>+</mo><mo>...</mo><mo>+</mo><msub><mover><mi>&mu;</mi><mo>^</mo></mover><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow><mo>)</mo><mo>/</mo><mi>N</mi><mo>;</mo></mrow>步骤六:构建二次约束二次优化问题对估计的期望信号导向矢量进行校正:<mfenced open = '' close = ''><mtable><mtr><mtd><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mrow><mi>e</mi><mo>&perp;</mo></mrow></munder></mtd><mtd><mrow><msup><mrow><mo>(</mo><msub><mover><mi>a</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mo>+</mo><msub><mi>e</mi><mo>&perp;</mo></msub><mo>)</mo></mrow><mi>H</mi></msup><msubsup><mi>R</mi><mrow><mi>i</mi><mo>+</mo><mi>n</mi></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><msub><mover><mi>a</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mo>+</mo><msub><mi>e</mi><mo>&perp;</mo></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mfenced open = '' close = ''><mtable><mtr><mtd><mrow><mi>s</mi><mi>u</mi><mi>b</mi><mi>j</mi><mi>e</mi><mi>c</mi><mi>t</mi><mi> </mi><mi>t</mi><mi>o</mi></mrow></mtd><mtd><mrow><mo>|</mo><msubsup><mover><mi>a</mi><mo>&OverBar;</mo></mover><mn>0</mn><mi>H</mi></msubsup><msub><mi>e</mi><mo>&perp;</mo></msub><mo>|</mo><mo>=</mo><mn>0</mn><mo>,</mo></mrow></mtd></mtr></mtable></mfenced><mrow><msup><mrow><mo>(</mo><msub><mover><mi>a</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mo>+</mo><msub><mi>e</mi><mo>&perp;</mo></msub><mo>)</mo></mrow><mi>H</mi></msup><msub><mi>R</mi><mrow><mi>i</mi><mo>+</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><msub><mover><mi>a</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mo>+</mo><msub><mi>e</mi><mo>&perp;</mo></msub><mo>)</mo></mrow><mo>&le;</mo><msubsup><mover><mi>a</mi><mo>&OverBar;</mo></mover><mn>0</mn><mi>H</mi></msubsup><msub><mi>R</mi><mrow><mi>i</mi><mo>+</mo><mi>n</mi></mrow></msub><msub><mover><mi>a</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub></mrow>得到校正后的期望信号导向矢量估计<mrow><msub><mover><mi>a</mi><mo>^</mo></mover><mn>0</mn></msub><mo>=</mo><msub><mover><mi>a</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mo>+</mo><msub><mi>e</mi><mo>&perp;</mo></msub><mo>;</mo></mrow>步骤七:得到自适应加权向量:利用w进行波束形成,得到天线阵列输出数据。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于哈尔滨工程大学,未经哈尔滨工程大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201510664765.X/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top