基于matlab实现稳健的自适应波束形成算法，常用mse准则算法.rar

%ASC旁瓣相消----MSE准则 clc;close all;clear all ima=sqrt(-1); M=32; %辅助天线的数目 d_lamda=.5; theta0=-30; %来波方向 theta1=60; %干扰方向 L=512; %采样单元数 s=zeros(1,512); %预划分一个区域 for ii=1:L; amp0=1*randn(1); %信号的幅度随机产生，保证信号之间是不相关的 amp1=200*randn(1); ampn=1; jam(:,ii)=amp1*exp(ima*2*pi*0.5*sin(theta1*pi/180)*[0:M-1]')+ampn*(randn(M,1)+ima*randn(M,1)); %干扰+噪声 s(ii)=amp0*exp(ima*2*pi*0.5*sin(theta0*pi/180))+amp1*exp(ima*2*pi*0.5*sin(theta1*pi/180))+ampn*(randn(1,1)+ima*randn(1,1));%接收信号（信号+干扰+噪声） s0(ii)=amp0*exp(ima*2*pi*0.5*sin(theta0*pi/180)); end Rx=1/L*jam*jam'; %噪声自相关矩阵，相当于X（t） r_xd=1/L*jam*s'; Wopt=pinv(Rx)*r_xd; delta=s0-(s-Wopt'*jam); delta1=abs(mean(delta.^2)-(mean(delta)).^2) %方差 theta=linspace(-pi/2,pi/2,200); for jj=1:length(theta) a=exp(ima*2*pi*.5*sin(theta(jj))*[0:M-1]'); f(jj)=Wopt'*a; end F=20*log10(abs(f)/(max(max(abs(f))))); figure(1) plot(theta*180/pi,F),grid on,hold on plot(theta0,-50:0,'.') plot(theta1,-50:0,'.') xlabel('theta/o'); ylabel('F/dB'); title('MSE准则下的方向图') axis([-90 90 -50 0]);%可为x轴和y轴设置一个极限范围