MVDR算法matlab程序.doc

MVDR算法matlab程序MVDR算法matlab程序MVDR算法matlab程序clcclearallcloseall％%常量定义Freqs=1。6e9；％工作频率c=3e8;％光速lamda=c/Freqs；%波长d=0。5*lamda;％单元间距M=16;%天线阵元数fs=2e6；％采样频率pd=10;%快拍数%%模型建立％—-—----—----—-第一个干扰模型—-——--—---------———-thetaJ1=20＊pi/180；％干扰方向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:pd）/fs）);%1*pd%---—----——----空域处理-—-—-—————-———-——-—-----—I1=zeros（M，1）；I2=zeros(M，1)；IS=zeros（M,1);forn=1:MI1(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）*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;forjj=1：1：length(p）％空间角度扫描w_scan=zeros（M,1）;scan=p(jj)；forn=1：1:Mw_scan(n）=exp（j＊2＊pi＊(n—1）*d*sin(scan）/lamda)；end%空间导向矢量1*MFmvdr(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-’)