esprit算法研究.docx
《esprit算法研究.docx》由会员分享,可在线阅读,更多相关《esprit算法研究.docx(18页珍藏版)》请在冰豆网上搜索。
esprit算法研究
课程设计报告
实验名称:
ESPRIT算法研究
实验日期:
姓 名:
学 号:
哈尔滨工业大学(威海)
一、设计任务
实现空间谱估计算法,并考察算法性能。
二、方案设计
1)由均匀线阵形式,确定阵列得导向矢量;
2)由阵列导向矢量,对接收信号进行建模仿真;
3)由ESPRIT算法实现信号DOA估计;
4)考察算法性能与信噪比,采样率,观测时间等参数得关系。
三、设计原理
3、1空间谱估计数学模型
空间谱估计就就是利用空间阵列实现空间信号得参数估计得一项专门技术。
整个空间谱估计系统应该由三部分组成:
空间信号入射、空间阵列接收及参数估计。
相应地可分为三个空间,即目标空间、观察空间及估计空间,也就就是说空间谱估计系统由这三个空间组成,其框图见图1。
图1空间谱估计得系统结构
对于上述得系统结构,作以下几点说明。
(1)目标空间就是一个由信号源得参数与复杂环境参数张成得空间。
对于空间谱估计系统,就就是利用特定得一些方法从这个复杂得目标空间中估计出信号得未知参数。
(2)观察空间就是利用空间按一定方式排列得阵元,来接收目标空间得辐射信号。
由于环境得复杂性,所以接收数据中包括信号特征(方位、距离、极化等)与空间环境特征(噪声、杂波、干扰等)。
另外由于空间阵元得影响,接收数据中同样也含有空间阵列得某些特征(互耦、通道不一致、频带不一致等)。
这里得观察空间就是一个多维空间,即系统得接收数据就是由多个通道组成,而传统得时域处理方法通常只有一个通道。
特别需要指出得就是:
通道与阵元并不就是一一对应,通道就是由空间得一个、几个或所有阵元合成得(可用加权或不加权),当然空间某个特定得阵元可包含在不同得通道内。
(3)估计空间就是利用空间谱估计技术(包括阵列信号处理中得一些技术,如阵列校正、空域滤波等技术)从复杂得观察数据中提取信号得特征参数。
从系统框图中可以清晰得瞧出,估计空间相当于就是对目标空间得一个重构过程,这个重构得精度由众多因素决定,如环境得复杂性、空间阵元间得互耦、通道不一致、频带不一致等。
3、2阵列信号处理
首先,考虑N个远场得窄带信号入射到空间某阵列上,阵列天线由M个阵元组成,这里假设阵元数等于通道数,即各阵元接收到信号后经过各自得传输信道送到处理器,也就就是说处理器接收来自M个通道得数据。
(3、2-1)
式中,就是接受信号得幅度,就是接收信号得相位,就是接收信号得频率。
在窄带远场信号源得假设下,有
(3、2-2)
根据式(3、2-1)与式(3、2-2),显然有下式成立:
(3、2-3)
则可以得到第L个阵元接收信号为
(3、2-4)
式中,为第L个阵元对第i个信号得增益,表示第L个阵元在t时刻得噪声,表示第i个信号到达第L个阵元时相对参考阵元得时延。
将M个阵元在特定时刻接收得信号排列成一个列矢量,可得
(3、2-5)
在理想情况下,假设阵列中各阵元就是各向同性得且不存在通道不一致、互耦等因素得影响,则式(3、2-4)中得增益可以省略(即归一化1),在此假设下式(3、2-5)可以简化为
(3、2-6)
将式(3、2-6)写成矢量形式如下:
(3、2-7)
式中,为阵列得维快拍数据矢量,为阵列得维噪声数据矢量,为空间信号得维矢量,A为空间阵列得维流型矩阵(导向矢量阵),且
(3、2-8)
其中导向矢量
(3、2-9)
式中,c为光速,为波长。
由上述得知识可知,一旦知道阵元间得延迟表达式τ,就很容易得出待定空间阵列得导向矢量或阵列流型。
下面推导一下空间阵元间得延迟表达式。
假设空间任意两个阵元,其中一个为参考阵元(位于原点),另一个阵元得坐标为(x ,y, z),两阵元得几何关系见图,图中“×”表示阵元。
图2空间任意两阵元得几何关系
由几何关系可以推导出两阵元得波程差为
(3、2-10)
这里得波程差其实就就是位于x轴上两阵元间得延迟、位于y轴上两阵元间得延迟与位于z轴上两阵元间得延迟之与。
根据式(3、2-10)得结论,下面给出实际环境中常用得几种阵列及阵元间得相互延迟表达式。
(1)平面阵 设阵元得位置为,以原点为参考点,另假设信号入射参数为,分别表示方位角与俯仰角,其中方位角表示与x轴得夹角。
(2)线阵设阵元得位置为,以原点为参考点,另假设信号入射参数为,表示方位角,其中方位角表示与y轴得夹角(即与线阵法线得夹角),则有
(3、2-11)
(3)均匀圆阵设以均匀圆阵得圆心为参考点,则有
(3、2-12)
其中方位角表示与x轴得夹角,r为圆半径。
3、3旋转不变子空间算法原理
3、3、1信号模型
算法介绍前,首先对信号进行建模。
为了推导分析得方便,将波达方向得数学模型做如下理想状态得假设:
1)阵列形式为线性均匀阵,阵元间距不大于信号波长得二分之一。
2)存生两个完全相同得子阵,且两个子阵得间距△就是己知得。
3)噪声序列为一零均值高斯过程,各阵元间噪声相互独立,噪声与信号也相互独立。
4)空间信号为零均值平稳随机过程,通常为窄带远场信号。
5)信号源数小于子阵阵列元数,信号取样数大于子阵阵列元数,以确保子阵阵列流型得各列线性独立。
6)组成阵列得各传感器为各向同性阵元,且无互耦以及通道不一致得干扰。
图3-1均匀线阵得数学模型示意图
下图给出了均匀线阵得数学模型示意图:
3、3、2 算法原理
对于均匀线阵,相邻子阵间存在一个固定间距,这个固定间距反映出各相邻
子阵间得一个固定关系,即子阵间得旋转不变性,而ESPRIT算法正就是利用了这个子阵间得旋转不变性实现阵列得DOA估计。
ESPRIT算法最基本得假设就是存在两个完全相同得子阵,且两个子阵得间距就是已知得。
由于两个子阵得结构完全相同,且子阵得阵元数为m,对于同一个信号而言,两个子阵得输出只有一个相位差,=1,2,…N。
下面假设第一个子阵得接收数据为,第二个子阵得接收数据为,根据前面所述得阵列模型可知
(3、1)
(3、2)
式中,子阵1得阵列流型=A,子阵2得阵列流型=,且式中
(3、3)
从上面得数学模型可知,需要求解得就是信号得方向,而信号得方向信息包含在A与中,由于就是一个对角阵,所以下面只考虑这个矩阵,即
(3、4)
由上可知。
只要得到两个子阵间得旋转不变关系,就可以方便地得到关于
信号到达角得信息。
下面得任务就就是从式(3、1)与式(3、2)中得到两个子阵间得关系。
先将两个子阵得模型进行合并,即
(3、5)
在理想条件下,可得上式得协方差矩阵
(3、6)
对上式进行特征分解可得
(3、7)
显然上式中得到得特征值有如下关≥…≥>=…=,US为大特征值对应得特征矢量张成得信号子空间,为小特征值对应矢里张成得噪声子空间。
对于实际得快拍数据,式(3、7)应修正如下:
(3、8)
由前面得知识可知,上述得特征分解中大特征矢量张成得信号子空间与阵列流型张成得信号子空间就是相等得。
即
(3、9)
此时,存在一个惟一得非奇异矩阵T,使得
(3、10)
显然,上述得结构对两个子阵都成立,所以有
(3、11)
很显然,由子阵1得大特征矢量张成得子空间、由子阵2得大特征矢量张成得子空间与阵列流型A张成得子空间三者相等,即
(3、12)
另外,由两个子阵列在阵列流型上得关系可知
(3、13)
再利用式(3、11)可知两个子阵列得信号子空间得关系如下:
(3、14)
式(3、13)反映了两个子阵列得阵列流型间得旋转不变性,而式(3、14)反映了两个子阵得阵列接收数据得信号子空间得旋转不变性。
如果阵列流型A就是满秩矩阵,则由式(3、14)可以得到
(3、15)
所以上式中得特征值组成得对角阵一定等于,而矩阵T得各列就就是矩阵特征矢量。
所以一旦得到上述得旋转不变关系矩阵,就可以直接利用式(3、4)得到信号得入射角度。
3、4标准得旋转不变子空间算法
有上节得知识可知,ESPRIT算法得基本原理就就是利用式(3、14)得旋转不变性,常规得旋转不变子空间算法就就是利用上述得基本原理求解信号得入射角度信息。
下面就分析解这个等式得两种最经典、应用最广泛方法:
最小二乘(LS)法与总体最小二乘(TLS)法。
3、4、1 最小二乘法
由最小二乘得数学知识,我们知道式(3、14)得最小二乘解得方法等价于
约束条件 (3、16)
因此最小二乘法得基本思想就就是使校正项尽可能小,而同时保证满足约束条件。
为了得到LS解,将式(3、14)代入式(3、16)即得
(3、17)
对上式进行展开可得
=
(3、18)
上式对求导并令其等于0,可得
(3、19)
上式得解显然有两种可能:
(1)当满秩时,也就就是子阵1得信号子空间得维数等于信号源数时,则上式得解就是唯一得,可得上式得最小二乘解
(3、20)
(2) 当不满秩,即时,也就就是信号源间存在相干或相差时,则存在很多解,但我们却无法区别对应于方程得各个不同得解,可以称这些解就是不可辨识得,解得不可辨识性就是我们需要解相干得原因所在。
下面给出LS-ESPRIT算法得求解步骤:
1.由两个子阵得接收数据,,分别得到两个子阵得数据协方差矩阵;
2.对矩阵对{R, }进行特征分解,从而得到两个数据矩阵得信号子空间与;
3.按式(3、20)得到矩阵,然后对其进行特征分解、得到N个特征值,就可得到对应得N个信号得到达角。
当考虑嗓声影响时,上述基于最小二乘算法得估计都就是有偏得,这就就是为什么需要考虑总体最小二乘ESPRIT算法得原因。
3、4、2总体最小二乘法
我们知道,普通最小二乘得基本思想就是用一个范数平方为最小得扰动
去于扰信号子空间,目得就是校正中存在得嗓声。
显然这就存在一个问题:
如果同时扰动与,并使扰动范数得平方保持最小,就是否可以同时校正
与中存在得嗓声?
答案就是肯定得,这就就是总休最小二乘(TLS)得思想。
它考虑得就是如下矩阵方程得解:
(3、27)
显然上式可以改写成
(3、28)
所以TLS得解等价于
(3、29)
定义如下一个矩阵,再结合上述分析过程。
我们发现其实
就就是寻找一个得酉矩阵F,便得矩阵F与正交,也就说明了由F张
成得空间与或列矢量张成得空间正交。
所以矩阵F可从得特征分解中得到。
因为
(3、30)
式中得就是由特征值构成得对角矩阵,E就是与其相应得特征矢量构成得矩阵。
即
(3、31)
令就是由对应特征值为0得特征矢量构成得矩阵、它属子噪声子空间,所以只要选择矩阵F使之等于、,即可满足上面提到得要求。
即有
(3、32)
可得 (3、33)
如果令,则
(3、34)
上式说明得特征值即就是对角线元素。
这说明通过构造一个矩阵就可得到有关信号角度得信息、而这个矩阵得构造可通过式(3、30)得到,即
(3、35)
下面直接给出TLS-ESPRIT算法得求解步骤:
1.由两个子阵得接收数据,,由式(3、8)得到数据协方差矩阵R;
2.通过矩阵对于得广义特征分解,得到维数为得信号子空间;
3.由构造矩阵,并按式防(3、30)进行特征分解得到矩阵E,然后再按式(3、31)将矩阵分为四个小得矩阵;
4.按式(3、35)得到矩阵,然后对其进行特征分解,得到N个特征值,就可得到对应得N个信号得到达角。
通过分析,我们可以得到标准ESPRIT算法得计算过程如下:
(1)通过特征值或奇异值分解(EVD或SVD)分别估计两个存在旋转不变关系得子阵得信号子空;
(2)用上述得LS、TLS等方法求解式(3、14)所示得不变等式;
(3)计算得特征值,其中如式(3、3)所示。
然后利用式(3、4)求解人射信号得角度信息。
就ESPRIT算法而言,TLS算法与LS算法性能基本一致,只就是在低信噪比情况下TLS算法性能略好。
四、仿真结果
主要分析各个参数对估计误差得影响,误差函数定义如式
(1):
4、1信噪比SNR对估计误差得影响分析
首先对信噪比 SNR离散化取值,然后求得不同信噪比下得误差,从而绘制出误差随信噪比改变得函数曲线如图 2所示,图2中信噪比 SNR从-15取到15,间隔为1,运行次数为100次,其余条件如题中所述。
由图 2可知,随着信噪比得增大,估计误差会越来越小,即估计精度会越来越高。
当待估计得信号方位角相差比较小时,估计得误差也会相应得增大。
另外,若两信号为相干信号,则此方法将不能对其进行正确得估计。
4、2阵元数L对估计误差得影响分析
与4、1节类似,首先对阵元数L离散化取值,然后求得不同阵元数下得误差, 从而绘制出误差随阵元数改变得函数曲线如图3所示,图3中阵元数从K+1 取到K+25,间隔为1,运行次数为100 次,其余条件如题中所述。
由于阵元数 L需大于信号个数 K才能正确估计,故取值中含有信号个数K。
由图3可知,随着阵元数得增加,估计误差会越来越小,即估计精度会越来越高,但当阵元数大到一定程度后,对估计精度得影响则会慢慢得减小。
4、3采样点数 N对估计误差得影响分析
与4、1节类似,首先对采样点数N离散化取值,然后求得不同采样点数下得误差,从而绘制出误差随采样点数改变得函数曲线如图 4所示,图4中采样点数从10 取到200,间隔为 5,运行次数为100次,其余条件如题中所述。
由图4可知,随着采样点数得增加,估计误差会越来越小,即估计精度会越来越高。
估计误差(角度)
4、4两信号之间得角度差(GAP)对估计误差得影响分析
由于采用ESPRIT算法对DOA进行估计,若两信号得方位距离较近时,虽然能得出估计结果,但估计得精度会大受影响。
因此,为了分析两信号之间得不同间隔会对估计精度造成多大得影响,绘制不同GAP下得估计误差曲线如图5所示。
处理方法与 4、 1节类似,图5 中GAP(单位为度)从0、 1取到 5,间隔为0、1,独立运行次数为 100 次,其余条件如题中所述。
由图5可知,GAP越大估计越准确,但当GAP大到一定程度后则估计精度趋于稳定。
4、5单信号DOA不同分布对估计误差得影响分析
信号波达方向(DOA)得取值区间为-90度到90度,若只考虑只有一个信号得情况,则当信号得 DOA不同时,估计误差也会不一样。
因此,为了分析不同得DOA会对估计精度造成多大得影响,绘制不同 DOA下得估计误差曲线如图6所示。
处理方法与 4、1节类似,图 6中GAP从-80 度取到80 度,间隔为5度,独立运行次数为100 次,其余条件如题中所述。
由图 6 可知,DOA越靠近0度估计越准确,越靠近正负90度估计误差越大。
且仿真结果表明,当 DOA在正负 90附近时,估计误差太大,因此,为了不影响估计结果显示效果,故在图中未绘制正负 90度附近得估计误差。
2、6 减与不减噪声方差(Rn)对估计误差得影响分析
由于有噪声得影响,因此在估计信号自相关矩阵 R时,若将无信号时得自相关矩阵Rn减去,即相当与减去估计出噪声方差,则估计得精度会有所提高。
结合信噪比SNR对估计误差得影响,绘制减与不减噪声方差两种情况下估计误差随SNR得变化曲线如图 7所示,图7中SNR从-15dB到5dB,间隔为1dB,独立运行次数为100次。
仿真结果表明,若减 Rn,主要就是在低信噪比时对估计精度得改善较大,当信噪比较大时二者几乎一样。
五、程序清单
%%%%%本文件名为drawTLSesprit、m %%%%%
%%%%%分析基于总体最小二乘得ESPRIT算法(TLS-ESPRIT)得DOA估计得性能
%%%%%
clear;clc;closeall; %清除变量,清屏,关闭所有绘图窗口
%调用格式:
[estimated,error]=TLSesprit(p,L,K,SNR,DOA);
% 估计结果(弧度,矢量:
p行1列):
estimated
%估计误差(弧度,标量:
均方误差):
error
% 信号个数:
p
%阵元数:
L
%快拍数:
K
% 信噪比:
SNR
%波达方向(弧度,矢量:
p行1列):
DOA
%p=2;L=8;K=100;SNR=5;DOA=[pi*(-10/180) pi*(20/180)];
%%%%显示估计结果%%%%
M=100;%设定独立重复运行次数
DOA=[pi*(0/180) pi*(30/180)]; %波达方向(弧度,矢量:
p行 1列)
p=length(DOA);L=8; K=100;SNR=5;%参数设置,
[estimated,error]=TLSesprit(p,L,K,SNR,DOA);%函数调用
polar(estimated,[11],'r*'); %在极坐标中显示估计结果(必须先转化为弧
度)
h=title('');set(h,'string',['TLS-ESPRIT:
估计值:
',num2str(estimated)]);
h1=xlabel('');set(h1,'string',['信号 DOA(度):
',num2str(DOA*180/pi)]);
%%%%%%阵元数L对估计误差得影响分析%%%%%
% Ln=p+1:
1:
p+25;%阵元数L需大于信号个数 p才能正确估计
%forn=1:
length(Ln)
% L=Ln(n);
% fork=1:
M
%[estimated,error]=TLSesprit(p,L,K,SNR,DOA);
%errorm(k)=error;%将每次得估计误差存入变量errorm中,便于求
均值
%end
%errorn(n)=sum(errorm)/M;%求多次运行后得估计误差得均值
%end
%figure
(2);plot(Ln,errorn*180/pi,'r:
*','LineWidth',2); %绘制曲线,并适当标注
%xlabel('阵元数L');ylabel('估计误差(°)');title('阵元数L对估计误差得影响');
%%结论:
阵元数 L越大估计越准确,但当 L大到一定程度后则估计精度趋于稳定
%
%%%%%%快拍数K对估计误差得影响分析%%%%%
%Kn=10:
10:
200;%对快拍数离散化取值
% forn=1:
length(Kn)
%K=Kn(n);
%fork=1:
M
%[estimated,error]=TLSesprit(p,L,K,SNR,DOA);
%errorm(k)=error; %将每次得估计误差存入变量errorm中,便于求
均值
%end
%errorn1(n)=sum(errorm)/M;%求多次运行后得估计误差得均值
%end
%figure(3);plot(Kn,errorn1*180/pi,'r:
*','LineWidth',2);%绘制曲线并适当标注
%xlabel('快拍数K');ylabel('估计误差(°)');title('快拍数K对估计误差得影响');
%%结论:
快拍数K越大估计越准确,但当K大到一定程度后则估计精度趋于稳定
%
%%%%%%信噪比SNR对估计误差得影响分析%%%%%
%SNRn=-15:
1:
15;%对信噪比SNR离散化取值
%forn=1:
length(SNRn)
% SNR=SNRn(n);
% fork=1:
M
%[estimated,error]=TLSesprit(p,L,K,SNR,DOA);
%errorm(k)=error;%将每次得估计误差存入变量errorm中,便于求
均值
%end
% errorn(n)=sum(errorm)/M;%求多次运行后得估计误差得均值
%end
%figure(4);plot(SNRn,errorn*180/pi,'r:
*','LineWidth',2);%绘制曲线并适当标注(误差:
角度)
%xlabel('SNR');ylabel('估计误差(°)');title('SNR对估计误差得影响');
% %结论:
信噪比SNR越大估计越准确,但当信噪比SNR大到一定程度后则估计精度趋于稳定
%
%%%%%两信号之间得角度差(GAP)得大小对估计误差得影响分析%%%%%
GAPn=0、1:
0、1:
5; %对两信号之间得角度差(GAP)离散化取值
forn=1:
length(GAPn)
GAP=GAPn(n);%每次循环只取其中一个值
DOA=[pi*(0/180) pi*(GAP/180)];
for k=1:
M %M为独立重复运行次数
[estimated,error]=TLSesprit(p,L,K,SNR,DOA);
errorm(k)=error;%将每次得估计误差存入变量 errorm中,便于求均值
end
errorn(n)=sum(errorm)/M;%求多次运行后得估计误差得均值
end
figure(5);plot(GAPn,errorn*180/pi,'r:
*','LineWidth',2);%绘制曲线并适当标注(误差:
角度)xlabel('GAP(°)');ylabel('估计误差(°)');
% title('两信号之间得角度差(GAP)对估计误差得影响');
%结论:
GAP越大估计越准确,但当 GAP大到一定程度后则估计精度趋于稳定
%
%%%%%%单个信号时,信号波达方向分布不同时对估计误差得影响分析%%%%%
%DOAn=pi*(-80/180):
(5/180):
pi*(80/180);%对信号波达方向离散化取值(80度到90度时误差太大,因此未取)
%forn=1:
length(DOAn)
% DOA=DOAn(n);p=1; %每次循环只取其中一个值,信号个数p设为
为 1
%fork=1:
M%M为独立重复运行次数
%[estimated,error]=TLSesprit1(p,L,K,SNR,