hallen3.m

用MOM求解HALLEN S INTEGRAL EQUATION

clc;clear all;close all;
e=8.854e-12;%介电常数
u=pi*4e-7;%磁导率
c=3e+8;%光速
N=501;%点数
lambda=1;%波长
L=lambda/2*3;%天线长度
f=c/lambda;%频率
w=2*pi*f;%角频率
k=2*pi/lambda;%波数
eta=120*pi;%波阻抗
a=5*10^(-3);%天线半径
v0=1%馈电电压
delta=L/(N-1);%每一段长度
for mm=1:N
    Zb(mm)=((N+1)/2-mm)*delta;%竖直方向上将天线离散化
end
z=zeros(N,N);%生成阻抗矩阵
for mm=1:N
    for nn=1:N-1
        if mm==nn%对应格林函数奇点
           z(mm,nn)=log(delta/a)/(2*pi)-(j*k*delta)/(4*pi);
        else
            R1=sqrt((Zb(mm)-Zb(nn))^2+a^2);
            z(mm,nn)=exp(-j*k*R1)*delta/(4*pi*R1);%基函数和权函数均取狄拉克函数
        end;
    end;
end;
for mm=1:N
    z(mm,N)=cos(k*Zb(mm));%把Z矩阵最后一列补足
end;
for mm=1:N
    v(mm,1)=v0*sin(k*abs(Zb(mm)))/(j*2*eta);
end;
I=inv(z)*v;%最后一个元素不是电流项!
current=zeros(N,0);
    current=I;
current(N)=I(1);
mm=-(N-1)/2:1:(N-1)/2;
plot(mm,abs(current),'r'),xlabel('z/delta（z表示天线在ｚ轴上的坐标，delta为天线分割长度）'),ylabel('电流强度'),title('电流分布图');