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

单元均匀线阵自适应波束形成图.doc4单元均匀线阵自适应波束形成图clearclcformatlong;v=1;M=4;N=1000;f0=21*10A3;f1=11*10A3;f2=15*10A3;omiga0=2*pi*f0;omiga1=2*pi*f1;omiga2=2*pi*f2;sita0=;%信号方向sita1=;%干扰方向1sita2=;%干扰方向2fort=1:Nadt(t)=sin(omigaO*t/(N*fO));a1t(t)=sin(omiga1*t/(N*f1));a2t(t)=sin(omiga2*t/(N*f2));endfori=1:Mad(i,1)=exp(j*(i-1)*pi*sin(sitaO));a1(i,1)=exp(j*(i-1)*pi*sin(sita1));a2(i,1)=exp(j*(i-1)*pi*sin(sita2));endR=zeros(M,M);fort=1:Nx=adt(t)*ad+a1t(t)*a1+a2t(t)*a2;%阵列对信号的完整响应R=R+xW;%信号的协方差矩阵endR=R/N;miu=1/(ad,*inv(R)*ad);w=miu*inv(R)*ad;forsita=0:pi/100:pifori=1:Mxji,1)=exp(j*(i-1)*pi*sin(sita)):endy(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单元均匀线阵自适应波束形成目标clearclcformatlong;v=1;M=4;阵元数N=100;f0=21*10A3;omiga0=2*pi*f0;sita0=;%信号方向fort=1:Nadt(t)=sin(omigaO*t/(N*fO));endfori=1:Mad(i,1)=exp(j*(i-1)*pi*sin(sitaO));endR=zeros(4,4);r=zeros(4,1);fort=1:Nx=adt(t)*ad;R=R+x*x.';endR=R/N;miu=1/(ad/*inv(R)*ad);w=miu*inv(R)*ad;forsita=0:pi/100:pi/2fori=1:Ma(i,1)=exp(j*(i-1)*pi*sin(sita));endy(1,v)=w/*a;v=v+1;endsita=0:pi/100:pi/2;plot(sita,y)xlabel('sita')ylabel('天线增益')这是程序全部,有的比较简单的部分就町以省事的,如果解释的好,再加分的,大大的。谢谢了我有更好的答案提问者采纳2010-05-0418:08你这里有两个程序,第二个程序与第一个实质上是一样的,区别就是信号与导向矢量的写法有点不同,这里我就不注释了。还有,我下面附了一•段我自己的写的程序,里面有SIM算法。G・S正交化算法等。是基于圆阵形式的,你的算法是基于线阵的,他们程序上的区别在于导向矢量的不同。我的算法是某项目中的,保证好使。建议学习波束形成技术,注意把程序分块,例如分成,求导向矢量;最优权值;形成波束等等。程序如下:4单元均匀线阵自适应波束形成图clearclcformatlong;v=1;M=4;N=1000;%%%%%%%快拍数f0=21*10A3;%%%%%%%%%%%信号与干扰的频率f1=11*10A3;f2=15*10A3;omiga0=2*pi*f0;%%%%%%%信号与干扰的角频率omiga1=2*pi*f1;omiga2=2*pi*f2;sita0=;%信号方向sita1=;%干扰方向1sita2=;%干扰方向2fort=1:N %%%%%%%%%%%%信号adt(t)=sin(omigaO*t/(N*fO));a1t(t)=sin(omiga1*t/(N*f1));a2t(t)=sin(omiga2*t/(N*f2));endfori=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);fort=1:Nx=adt(t)*ad+a1t(t)*a1+a2t(t)*a2;%阵列对信号的完整响应R=R+x*x,;%信号的协方差矩阵endR=R/N;%%%%%%%%%协方差矩阵,所有快拍数的平均miu=1/(ad,*inv(R)*ad);%%