对称阵子天线matlab
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
方向图
clear
lamda=10;%自由空间的波长
n=1.125;
l=n*lamda;
k=2*pi/lamda;%自由空间的相移常数
theta0=[0.0001:0.1:360];
theta=theta0*pi/180;
for i=1:length(theta0)
fe(i)=abs((cos(k*l*cos(theta(i)))-cos(k*l))/sin(theta(i)));
% fh(i)=1-cos(k*l);
end
% figure
polar(theta,fe/max(fe));
xlabel('l=1.125λ');
% figure
% polar(theta,fh/max(fh));
辐射阻抗
clear
lamda=10;%自由空间的波长
n=0.1:0.01:0.9;
l=n*lamda;
k=2*pi/lamda;%自由空间的相移常数
gama=0.5772156;
for i=1:length(n)
Rr1(i)=sin(2*k*l(i))*(Si(4*k*l(i))-2*Si(2*k*l(i)));
Rr2(i)=cos(2*k*l(i))*(gama+log(k*l(i))+Ci(4*k*l(i))-2*Ci(2*k*l(i)));
Rr3(i)=2*(gama+log(2*k*l(i))-Ci(2*k));
Rr(i)=(Rr1(i)+Rr2(i)+Rr3(i))*30;
end
plot(n,Rr)
阻抗特性
clc
clear
lamda=100;%自由空间的波长
n=0.1:0.01:0.9;
l=n*lamda;
a=[30,300,3000];%l和a的比值
k=2*pi/lamda;%自由空间的相移常数
gama=0.5772156;
for j=1:length(a)
for i=1:length(n)
Rr1(i)=sin(2*k*l(i))*(Si(4*k*l(i))-2*Si(2*k*l(i)));
Rr2(i)=cos(2*k*l(i))*(gama+log(k*l(i))+Ci(4*k*l(i))-2*Ci(2*k*l(i)));
Rr3(i)=2*(gama+log(2*k*l(i))-Ci(2*k));
Rr(i)=(Rr1(i)+Rr2(i)+Rr3(i))*30;
Xr1(i)=2*Si(2*k*l(i));
Xr2(i)=sin(2*k*l(i))*(gama+log(k*l(i))+Ci(4*k*l(i))-2*Ci(2*k*l(i))-2*log(a(j)));
Xr3(i)=cos(2*k*l(i))*(2*Si(2*k*l(i))-Si(4*k*l(i)));
Xr(i,j)=(Xr1(i)+Xr2(i)+Xr3(i))*30;
Zca(i)=120*(log(2*(a(j)))-1);
afa(i)=Rr(i)/(Zca(i)*l(i)*(1-sin(2*k*l(i))/(2*k*l(i))));
Rin(i,j)=Zca(i)*(sinh(2*afa(i)*l(i))-afa(i)/k*sin(2*k*l(i)))/(cosh(2*afa(i)*l(i))-cos(2* k*l(i)));
Xin(i,j)=-Zca(i)*(afa(i)/k*sinh(2*afa(i)*l(i))+sin(2*k*l(i)))/(cosh(2*afa(i)*l(i))-cos(2 *k*l(i)));
end
end
figure
plot(n,Rr)%辐射电阻
title('辐射电阻随l/λ的变化曲线');
xlabel('l/λ');
ylabel('Rr (Ω)');
figure
plot(n,Xr)%辐射电抗
title('l/a不同时辐射电抗随l/λ的变化曲线');
xlabel('l/λ');
ylabel('Xr (Ω)');
figure
plot(n,Rin)%输入电阻
title('不同平均特性阻抗时输入电阻随l/λ的变化曲线');
xlabel('l/λ');
ylabel('Rin (Ω)');
figure
plot(n,Xin)%输入电抗
title('不同平均特性阻抗时电抗随l/λ的变化曲线');
xlabel('l/λ');
ylabel('Xin (Ω)');
对称阵子
clc
clear
lamda=10;%自由空间的波长
n=0.25;
l=n*lamda;
k=2*pi/lamda;%自由空间的相移常数
theta0=[0.0001:0.1:360];
theta=theta0*pi/180;
for i=1:length(theta0)
fe(i)=abs((cos(k*l*cos(theta(i)))-cos(k*l))/sin(theta(i)));
fh(i)=1-cos(k*l);
end
figure
polar(theta,fe/max(fe));
figure
polar(theta,fh);
用MATLAB语言编程计算半波对称振子场图，电流分布和输入阻抗n=2*pi;
d=n/50+0.001;
[x,y] = meshgrid(-n:d:n,-n:d:n);
nn=30
j=0;
M = moviein(nn);
for j=1:nn
z=(x./sqrt(x.^2+y.^2)).^2.*(...
-sin(sqrt(x.^2+y.^2)-j*2*pi/nn)...
-cos(sqrt(x.^2+y.^2)-j*2*pi/nn)./sqrt(x.^2+y.^2)...
);
v=[-1 -0.9 -0.8 -0.75 -0.5 -0.25 -0.1 -0.05 0 0.05 0.1 0.25 0.5 ...
0.75 0.8 0.9 1];
%v=[ -0.75 -0.5 -0.25 0 0.25 0.5 0.75 ];
contour(z,v);
%contour(z,30);
%colorbar
axis equal
%grid on
M(:,j) = getframe;
end