偶数平面阵切比雪夫矩阵及波束图形.m

相控阵雷达的重要性就不用多说了。本人仔细编写了相控阵雷达的仿真波束图

%%%%%%% 偶数平面阵切比雪夫矩阵及波束图形
clear,clc
N=5;PI=pi/180;
p1=zeros(180,180); 
p2=zeros(180,180);
sita=-90:89;DLLR=0.5;
fai=-90:89;
[sita,fai]=meshgrid(sita,fai);
u=pi*DLLR*(sin(sita*PI).*cos(sita*PI));
v=pi*DLLR*(sin(sita*PI).*sin(sita*PI));
R=0;
sll=20;
SLL=cosh(acosh(10^(sll/20))/(2*N-1));
% for p=1:N
%     R=R+configure(2*N-1,N-p);
% end
I=zeros(N,N);
for m=1:N
    for n=1:N
        if m>=n
            p=m;
        else 
           p=n;
        end
        for r=p:N
            I(m,n)=I(m,n)+(-1)^(N-r)*(2*N-1)*configure(N+r-1,2*r-1)*configure(2*r-1,r-m)*configure(2*r-1,r-n)*(SLL/2)^(2*r-1)/(2*(N+r-1));
        end
      a=cos(pi*DLLR*sin(sita*PI)*(2*m-1));
      b=cos(pi*DLLR*sin(sita*PI)*(2*n-1));
 %a=cos((2*m-1)*cos(u));
 %b=cos((2*n-1)*cos(v));
       p1=p1+I(m,n).*a'*b;%
    end
end
p1=abs(p1);C=max(max(p1));
p1=p1/C;  
p1=20*log10(p1); 
for ele=1:180,
    for aza=1:180,
       if p1(ele,aza)<-40,
           p1(ele,aza)=-40;
       end;
   end;
end;
% mesh(I);

% figure (2) 
mesh(-90:89,-90:89,p1);%view(90,0)
title('10x10阵元切比雪夫矩形阵加权')