yuanzhenzhixiangxing2.m

矢量阵指向性程序

clear;
f0=1000;
fs=80*f0;
Ns=2048;
N1=8;
r1=1;
r2=0.3
c=1500;
fai=0;
w=2*pi*f0;
ph1=-180;                   
ph2=180;                   
phn=1;                  
ax=(ph2-ph1)/phn+1;

for n=1:N1;
       a=2;
    if mod(n,a)~=0
       for k=1:Ns;
        p1(n,k)=exp(j*w*(k-1)/fs-j*w*r1*cos(fai*pi/180-(n-1)*2*pi/N1)/c);
        vx1(n,k)=p1(n,k)*cos(fai*pi/180);
        vy1(n,k)=p1(n,k)*sin(fai*pi/180);  
       end
    else
       for k=1:Ns;
        p2(n,k)=exp(j*w*(k-1)/fs-j*w*r2*cos(fai*pi/180-(n-1)*2*pi/N1)/c);
        vx2(n,k)=p2(n,k)*cos(fai*pi/180);
        vy2(n,k)=p2(n,k)*sin(fai*pi/180);  
       end
    end
%     plot((1:Ns),s)
%     hold on
end

for i=1:ax;     
    num=ph1+(i-1)*phn;
    fain=num*pi/180;
     Xp1=0;
    Xv1=0;
    for n=1:N1;
        a=2;
      if mod(n,a)~=0
        Xp01(n,1,i)=p1(n,1)*exp(j*w*r1*cos(fain-(n-1)*2*pi/N1)/c); 
        Xp1=Xp1+Xp01(n,1,i);
        Xv01(n,1,i)=Xp01(n,1,i)*(cos(fai*pi/180)*cos(fain)+sin(fai*pi/180)*sin(fain));
        Xv1=Xv1+Xv01(n,1,i);
      else
        Xp02(n,1,i)=p2(n,1)*exp(j*w*r2*cos(fain-(n-1)*2*pi/N1)/c); 
        Xp1=Xp1+Xp02(n,1,i);
        Xv02(n,1,i)=Xp02(n,1,i)*(cos(fai*pi/180)*cos(fain)+sin(fai*pi/180)*sin(fain));
        Xv1=Xv1+Xv02(n,1,i);
      end
    end
    Rp(i)=abs(Xp1*Xv1);
%       Rp(i)=abs((Xp1+Xv1)*Xv1);
%         Rp(i)=abs((Xp1+Xv1).^2);
%           Rp(i)=abs((Xp1+2*Xv1)*Xv1);
%             Rp(i)=abs((Xp1+2*Xv1).^2);
%               Rp(i)=abs((Xp1+2*Xv1)*(Xp1+Xv1));
end 
i=1:ax;
maxRp=max(Rp);
Rpp=20*log10(Rp(i)/maxRp);
figure
plot(ph1+(i-1)*phn,Rpp);
grid on
hold on
axis([-180 180 -100 0])
xlabel('角度')
ylabel('dB')