
esprit算法研究.doc
26页课程设计报告实验名称:ESPRIT算法研究实验日期:姓 名:学 号:哈尔滨工业大学(威海)设计任务实现空间谱估计算法,并考察算法性能方案设计1) 由均匀线阵形式,确定阵列的导向矢量;2) 由阵列导向矢量,对接收信号进行建模仿真;3) 由ESPRIT算法实现信号DOA古计;4) 考察算法性能与信噪比,采样率,观测时间等参数的关系三、设计原理3.1空间谱估计数学模型空间谱估计就是利用空间阵列实现空间信号的参数估计的一项专门技术 整个空间谱估计系统应该由三部分组成:空间信号入射、空间阵列接收及参数估计相应地可分为三个空间, 即目标空间、观察空间及估计空间, 也就是说空间谱估计系统由这三个空间组成, 其框图见图1图1空间谱估计的系统结构对于上述的系统结构,作以下几点说明1) 目标空间是一个由信号源的参数与复杂环境参数张成的空间 对于空间谱估计系统,就是利用特定的一些方法从这个复杂的目标空间中估计出信号的未知参数2) 观察空间是利用空间按一定方式排列的阵元,来接收目标空间的辐射信号由于环境的复杂性,所以接收数据中包括信号特征 (方位、距离、极化等)和空间环境特征(噪声、杂波、干扰等)另外由于空间阵元的影响, 接收数据中同样也含有空间阵列的某些特征 (互耦、通道不一致、频带不一致等 )。
这里的观察空间是一个多维空间,即系统的接收数据是 由多个通道组成,而传统的时域处理方法通常只有一个通道 特别需要指出的是: 通道与阵元并不是一一对应, 通道是由空间的一个、 几个或所有阵元合成的(可用加权或不加权),当然空间某个特定的阵元可包含在不同的通道内3)估计空间是利用空间谱估计技术 (包括阵列信号处理中的一些技术,如阵列校正、空 域滤波等技术)从复杂的观察数据中提取信号的特征参数从系统框图中可以清晰的看出,估计空间相当于是对目标空间的一个重构过程,这个重构的精度由众多因素决定,如环境的复杂性、空间阵元间的互耦、通道不一致、 频带不一致等3.2 阵列信号处理首先,考虑N个远场的窄带信号入射到空间某阵列上,阵列天线由 M个阵元组成,这里假设阵元数等于通道数, 即各阵元接收到信号后经过各自的传输信道送到处理器, 也就是说 处理器接收来自M个通道的数据si(t)ui (t)ej( 0t(t))si (tj( 0 (t ) (t) ui (t )e 0))(3.2-1)式中, ui (t) 是接受信号的幅度,窄带远场信号源的假设下,有(t) 是接收信号的相位,0是接收信号的频率在ui (t(tui (t)(t)(3.2-2)根据式 (3.2-1) 和式 (3.2-2) ,显然有下式成立:si (tsi (t)ej 0(3.2-3)则可以得到第L个阵元接收信号为Nxl(t) glisi(t li ) nl(t)i1l 1,2,,M (3.2-4)式中,go为第L个阵兀对第i个信号的增益,nl (t)表示第L个阵元在t 时刻的噪声, li表示第i个信号到达第L个阵兀时相对参考阵兀的时延。
将M个阵兀在特定时刻接收的信号排列成一个列矢量,可得x1(t) g11ej 011 g12ej 0 12j 0 1N g1Nes1(t)n1(t)x2(t) g21ej 021 g22ej 022j 0 2 N g2Nes2 (t)n2(t)2 (3.2-5)xM (t) gM1ej 0 M1 gM2ej 0 M2g ej 0 MNgMN esN (t)nM (t)在理想情况下,假设阵列中各阵元是各向同性的且不存在通道不一致、互耦等因素的影响,则式 (3.2-4) 中的增益 gli 可以省略 (即归一化 1) ,在此假设下式 (3.2-5) 可以简化为x1(t)j 0 11 ej 0 12 ej 0 1N es1(t)n1(t)x2 (t)j 0 21 ej 0 22 ej 0 2N es2(t)n2(t)(3.2-6)xM (t)e j 0 M1ej 0 M 2e j 0 MNsN (t)nM (t)将式 (3.2-6) 写成矢量形式如下:x(t)As(t)n(t)(3.2-7)式中, x(t) 为阵列的M 1 维快拍数据矢量,n(t) 为阵列的 M1维噪声数据矢量,s(t)为空间信号的N 1维矢量,A为空间阵列的MN维流型矩阵(导向矢量阵),且A ai ( 0 ) a2(0) aN(0) (3.2-8)其中导向矢量exp(j o ii)ai( o) eXP(j 0 2i) i 1,2, , N (3.2-9)exp(j 0 Mi)式中o 2 f 2 c , c为光速,为波长。
由上述的知识可知,一旦知道阵元间的延迟表达式T,就很容易得出待定空间阵列的导向矢量或阵列流型 下面推导一下空间阵元间的延迟表达式 假设空间任意两个阵元, 其中一个为参考阵元(位于原点),另一个阵元的坐标为(x ,y, z),两阵元的几何关系见图, 图中“X”表示阵元'7方惊角ML八...............■二图2空间任意两阵元的几何关系由几何关系可以推导出两阵元的波程差为1 (x cos cos ysin cos zsin ) (3.2-10)c这里的波程差其实就是位于 x轴上两阵元间的延迟、位于y轴上两阵元间的延迟和位于 z轴上两阵元间的延迟之和根据式(3.2-10)的结论,下面给出实际环境中常用的几种阵列及阵元间的相互延迟表 达式1)平面阵 设阵元的位置为(Xk,yQ(k 1,2, ,M),以原点为参考点,另假设信号 入射参数为(i, i)(i 1,2, , N),分别表示方位角与俯仰角,其中方位角表示与 x轴的夹角⑵线阵设 阵元的位置为Xk(k 1,2, ,M ),以原点为参考点,另假设信号入射参XZ'XX……X/ L—I—d―I »Xjt) X2(t) 1 ——X3(t)图3-1均匀线阵的数学模型示意图 Xm (t)(3.2-11)数为 ,i 1,2, , N ),表示方位角,其中方位角表示与y轴的夹角(即与线阵法线的夹角) 则有ki — (XkSin i)c(3)均匀圆阵 设以均匀圆阵的圆心为参考点,则有kix cosc Mcos(3.2-12)其中方位角表示与x轴的夹角,r为圆半径。
3.3旋转不变子空间算法原理3.3.1信号模型算法介绍前,首先对信号进行建模为了推导分析的方便,将波达方向的 数学模型做如下理想状态的假设:1) 阵列形式为线性均匀阵,阵元间距不大于信号波长的二分之一2) 存生两个完全相同的子阵,且两个子阵的间距△是己知的3) 噪声序列为一零均值高斯过程,各阵元间噪声相互独立,噪声与信号也 相互独立4) 空间信号为零均值平稳随机过程,通常为窄带远场信号5) 信号源数小于子阵阵列元数,信号取样数大于子阵阵列元数,以确保子 阵阵列流型的各列线性独立6) 组成阵列的各传感器为各向同性阵元,且无互耦以及通道不一致的干扰 下图给出了均匀线阵的数学模型示意图:332算法原理对于均匀线阵,相邻子阵间存在一个固定间距,这个固定间距反映出各相邻 子阵间的一个固定关系,即子阵间的旋转不变性,而 ESPRIT算法正是利用了这 个子阵间的旋转不变性实现阵列的 DOA估计ESPRIT算法最基本的假设是存在两个完全相同的子阵,且两个子阵的间距 是已知的由于两个子阵的结构完全相同,且子阵的阵元数为m,对于同一个 信号而言,两个子阵的输出只有一个相位差 i , i=1, 2,…N下面假设第一个子阵的接收数据为 X!,第二个子阵的接收数据为X2,根据 前面所述的阵列模型可知X! [a( !)L a( n)]S N! AS 2(3.1)(3.2)X2 [a( i)ej !L a( N)ej n]S N2 A S N2式中,子阵1的阵列流型A1=A,子阵2的阵列流型A= A ,且式中diag[ej1…ejN] ( 3.3)从上面的数学模型可知,需要求解的是信号的方向,而信号的方向信息包含在A和 中,由于 是一个对角阵,所以下面只考虑这个矩阵,即k (2 sin Q/ (3.4)由上可知。
只要得到两个子阵间的旋转不变关系 ,就可以方便地得到关于 信号到达角的信息下面的任务就是从式(3.1)和式(3.2)中得到两个子阵间的关系先将两个子阵的模型进行合并,即X1X2S N AS N在理想条件下,可得上式的协方差矩阵H — —HR E[XXh] ARSA Rn对上式进行特征分解可得2mR ieeHUs sUsHUn nUnH(3.5)(3.6)(3.7)显然上式中得到的特征值有如下关 i >-> n > N1=-=2m , US为大特征 值对应的特征矢量张成的信号子空间,U n为小特征值对应矢里张成的噪声子空间对于实际的快拍数据,式(3.7)应修正如下:R L?s?sU?s L?n?nU?H (3.8)由前面的知识可知,上述的特征分解中大特征矢量张成的信号子空间与阵 列流型张成的信号子空间是相等的即spanUS span(A( )) (3.9)(3.10)此时,存在一个惟一的非奇异矩阵 T,使得Us A( )T显然,上述的结构对两个子阵都成立,所以有UsUsiU S2ATA T(3.11)很显然,由子阵1的大特征矢量张成的子空间Usi、由子阵2的大特征矢量张成 的子空间Us2与阵列流型A张成的子空间三者相等,即span U S1 span A( ) span U S2 ( 3.12)另外,由两个子阵列在阵列流型上的关系可知A2 A1 ( 3.13)再利用式(3.11)可知两个子阵列的信号子空间的关系如下:Us2 US1T 1 T Us1 ( 3.14)式(3.13)反映了两个子阵列的阵列流型间的旋转不变性,而式( 3.14)反映了两个子阵的阵列接收数据的信号子空间的旋转不变性。
如果阵列流型A是满秩矩阵,则由式(3.14)可以得到T T 1 ( 3.15)所以上式中 的特征值组成的对角阵一定等于 ,而矩阵T的各列就是矩阵特征矢量所以一旦得到上述的旋转不变关系矩阵 ,就可以直接利用式(3.4)得到信号的入射角度3.4标准的旋转不变子空间算法有上节的知识可知,ESPRIT算法的基本原理就是利用式(3.14)的旋转不 变性,常规的旋转不变子空间算法就是利用上述的基本原理求解信号的入射角度 信息下面就分析解这个等式的两种最经典、应用最广泛方法 :最小二乘(LS)法和总体最小二乘(TLS )法3.4.1 最小二乘法由最小二乘的数学知识,我们知道式(3.14 )的最小二乘解的方法等价于2min| Us2〔,约束条件 U si U S2 U S2 ( 3.16)因此最小二乘法的基本思想就是使校正项 U S2尽可能小,而同时保证满足约束条件为了得到LS解,将式(3.14)代入式(3。












