LMS.rar_lms波束_lms波束形成_lms算法_自适应波束形成算法

clear all close all clc M=9; %天线数 K=3; %信源数 theta=[15 -45 75]; %信号入射角度 d=0.3; %天线间距 N=500; %采样数 Meann=0; varn=1; %噪声均值、方差 SNR=10; %信噪比 INR=3; %干噪比 pp=zeros(100,N); pp1=zeros(100,N); rvar1=sqrt(varn)*10^(SNR/20); %信号功率 rvar2=sqrt(varn)*10^(INR/20); %干扰功率 rvar3=sqrt(varn)*10^(INR/20); for q=1:100 s=[rvar1*exp(j*2*pi*(50*0.001*[0:N-1]));rvar2*exp(j*2*pi*(100*0.001*[0:N-1]+rand));rvar3*exp(j*2*pi*(100*0.001*[0:N-1]+rand))]; %生成源信号 A=exp(-j*2*pi*d*[0:M-1].'*sin(theta*pi/180)); %方向矢量 e=sqrt(varn/2)*(randn(M,N)+j*randn(M,N)); %噪声 Y=A*s+e; %接收信号 %LMS算法 L=200; de=s(1,:); mu=0.001; w=zeros(M,1); for k=1:50 y(k)=w'*Y(:,k); %预测下一个采样和误差 e(k)=de(k)-y(k); %误差 w=w+mu*Y(:,k)*conj(e(k)); %调整权矢量 end end %波束形成 beam=zeros(1,L); for i=1:L a=exp(-j*2*pi*d*[0:M-1].'*sin(-pi/2+pi*(i-1)/L)); beam(i)=20*log10(abs(w'*a)); end %绘图 figure angle=-90:180/200:(90-180/200); plot(angle,beam); grid on xlabel('方向角/degree'); ylabel('幅度响应/dB'); figure for k=1:N en(k)=(abs(e(k)))^2; end %num=1:1:k %plot(num,en); semilogy(en); hold on xlabel('迭代次数'); ylabel('MSE');