幅相误差下的MUSIC谱峰比较

通过存在幅相误差与未存在幅相误差时的MUSIC谱峰比较，可看出幅相误差对DOA估计有很大的影响。此程序中加入的幅相误差参数服从高斯分布。

clear all
N=200;%%快拍数
doa=[20 40]/180*pi; %%信号到达角，
w=[pi/4 pi/3]‘;%%信号频率
M=8;%%阵元数
P=length（doa）; %%信号个数
lma=150;%波长
d=lma/2;%阵元间距
snr=15;%信噪比
%%%导向向量
B=zeros（PM）;
for k=1:P
B（k:）=exp（-j*2*pi*d*sin（doa（k））/lma*[0:M-1]）;
end
B=B‘;
 
xx=2*exp（j*（w*[1:N]））;
x=B*xx;
 
%%%增益和相位误差，

DeltaAm=0.2;%增益方差
DeltaPh=40;%相位方差度
Ame=[1 1+（rand（1M-1）-0.5）*DeltaAm*sqrt（12）]; %幅度误差
Phe=[0 （rand（1M-1）-0.5）*DeltaPh*sqrt（12）];  %相位误差，度
G=diag（Ame.*exp（j*Phe*pi/180））;   %幅相误差矩阵 

Noise=sqrt（0.5）*（randn（MN）+j*randn（MN））; %噪声

Am=10.^（snr/20）;%信号幅度
S=diag（Am）; 
x=G*x*S+Noise;%快拍数据矢量
R=x*x‘/N; %数据协方差矩阵

[UV]=eig（R）;

UU=U（:1:M-P）; 

theta=-90:0.5:90;
for ii=1:length（theta）
    AA=zeros（1length（M））;
        for jj=0:M-1
            AA（1+jj）=exp（-j*2*jj*pi*d*sin（theta（ii）/180*pi）/lma）;
        end
    PP=AA*UU*UU‘*AA‘;
    Pmusic（ii）=abs（1/PP）;
end
Pmusic=10*log10（Pmusic）;
figure（1）
plot（thetaPmusic）
xlabel（‘ \theta/°‘）
ylabel（‘P（\theta） /dB‘）
title（‘幅相误差影响‘）
grid on

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1108  2014-04-08 20:54  幅相误差影响\fuxiangwucha.m

     文件       1484  2014-04-08 15:52  幅相误差影响\fuxiangwucha_compare.asv

     文件       1511  2014-04-08 15:55  幅相误差影响\fuxiangwucha_compare.m

     文件       1483  2014-04-10 21:39  幅相误差影响\fuxiangwucha_gaosi.m

     目录          0  2014-04-12 16:50  幅相误差影响

----------- ---------  ---------- -----  ----

                 5586                    5