电磁场与电磁波实验有限差分法
电磁场与电磁波实验报告
实验项目：有限差分法
一、实验目的及要求
1、学习有限差分法的原理与计算步骤；
2、学习用有限差分法解静电场中简单的二维静电场边值问题；
3、学习用Matlab语言描述电磁场与电磁波中内容，用matlab求解问题并用图形表示出了，学习matlab语言在电磁波与电磁场中的编程思路。
二、实验内容
理论学习：学习静电场中边值问题的数值法中的优先差分法的求解知识；
实践学习：学习用matlab语言编写有限差分法计算二维静电场边值问题；
三、实验仪器或软件
Matlab7.0
电脑
四、实验原理
有限差分法的基本思想
将计算场域划分成网格，把求解场域内连续的场分布用求解网格节点上的离散数值解来代替；即用网格节点的差分方程近似代替场域内的偏微分方程来求解。
简单迭代法
先对场域内的节点赋予初始值
）（0,j i Φ ，这里上标（0）表示第0次近似值，即初
始值。
然后再按照： ][4
1k 1,k ,1k 1,k ,11k ,）（）（）（）（）（++--+Φ+Φ+Φ+Φ=Φ
j i j i j i j i j i
进行反复迭代。
若当第N 次迭代结束后，所有内节点相邻两次迭代值之间的绝对误差小于事先给定的精度，则迭代停止。
W MAX N j
i N j i 〈Φ-Φ-）（）（1,, 注意：
初始值的赋予是任意的；
赋予初始值后，请按“从左到右、从下到上”的固定顺序依次计算各节点值； 当所有节点都算完一遍后，再用它们的新值代替旧值，即完成一次迭代。
五、实验步骤
复习理论知识；
编写matlab 程序；
六、结果分析与问题讨论
1、程序：
clear
X=[0,0,0,0,0;0,25,25,25,0;0,50,50,50,0;0,75,75,75,0;100,100,100,100,100]
Pot=[0,0];
for i=2:4
for j=2:4
PotX(i-1,j-1)=(X(i-1,j)+X(i,j-1)+X(i+1,j)+X(i,j+1))/4 Pot(1)=abs(PotX(i-1,j-1)-X(i,j));
Pot(2)=max(Pot)
end
end
X(2:4,2:4)=PotX
num=1;
while(max(1000.*Pot)>1)
Pot(2)=0; (,1,2,......) (0,1,2,......)
i j k ==
for i=2:4
for j=2:4
PotX(i-1,j-1)=(X(i-1,j)+X(i,j-1)+X(i+1,j)+X(i,j+1))/4 Pot(1)=abs(PotX(i-1,j-1)-X(i,j));
Pot(2)=max(Pot)
end
end
X(2:4,2:4)=PotX
num=num+1
end
surf([0:4],[0:4],X);
shading interp
colorbar('horiz')
title('有限差分法计算电位图');
2、运行结果
X =
0 0 0 0 0
0 25 25 25 0
0 50 50 50 0
0 75 75 75 0
100 100 100 100 100
%%第一次迭代
PotX =
18.7500
Pot =
6.2500 6.2500
PotX =
7.1440 9.8230 7.1440
18.7515 25.0023 18.7515
42.8583 52.6801 42.8583
Pot =
1.0e-003 *
0.3815 0.7629
%%第28次迭代
X =
0 0 0 0 0
0 7.1440 9.8230 7.1440 0
0 18.7515 25.0023 18.7515 0
0 42.8583 52.6801 42.8583 0
100.0000 100.0000 100.0000 100.0000 100.0000
num =
28
3、波形图
matlab软件在使用有限差分法研究静电场边值问题中有着重要的作用，它能够快捷有效并且准确的解决边值问题，是解决计算相对复杂问题的有效工具。
由于很久没有用过matlab
程序，忘记了matlab的函数使用和初值定义方法等，导致实验的进度十分缓慢。
对于编程思路，理清了有限差分法的概念就使之变得容易很多，所以我复习了优先差分法的求解方法，并自己做了一遍题，思路就有了，后期编程就是要知道运用什么函数就可以了。
有限差分法在电磁场理论教学中的应用研究
分布 ，而 电磁场与 电磁波又都具 有不可 见和 不可触摸
的特性 ，只能进行抽象 的想象或通过仪器进行 数据测
量 ；三是 电磁波是动态 的 ，电磁波是 电磁场相 互激 发 的结果 ，它在空间 的传播每时每刻 它的位置和状 态都 在发生变化。
传 统 的 电磁 场 与 电磁 波 课 堂 理 论 教 学 方式 ，一是
１静 态 场 模 拟
在 均 匀 介质 内 ，静 电势 （满 足泊 松 方程 ｐ
Ｖ ＝
如 各种复 杂的边 界条件 等 ，这种枯燥 的教学方法无法 引起学生 的学 习兴趣 ，找到一种方便 快捷的方法计算
和 模 拟 电磁 场 ， 并 以形 象 化 的 图形 演 示 电磁 场 是 一 种
行 之有效 的教学手段…。二是偏重理论教学 ，而忽视
为求解 由偏微 分方程定解 问题所构 造的数学模型 ，有
限差分法是将 定解区域 （ 区 ）离散化为 网格离散 节 场 点的集合 。并 以各离散点上 函数的差商来近似该 点的 偏导数 ，使待 求的偏 微分 方程定解 问题转化 为一组 相 应 的差 分方程 。根据 差分 方程组解 出各离散点 处的待
解不仅耗 时费力 ，容易出差错 ，并且求解 的电磁场和
程 中，学生就会思考 ：为什 么可以这么求解 ；如何 求
解 ；解决一 些什 么问题？ 带着 这些 问题学 习，学 习效
果 明 显 提 高。
电磁 波问题均 为设计 的理想化模型 ，只对一些特殊对
称 的边界 才能求解 ，而实际工程中的问题是变化 的，
电磁场 与 电磁 波理论作 为电子信 息类专业 的一门
的电位 分布问题 ，利用时域有 限差分法数值计 算波导
中的 电磁 波传 播问题 ，实现 电磁场和 电磁波可视化教 学简 化 了繁 琐的数学推 导 ，能够形象而直观地输 出可
电磁场与电磁波教案4(hao)
边值问题的求解是偏微分方程的求解，同 时要考虑解得存在性、唯一性和稳定性。 已知整个边界上的电位函数（第一类边界 条件），则场域的解答是唯一的。 已知整个边界上的电位法向导数（第二类 边界条件），则场域的解答是唯一的。 已知一部分边界上的电位函数和另一部分 边界上的电位法向导数（第三类边界条 件），则场域的解答是唯一的。
使用镜像法时必须注意：
边界必须是无限大，或者是具有球或圆柱 对称性。 为了保证待求场区域内原方程成立，镜像 电荷不能出现在待求场区域内。 镜像电荷的确定必须保证原有的边界条件 全部满足。
1 q q '' q q '' ( z 0) 2 ( x, y, z ) ( ) 4 2 R 4 2 x2 y 2 ( z h)2
1
q
z q R1 O h R1′ P z q q h
1
h
2
O R2 P
q′
在z=0面上应用电位边界条件
1 2 q ' q 1 2 q '' q ' 1 2 q 1 2
z
( x, y, z)
q 4 0
[
1 x y ( zபைடு நூலகம் 1)
2 2 2
1 x y ( z 1)
2 2 2
P
]
O x
106 1 1 (0,0, z) [ ] 104 4 0 z 1 z 1
z1 1.67 z2 0.45
4.3.2 导体球面的镜像
1 z 0 2 z 0 1 2 2 1 z z z 0
q q ' q q '' 2 1 q q ' q q '' z 0
电磁场数值计算
电磁场数值计算引言：电磁场是电荷和电流产生的物理现象，它在现代科技和工程中起着至关重要的作用。
对电磁场的数值计算是研究和应用电磁学的基础。
本文将介绍电磁场数值计算的原理和方法，并探讨其在实际问题中的应用。
一、电磁场的数值计算方法：电磁场的数值计算可以通过求解麦克斯韦方程组来实现，这是描述电磁场的基本方程。
麦克斯韦方程组包括四个方程，分别是电场的高斯定律、磁场的高斯定律、法拉第电磁感应定律和安培环路定律。
通过数值方法求解这些方程，可以得到电磁场在空间中的分布情况。
1. 有限差分法：有限差分法是一种常用的数值计算方法，通过将空间离散化为有限个点，时间离散化为有限个步骤，将偏微分方程转化为差分方程进行求解。
在电磁场计算中，可以将空间划分为网格，通过有限差分法计算电场和磁场在网格节点上的数值。
2. 有限元法：有限元法是一种广泛应用于工程领域的数值计算方法，它通过将计算域划分为许多小的有限元，将偏微分方程转化为代数方程组进行求解。
在电磁场计算中，可以将计算域划分为三角形或四边形网格，通过有限元法计算电场和磁场在每个有限元上的数值。
3. 边界元法：边界元法是一种适用于边界值问题的数值计算方法，它将偏微分方程转化为积分方程进行求解。
在电磁场计算中，可以通过边界元法计算电场和磁场在边界上的数值，然后利用边界条件求解整个计算域内的电磁场分布。
二、电磁场数值计算的应用：电磁场数值计算在科学研究和工程应用中具有广泛的应用价值，以下是一些常见的应用领域：1. 电磁场仿真：电磁场数值计算可以用于电磁场仿真，模拟和预测电磁场在不同结构和材料中的分布情况。
例如，可以通过数值计算预测电磁波在天线中的传播情况，从而优化天线设计和布局。
2. 电磁场辐射：电磁场数值计算可以用于估计电磁场辐射对人体和环境的影响。
例如，可以通过数值计算评估电磁辐射对人体健康的潜在风险，从而制定相应的防护措施。
3. 电磁场感应：电磁场数值计算可以用于分析电磁感应现象，研究电磁场对电路和设备的影响。
电磁学的数值计算方法
电磁学的数值计算方法电磁学是研究电场和磁场相互作用的学科，它在日常生活和科学研究中起着重要的作用。
随着计算机技术的快速发展，数值计算方法在电磁学中的应用也越来越广泛。
本文将介绍几种常用的电磁学数值计算方法，并探讨其原理和应用。
一、有限差分法（Finite Difference Method）有限差分法是一种基于离散化空间和时间的数值计算方法，常用于求解求解具有边值条件的偏微分方程。
在电磁学中，有限差分法可以用来求解电磁场的静电场、静磁场以及时变电磁场等问题。
该方法通过将空间和时间进行网格离散化，将偏微分方程转化为差分方程，并用迭代方法求解得到数值解。
二、有限元法（Finite Element Method）有限元法是一种广泛应用于各种物理问题求解的数值计算方法，电磁学也不例外。
该方法通过将求解区域划分为有限的小元素，并在局部内部逼近真实场量的变化。
在电磁学中，有限元法可以用来求解电场、磁场以及电磁波传播等问题。
通过选择合适的元素类型和插值函数，以及建立元素之间的边界条件，可以得到电磁场的数值解。
三、时域积分法（Time Domain Integral Method）时域积分法是一种基于格林函数的数值计算方法，通过积分形式表示电磁场的边界条件和过渡条件，进而求解电磁场。
时域积分法广泛应用于求解电磁波的辐射和散射问题，如天线辐射和散射、电磁波在介质中的传播等。
该方法通过离散化电磁场的源和观测点，并利用格林函数的性质进行数值积分，得到电磁场的数值解。
四、有限时域差分法（Finite-Difference Time-Domain Method）有限时域差分法是一种基于电磁场的离散化网格和时间的有限差分法，是求解各种电磁问题最常用的数值计算方法之一。
有限时域差分法通过离散化时空域，将麦克斯韦方程组转化为差分方程组，并通过时间步进的方式求解得到电磁场的数值解。
该方法适用于求解各种电磁波传播、辐射和散射等问题。
电磁场数值分析,有限差分法
可得：
(K ) x (x x 0 )K 0 x x 0 )n ) (( K 0 K ! 1 2 2 1 3 3 1 0 h( ) 0 h ( 2 ) 0 h ( 3 ) 0 x 2! x 3! x 1 2 2 1 3 3 3 0 h( ) 0 h ( 2 ) 0 h ( 3 ) 0 x 2! x 3! x
有限差分法的基本步骤 （1）剖分场区，确定离散点。将所研究的电位分布 按某种几何形状(如矩形、任意多边形等)剖分成网络系统。 （2）建立电位分布问题的差分方程组。
（3）求解差分方程组。可以采用各种迭代法，如简 单迭代法，塞德尔迭代法，超松弛迭代法等等。
100 V 1 4 OV 7 2 5 8 3 6 9
在xoy 平面内把所求解区域 划分为若干个相同的正方形格 子，边长均为h，假设0点点位 为φ0 ，其余各点为φ1，φ2，φ3，
φ4，φ5。
将这几个点的点位用泰勒级数展开：
f 1 2f 1 3f f f0 (x x 0 ) 2 (x x 0 )2 3 (x x 0 )3 ... x 2! x 0 3! x 0 0
四，计算实例
1V
如图，一正方形区域，四个边的电位分 别是0V,0V,1V,100V，求解该区域内部 的电位分布。
0V
解： 将该正方形区域分割成4X4的小正 方形区域，则一共有9个内点。按照前 面得出的结论，任意一点的电位等于他 周围四个点电位的平均值。可以得到方 程组：
1 2 3 4 5 6 7 8 9 二源自差分方程的数值解法简单迭代法
《电磁场与电磁波》仿真实验
年《电磁场与电磁波》仿真实验————————————————————————————————作者：————————————————————————————————日期：《电磁场与电磁波》仿真实验2016年11月《电磁场与电磁波》仿真实验介绍《电磁场与电磁波》课程属于电子信息工程专业基础课之一，仿真实验主要目的在于使学生更加深刻的理解电磁场理论的基本数学分析过程，通过仿真环节将课程中所学习到的理论加以应用。
受目前实验室设备条件的限制，目前主要利用MATLAB 仿真软件进行，通过仿真将理论分析与实际编程仿真相结合，以理论指导实践，提高学生的分析问题、解决问题等能力以及通过有目的的选择完成实验或示教项目，使学生进一步巩固理论基本知识，建立电磁场与电磁波理论完整的概念。
本课程仿真实验包含五个内容：一、电磁场仿真软件——Matlab的使用入门二、单电荷的场分布三、点电荷电场线的图像四、线电荷产生的电位五、有限差分法处理电磁场问题目录一、电磁场仿真软件——Matlab的使用入门 (4)二、单电荷的场分布 (10)三、点电荷电场线的图像 (12)四、线电荷产生的电位 (14)五、有限差分法处理电磁场问题 (17)实验一电磁场仿真软件——Matlab的使用入门一、实验目的1. 掌握Matlab仿真的基本流程与步骤；2. 掌握Matlab中帮助命令的使用。
二、实验原理（一）MATLAB运算1.算术运算(1)．基本算术运算MATLAB的基本算术运算有：＋(加)、－(减)、*(乘)、/(右除)、\(左除)、^(乘方)。
注意，运算是在矩阵意义下进行的，单个数据的算术运算只是一种特例。
(2)．点运算在MATLAB中，有一种特殊的运算，因为其运算符是在有关算术运算符前面加点，所以叫点运算。
点运算符有.*、./、.\和.^。
两矩阵进行点运算是指它们的对应元素进行相关运算，要求两矩阵的维参数相同。
例1：用简短命令计算并绘制在0≤x≦6范围内的sin(2x)、sinx2、sin2x。
项目训练二 求二维金属槽内的点位分布
将这几点的电位用泰勒级数展开，化简，近似可得：
1
0
1 ( ) 2 3 4 4 1
上式表明，任一点的电位等于它周围四个点电位的平均值。
图 3-1 差分网格
差分方程的数值解法: 平面内有多少个节点，就能得到多少个差分方程，当这些节点数目较大时，使用迭代 法求解差分方程组比较方便。 a. 简单迭代法： 用迭代法解二维电位分布时，将包含边界在内的节点均以下标（i，j）表示，i，j 分 别表示沿 x，y 方向的标点。次序是 x 方向从左到右，y 方向从上到下，我们用上标 n 表示 某点电位的第 n 次的迭代值。下式得出点（i，j）的第 n+1 次电位的计算公式：
n n 1 n 1 1 n ( ) i , j 1 i 1, j i , j 1 4 i 1, j
c. 超松弛迭代法：
n 1 i, j
此式也称为异步迭代法，异步迭代法比简单迭代法收敛速度加快一倍左右。
为了加快收敛速度，常采用超松弛迭代法。计算时，将某点的新老电位值之差乘以一
图 4-3 三维电位分布图 等位线分布图
4.6 中心点电位的数值解和精确解的比较
4.6.1 中心点 p ( , ) 的数值解：
因为电位是从上至下逐渐呈递减状，中间点的所在等位线的 a b 为 50V。
, 2 2
a b 2 2 a b 2 2
4.6.2 中心点 p ( , ) 的精确解：
在不同的 m 取值下， a b 的取值会不同。
, 2 2
m=1.2， a b =49.5133V
, 2 2 , 2 2
m=1.4， a b =49.5149V M=1.6， a b =49.5161V
电磁场与电磁波实验报告
工程电磁场实验报告实验1 :熟悉Matlab、矢量运算要求：学习矢量的定义方法（例A=[1,2,3]），加减运算，以及点积积cross（A,B）、求模运算norm（A）。
1）通过调用函数，完成下面计算【p31，习题1.1】。
给定三个矢量A、B和C如下：求（1）e；（2）|A-B| ；（3）AB ；（4）AB（5）A在B上的投影（6）A C ；(7)A (B C)和C (A B) ；( 8)(A B) C 和A (B C)答案：(1)e A =[0.2673,0.5345,—0.8018]； (2)| A —B|= 7.2801 ； (3)=；(4)% =2.3646(135.4815) ；( 5)-2.6679 ；(6)A C =[-4,-13,-10]；(7) A (B C) =C (A B) =-42(8) (A B) C =[2, -40,5] ；A (B C)珂55,-44,-11]程序代码：A=[1,2,-3]B=[0,-4,1]C=[5,0,-2](1) eA=A/norm(A) ( 2)X=norm(A-B)(4) D=acos(dot(A,B)/(norm(A)*norm(B)))(6)G=cross(A,C)(7)I=dot(A,cross(B,C))(8)M=cross(cross(A,B),C) 2、三角形的三个顶点位于A(6,-1,2), B(-2,3,-4), C(-3, 1,5)点，求(1)该三角形的面积；(2)与该三角形所在平面垂直的单位矢量。
(答案S=42.0119, n =[0.2856,0.9283,0.238])程序代码：A=[6，-1,2]B=[-2,3,-4]C=[-3, 1,5]D=acos(dot(A-B,C-A)/( norm(A-B)* norm(B-C)))A = e x 2e y _ 3e zC = 5e x- 2e zdot(A,B)、叉(3)Y=dot(A,B)(5) H=norm(A)*cos(F)K=dot(C,cross(A,B))P=cross(A,cross(B,C))E=norm(A-B)*si n(D)S=norm(B-C)*E/2n=cross(A-B,C-A)/( norm(C-A)* norm(A-B)*(sqrt(1-(dot(C-A,A-B)/( no rm(C-A- A)* norm(A-B))F2)))3、在直角坐标系中，在点P(3,4,2 )处的电场强度为4e x 2e y 3e z。
通信电子中的电磁场计算方法
通信电子中的电磁场计算方法电磁场在通信电子领域中扮演着至关重要的角色。
无论是无线电通信还是光通信，都需要计算电磁场的传输效应和衰减规律。
为了得到准确的计算结果，研究人员需要将电磁场模型转化为数学模型进行计算，这就涉及到了电磁场计算方法。
本文将介绍通信电子中常用的电磁场计算方法。
有限差分法有限差分法（Finite Difference Method，FDM）是求解偏微分方程的常用数值方法。
在电磁场计算中，有限差分法主要用于求解Maxwell方程组（即电磁场方程），通过数值模拟的方式来得到电磁场分布的规律。
该方法的关键是将空间离散化为网格，电磁场在各个节点上的取值通过差分来表示。
有限差分法的优点是易于理解和实现，可以处理复杂的几何形状和不规则介质。
不过，该方法对于某些边界条件的处理比较困难，且计算过程中需要大量存储和计算空间，其计算效率并不高。
有限元法有限元法（Finite Element Method，FEM）是一种求解偏微分方程的数值方法。
与有限差分法不同，有限元法将连续介质离散为若干个简单的几何单元，通过组建形函数来表示电磁场在每一个单元中的分布情况。
通过将所有单元的贡献组合起来，可以得到整个电磁场的分布情况。
有限元法的优点是能够完美地处理复杂的几何形状和不规则介质，并且可以方便地处理不同类型的边界条件。
但是需要高度复杂的数值计算，且对计算机性能的要求较高。
时域积分法时域积分法（Time Domain Integral Method，TDIM）是一种计算电磁波传输的数值方法。
与有限差分法和有限元法不同，时域积分法并不对电磁场进行离散化处理，而是通过求解Maxwell方程组的积分形式来得到电磁场随时间的演变规律。
时域积分法的优点是可以考虑到时域的变化，从而更加真实地模拟电磁场传播的过程。
但是由于每一时刻都需要计算全局的积分，其计算效率和精度都较低。
频域方法频域方法是一种将电磁波转化为频域信号进行计算的方法。
计算电磁场理论中的有限差分法与有限元法
计算电磁场理论中的有限差分法与有限元法电磁场理论是电磁学的重要组成部分，研究电磁场的分布和变化规律对于解决实际问题具有重要意义。
在计算电磁场中，有限差分法和有限元法是两种常用的数值计算方法。
本文将从理论原理、应用范围和优缺点等方面对这两种方法进行探讨。
有限差分法是一种将连续问题离散化的方法，通过将连续的电磁场分割成网格，然后在每个网格上进行离散计算。
这种方法的基本思想是将微分方程转化为差分方程，然后利用差分方程进行求解。
有限差分法的优点是简单易懂，计算过程直观，适用于各种电磁场问题的求解。
然而，由于差分法中的网格离散化会引入一定的误差，所以在计算精度上存在一定的限制。
与有限差分法相比，有限元法是一种更加精确的数值计算方法。
有限元法将电磁场问题的求解区域划分为有限个小单元，然后在每个小单元上建立适当的插值函数，通过求解代数方程组得到电磁场的近似解。
有限元法的优点是可以处理复杂的几何形状和材料特性，适用于各种边界条件和非线性问题。
然而，有限元法的计算过程相对较为复杂，需要对问题进行合理的离散化和网格划分，同时对于大规模问题，计算量也较大。
在实际应用中，根据具体问题的特点和求解要求，选择合适的数值计算方法是十分重要的。
对于简单的电磁场问题，如一维导线的电流分布，可以选择有限差分法进行求解。
而对于复杂的电磁场问题，如三维空间中的电磁波传播，有限元法更适合。
此外，有限差分法和有限元法还可以结合使用，通过将两种方法的优点相结合，提高计算精度和效率。
除了理论原理和应用范围，有限差分法和有限元法的优缺点也值得关注。
有限差分法的优点是简单易懂，计算过程直观，而且对于一些简单问题可以得到较为准确的结果。
然而，由于差分法中的网格离散化会引入一定的误差，对于复杂问题的求解精度有限。
相比之下，有限元法可以处理复杂的几何形状和材料特性，适用于各种边界条件和非线性问题，计算精度较高。
然而，有限元法的计算过程相对复杂，需要对问题进行合理的离散化和网格划分，同时对于大规模问题计算量较大。