2022年阵列信号处理中DOA算法分类总结 .pdf





《2022年阵列信号处理中DOA算法分类总结 .pdf》由会员分享,可在线阅读,更多相关《2022年阵列信号处理中DOA算法分类总结 .pdf(12页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、阵列信号处理中的DOA (窄带)空域滤波波束形成:主要研究信号发射/接收过程中的信号 增强。空间谱估计空域参数估计:从而对目标进行定位/ 给空域滤波提供空域参数。测向波达方向估计 (DOA) 空间谱: 输出功率 P 关于波达角 的函数, P().延迟相加法 /经典波束形成器注,延迟相加法和CBF 法本质相同,仅仅是 CBF 法的最优权向量是归一化 了的。1、传统法常规波束形成 CBF / Bartlett 波束形成器常规波束形成(CBF: Conventional Beam Former )Capon最小方差法 /Capon 波束形成器 / MVDR 波束形成器最小方差无畸变响应(MVDR :
2、minimum variance distortionless response )Root-MUSIC 算法多重信号分类法解相干的 MUSIC 算法(MUSIC)基于波束空间的 MUSIC 算法TAM旋转不变子空间法LS-ESPRIT(ESPRIT)TLS-ESPRIT确定性最大似然法( DML :deterministic ML )3、最大似然法随机性最大似然法( SML:stochastic ML)4、综合法: 特性恢复与子空间法相结合的综合法,首先利用特征恢复方案区分多个信号,估计空间特征,进而采用子空间法确定波达方向最大似然估计法是 最优的方法,即便是在信噪比很低的环境下仍然具有良好
3、的性能,但是通常计大 特征值对应的特征矢量:信号子空间小 特征值对应的特征矢量 :噪声子空间无法超过瑞利 限的 制约, 分辨率上有 本质 的局限性。2、 子空间法阵列信号处理波达方向估计的算法子空间算法可以 突 破瑞 利限,达 到 较 高的分辨率计算量小, 不需进行谱峰搜索精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页,共 12 页算量很大。同子空间方法不同的是,最大似然法在原信号为相关信号 的情况下也能保持 良好的性能 。阵列流形矩阵 (导向矢量矩阵) 只要确定了阵列各阵元之间的延迟,就可以很容易地得出一个特定阵列天线的阵列流形矩阵A 。
4、传统的波达方向估计方法是基于波束形成 和零波导引 概念的,并没有利用接收信号向量的模型(或信号和噪声的统计特性) 。知道阵列流形A 以后,可以对阵列进行 电子导引 ,利用电子导引可以把波束调整到任意方向上,从而寻 找输出功率的峰值 。常规波束形成 (CBF)法CBF 法,也称延迟 相加法 /经典波束形成器法 /傅里叶法 / Bartlett 波束形成法 ,是最简单的DOA 估计方法之一。这种算法是使波束形成器的输出功率相对于某个信号为最大。(参考自:阵列信号处理中DOA估计及 DBF技术研究 _赵娜)注意:理解 信号模型注意:上式中,导向矩阵A 的行向量 表示第 K 个天线阵元对N个不同的信号
5、s(i) 的附加权值, 列向量 表示第 i 个信号 s(i) 在 M个不同的天线上的附加权值。将式( 2.6 )的阵元接收信号,写成矢量形式为:X(t)=AS(t)+N(t)精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页,共 12 页其中, X(t)为阵列的 M 1 维 快拍数据矢量, N(t)为阵列的 M 1 维噪声数据矢量, S(t)为信号空间的 N 1维矢量, A 为空间阵列的 M N 维阵列流型矩阵(导向矢量矩阵) ,且A = a1(0)a2(0)aN(0)其中,导向矢量 ai(0)为列矢量 ,表示第 i 个信号在 M个天线上的附加
6、权值ai(0) = ? (-?0?1?)? (-?0?2?)? (-?0?) ,? = 1,2,?式中,?0= 2? = 2? , 其中, c 为光速, 为入射信号的波长。对于均匀线阵,第k 个天线阵元对接收到的第i 个信号 s(i)的时间延时为 ?,则有:?= (? - 1)? ?,k = 1,2,M,其中, d 为阵元间距,一般取d=? /2。第 i 个信号在天线阵元上的入射角为?i。由上述的知识可知 , 一旦知道阵元间的延迟表达式 , 就很容易得出特定空间阵列的导向矢量或阵列流型。波束形成 技术的 基本思想 :通过将各阵元输出进行加权求和, 在一时间内将阵列波束“导向”到一个方向上 ,
7、对期望信号得到最大输出功率的导向位置即是波达方向估计值, 如图 1 所示。假设空间存在 M个阵元组成的阵列 ,N 个信号源,各阵元的 权矢量 为?= ?1?2?阵列的输出为? (? ) = ? (? ) = ?=1?(?)则整个阵列输出的平均功率为? (? ) =1? |?(?)|2= ? ? (? )? (? )?= ?=1其中, R为接收信号矢量x(t)的自相关矩阵图 1 阵列信号处理示意图假设来自 方向的输出功率最大 , 则该最大化问题可表述为:?= ? (? )精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页,共 12 页= ? ?
8、(? )? (? )? = ? ? (? )? (? )? = ? |? (? )|2|? (? )|2+ ?2 ? 2为了使加权向量 w的权值不影响输出信噪比, 在白化噪声方差 ?2一定的情况下, 取 ? 2= 1,此时求解为:?=?(?) ?(? )?(?)此时 Bartlett 波束形成器的空间谱为:?(? ) = ?=?(? )?(?)?(?)?(?)延迟相加法(本质和 Bartlett算法相同,仅最优权向量不同,后者的最优权是归一化了的。)(参考自:阵列信号处理中的DOA 估计技术研究_白玉)k 时刻,令 x(t)=u(k),s(t)=s(k),n(t)=n(k) ,上面公式中: ?
9、(?) = ? (? ),? (? ) = ? (? ),令u(k)=a( )s(k)+n(k),波束形成器输出信号y(k)是传感器阵元输出的线性加权之和,即y(k)=wHu(k)(2-1)传统的波束形成器总的输出功率可以表示为:Pcbf =E|y(k)|2=E|wHu(k)|2=wHEu(k)uH(k)w=wHRuuw(2-2)式中,Ruu定义为阵列输入数据的自相关矩阵。式(2-2)在传统 DOA 估计算法中的地位举足轻重。自相关矩阵 Ruu包含了阵列响应向量和信号自身的有用信息,仔细分析 Ruu,可以估计出信号的参数。考察一个以角度 入射到阵列上的信号s(k),则有 u(k)=a( )s(
10、k)+n(k)。根据窄带输入数据模型 ,波束形成器的输出功率可以表示成:Pcbf( )=E|wHu(k)|2=E|wH(a( )s(k)+n(k) )|2 =|wHa( )|2s2+|wH|2n2(2-3)式中, s2= Es(k)2, a( )是关于 DOA 角 的导引向量,n(k)是阵列输入端的噪声向量。 当 w=a( )时,系统的输出 (信号)功率达到最大。这是因为,权值向量w 在传感器阵元处和来自方向的信号分量相位对齐,使得它们能够同相相加,从而使系统的输出功率相对于某个信号为最大。在 DOA 估计的经典波束形成方法中,波束形成器产生的波束在感兴趣的区域中离散地扫描,对应不同的 可以产
11、生不同的权向量:?= ?(?)从而得到的输出功率也不相同。利用式(2-3),经典波束形成器的输出功率与波达方向的关系由下式给出 :Pcbf( )=wHRuuw=aH(? )Ruua(? ) (2-4)因此,如果我们对输入自相关矩阵进行估计,知道对所有感兴趣的导引向量(通过校准或分析计算) ,就可能估计出输出功率关于波达角 的函数。输出功率关于波达角的函数通常称为空间谱(spatial spectrum) 。很明显 ,通过锁定式 (2-4)定义的空间谱的峰值就可以估计出波达方向。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页,共 12 页最大
12、的功率对应着 最大的峰值 ,而最大的空间谱峰所对应的 角度方向 即为信号的 波达方向。延迟相加法 (常规波束形成器法 ),CBF 法(Bartlett 波束形成器法)具有一定的局限性,可以很好的识别 单个信号, 但是当存在着来自多个方向的信号时, 该方法要受到 波束宽度 和旁瓣高度 的制约,因而这种方法的分辨率较低, 只能大致分辨出信号所处的角度范围。这是因为, 延迟相加法是把阵列形成的波束指向 某个方向 ,由此可以获得来自于这个方向的信号的最大功率。就单个信号而言, 延迟相加法可以很好地估计出它的波达方向。但是当信号空间中存在多个信号的时侯,因为波束宽度的限制,受到同一个波束内信号之间的相互
13、干扰,延迟相加法的估计性能就会急剧的下降。增加阵列的 阵元数(M ) 可以改善延迟 相加法的性能,提高分辨率,但是这会使系统更加复杂,还会增加算法的计算量和数据存储空间。Capon 最小方差法(Capon 波束形成器,也称MVDR 波束形成器 )最小方差无畸变响应 (MVDR )波束形成器解决了延迟 相加法分辨率差的缺点, 用一部分自由度在期望方向 上形成一个波束, 利用剩余的一部分自由度在干扰方向 形成零陷。这种方法使得输出功率 和 信号方差达到最小,使得非期望干扰信号的贡献为最小,同时使观测方向上的增益达到最大,约束条件为 wHa( )=1 ,使得来自期望方向的信号功率不变。其优化问题表述
14、为:?= ? (? )约束条件为:? (? ) = 1综合上式求解 w 为:?=?-1?(?)?-1?(?)此时 Capon 波束形成器的空间谱为:?= ?=1?-1?(?)Capon算法比延迟 相加法有了一定程度的改进, 可以对多个信号进行DOA 估计。但是 Capon 算法只能分辨 非相干 信号,当存在与感兴趣信号相关的其它信号时,它就不能起作用了。这是因为Capon 算法在运算的过程中使用到了信号的自相关矩阵 , 因而不能对干扰信号形成 零陷。也就是说,在使得输出功率为最小的过程当中,相关分量可能会恶性合并。此外,Capon算法运算时需要对信号的自相关矩阵求逆,当阵列加大时会有巨大的运算
15、量。对于任意的 ,PCapon( )是来自方向 的信号功率的最大似然估计。多重信号分类( MUSIC )算法为代表的 子空间分解 类算法开始兴起。这一类算法有一个共同的特点,就是需要对阵列的 接收数据矩阵进行数学分解(如奇异值分解、特征值分解和QR 分解等),将数据分解成两个互相 正交的特征子空间:一个是信号子空间,另一个是噪声子空间。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页,共 12 页子空间类算法按照处理方式的不同可以分成两类:一种是以 MUSIC 算法为代表的 噪声子空间类算法另一种是以ESPRIT 算法为代表的 信号子空间类算
16、法。式中, Rs是信号相关矩阵 ( signal correlation matrix ),EssH 。R的特征值为 0, 1, 2, ., M-1 , 使得|R- iI|=0 (2-12)利用式 (2-11),我们可以把它改写为|ARsAH+?2I-?iI|=|ARsAH-(?i-?2)I|=0 (2-13) 因此 ARsAH的特征值 (eigenvalues)?为?=?i- n2(2-14)因此 A 是由线性独立 的导引向量构成的 ,因此是 列满秩 的,信号相关矩阵 Rs也是非奇异 的,只要入射信号不是高度相关的。列满秩的A 和非奇异的Rs可以保证 ,在入射信号数L 小于阵元数M 时,M
17、M的矩阵 ARsAH是半正定 的,且秩为D。这意味着ARsAH的特征值?中,有 M-L 个为 零。由式 (2-14)可知 ,R 的特征值?i中有 M-L 个等于噪声方差n2。该M-L 个最小特征值?i相关的 特征向量 ,和构成 A 的 L 个导引向量 正交 。噪声子空间和信号子空间是相互正交的,而由导向矢量 所张成的空间与 信号子空间 是一致的。应当指出,与传统方法不同, MUSIC算法在估计信号功率时并没有考虑波达角。在噪声与信号源非相关 的环境下,可以确保PMUSIC() 的谱峰对应着信号的真实方向。由于PMUSIC()的峰值是可以分辨的, 并且与信号之间的真实角度间隔没有关系,因此从理论
18、上来讲,只要阵元位置校准的足够准确,MUSIC算法就可以分辨出两个邻近的信号。但是当入射信号之间彼此高度相关 时,自相关矩阵 Rxx会变成奇异矩阵 ,这将导致 MUSIC 算法失效。旋转不变子空间算法( ESPRIT )是空间谱估计算法中的典型算法之一,它和前面介绍的 MUSIC 算法一样,也需要对阵列接收数据的协方差矩阵进行特征分解。但是两者也存在着明显的不同点,即MUSIC 算法利用了阵列接收数据的协方差矩阵的噪声子空间和导向矢量之间的正交特性 ,而ESPRIT 算法则利用了阵列接收数据的协方差矩阵信号子空间的旋转不变性 ,所以 MUSIC 算法与 ESPRIT 算法可以看成为是一种互补的
19、关系。和 MUSIC 算法相比, ESPRIT 算法直接给出了待估角的闭式解,不需谱峰搜索,计算复杂度更小总的来说 ESPRIT 算法的性能要 差于 MUSIC 算法。 ESPRIT 算法的优势在于它的实时性 ,一般的情况下,只要有两个子阵列满足旋转不变性,就可以用ESPRIT 算法来实现,且其实现速度要优于MUSIC 算法。与精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页,共 12 页MUSIC 算法不同的是, ESPRIT 算法不需要 知道精确的导向矢量,只需要各 子阵之间的 阵元保持一致 。这种算法思想是将接收阵列在 几何结构 上分为
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 2022年阵列信号处理中DOA算法分类总结 2022 阵列 信号 处理 DOA 算法 分类 总结

限制150内