資源簡(jiǎn)介
通過存在幅相誤差與未存在幅相誤差時(shí)的MUSIC譜峰比較,可看出幅相誤差對(duì)DOA估計(jì)有很大的影響。此程序中加入的幅相誤差參數(shù)服從高斯分布。
代碼片段和文件信息
clear?all
N=200;%%快拍數(shù)
doa=[20?40]/180*pi;?%%信號(hào)到達(dá)角,
w=[pi/4?pi/3]‘;%%信號(hào)頻率
M=8;%%陣元數(shù)
P=length(doa);?%%信號(hào)個(gè)數(shù)
lma=150;%波長
d=lma/2;%陣元間距
snr=15;%信噪比
%%%導(dǎo)向向量
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);%信號(hào)幅度
S=diag(Am);?
x=G*x*S+Noise;%快拍數(shù)據(jù)矢量
R=x*x‘/N;?%數(shù)據(jù)協(xié)方差矩陣
[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
?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件???????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
評(píng)論
共有 條評(píng)論