电磁场与电磁波实验报告

  • 格式:doc
  • 大小:2.99 MB
  • 文档页数:12

. .. . . . . . w 电磁场与电磁波 实 验 报 告

哈尔滨工业大学(威海) 实验名称: 有限差分法解电场边值问题 实验日期: 2012年12月8日 姓 名: 赵文强 学 号: 100240333 . .. . .

. . . w 问题陈述 如下图无限长的矩形金属导体槽上有一盖板,盖板与金属槽绝缘,盖板电位为U0,金属槽接地,横截面如图所示,试计算此导体槽内的电位分布。

参数说明:a=b=10m, 0U=100v

实验要求 1) 使用分离变量法求解解析解; 2) 使用简单迭代发求解,设-10=100.1,1xy,两种情况分别求解数值解; 3) 使用超松弛迭代法求解,设-10=100.1xy,确定(松弛因子)。 求解过程 一、 分离变量法求解 因为矩形导体槽在z方向为无限长,所以槽内电位函数满足直角坐标系中的二维拉普拉斯方程。 2222

00(0,)0,(,)0(0)(,0)0,(,)(0)xyyayybxxbUxa

. .. . . . . . w 根据边界条件可以确定解的形式:

1ππ(,)sin()sinh()nnnxnyxyAaa

利用边界条件0(,)xbU求解系数。 01ππsin()sinh()nnnxnbAUaa



01πsin()nnnxUfa



00041,3,5,2πsin()dπ02,4,6,an

UnnxfUxnaan





011πππsin()sinh()sin()nnnnnxnbnxAUfaaa





041,3,5,πsinh(π/)'πsinh()02,4,6,nnUnfnnbaAnbna





01,3,5,4ππ(,)sin()sinh()πsinh(π/)nUnxnyxynnbaaa

简单迭代法求解

二、 有限差分法 有限差分法(Finite Differential Method)是基于差分原理的一种数值计算法。其基本思想:将场域离散为许多小网格,应用差分原理,将求解连续函数的泊松方程的问题转换为求解网格节点上

的差分方程组的问题。 泊松方程的五点差分格式 )(414243210204321Fh Fh 当场域中,0得到拉普拉斯方程的五点差分格式 . .. . . . . . w )(41044321004321 差分方程组的求解方法(1) 高斯——赛德尔迭代法 ][)(,)(,)(,)(,)(,2k1jikj1i1k1ji1kj1i1kjiFh41

(1-14)

式中:,2,1,0,2,1,kji, • 迭代顺序可按先行后列,或先列后行进行。 • 迭代过程遇到边界节点时,代入边界值或边界差分 格式,直到所有节点电位满足)(,)(,kjilkji为止。 (2)超松弛迭代法 ][)(,)(,)(,)(,)(,)(,)(,kji2k1jikj1i1k1ji1kj1ikji1kji4Fh4 (1-15)

式中:——加速收敛因子)21( 可见:迭代收敛的速度与有明显关系 (一) 简单迭代法 简单迭代法程序: 1) 步长=1 clear all;clc;close all; %设置节点数,步长1 hx=11; hy=11; v1=ones(hy,hx); %% %% %设置边界条件 v1(hy,:)=ones(1,hx)*100; v1(1,:)=zeros(1,hx); v1(1:hy,1)=0; v1(1:hy,hx)=0; %% %% %初始化 v2=v1;

图1-4 高斯——赛德尔迭代法 . .. . . . . . w maxt=1; t=0; k=0; %% %% while(maxt>1e-10) k=k+1; %计算迭代次数 maxt=0; for i=2:hy-1 for j=2:hx-1 v2(i,j)=(v1(i,j+1)+v1(i+1,j)+v2(i-1,j)+v2(i,j-1))/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]); title('步长=1,各点电位'); subplot(1,2,2),contour(v2); %画等势线 title('等位线');

实验结果:

051005

10

020406080100步长=1,各点电位等位线

2468101234567891011

图1,简单迭代法结果,步长1 . .. . .

. . . w 步长1,迭代次数 k = 246 各节点电位数据:

2) 步长=0.1 实验结果:

050100050100020406080100步长=0.1,各点电位等位线2040608010010

20

30405060708090100

图2,简单迭代法步长0.1

0 0 0 0 0 0 0 0 0 0 0 0 1.107499 2.099344 2.877502 3.371569 3.540667 3.371569 2.877502 2.099344 1.107499 0 0 2.330652 4.412375 6.039095 7.068108 7.419529 7.068108 6.039095 4.412375 2.330652 0 0 3.802735 7.180408 9.798395 11.44224 12.00123 11.44224 9.798395 7.180408 3.802735 0 0 5.699881 10.70813 14.53184 16.90122 17.70092 16.90122 14.53184 10.70813 5.699881 0 0 8.28866 15.42038 20.7196 23.9299 25 23.9299 20.7196 15.42038 8.28866 0 0 12.03438 21.96514 28.99628 33.09878 34.43928 33.09878 28.99628 21.96514 12.03438 0 0 17.88372 31.40952 40.20161 45.02964 46.55957 45.02964 40.20161 31.40952 17.88372 0 0 28.09096 45.58763 55.37098 60.25862 61.73971 60.25862 55.37098 45.58763 28.09096 0 0 48.8925 67.47904 75.43605 78.89417 79.88201 78.89417 75.43605 67.47904 48.8925 0 0 100 100 100 100 100 100 100 100 100 0 . .. . .

. . . w 步长0.1,迭代次数 k = 20051 部分实验结果数据截图:

图3,简单迭代法步长0.1部分数据 (二) 超松驰迭代法 1. 理论最佳松弛因子实验结果 实验程序: clear all;clc;close all; %设置节点数,步长0.1 hx=101; hy=101; m=100; n=100; v1=ones(hy,hx); %% %% %设置边界条件 v1(hy,:)=ones(1,hx)*100; v1(1,:)=zeros(1,hx); v1(1:hy,1)=0; v1(1:hy,hx)=0; %% %%