用矩量法求解双位方程的matlab源程序.txt

MOM_for_Antenna matlab

[file head]:
clc;clear all;close all; 
e=8.854e-12;>介电常数 
u=pi*4e-7;>磁导率 
c=3e+8;>光速 
N=101;>点数 
lambda=1;>波长 
L=lambda/2;>天线长度 
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为天
... ...
[file tail]:
... ...
线分割长度）'),ylabel('电流强度'),title('电流分布图');