同轴线_matlab_电磁场数值分析_源程序.m

clcclear all close all
ua = 10;%赋初值
i = 131;
j = 51;
A_chafen（ij） = 0;
Temp（ij） = 0;
derta = 0.01;%精度
womiga = 1.0;%超松弛法系数
epise = 8.85E-12;%介电常数
R1 = 10;%赋初值
R2 = 35;
R3 = 50;
C1 = 2*3.14*70*epise/log（R2/R1）;%计算C1
C2 = 2*3.14*60*epise/log（R3/R1）;%计算C2
for j = 12:35                   	%B边界
    A_chafen（1j） = 10 - （ua/log（R2/R1）*log（（j-1）/R1））;
end
for j = 12:50                       %C边界
    A_chafen（131j） = 10 - （ua/log（R3/R1）*log（（j-1）/R1））;
end
A_chafen（1:1311:11） = 10;  %内边界
Temp = A_chafen;
 for n=1:50000
    for i = 2:71
        for j = 12:35
            A_chafen（ij） = 0.25*（A_chafen（i+1j） + A_chafen（i-1j） + （1+1/（2*（j-1）））*A_chafen（ij+1） + （1-1/（2*（j-1）））*A_chafen（