天线阵代码
http://www.pudn.com/downloads164/
sourcecode/math/detail750575.html
一、
clc
clear all
f=3e9;
N1=4;N2=8;N3=12;
a=pi/2; %馈电相位差
i=1; %天线电流值
lambda=(3e8)/f; %lambda=c/f 波长
d=lambda/2;
beta=2.*pi/lambda;
W=-2*pi:0.001:2*pi;
y1=sin((N1.*W./2))./(N1.*(sin(W./2))); %归一化阵因子
y1=abs(y1);
r1=max(y1);
y2=sin((N2.*W./2))./(N2.*(sin(W./2))); %归一化阵因子
y2=abs(y2);
r2=max(y2);
y3=sin((N3.*W./2))./(N3.*(sin(W./2))); %归一化阵因子
y3=abs(y3);
r3=max(y3);
%归一化阵因子绘图程序,
&gure(1)
subplot(311);plot(W,y1) ; grid on; %绘出 N=4 等幅等矩阵列的归一化
阵因子
xlabel('f=3GHz,N=4,d=1/2 波长,a=π/2')
subplot(312);plot(W,y2) ; grid on; %绘出 N=8 等幅等矩阵列的归一化
阵因子
xlabel('f=3GHz,N=8,d=1/2 波长,a=π/2')
subplot(313);plot(W,y3) ; grid on; %绘出 N=12 等幅等矩阵列的归一
化阵因子
xlabel('f=3GHz,N=12,d=1/2 波长,a=π/2')
%---------------------
%只有参数 N 改变的天线方向图
t=0:0.01:2*pi;
W=a+(beta.*d.*cos(t));
z1=(N1/2).*(W);
z2=(1/2).*(W);
W1=sin(z1)./(N1.*sin(z2)); %非归一化的阵因子 K1