MVDR算法.docx

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）'）;