单元均匀线阵自适应波束形成图.doc

------------------------------------------------------------------------------------------------ ——————————————————————————————————————单元均匀线阵自适应波束形成图 clear clc format long; v=1; M=4; N=1000; f0=21*10^3; f1=11*10^3; f2=15*10^3; omiga0=2*pi*f0; omiga1=2*pi*f1; omiga2=2*pi*f2; sita0=; % 信号方向 sita1=; % 干扰方向 1 sita2=; % 干扰方向 2 for t=1:N adt(t)=sin(omiga0*t/(N*f0)); a1t(t)=sin(omiga1*t/(N*f1)); a2t(t)=sin(omiga2*t/(N*f2)); end ------------------------------------------------------------------------------------------------ —————————————————————————————————————— for i=1:M ad(i,1)=exp(j*(i-1)*pi*sin(sita0)); a1(i,1)=exp(j*(i-1)*pi*sin(sita1)); a2(i,1)=exp(j*(i-1)*pi*sin(sita2)); end R=zeros(M,M); for t=1:N x=adt(t)*ad+a1t(t)*a1+a2t(t)*a2; % 阵列对信号的完整响应 R=R+x*x';% 信号的协方差矩阵 end R=R/N; miu=1/(ad'*inv(R)*ad); w=miu*inv(R)*ad; for sita=0:pi/100:pi for i=1:M x_(i,1)=exp(j*(i-1)*pi*sin(sita)); end y(1,v)=w'*x_; v=v+1; end y_max=max(y(:)); y_1=y/y_max; ------------------------------------------------------------------------------------------------ —————————————————————————————————————— y_db=20*log(y_1); sita=0:pi/100:pi; plot(sita,y) Xlabel (‘ sitaa ’) Ylabel (‘天线增益 db’) 4 单元均匀线阵自适应波束形成目标 clear clc format long; v=1; M= 4; 阵元数 N=100; f0=21*10^3; omiga0=2*pi*f0; sita0=;% 信号方向 for t=1:N adt(t)=sin(omiga0*t/(N*f0)); end for i=1:M ad(i,1)=exp(j*(i-1)*pi*sin(sita0)); end R=zeros(4,4); r=zeros(4,1); for t=1:N ------------------------------------------------------------------------------------------------ —————————————————————————————————————— x=adt(t)*ad; R=R+x*x.'; end R=R/N; miu=1/(ad.'*inv(R)*ad); w=miu*inv(R)*ad; for sita=0:pi/100:pi/2 for i=1:M a(i,1)=exp(j*(i-1)*pi*sin(sita)); end y(1,v)=w.'*a; v=v+1; end sita=0:pi/100:pi/2; plot(sita,y) xlabel('sita') ylabel(' 天线增益’) 这是程序全部,有的比较简单的部分就可以省事的, 如果解释的好,再加分的,大大的。谢谢了我有更好的答案提问者采纳------------------------------------------------------------------------------------------------ ————————————————————