空间谱专题11:子阵平滑与秩亏缺

作者:桂。

时间:2017-09-29  21:20:18

链接:http://www.cnblogs.com/xingshansi/p/7612984.html 


前言

主要分析在解相干算法中,子阵平滑的有效性。

一、前向平滑

以均匀线阵(ULA)为例,第l个接收阵元的信号为:

其中,M为阵元数,N为信号个数。

以前向平滑为例:

令β = ,对于第k个子阵有

其中D为:

假设信号与噪声不相关,且噪声为白噪声,计算相关矩阵:

前向平滑修正的协方差矩阵:

证明:当子阵阵元数m>=N,且当p>=N时,Rf为满秩矩阵。

对于相干情况,Rs的秩为1,故可以用矢量相乘的形式表示:

显然C与Rsf的秩是一致的。式中

C可进一步表示为

其中,α为对应信号的幅度,因此该矩阵rank = K

AL是Vandermonde矩阵:

从而有:

进一步得出结论:

得证。

 上面的推导是教科书中的推导,以两个子阵、两个信号为例,假设两个子阵的间距为D,对应的A可写为:

虽然具有Vandermonde结构,但即使入射角度不同,上面的矩阵仍然可能出现rank = 1的情况,只要满足(k为整数):

任意取D = 4*λ,theta1、2分别取-60°~60°,仿真结果:

D = 2*λ,进行仿真:

当两个相位差2pi整数倍,如图蓝线所示,此时相关矩阵的秩仍然是亏缺的,此时谱估计解相干失效,这里称蓝色区域为盲区。以MUSIC算法为例,自然也无从谈起矩阵分解以及谱估计。另一方面,D与λ比值的取值越小,则盲区越小

这个盲区类似:子阵间距(孔径)减小,则相位模糊较少(盲区减少),但分辨率下降;子阵间距增大,则相位模糊增多(盲区增加),对应分辨率提高。

 

二、后向平滑算法

后向平滑:

子阵信号为:

其中J为交换矩阵。

以第p-k+1个子阵为例:

后向空间平滑修正的相关矩阵为:

得出的性质与前向平滑类似。当子阵阵元数m>=N,且当p>=N时,Rb为满秩矩阵。

 

三、前后向平滑算法

前后向平滑,即是将前向平滑、后向平滑的相关矩阵累加。

当子阵阵元数m>=N,且当2p>=N时,R为满秩矩阵。这一点仅当子阵结构为中心对称结构才成立

 

四、前向、后向、前后向平滑与阵列的关系

前向平滑可以理解成原始子阵的平滑:

后向平滑则可以理解为镜像对称的虚拟(子阵):

因为到原点的距离由d变为-d,在形式上体现的就是相位相反,也就是接收数据取共轭。顺序就是翻转,因此后向平滑就是前向平滑的共轭+翻转

假设两个子阵的原始子阵为1、 2,翻转后为1‘,2’,前向平滑就是利用1,2;后向平滑为1‘,2’。前后向平滑为1,1‘,2,2’,由于形式上利用了更多的子阵,因此前后向平滑需要的子阵个数更少。但子阵必须中心对称,即镜像之后与原始子阵结构一致

 

五、前后向与非均匀的关系

镜像的思路(即后向平滑),也给了非均匀布阵一些启发。

在构造子阵时,如果是一般的子阵,进行copy之后,即可利用空间平滑的思想:

但可以借助镜像特性,进一步节省资源,例如一个5阵元的非均匀线阵:

利用镜像的关系,可以得出两个4元的非均匀线阵(红框所示),这就给布阵拓展了思路。

 

六、非均匀镜像中坐标点的影响

该章节分析,从理论到结论,均错误。

根据上文第一部分的分析,子阵的间距D会造成盲区,且D/λ越大,盲区越大

对于直接copy的空间平滑子阵思路,一旦子阵确定,子阵间距D也随之确定,因此盲区变得不可控

而对于镜像的思路,子阵间距D是可控的,这就使得盲区变得相对可控。r对于上面提到的非均匀线阵,取不同的坐标原点:

子阵之间的距离也随之改变,对于两个子阵,如果借助镜像的思想,理论上间距D可以取的范围是:0 - +∞。

这里就难免带来一个问题:改变坐标原点,即对应改变导向矢量,近而改变原始的阵列接收数据,可对于阵列来讲,接收的数据是固定不变的。这个可以换一个角度理解:X = AS,改变参考点,调节的是A,在X给定不变的情况下,S是一个相对的信息,尽管改变参考系可以使得盲区减小,但S可能得到的是小信号。一个解决办法是在谱估计(MUSIC为例)时取不同的参考系,假设X给定,从何借助两组数据进行DOA估计(镜像阵列,对于不同参考系,盲区不同。)

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/242746.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

空间谱专题12:二维测向的基本方法

作者:桂。 时间:2017-10-01 16:50:08 链接:http://www.cnblogs.com/xingshansi/p/7617442.html 前言 圆阵、面阵、L阵、十字阵、V形阵,都是直接地二维求解,复杂度较高,这里梳理复杂度较低的求解算法。 一…

使用有限状态机(FSM)解释shell 命令

一、有限状态机(Finite State Machine,FSM)是表示有限个状态及在这些状态之间的转移和动作等行为的数学模型,在计算机领域有着广泛的应用。FSM一个常见的应用就是用来负责Parser复杂的数据结构,比如解释URI协议&#x…

搭建JAVAEE环境

搭建JAVAEE环境 ./表示在当前目录下查找文件所以需要配置环境变量 在Linux下下载 或者在window下传输

小数延迟滤波器

作者:桂。 时间:2017-10-10 22:38:46 链接:http://www.cnblogs.com/xingshansi/p/7648274.html 前言 阵列信号处理中,经常用到小数延迟(fractional delay,FD)的思路,例如Beamformi…

空间谱专题13:联合解算DOA(ML/AP)

其中作者:桂。 时间:2017-10-16 07:51:40 链接:http://www.cnblogs.com/xingshansi/p/7675380.html 前言 主要记录二维测向中,分别利用两个一维阵联合解算的思路。 一、AP算法思想 信号模型: 对应相关矩阵 假设噪声为…

机器学习之正则化

根据上一篇博客《统计学习概论》可以知道,正则化的作用是选择经验风险和模型复杂度同时较小的模型。下面从过拟合的角度来理解正则化。 #过拟合问题 例子说明,线性回归问题(房价) 分析: 1)左边第一幅图&am…

python开发平台Ubuntu

python开发平台Ubuntu APT软件管理和远程登录

投影矩阵

作者:桂。 时间:2017-10-19 06:02:00 链接:http://www.cnblogs.com/xingshansi/p/7690292.html 前言 最近在交替投影算法中,用到投影矩阵,简单记录。 一、投影矩阵定义 此处以列满秩为例,行满秩可依次类推…

不同维度极值点查找

作者:桂。 时间:2017-10-19 17:00:12 链接:http://www.cnblogs.com/xingshansi/p/7693557.html 前言 主要梳理不同维度信号极值点的查找思路。 一、思想 思想都是一个,通过极值点的定义:数值高于相邻的点&#xff0c…

俯仰角/偏航角的转化

作者:桂。 时间:2017-10-20 10:29:52 链接:http://www.cnblogs.com/xingshansi/p/7698237.html 前言 主要记录坐标系的转化。 一、坐标转化 对于坐标系: 设方位角为φ,俯仰角为theta,仰角β与偏航角α&am…

安装Centos8.1

安装Centos8.1 按回车

逾期后,如何修复个人征信?

个人征信大家都是知道很宝贵,但是有些朋友会在有意无意之间造成逾期,结果给申请贷款、申请信用卡带来了诸多不便。逾期是谁也不想看到的,但是如果逾期已经发生了,我们就无法改变,唯一能做的就是努力去修复,…

空间谱专题16:信号个数估计

作者:桂。 时间:2017-10-24 21:50:16 链接:http://www.cnblogs.com/xingshansi/p/7726082.html 前言 记录阵列信号在DOA估计中,信源个数估计的基本方法。 一、基本估计方法 参考:王永良《空间谱估计》p42: 以MDL为例&#xff1a…

根据verilog代码画电路图

根据verilog代码画电路图 FPGA设计的本质是硬件设计,而且verilog是描述硬件设计的语言(也就是描述电路),一个标准的工程师需要学会建立电路和Verilog对应的关系,学会看到电路图,就能写出相应的Verilog代码…

复数矩阵分解的拆解思路(矩阵求逆/特征值分解)

作者:桂。 时间:2017-10-26 07:11:02 链接:http://www.cnblogs.com/xingshansi/p/7735016.html 前言 主要记录特征值分解的硬件实现思路。 一、实数矩阵转化 在FPGA运算中,对实数运算通常优于对复数运算。假设C为复数矩阵&#…

贷款机构如何审核个人征信?

贷款是要查看征信的, 但是大家知道贷款机构是如何审查借款人的征信吗?下面我们一起来看下。1看征信报告打印时间一般贷款机构要求提供的是最新的征信报告,当然每个机构要的具体期限不一样,有的要求是最近一个月的,有的…

COMS技术

COMS技术 n型MOS(NMOS)三极管的结构如图所示,该图不是按照实际比例绘制的。三极管的衬底是被掺杂后成为p型半导体材料的硅晶片。NMOS衬底的厚度远比其他三极管要厚。在每个三极管上有两个区域,被掺入大量的杂质,成为n型半导体区域。这两个区域分别形成了三极管的源极和漏极…

VS2015 + CUDA 8.0 配置GTX1070的OpenCL 开发环境

一、查看计算机对OpenCL异构计算的支持情况 使用 GPU Caps Viewer 查看计算机对OpenCL的支持情况,目前最新的版本是 gpu-caps-viewer-1-32-0, 下载地址:http://www.geeks3d.com/20161107/gpu-caps-viewer-1-32-0-released/ 从上面两个图可以…