MVDR算法matlab程序.doc

《MVDR算法matlab程序.doc》由会员分享，可在线阅读，更多相关《MVDR算法matlab程序.doc（9页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、Four short words sum up what has lifted most successful individuals above the crowd: a little bit more.-author-dateMVDR算法matlab程序clcclcclear allclose all % 常量定义Freqs=1.6e9; %工作频率c=3e8; %光速lamda=c/Freqs; %波长d=0.5*lamda; %单元间距M=16; %天线阵元数fs=2e6; %采样频率pd=10; %快拍数 % 模型建立%-第一个干扰模型-thetaJ1=20*pi/180; %干扰方

2、向FreqJ1=5e5; %第一个干扰的频率J1NR=sqrt(10(60/10);J1one=J1NR*exp(j*(2*pi*FreqJ1*(1:1:pd)/fs); %1*pd%-第二个干扰模型-ThetaJ2=60*pi/180; %干扰方向FreqJ2=6e5; %第二个干扰的频率J2NR=sqrt(10(60/10);J2one=J2NR*exp(j*(2*pi*FreqJ2*(1:1:pd)/fs); %1*pd%-信号模型-ThetaS=30*pi/180; FreqS=7e5;SNR=sqrt(10(40/10);Sone=SNR*exp(j*(2*pi*FreqS*(1:1

3、:pd)/fs); %1*pd%-空域处理-I1=zeros(M,1);I2=zeros(M,1);IS=zeros(M,1); for n=1:M I1(n)=exp(j*2*pi*(n-1)*d*sin(thetaJ1)/lamda); I2(n)=exp(j*2*pi*(n-1)*d*sin(ThetaJ2)/lamda); IS(n)=exp(j*2*pi*(n-1)*d*sin(ThetaS)/lamda); endJ1=I1*J1one;J1=J1.;J2=I2*J2one;J2=J2.;%-产生噪声-noise=sqrt(1/2)*randn(pd,M)+j*sqrt(1/2)*

4、randn(pd,M);%-产生接收信号-XX=J1+noise+J2; % MVDR算法Rxx=zeros(M,M);Rxx=(XX.*(XX.)/pd;%M*M维detaS=IS;Wmvdr=inv(Rxx)*detaS*inv(detaS*inv(Rxx)*detaS); %MVDR 算法 % 直线阵布阵p=-pi/2:pi/180:pi/2; for jj=1:1:length(p) %空间角度扫描 w_scan=zeros(M,1); scan=p(jj); for n=1:1:M w_scan(n)=exp(j*2*pi*(n-1)*d*sin(scan)/lamda); end %空间导向矢量 1*M Fmvdr(jj)=abs(Wmvdr*w_scan); end % 画出方向图Fmax=max(max(Fmvdr);Fa=Fmvdr./Fmax;Fa=20*log10(Fa);figure(1);plot(180*p/pi,Fa,b-)-