MVDR算法matlab程序mvdr算法代码_阵列干扰

clcclear allcloseall%%常量定义・Freqs=l6e9;c=3e8;%工作频率lamda=c/Freqs;d=

0.5*lamda;%光速M=16;fs=2e6;%波长pd=10;%单元间距%天线阵元数%%模型建立%采样频率%--------------------------第一个干扰模型---------------------------------------%快拍数%干扰方向theta J1=20*pi/180;%第一个干扰的频率FreqJl=5e5;AJlNR=sqrt1060/10;Jlone=JlNR*expj*2*pi*FreqJl*1:1:pd/fs;%l*pd%--------------------------第二个干扰模型---------------------------------------%干扰方向ThetaJ2=60*pi/180;%第二个干扰的频率FreqJ2=6e5;八J2NR=sqrt1060/10;/fs;J2one=J2NR*expj*2*pi*FreqJ2*1:1:pd%l*pd%--------------------------信号模型--------------------------------------ThetaS=30*pi/180;FreqS=7e5;人SNR=sqrt1040/10;Sone=SNR*expj*2*pi*FreqS*1:1:pd/fs;%l*pd%--------------------------空域处理------------------------------------------------Il=zerosM,1;I2=zerosM,1;IS=zerosM,1;for n=l:MIl n=exp j*2*pi*n-1*d*sinthetaJI/lamda;12n=exp j*2*pi*n-1*d*sinThetaJ2/lamda;IS n=exp jn-1*d*sinThetaS/lamda;endJl=Il*Jlone;1J1=J

1.;J2=I2*J2one;J2=J

2.;%----------------------产生噪声------------------------------------------------noise=sqrt1/2*randnpd,M+j*sqrt1/2*randnpd,M;%----------------------产生接收信号------------------------------------------------XX=Jl+noise+J2;算法%%MVDRRxx=zerosM,M;维1f TRxx=XX.*XX./pd;%M*MdetaS=IS;算法1Wmvdr=inv Rxx*detaS*inv detaS*inv Rxx*detaS;%MVDR%%直线阵布阵p=-pi/2:pi/180:pi/2;%空间角度扫描for jj=l:1:length pw_scan=zerosM,1;scan=p jj;for n=l:1:Mw_scann=expj*2*pi*n-1*d*sinscan/lamda;%空间导向矢量end1Fmvdrj j=absWmvdr*w_scan;end%%画出方向图Fmax=maxmaxFmvdr;Fa=Fmvdr./Fmax;Fa=20*logl0Fa;figure1;plot180*p/pi,Fa,。