同轴电缆的特性阻抗计算.doc

同轴电缆的特性阻抗计算
同轴电缆特性阻抗拉普拉斯方程矩形网格
同轴电缆的横截面可以看做是两个同心圆。外圆半径为2,内圆半径为1。外圆上的电势为1,内圆上的电势为0。我们依据这些条件,通过编写matlab程序来计算出同轴缆线的特性阻抗。
首先介绍一下计算中所用到的物理学公式。特性阻抗的公式为如下所示,C为电容,C0为光速。

由这两个公式,我们可将求解阻抗的问题转化为求解电量的问题。此时我们可以使用高斯公式。

为了处理截面上的问题,我们将面积分化为线积分。

本次计算过程中编程采用的方法是逐次超松弛迭代法。先将同轴电缆的截面按矩形网格进行划分。由于同轴电缆截面具有对称性,为了缩短程序运行时间,我们可以先计算四分之一截面内的电位分布。电位的迭代公式如下。
由于这个程序采用矩形网格来处理圆的问题,所以处理精度和处理速度都没有采用极坐标处理理想。如果希望得到跟极坐标情况下同样误差的结果,则需要耗费更多的计算时间。
图一为基本算法。图二、图三、图四分别是将代误差率为百万分之一时的特性阻抗、电势分布图和电场分布图。在文章的最后附有程序的代码。
图一图二
图三
图四
程序代码:
clc
clear all;
tic
r1=2;
r2=1;
n=.01;
c=299792458;%
err=-12;
wuchalv=.0001;
x=-r1:n:r1;
y=r1:-n:-r1;
l=length(x);
dones=ones((l+1)/2);
dlens=n*dones;
dianwei_1=NaN((l+1)/2);
[X,Y]=meshgrid(x,y);
for i=1:(l+1)/2
for j=1:(l+1)/2
if X(i,j)^2+Y(i,j)^2<=4&&X(i,j)^2+Y(i,j)^2>=1
dianwei_1(i,j)=0;
else
end
end
end
dianwei_2=isnan(dianwei_1);
len3=dlens;
for i=1:(l+1)/2
for j=1:(l+1)/2-1
if dianwei_2(i,j)==1&&dianwei_2(i,j+1)==0
len3(i,j+1)=abs(abs(sqrt(r1^2-Y(i,j+1)^2))-abs(X(i,j+1)));
else
end
end
end
len3((l+1)/2,1)=0;
len2=len3';
len1=dlens;
for i=1:(l+1)/2
for j=1:(l+1)/2-1
if dianwei_2(i,j)==0&&dianwei_2(i,j+1)==1
len1(i,j)=abs(abs(sqrt(r2^2-Y(i,j)^2))-abs(X(i,j)));
else
end
end
end

len4=len1';
c1=len1./n;
c2=len2./n;
c3=len3./n;
c4=len4./n;
dianwei_3=[dianwei_1 dianwei_1(:,(l+1)/2);dianwei_1((l+1)