matlab矩量法（Hallen公式）编的半波振子，但是电流总是不对

程序如下，希望大家能帮忙。只要成功，立刻给分。
%%%%矩量法解Hallen方程，求直线天线表面电流分布%%%%%
clear all
clc
ep=8.854e-12;%介电常数
c=3e8;%光速
ft=3e8;%频率
wt=2*pi*ft;%角频率
lamda=c/ft;%波长
a=0.005*lamda;%天线半径
L=0.5*lamda;%天线长度
k=2*pi/lamda;%波数
V=1;%激励电压
no=120*pi;%波阻抗
ct=V/(j*2*no);%系数

N=5;%分段数
pmn=zeros(N,N);
tm=zeros(N,1);
In=zeros(N,1);%电流
dL=L/N;%小段长

Q=linspace(-L/2,L/2,N+1);
for i=1:N
W(i)=(Q(i)+Q(i+1))/2; %每小段中点坐标
end

for i=1:N
pmn(i,N-1)=cos(k*W(i));%公式中的qm 列
pmn(i,N)=sin(k*W(i));%公式中的sm 列
tm(i,1)=ct*sin(k*abs(W(i))); %公式中的tm列
end

for i=1:N
for j=1:N-2
if i-j==0
pmn(i,j)=1/(2*pi)*log(dL/a)-j*k*dL/(4*pi);
else
Rmn=sqrt((W(i)-W(j))^2+a^2);
pmn(i,j)=exp(-j*k*Rmn)*dL/(4*pi*Rmn);
end
end
end

In=inv(pmn)*tm;
abs(In)

matlab？？ matlab下载 matlab 动画 MATLAB是什么？ 关于matlab matlab 编程 Matlab 使用 matlab电子书 关于Matlab matlab问题