用hallen积分方程求偶极子天线电流分布的matlab程序

%***************start*******************%
clear;
clc;
lamda=1;%波长归一化为1
l=0.5*lamda;%半波长偶极子天线的长度，通过更改l的值，可
            %以随意求出其他长度的偶极子天线的电流分布图
k=2*pi/lamda;%k为波数
a=0.001*lamda;%偶极子天线半径
omega=2*log（l/a）;
V0=1;%V0为常数，这里取为1
step1=0.01;
step2=0.001;
x=1;
Iz=zeros（l/step1+11）;%预先为Iz矩阵分配空间

%****************求Iz*********************************************%
for z=-0.5*l:step1:0.5*l

%***************求积分****************************%
   M=0;I=0;%赋初值   
   for z1=-0.5*l:step2:0.5*l   
      r=sqrt（a^2+（z-z1）^2）;
      r1=sqrt（a^2+（0.5*l-z1）^2）;
      F0z=cos（k*z）-cos（k*l*0.5）;
      F0z1=cos（k*z1）-cos（k*l*0.5）;
      I=I+（（F0z1*exp（-1i*k*r）-F0z）/r）*step2; %化积分为累加求和
      M=M+（（F0z1*exp（-1i*k*r1））/r1）*step2;   %化积分为累加求和   
   end
 %****************求Iz**********************