空间谱专题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…

从今天起,逃离迷茫与枯燥,做一个与快乐常伴的人

-0-又是一个下雨天,郑州连着下了一个星期的雨了。国庆长假,街道上褪去了往日的繁华,空荡荡的,连几家最爱吃的商铺也关门回家团圆了。今天中秋节,你在哪里?又在做着什么呢?-1-一个人猫在家里&…

搭建JAVAEE环境

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

小数延迟滤波器

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

387. First Unique Character in a String QuestionEditorial Solution

Given a string, find the first non-repeating character in it and return it’s index. If it doesn’t exist, return -1. Examples: s “leetcode” return 0. s “loveleetcode”, return 2. Note: You may assume the string contain only lowercase letters. …

AD提高动态的方法(附SNR计算)

作者:桂。 时间:2017-10-10 23:36:40 链接:http://www.cnblogs.com/xingshansi/p/7648475.html 前言 主要分析AD的基本特性,以及动态的计算方法。并总结量化位数确定的前提下,提高信号精度的一般性方法。 未完待续..…

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

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

现代C++函数式编程

链接 http://geek.csdn.net/news/detail/96636 概述 函数式编程是一种编程范式,它有下面的一些特征: 函数是一等公民,可以像数据一样传来传去。 高阶函数 递归 pipeline 惰性求值 柯里化 偏应用函数 C98/03中的函数对象&#xff0c…

机器学习之正则化

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

python开发平台Ubuntu

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

a.out、coff、elf三种文件格式

转自:http://blog.chinaunix.net/uid-11469366-id-1747286.html 补充:a.out早期并不是elf格式的,而是unix下另一种可执行格式,新的a.out是 本文讨论了 UNIX/LINUX 平台下三种主要的可执行文件格式:a.out(…

投影矩阵

作者:桂。 时间: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 按回车

Ubuntu 14.04 LTS 下升级 gcc 到 gcc-4.9、gcc-5 版本

转载: http://www.cnblogs.com/BlackStorm/p/5183490.html Ubuntu 14.04 LTS 下升级 gcc 到 gcc-4.9、gcc-5 版本 如果没记错的话,阿里云ECS上的Ubuntu也是LTS版本。 如果还在使用较旧版本的Ubuntu,或者是Ubuntu LTS,那么我们是…

SVM学习——在matlab上安装libsvm库(一)

环境搭建平台: Windows PCMATLAB 软件libsvm库(SVM工具箱)安装步骤 准备工作 安装对应的编译器,在网站上查看当前matlab版本支持的编译器版本。本文使用的matlab 2015a版本(支持编译器详情)。通过查询&…