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

4单元均匀线阵自适应波束形成图clearclcformat 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=0.8; %信号方向sita1=0.4; %干扰方向1sita2=2.1; %干扰方向2for t=1:N adt(t)=sin(omiga0*t/(N*f0)); a1t(t)=sin(omiga1*t/(N*f1)); a2t(t)=sin(omiga2*t/(N*f2));endfor 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));endR=zeros(M,M);for t=1:N x=adt(t)*ad+a1t(t)*a1+a2t(t)*a2; %阵列对信号的完整响应 R=R+x*x;%信号的协方差矩阵endR=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;endy_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单元均匀线阵自适应波束形成目标clearclcformat long;v=1;M=4;阵元数N=100;f0=21*10^3;omiga0=2*pi*f0;sita0=0.6;%信号方向 for t=1:N adt(t)=sin(omiga0*t/(N*f0));endfor i=1:M ad(i,1)=exp(j*(i-1)*pi*sin(sita0));endR=zeros(4,4);r=zeros(4,1);for t=1:N x=adt(t)*ad; R=R+x*x.;endR=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;endsita=0:pi/100:pi/2;plot(sita,y)xlabel(sita)ylabel(天线增益’)这是程序全部，有的比较简单的部分就可以省事的，如果解释的好，再加分的，大大的。谢谢了我有更好的答案提问者采纳2010-05-04 18:08你这里有两个程序，第二个程序与第一个实质上是一样的，区别就是信号与导向矢量的写法有点不同，这里我就不注释了。还有，我下面附了一段我自己的写的程序，里面有SIM算法。G-S正交化算法等。是基于圆阵形式的，你的算法是基于线阵的，他们程序上的区别在于导向矢量的不同。我的算法是某项目中的，保证好使。建议学习波束形成技术，注意把程序分块，例如分成，求导向矢量；最优权值；形成波束等等。程序如下：4单元均匀线阵自适应波束形成图clearclcformat 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=0.8; %信号方向sita1=0.4; %干扰方向1sita2=2.1; %干扰方向2for t=1:N %%%%%%%%%%%%信号 adt(t)=sin(omiga0*t/(N*f0)); a1t(t)=sin(omiga1*t/(N*f1)); a2t(t)=sin(omiga2*t/(N*f2));endfor 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));endR=zeros(M,M);for t=1:N x=adt(t)*ad+a1t(t)*a1+a2t(t)*a2; %阵列对信号的完整响应 R=R+x*x;%