clear all
close all
derad=pi/180; %角度->弧度
radeg=180/pi; %弧度->角度
twpi=2*pi;
kelm=8; %阵元数
dd=0.5; %阵元间距
d=0:dd:(kelm-1)*dd;
iwave=3; %信源数
theta=[-50 20 60]; %波达方向
snr= 10; %信噪比
n=500; %采样数
A=exp(-j*twpi*d.'*sin(theta*derad)); %方向矢量
S=randn(iwave,n); %信源信号
X=A*S; %接收信号
Xl=awgn(X,snr, 'measured'); %添加信号
Rxx=Xl*Xl'/n; %计算协方差矩阵
InvS=inv(Rxx);
[EV,D]=eig(Rxx); %特征值分解
EVA=diag(D)';
[EVA,I]=sort(EVA); %特征值从小到大排序
EVA=fliplr(EVA); %左右翻转,从大到小排序
EV=fliplr(EV(:,I)); %对应特征矢量排序
%构造MUSIC谱函数
for iang=1:361
angle(iang)=(iang-181)/2;
phim=derad*angle(iang);
a=exp(-j*twpi*d*sin(phim)).';
本内容试读结束,登录后可阅读更多
下载后可阅读完整内容,剩余1页未读,立即下载