clc;
clear all
clf;
tic; %计时
lambda=1;
N=31;a=0.0000001;%已知天线和半径
ii=1;
for h=0.2:0.1:0.9
L=h*lambda;
len=L/N;%将线分成奇数段,注意首末两端的电流为0
e0=8.854e-012;u0=4*pi*10^(-7);k=2*pi/lambda;
c=3e+008;w=2*pi*c;%光速,角频率
ata=sqrt(u0/e0);
z(1)=-L/2+len/2;
for n=2:N
z(n)=z(n-1)+len;
end
for m=1:N
for n=1:N
if (m==n)
p(m,n)=log(len/a)/(2*pi)-j*k*len/4/pi;
else
r(m,n)=sqrt((z(m)-z(n))^2+a^2);
p(m,n)=len*exp(-j*k*r(m,n))/(4*pi*r(m,n));
end
end
end
for m=1:N
q(m)=cos(k*z(m));
s(m)=sin(k*z(m));
t(m)=sin(k*abs(z(m)))/(j*2*ata);
end
pp=p(N+1:N^2-N);
pp=reshape(pp,N,N-2);