文档介绍:电磁场与电磁波
实验报告
实验名称: 有限差分法解电场边值问题
实验日期: 2012 年 12 月 8 日
名:
学 号:
哈尔滨工业大学(威海)
问题陈述
如下图无限长的矩形金属导体槽上有一盖板,盖板与金属槽绝缘, 盖hy-1
for j=2:hx-1
v2(i,j) = (vl(i,j+l)+vl(i+l,j)+v2(i-l,j)+v2(i,j-l))/4;%拉普拉斯方程差分形式 t=abs(v2(i,j)-v1(i,j));
if(t>maxt) maxt=t;end
end
end
v1=v2;
end
%%
%%
%可视化显示
subplot(1,2,1),mesh(v2); %画电势的三维曲面图
axis([0 ,11,0,11,0,100]);
tit le ('步长=1,各点电位');
subplot(1,2,2),contour(v2); %画等势线
title('等位线');
实验结果:
步 长 =1 , 各 点 电 位
5
5
100
等位线
11
10
9
8
7
6
5
4
3
2
图1,简单迭代法结果,步长1
步长1,迭代次数
k =
246
各节点电位数据:
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
25
0
0
0
0
0
0
0
0
0
0
100
100
100
100
100
100
100
100
100
0
2) 步长=
实验结果
步长=,各点 电 位
100
0
100
100
00
等位线
100
90
80
70
60
50
40
30
20
10
20 40 60 80 100
图2,
,迭代次数
k 二
20051
部分实验结果数据截图:
0
2
D
OLGtOS
DjOSID
ajtJjt37
C.^1545
0J376D
£
010071
-n?'
=
u
烬紬
0JJ87^
ClU第
01520
Gd型
0252
4
c
axjisjT
?
aiTCL
OJ*2EZ
C-.2K2 Q29LS
M2
<4
S
&
0U4J9
M&77
i/in-j
01751