siyuanshiliangzhen.m

矢量阵指向性程序

clear;
f0=1000;   %信号频率
fs=80*f0;   %采样频率 
Ns=2048;   %采样点数 
N1=5;      %阵元数
c=1500;    %声速
r=1;      %阵元半径 
fai=0;     %信号入射方向 
w=2*pi*f0;
ph1=-180;                   
ph2=180;                   
phn=1;                  
ax=(ph2-ph1)/phn+1;

for n=1:N1;
    for k=1:Ns;
        p(n,k)=exp(j*w*(k-1)/fs-j*w*r*cos(fai*pi/180-(n-1)*2*pi/N1)/c);
        vx(n,k)=p(n,k)*cos(fai*pi/180);
        vy(n,k)=p(n,k)*sin(fai*pi/180);  
    end
end

for i=1:ax;     
    num=ph1+(i-1)*phn;
    fain=num*pi/180;
     Xp1=0;
    Xv1=0;
    for n=1:N1;
        Xp(n,1,i)=p(n,1)*exp(j*w*r*cos(fain-(n-1)*2*pi/N1)/c); 
        Xp1=Xp1+Xp(n,1,i);
        Xv(n,1,i)=Xp(n,1,i)*(cos(fai*pi/180)*cos(fain)+sin(fai*pi/180)*sin(fain));
        Xv1=Xv1+Xv(n,1,i);
    end
    Rp1(i)=abs(Xp1*Xv1);
  Rp2(i)=abs((Xp1+Xv1)*Xv1);
  Rp3(i)=abs((Xp1+Xv1).^2);
  Rp4(i)=abs((Xp1+2*Xv1)*Xv1);
  Rp5(i)=abs((Xp1+2*Xv1).^2);
  Rp6(i)=abs((Xp1+2*Xv1)*(Xp1+Xv1));
end 
figure
subplot(321)
i=1:ax;
maxRp1=max(Rp1);
Rpp1=20*log10(Rp1(i)/maxRp1);
plot(ph1+(i-1)*phn,Rpp1);
grid on
hold on
axis([-180 180 -100 0])
title('四元矢量阵(pv)')
xlabel('角度')
ylabel('dB')
subplot(322)
i=1:ax;
maxRp2=max(Rp2);
Rpp2=20*log10(Rp2(i)/maxRp2);
plot(ph1+(i-1)*phn,Rpp2);
grid on
hold on
axis([-180 180 -100 0])
title('四元矢量阵((p+v)v)')
xlabel('角度')
ylabel('dB')
subplot(323)
i=1:ax;
maxRp3=max(Rp3);
Rpp3=20*log10(Rp3(i)/maxRp3);
plot(ph1+(i-1)*phn,Rpp3);
grid on
hold on
axis([-180 180 -100 0])
title('四元矢量阵((p+v)^2)')
xlabel('角度')
ylabel('dB')
subplot(324)
i=1:ax;
maxRp4=max(Rp4);
Rpp4=20*log10(Rp4(i)/maxRp4);
plot(ph1+(i-1)*phn,Rpp4);
grid on
hold on
axis([-180 180 -100 0])
title('四元矢量阵((p+2v)v)')
xlabel('角度')
ylabel('dB')
subplot(325)
i=1:ax;
maxRp5=max(Rp5);
Rpp5=20*log10(Rp5(i)/maxRp5);
plot(ph1+(i-1)*phn,Rpp3);
grid on
hold on
axis([-180 180 -100 0])
title('四元矢量阵((p+2v)^2)')
xlabel('角度')
ylabel('dB')
subplot(326)
i=1:ax;
maxRp6=max(Rp6);
Rpp6=20*log10(Rp6(i)/maxRp6);
plot(ph1+(i-1)*phn,Rpp6);
grid on
hold on
axis([-180 180 -100 0])
title('四元矢量阵((p+2v)(p+v)')
xlabel('角度')
ylabel('dB')