MVDR算法.docx
- 文档编号:7394736
- 上传时间:2023-05-11
- 格式:DOCX
- 页数:2
- 大小:12.51KB
MVDR算法.docx
《MVDR算法.docx》由会员分享,可在线阅读,更多相关《MVDR算法.docx(2页珍藏版)》请在冰点文库上搜索。
%%%%%%非相干信号,矢量阵,
clc;
closeall;
clearall;
M=10;%阵元数
m=[0:
M-1]';
N=1000;%采样点数
n=0:
N-1;
f=1000;%频率
fs=10*f;%采样频率
c=1500;%声速
lamd=c/f;%波长
%d=lamd/2;%阵元间距
d=2*lamd;
sita0=80*pi/180;%期望方向
sita1=80*pi/180;%目标1方位
sita2=100*pi/180;%目标2方位
sita3=225*pi/180;%目标3方位
snr1=20;%信噪比
snr2=20;
snr3=-10;
s1=exp(-j*(2*pi*f*n/fs+2*pi*rand(1,N)));
s2=exp(-j*(2*pi*f*n/fs+2*pi*rand(1,N)));
s3=exp(-j*(2*pi*f*n/fs+2*pi*rand(1,N)));
s=[s1;s2;s3];
a0=exp(-j*2*pi*f*m*d*cos(sita0)/c);
u0=[1cos(sita0)sin(sita0)].';
Av0=kron(a0,u0);
a1=exp(-j*2*pi*f*m*d*cos(sita1)/c);
a2=exp(-j*2*pi*f*m*d*cos(sita2)/c);
a3=exp(-j*2*pi*f*m*d*cos(sita3)/c);
u1=[1cos(sita1)sin(sita1)].';
u2=[1cos(sita2)sin(sita2)].';
u3=[1cos(sita3)sin(sita3)].';
Av=[kron(a1,u1)kron(a2,u2)kron(a3,u3)];
%%%有噪声
S1=awgn(Av(:
1)*s1,snr1,'measured','db');
S2=awgn(Av(:
2)*s2,snr2,'measured','db');
S3=awgn(Av(:
3)*s3,snr3,'measured','db');
S=S1+S2;%+S3;
%%%无噪声
%S=Av*s;%无噪声
R=S*S'/N;
[V,D]=eig(R);
Vm=V(:
1:
25);
IR=inv(R);
Wmvdr=IR*Av0/(Av0'*IR*Av0);
step=0.1;
Ang=0:
step:
360;
A=zeros(3*M,0);
forkk=1:
length(Ang)
u=[1cos(Ang(kk)*pi/180)sin(Ang(kk)*pi/180)]';
a=kron(exp(-i*(2*m*pi*f*d*cos(Ang(kk)*pi/180)/c)),u);
PP(kk)=1/(a'*IR*a);%MVDR方位谱
end
figure
(1);
plot(Ang,10*log10(abs(PP)/max(abs(PP))));%,'LineWidth',2);
gridon;
axis([0360-350]);
xlabel('方位(^o)');
ylabel('归一化功率(dB)');
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MVDR 算法