xianzhen.m

来自「线列阵加权程序」· M 代码 · 共 27 行

N=12;zpb=40.5%主旁瓣
q=chbe(N,zpb)
c=1500;
f0=1500;
fs=50000;
nfft=512;
A=2;
d=0.4;
n=1:nfft;
w=2*pi*f0;
s=zeros(5,nfft);
y=[];
for x=-pi/2:pi/180:pi/2
    tao=d*sin(x)/c;
   for i=1:N
        s(i,:)=q(i)*A*cos(w*(n/fs+tao*(i-1)));
    end
     p=sum(s);
b =sqrt(sum(p.^2)) ;
y=[y,b];
end
x=-pi/2:pi/180:pi/2;
plot(x,y);
title('自然指向性');xlabel('入射信号角度');ylabel('波束输出');