窄带波束成形最小二乘方法matlab实现

clc
clear
f=2.4e9;
w=2*pi*f;
v=3e8;
lamda=v/f;
d=0.5*lamda;
%设置初始参数
doas=[90,40]*pi/180;%设置来波方向
k=length(doas);
m=10;
HC=[1,0]';%设置来波方向的绝对限制条件
HD=zeros(1,181)';%设置期望的方向图
tao=d*cos(doas)/v;%行向量1*2
tao2=d*cos((0:180)*pi/180)/v;%时间常数tao
ac1=zeros(1,10);%信号导向矢量
ac2=zeros(1,10);%信号导向矢量
for i=1:10
    ac1(i)=exp(1i*w*tao(1)*(i-1));
    ac2(i)=exp(1i*w*tao(2)*(i-1));
end
C=[ac1.',ac2.'];
ad = zeros(181,10);%信号导向矢量
for i=1:181
    for j=1:10
        ad(i,j)=exp(1i*w*tao2(i)*(j-1));
    end
end
D=ad.';
weight=inv(D*D')*(D*HD-C*inv(C'*inv(D*D')*C)*(C'*inv(D*D')*D*HD-HC));
Hwtheta=10*log(abs(weight'*D));
angles=0:1:180;
plot(angles,Hwtheta,'k');

来波方向为90度与40度，其中90度为信号方向，40度为噪声方向
结果如图所示