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