现代信号处理课件之空间谱估计(空域MUSIC算法)概要.ppt

1、11.1.阵列信号处理简介阵列信号处理简介（array signal processing)九、九、空间谱估计空间谱估计n虽然早在1978年就已在军用通信系统中使用了自适应天线，但在民用的蜂窝式通信中使用天线阵列却只是在1990年才开始的。n近20年来，阵列信号处理作为信号处理的一个重要分支，在通信、雷达、声呐、地震、勘探、射电天文等领域获得了广泛应用和迅速发展。阵列信号处理的最重要应用包括阵列信号处理的最重要应用包括：信（号）源定位确定阵列到信源的仰角和方位角，甚至距离（若信源位于近场）；信源分离确定各个信源发射的信号波形。各个信源从不同方向到达阵列，这一事实使得这些信号波形得以分离，即使它

2、们在时域和频域是叠加的。信道估计确定信源与阵列之间的传输信道的参数（多径参数）。n阵列信号处理的主要问题包括：波束形成技术-使阵列方向图的主瓣指向所需的方向；零点形成技术-使天线的零点对准干扰方向；空间谱估计对空间信号波达方向的分布进行超分辨估计。22.2.阵列信号处理问题阵列信号处理问题（array signal processing)阵列：多个天线的组合阵列：多个天线的组合每个天线阵元：天线、传感器每个天线阵元：天线、传感器假设假设:(1)窄带信号窄带信号 ：点信源：点信源 (2)远场远场(far field):波前平面波波前平面波()ins九、九、空间谱估计空间谱估计*()in 1d23

3、4mi波达方向 (DOA:direction of arrival)：入射线与法线之间的夹 角，可以有正有负。2sin/2iidd波长（半波长条件）：若不满足该条件，会出现（半波长条件）：若不满足该条件，会出现DOA估计估计的模糊的模糊12m()is n()ijis n e(1)()ij mis n eP个信号个信号()is n信号信号 的方向向量，的方向向量，(阵列响应阵列响应)向量向量：1212(1)1(1)(1)(1)()1,.,()(),.,()111iippTjj mipjjjj mj mj meeeeeeeeaAaa方向矩阵方向矩阵满列秩满列秩12.pVandermonde矩阵矩阵

4、4信号模型信号模型1()()()(),1,.,pkkiikix nas ne nkm阵元阵元k上的观测数据上的观测数据11()()()()()()TTmmnx nxnne nenxe，.，.，111()()()()()()()TTpmpmxpns nsnsAaa，.，.，阵列信号处理的数学模型阵列信号处理的数学模型：()()()()nnnxAse阵列信号处理的问题阵列信号处理的问题：已知数据向量已知数据向量 ，求空间参数，求空间参数 (1)()Nxx，.,1p，.,N个快拍个快拍波达方向波达方向空间谱估计空间谱估计53.最优波束形成器（最小方差谱估计最优波束形成器（最小方差谱估计）DOA 估计

5、：波束形成器估计：波束形成器加权求和)()(1*nxwnzimii干扰信号拒绝其它信号期望信号只包含输出信号)()(nxnzdmww,1设计一个滤波器抽头(权系数)空间谱估计空间谱估计最小输出能量最小输出能量(MOE:minimum output energy)准则：准则：5*1()(),()()mHiiiz nnnw x nw xw x22111211111()()()()1min()min1()()NNNHHHnnnNHxxnNHxxnz nnnnNNNz nNnnNw xwxxww R wRxx则其中空间谱估计空间谱估计211()minNnz nN6211lim min()minNHxx

6、Nnz nNw R w1,()()()()()()()()()pkkiiii knnns ns nnxAseaae2212222221,1()lim()()()()()()()NHHNnpHHkkiiii kE z nz nEnnNE s nE s nwxxww aw aw()1()0,HkHiikw aw a波束形成条件干扰拒绝条件，零点形成条件(期望信号期望信号)(干扰信号干扰信号)(加性噪声加性噪声)空间谱估计空间谱估计722222()()()1min()kHkEz nEsnEz nww a则，在约束条件下，使Largange乘子法乘子法：221*1()1()()0()()1()1()(

7、)HkHxxoptxxkHHoptkkoptHkxxkJE z nE z nJ ww aw R wwwR awwaawaR a()=其中()又代入上式空间谱估计空间谱估计8()1xxkoptH1kxxkR a()wa()R a()最佳滤波器由由capon提出，称为最小方差无畸变提出，称为最小方差无畸变(MVDR:minimum variance distortionless response)波束形成器波束形成器即为所求。最大幅值对应的空间谱：求关键：kk1xxkHkkk)a()R(a1)P()a(空间谱估计空间谱估计93.子空间方法子空间方法0)()()()(02nnEnnE(n)ETHee

8、Ieee假设1：对于不同的 值，向量 线性独立假设2：各阵元上复加性噪声具有零均值、相同方差，且不相关i)(iajineneEnjynxnejiiii,0)()(),()()(则令复白噪声分量0)()(2)()()(222nynxEjnyEnxEneEiiiii（实部和虚部不相关，具有相同方差）假设3：满秩矩阵（非奇异）)()(nnEHssP 空间谱估计空间谱估计n虽然早在1978年就已在军用通信系统中使用了自适应天线，但在民用的蜂窝式通信中使用天线阵列却只是在1990年才开始的。n近20年来，阵列信号处理作为信号处理的一个重要分支，在通信、雷达、声呐、地震、勘探、射电天文等领域获得了广泛应用

9、和迅速发展。阵列信号处理的最重要应用包括阵列信号处理的最重要应用包括：信（号）源定位确定阵列到信源的仰角和方位角，甚至距离（若信源位于近场）；信源分离确定各个信源发射的信号波形。各个信源从不同方向到达阵列，这一事实使得这些信号波形得以分离，即使它们在时域和频域是叠加的。信道估计确定信源与阵列之间的传输信道的参数（多径参数）。n阵列信号处理的主要问题包括：波束形成技术-使阵列方向图的主瓣指向所需的方向；零点形成技术-使天线的零点对准干扰方向；空间谱估计对空间信号波达方向的分布进行超分辨估计。112()()()()()()()()()()()()()()HxxHHHHHEnnEnns nnEnnE

10、nnRxxAseAeAssAeeAPAImp pp pm特征值分解特征值分解：2112200ppI2HHU APA UIURUxxH即即空间谱估计空间谱估计12的特征值：xxR222,1,1,iiiaipipm按特征值大小，将U U分为两部分：mppuuuuGSU,|,|11信号S噪声G子空间子空间：向量组 线性组合的集合，称为 张成的空间。paa,1paa,111,ppjjjjCspanaaa空间谱估计空间谱估计13信号子空间：噪声子空间：观测空间：ppuuspanssspan,11mppuuspanggspan,11mNuuspanxxspan,)(,),1(1观测空间信号子空间噪声子空间

11、特征值分解后，与大特征值对应 与小特征值对应空间谱估计空间谱估计14,US G子空间的几何意义：子空间的几何意义：,HHHHHHHSS SS GU US GIGG SG G,0,0HHHHpm pS SIG GIG SS G11,(,(HHHsHHHGPS S SSSSSSPG G GGGGGG称为信号子空间)称为噪声子空间)投影矩阵投影矩阵空间谱估计空间谱估计15ssPIP正交投影矩阵正交投影矩阵,HHHHHSUUS GSSGGIG1,HHHs GGISSIS S SSP噪声子空间是信号子空间的正交补噪声子空间是信号子空间的正交补，几何意义：信号子空间和噪声子空间正交几何意义：信号子空间和噪

12、声子空间正交空间谱估计空间谱估计164.MUSIC 方法方法HxxRUUUS G121222,0,HHxxHxx SR G U U GSGGI GSGGI IR GG即空间谱估计空间谱估计222000(0iff0)()()()()0HxxHxxHHHHHHTiHHiiRAPAIR GAPA GGGAPA GG APA GA Gt RttaGoaGG a行向量标量17MUSIC空间谱空间谱：)()(1)()(1)(aSSIaaGGaHHHHP噪声子空间方法信号子空间方法（需一维搜索），就给出个取峰值的ppp1)(波束形成器波束形成器：)()()(11dHxxddHxxoptaRaaRw空间谱估计

13、空间谱估计作为比较作为比较,这里重新给出这里重新给出MV谱估计结果：谱估计结果：)()(1)(1aRapHMVP（最小方差谱最小方差谱)(最佳滤波器系数最佳滤波器系数)18:MUSIC方法恢复谐波功率谱应用于雷达动目标检测功率谱估计应用于3G移动通信v雷达回波信号功率谱雷达回波信号功率谱v我们想检测一个运动目标，但是该目标到雷达的路径上的建筑我们想检测一个运动目标，但是该目标到雷达的路径上的建筑物、海浪、云雾等引起杂波信号，因此雷达接受的回波信号，包物、海浪、云雾等引起杂波信号，因此雷达接受的回波信号，包括来自运动目标的信号和杂波的信号。括来自运动目标的信号和杂波的信号。v 而此时，杂波电平往

14、往比运动目标的回波电平高很多，运动而此时，杂波电平往往比运动目标的回波电平高很多，运动目标回波目标回波“淹没淹没”在杂波中，从时域分开已经不可能了。在杂波中，从时域分开已经不可能了。v功率谱估计应用于功率谱估计应用于3G移动通信移动通信v电信科学技术研究院代表我国电信主管部门向国际电联提交的电信科学技术研究院代表我国电信主管部门向国际电联提交的TD-SCDMA中采用的智能天线技术中采用的智能天线技术v 智能天线原名自适应天线阵列（智能天线原名自适应天线阵列（AAA,Adaptive Antenna Array),自自适应天线阵由多个天线单元组成。适应天线阵由多个天线单元组成。v利用基带数字信号

15、处理技术识别用户信号到达方向（利用基带数字信号处理技术识别用户信号到达方向（DOA），并在此），并在此方向形成空间定向波束，使天线主波束对准用户信号到达方向，同时使方向形成空间定向波束，使天线主波束对准用户信号到达方向，同时使旁瓣或零陷对准干扰信号到达方向，从而给有用信号带来最大增益，有旁瓣或零陷对准干扰信号到达方向，从而给有用信号带来最大增益，有效减少多径影响，同时达到对干扰信号删除和抑制的目的。效减少多径影响，同时达到对干扰信号删除和抑制的目的。v智能天线另一个可能的用途是进行紧急呼叫定位，并提供更高的定位智能天线另一个可能的用途是进行紧急呼叫定位，并提供更高的定位精度，因为在获得可用于定位的时延、强度等信息的同时它还可以获得精度，因为在获得可用于定位的时延、强度等信息的同时它还可以获得波达角信息。波达角信息。v 采用智能天线的基站可以获得接收信号的空间特征矩阵，由此获得信采用智能天线的基站可以获得接收信号的空间特征矩阵，由此获得信号的功率估值和到达方向（号的功率估值和到达方向（DOA)。v