矩量法求解线天线传输模型积分方程（含代码）

一.实验原理

对于线天线的矩量法分析，我们得知Hallen积分方程后，将边界条件带入后即可进行求解

1.对整个线天线的长度进行分割

将z=-l~z=l区间分成N个节点，并将每个节点进行编号；

2.矩形基函数展开

将转化为，其中z’为源点位置（如图一所示）；

图一

3.匹配点选择与匹配方程

由Hallen积分方程可得=，将第一步所得节点分别带入，代入可得M×M的矩阵，其中具体方程如下；

4.矩阵方程与矩阵中系数求解

将第三步所得方程得出矩阵

其中矩阵方程与矩阵元素表示为：

通过格林公式与第一步分隔区域后的节点可以解得输入阻抗矩阵元素，但是为解决对格林公式积分方程中源点与场点重合的问题以及编程时处理积分问题，需对积分进行近似处理，

在m与n相离较远时，将近似处理为：

在m=n时，会出现奇异点，对分子在零点附近做泰勒展开后进行积分，可得：

5.求解矩阵方程

得到形如的方程，将边界条件带入，即可求解a。求得表面电流分布后可得到输入阻抗矩阵，其中远区辐射远程、方向系数等亦可求得并直观表示。

二、实验步骤及代码讲解

1.首先设置线天线参数

c=3.0e8;
f=3.0e8;
k=2*pi*f/c;
a=0.001;
N=21; %分段个数
len=0.625; %线天线长度
%len=0.5;
%len=0.625;
delta=2*len/(N-1);
Z=120*pi;

2.分割线天线并初始化匹配点

zn=zeros(N,2);
zn(1,1)=-len;
zn(1,2)=zn(1,1)+delta/2;
for i=2:N-1
    zn(i,1)=zn(i-1,2);
    zn(i,2)=zn(i,1)+delta;
end
zn(N,1)=zn(N-1,2);
zn(N,2)=zn(N,1)+delta/2;
zm=-len:delta:len;

在该代码片段中，匹配点从一一被赋值，而积分区域与分别为相近两匹配点之间的中点，由此得到了匹配点具体坐标；

3.计算各个矩阵元素

将输入阻抗矩阵元素通过穷举法按照书中给出的公式直接计算即可，代码如下：

Zn=zeros(N,N);
znn_1=(1/(4*pi))*log((sqrt(delta^2+4*a^2)+delta)/(sqrt(delta^2+4*a^2)-delta));
znn_2=-(1j*k/(4*pi))*delta;
znn_3=-(k^2/2)*(1/(8*pi))*(delta/2*sqrt(4*a^2+delta^2)+a^2*log((sqrt(delta^2+4*a^2)+delta)/(sqrt(delta^2+4*a^2)-delta)));
for i=1:N
    for j=1:N
        if(j==1)
            Zn(i,j)=-1*cos(k*zm(i));
        else
            if(i==j)
                Zn(i,j)=znn_1+znn_2+znn_3;
            else 
                Zn(i,j)=exp(-1j*k*sqrt(a^2+(zm(j)-zm(i))^2))*(zn(j,2)-zn(j,1))/(4*pi*sqrt(a^2+(zm(j)-zm(i))^2));
            end
        end
    end
end

之后对B矩阵中元素进行计算，其中表面电流，其中；

Bm=-1j*(sin(k*abs(zm)))/(2*Z);
IA=Zn\Bm';
IA(1)=0;

4. 计算电场幅值

r=10;
%r=100;
theta=0:0.01:pi;
for i=1:length(theta)
    Er(i)=1j*k*Z*exp(-1j*k*r)*sin(theta(i))*sum(IA'.*exp(1j*k*zm*cos(theta(i))),'all')*delta/(4*pi*r);
end

5.输出线天线电流分布、电场幅值以及三维方向图

figure(1)
plot(zm,abs(IA));
xlabel('位置z/m'),ylabel('电流I幅值/A');
title('对称线天线的电流分布');
Z_in=1/IA((N+1)/2);   %输入阻抗Z_i
disp(Z_in)

figure(2)
polarplot(theta,abs(Er))
title('对称线天线的电场幅值');

theta1=0:0.01:2*pi;
F_theta=(cos(k*len*cos(theta1))-cos(k*len))./sin(theta1);
figure(3)
polarplot(theta1,abs(F_theta))
title('对称线天线电场方向系数')

theta2=meshgrid(0:pi/180:pi);
phi=meshgrid(0:2*pi/180:2*pi)';
F_theta_3=(cos(k*len*cos(theta2))-cos(k*len))./sin(theta2);
[x,y,z]=sph2cart(phi,pi/2-theta2,abs(F_theta_3));
figure(4);
mesh(x,y,z);
xlabel('x轴');
ylabel('y轴');
zlabel('F轴');
title('对称线天线的三维方向图')

三、实验结果

①当len=0.25时

N=11时，电流分布为：

电场幅值

三维方向图

线天线方向系数图

其中输入阻抗为

       

当N=21时

电流分布为：

电场幅值为：

      

电场方向系数为：

三维方向图为：

输入阻抗为：

②当len=0.5

N=11时，电流分布为：

电场幅值为：

电场方向系数图为：

三维方向系数图为：

输入阻抗为：

当N=21时

电流分布为：

电场幅值图为：

电流方向系数图为：

电场三维方向图为：

输入阻抗大小为：

③当len=0.625时

N=11时电流分布为：

电场幅值分布为：

电场方向系数为：

三维方向图为：

当N=21时

电流分布为：

电场幅值图为：

电场方向系数图为：

三维方向图为：