工程电磁场数值分析试题
一、一同轴电缆，内导体（铜）外半径为0.01m，外导体（铜）内半径为0.03m，
外导体厚度为0.003m，内外导体间有两层电介质，一层电介质为聚乙烯(
13
r
ε=)，
另外一层电介质为聚氯乙烯(
26
r
ε=)，两层电介质厚度均为0.01m，内导体电位为5kV，外导体电位为0V。
（1）试用有限元法求内外导体间的电位和电场强度分布，
（2）求此电缆中最大场强的位置和最大值，能否击穿电介质或发生局部击穿，（3）在不击穿的前提下，此电缆能承载的最大电压为多少？
分析：
参数设置：铜相对介电常数ε=1，电阻率ρ=1e-7Ω/m
聚乙烯相对介电常数ε=3，电阻率ρ=1e+13Ω/m
聚氯乙烯相对介电常数ε=6，电阻率ρ=1e+14Ω/m
（1）其中电位分布及场强分布如下：
通过定义路径（-0.033,0）到（0.033,0）分析其场强分布如下图：
电压分布如下图：
（2）从图中可以看出其场强最大值位于内导体外半径附近处取得距离内导体圆心0.0132m处取最大值422728V/m。聚乙烯击穿场强为35-50MV/m，聚氯乙烯击穿场强为20-35MV/m，计算可知无法击穿电介质。
（3）聚乙烯击穿场强为35-50MV/m，聚氯乙烯击穿场强为20-35MV/m，按最小值计算理论上0.01m距离上其耐压分别为350kv和200kv，所以电介质不会被击穿。如不击穿理论上应能够承压200kv。
二、一同轴电缆，内导体（铜）外半径为0.01m，外导体（铜）内半径为0.03m，
外导体厚度为0.003m，内外导体间有一层电介质，电介质为聚氯乙烯(
26
ε=)，
电介质厚度均为0.02m，内导体电位为5kV，外导体电位为0V。
（1）试用有限元法求内外导体间的电位和电场强度分布，
（2）求此电缆中最大场强的位置和最大值，能否击穿电介质或发生局部击穿，（3）在不击穿的前提下，此电缆能承载的最大电压为多少？
（4）通过一题和二题的对比，说明同轴电缆的内外导体间用一层还是二层电介质比较好，为什么？
分析：
参数设置：铜相对介电常数ε=1，电阻率ρ=1e-7Ω/m
聚氯乙烯相对介电常数ε=6，电阻率ρ=1e+14Ω/m （1）其中电位分布及场强分布如下：
通过定义路径（-0.033,0）到（0.033,0）分析其场强分布如下图：
电压分布如下图：
（2）从图中可以看出其场强最大值位于内导体外半径附近处取得距离内导体圆心0.0132m处，最大值343991.6V/m。聚氯乙烯击穿场强为20-35MV/m，计算可知无法击穿电介质。
（3）聚氯乙烯击穿场强为20-35MV/m，按最小值计算理论上0.01m距离上其耐压分别200kv，所以电介质被击穿。即能够承压200kv。
（4）通过一、二两题的对比，可知同轴电缆的内外导体间用两层电介质比较好，使用一层电介质时场强的峰值要小于使用两层电介质，但使用两层电介质时场强下降更为迅速，到达外层时场强也相对更低，则发生击穿的可能性相对更低，所以使用两层电介质使同轴电缆工作性能更强。
三、一同轴电缆，内导体（铜）外半径为0.01m，外导体（铜）内半径为0.03m，
外导体厚度为0.003m，内外导体间有两层电介质，一层电介质为聚乙烯(
13
ε=)，
另外一层电介质为聚氯乙烯(
26
ε=)，两层电介质厚度均为0.01m，若内导体通电
流密度J1=500A/m2，外导体通电流密度为J2=-500A/m2，求内外导体间的磁场分布，画出一路经找出最大磁场和最小磁场位置。
分析：
参数设置:铜相对磁导率μ=1，电阻率ρ=1e-7Ω/m
聚乙烯相对磁导率μ=1，电阻率ρ=1e+13Ω/m
聚氯乙烯相对磁导率μ=1，电阻率ρ=1e+14Ω/m （1）其磁场分布如下图所示：
（2）定义路径（-0.033，0）到（0.033，0）可得如下磁场分布图：
数值分析试题 一、 填空题（2 0×2′） 1. ?? ????-=? ?????-=32,1223X A 设x =是精确值x *=的近似值，则x 有 2 位 有效数字。 2. 若f (x )=x 7－x 3＋1，则f [20,21,22,23,24,25,26,27]= 1 ， f [20,21,22,23,24,25,26,27,28]= 0 。 3. 设，‖A ‖∞＝___5 ____，‖X ‖∞＝__ 3_____， ‖AX ‖∞≤_15_ __。 4. 非线性方程f (x )=0的迭代函数x =?(x )在有解区间满足 |?’(x )| <1 ，则使用该迭代 函数的迭代解法一定是局部收敛的。 5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。 6. 当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差商 公式的 前插公式 ，若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的 后插公式 ；如果要估计结果的舍入误差，应该选用插值公式中的 拉格朗日插值公式 。 7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是：=∑=n i i x a 0)( 1 ；所以当 系数a i (x )满足 a i (x )>1 ，计算时不会放大f (x i )的误差。 8. 要使 20的近似值的相对误差小于%，至少要取 4 位有效数字。 9. 对任意初始向量X (0)及任意向量g ，线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…)收 敛于方程组的精确解x *的充分必要条件是 ?(B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。 11. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。 12. 线性方程组的松弛迭代法是通过逐渐减少残差r i (i =0,1,…,n )来实现的，其中的残差 r i ＝ (b i -a i1x 1-a i2x 2-…-a in x n )/a ii ，(i =0,1,…,n )。 13. 在非线性方程f (x )=0使用各种切线法迭代求解时，若在迭代区间存在唯一解，且f (x )
《电磁场数值分析》（期末作业） --- 2019学年 --- 学院：电子工程学院 学号： 姓名： 联系方式： 任课教师： 2019年5月
作业1 模拟真空中二维TM 电磁波的传播，边界设置为一阶Mur 吸收边界，观察电磁波的传播过程。波源为正弦函数： sin()sin(2)25 z t c E t n t ωπ ==? 代码： clc clear close all xmesh =150; ymesh =150; mu0=4*pi*1.0E-7; eps0=8.85E-12;C= 3.0E8; dx=1.0; dt=0.7*dx/C; timestep=150; ez( 1:xmesh+1,1:ymesh+1 ) = 0.0; hx( 1:xmesh+1,1:ymesh ) = 0.0; hy( 1:xmesh,1:ymesh+1 ) = 0.0; coef1 = dt/( mu0 * dx ); coef2 =dt/( eps0 * dx );coef3=(C*dt-dx)/(C*dt+dx); ez1=ez; for now = 1 : timestep hx = hx - coef1 * ( ez( :, 2 : ymesh+1 ) - ez( :, 1 : ymesh ) ); hy = hy + coef1 * ( ez(2 : xmesh+1, : ) - ez(1 : xmesh, : )); ez( 2 : xmesh , 2 : ymesh ) = ez( 2 : xmesh , 2 : ymesh ) - ... coef2 * ( hx( 2 : xmesh, 2 : ymesh ) - hx( 2 : xmesh , 1 : ymesh - 1) ) + ... coef2 * ( hy( 2 : xmesh ,2 : ymesh ) - hy( 1 : xmesh - 1,2 : ymesh) ); ez(1,:)=ez1(2,:)+coef3*(ez(2,:)-ez1(1,:)); ez(xmesh+1,:)=ez1(xmesh,:)+coef3*(ez(xmesh,:)-ez1(xme sh+1,:)); ez(:,1)=ez1(:,2)+coef3*(ez(:,2)-ez1(:,1));
注：考生属哪种类别请划“√” （博士、在校硕士、工程硕士、师资硕士、同等学力、研究生班） 辽宁工程技术大学 研究生考试试卷 考试科目：电磁场数值分析 考生班级：电控研 考生姓名： 学号： 考试分数： 注意事项 1、考前研究生将上述项目填写清楚 2、字迹要清楚，保持卷面清洁 3、试题、试卷一齐交监考老师 4、教师将试题、试卷、成绩单，一起送研究生学院； 专业课报所在院、系
直流无刷电机的内部电磁分析 1提出问题 在电磁学里，电磁场是一种由带电物体产生的一种物理场。处于电磁场的带电物体会感受到电磁场的作用力。电磁场与带电物体之间的相互作用可以用麦克斯韦方程和洛伦兹力定律来描述。电磁场是有内在联系、相互依存的电场和磁场的统一体的总称。随时间变化的电场产生磁场，随时间变化的磁场产生电场，两者互为因果，形成电磁场。 电磁场可由变速运动的带电粒子引起，也可由强弱变化的电流引起，不论原因如何，电磁场总以光速向四周传播，形成电磁波。电磁场是电磁作用的媒介，具有能量和动量，是物质存在的一种形式，电磁场的性质、特征及其运动变化规律由麦克斯韦方程确定。 ANSYS软件提供了图形用户界面与命令流两种方式来分析电机电磁场问题。在电机电磁场计算中,命令流方式和图形用户界面方式相比,具有以下优点:通用性好,对于同系列、同型号的电机电磁场计算只要对电机的尺寸参数进行修改即可,而采用ANSYS的图形用户界面方式进行电机电磁场计算,每次计算都要重新输入图形,没有通用性; 通过合理应用ANSYS的APDL语言编写一个两重循环程序就可实现转子自动旋转和自动施加励磁电流的功能,与ANSYS的图形用户界面方式相比,减少了人机交互的次数,缩短了计算时间。 电机的电磁分析，常用的软件是Maxwell，他是一个功能强大、灵活的,融结构、热、流体、电磁、声学于一体的大型通用有限元分析软件。广泛用于核工业、石油化工、航空航天、国防军工、机械制造、土木工程等一般工业及科学研究领域的设计分析。 本次作业中，将对直流无刷电机的内部电磁进行分析，采用Maxwell3D来建模，并进行磁场分析。 2直流无刷电机 直流无刷电机被广泛的用于日常生活用具、汽车工业、航空、消费电子、医学电子、工业自动化等装置和仪表。顾名思义，直流无刷电机不使用机械结构的换向电刷而直接使用电子换向器，在使用中直流无刷电机相比有刷电机有许多
1、（本题5分）试确定7 22 作为π的近似值具有几位有效数字，并确定其相对误差限。 解 因为 7 22 =3.142857…=1103142857 .0-? π=3.141592… 所以 312102 11021005.0001264.0722--?=?=<=- π （2分） 这里，3,21,0=-=+-=n n m m 由有效数字的定义可知7 22 作为π的近似值具有3位有效数字。 （1分） 而相对误差限 3102 1 0005.00004138.0001264.07 22-?= <≈= -= π π πε r （2分） 2、（本题6分）用改进平方根法解方程组：??? ?? ??=????? ??????? ??--654131*********x x x ； 解 设???? ? ??????? ? ?????? ??===????? ??--11111 1 131321112323121 32 132 31 21 l l l d d d l l l LDL A T 由矩阵乘法得： 5 7,21,215 27 ,25,2323121321- ==-== -==l l l d d d （3分） 由y D x L b Ly T 1 ,-==解得 T T x y )9 23 ,97,910(,)563, 7,4(== （3分） 3、（本题6分）给定线性方程组???????=++-=+-+=-+-=-+17 7222382311387 510432143213 21431x x x x x x x x x x x x x x 1）写出Jacoib 迭代格式和Gauss-Seidel 迭代格式； 2）考查Jacoib 迭代格式和Gauss-Seidel 迭代格式的敛散性； 解 1）Jacoib 迭代格式为
电磁场数值分析方法的若干研究张晨颜 发表时间：2018-05-23T17:08:53.507Z 来源：《基层建设》2018年第8期作者：张晨颜 [导读] 摘要：文章主要阐述电磁场分析的重要方法，即电磁场数值计算方法。 中国矿业大学孙越崎学院电气系江苏省徐州市 221000 摘要：文章主要阐述电磁场分析的重要方法，即电磁场数值计算方法。具体是对直接积分法、有限差分法以及有限元法三种常见分析方法的原理以及优势特征进行探究。同时采用不同方法技能型求解，并对结果进行简要分析，发现只有在合理应用前处理技术基础上，各种数值分析方法的计算精确性才会有所保障。 关键词：电磁场；数值分析方法；直接积分法；有限差分法；有限元法 测算与处理电磁场边值问题的方法主要有模拟法、图解法、解析法、数值法四种类型。前处理、计算和后处理是所有工程电磁场数值分析的三大要素。在静态条件下，电磁场分布均可归纳在一定边界条件下求解Poisson或Laplace方程。阐述电磁场的麦克斯韦方程组有有微分与积分两种类型，不同数值计算方法的计算量离散化所参照的基本方程形式，可以细化为积分方程法与微分方程法。本文对相关计算过程实施简化措施，并对不同电磁场数值计算的方法优势与弊端进行归纳。 1直接积分法 ③选择一定的代数解法(通常应用迭代法)，编写相关计算流程，以获得相应待求边值问题的差分方程组，得到边值问题的数值解. 有限差分法的主要内容通常涵盖三个方面：①差分方程的形成；②边界条件的处理；③方程的求解。差分方程的推导通常采用泰勒级数法。将电磁场的微分方程形式——泊松方程或拉普拉斯方程设为初始点，借助展开泰勒绿数的方式，列算差分方程。结合现存的边界条件，结合具体情况修整边界上的节点的差分方程形式。最后是对代数方程组———差分方程进行计算以获得最终结果。同步迭代法、异步迭代法和超松弛迭代法石常见解题方法。通常采用点超松弛迭代法和线迭超松弛迭代法。但应用过程中的重点是合理选择松弛因子，只有在选择得当时迭代加速进程才会得到有效管控。 有限差分法的优点是能够较为快速的找出差分方程组，同时差分方程组自体也体现出简洁化特征，网格的剖分过程也没有太大技术含量，数据信息准备工作不会耗用太多时间，计算流程制定相对简易。但是对于曲线边界等不规则的边界，处理难度会相应增加。若区域的边缘线以及内部媒介分界线形体较为繁杂，并且场域布设形式多变时，因为差分法的网格剖分灵敏性较差，故此计算结果的测算过程将会受到层层阻碍。有限差分法适用于对象有如下几种类型：①边界形状规则的第一类边界，第二类齐次边界；②静态场，时变场；③线性场，非线性场等。 3有限元法 有限元法是采用变分原理和离散化去获得近似解的方法。电磁场的问题通常都可总结为求解的偏微分方程的边值问题。有限元法不是采用直接偏微分方程去求解电磁场的，其将偏微分方程边值问题设为始发点，探寻一个能量泛函的积分式，并促使其其在满足第一类边界条件的前提提取取极值，等同于构建条件变分问题。这个条件变分问题等同于偏微分方程边值问题。在求解过程中，将场的求解区域细化陈可以量化的单元，在每一单元中，策略的认为对每一点的求解函数是在单元节点的函数值间随坐标变化而产生相应变化的。故此插值函数在单元格中产生，把插值函数整合到能量泛函的积分式，继而将泛函离散化转型为数个多元函数。继而求解极值。借此方式获得一个代数方程组。最后由此方程组求解得到数值解。对第二有限元法是结合变分原理和离散化而获得相似值解的方法。若场域中存有不同的
数值分析整理版试题及答案
例1、 已知函数表 x -1 1 2 ()f x -3 0 4 求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。 解： （1）k x -1 1 2 k y -3 0 4 插值基函数分别为 ()()()()()()()()()() 1200102121()1211126 x x x x x x l x x x x x x x ----= ==-------- ()()()()()()()() ()()021******* ()1211122x x x x x x l x x x x x x x --+-= ==-+---+- ()()()()()()()()()()0122021111 ()1121213 x x x x x x l x x x x x x x --+-= ==-+--+- 故所求二次拉格朗日插值多项式为 () ()()()()()()()()()()2 20 2()11131201241162314 121123537623k k k L x y l x x x x x x x x x x x x x ==?? =-? --+?-+-+?+-????=---++-=+-∑ （2）一阶均差、二阶均差分别为
[]()()[]()()[][][]010********* 011201202303 ,11204 ,412 3 4,,5 2,,126 f x f x f x x x x f x f x f x x x x f x x f x x f x x x x x ---===-----= = =----=== --- k x ()k f x 一阶 二阶 -1 -3 1 0 3/ 2 2 4 4 5/6 故所求Newton 二次插值多项式为 ()()[]()[]()() ()()()20010012012,,,35 311126537623P x f x f x x x x f x x x x x x x x x x x x =+-+--=-+ +++-=+- 例2、 设2 ()32f x x x =++，[0,1]x ∈，试求()f x 在[0, 1]上关于()1x ρ=，{} span 1,x Φ=的最佳平方逼近多项式。 解： 若{}span 1,x Φ=，则0()1x ?=，1()x x ?=，且()1x ρ=，这样，有
第一章 一、矢量代数 A ?B =AB cos θ A B ?= AB e AB sin θ A ?(B ?C ) = B ?(C ?A ) = C ?(A ?B ) ()()()C A C C A B C B A ?-?=?? 二、三种正交坐标系 1. 直角坐标系 矢量线元x y z =++l e e e d x y z 矢量面元=++S e e e x y z d dxdy dzdx dxdy 体积元d V = dx dy dz 单位矢量的关系?=e e e x y z ?=e e e y z x ?=e e e z x y 2. 圆柱形坐标系 矢量线元=++l e e e z d d d dz ρ?ρρ?l 矢量面元=+e e z dS d dz d d ρρ?ρρ? 体积元dz d d dV ?ρρ= 单位矢量的关系?=??=e e e e e =e e e e z z z ρ??ρ ρ? 3. 球坐标系 矢量线元d l = e r d r + e θ r d θ + e ? r sin θ d ? 矢量面元d S = e r r 2sin θ d θ d ? 体积元 ?θθd d r r dV sin 2= 单位矢量的关系?=??=e e e e e =e e e e r r r θ? θ??θ 三、矢量场的散度和旋度 1. 通量与散度 =?? A S S d Φ 0 lim ?→?=??=??A S A A S v d div v 2. 环流量与旋度 =??A l l d Γ max n rot =lim ?→???A l A e l S d S 3. 计算公式 ????= ++????A y x z A A A x y z 11()z A A A z ?ρρρρρ?????= ++????A 22111()(s i n )s i n s i n ????= ++????A r A r A A r r r r ? θ θθθθ? x y z ? ????= ???e e e A x y z x y z A A A 1z z z A A A ρ?ρ?ρρ?ρ? ?? ??= ???e e e A
电磁场数值计算方法地发展及应用 专业：电气工程 姓名：毛煜杰 学号： 一、电磁场数值计算方法产生和发展地必然性 麦克斯韦尔通过对以往科学家们对电磁现象研究地总结，认为原来地研究工作缺乏严格地数学形式，并认为应把电流地规律与电场和磁场地规律统一起来.为此，他引入了位移电流和涡旋场地概念，于年提出了电磁场普遍规律地数学描述—电磁场基本方程组，即麦克斯韦尔方程组.它定量地刻画了电磁场地转化和电磁波地传播规律.麦克斯韦尔地理论奠定了经典地电磁场理论，揭示了电、磁和光地统一性.资料个人收集整理，勿做商业用途 但是，在电磁场计算地方法中，诸如直接求解场地基本方程—拉普拉斯方程和泊松方程地方法、镜象法、复变函数法以及其它种种解析方法，其应用甚为局限，基本上不能用于求解边界情况复杂地、三维空间地实际问题.至于图解法又欠准确.因此，这些电磁场地计算方法在较复杂地电磁系统地设计计算中，实际上长期未能得到有效地采用.于是，人们开始采用磁路地计算方法，在相当长地时期内它可以说是唯一实用地方法.它地依据是磁系统中磁通绝大部分是沿着以铁磁材料为主体地“路径”—磁路“流通”.这种计算方法与电路地解法极其相似，易于掌握和理解，并得以沿用至今.然而，众所周知，对于磁通是无绝缘体可言地，所以磁路实际上是一种分布参数性质地“路”.为了将磁路逼近实际情况，当磁系统结构复杂、铁磁材料饱和时，其计算十分复杂.资料个人收集整理，勿做商业用途 现代工业地飞速发展使得电器产品地结构越来越复杂，特殊使用场合越来趁多.电机和变压器地单机容量越来越大，现代超导电机和磁流体发电机必须用场地观点和方法去解决设计问题.由于现代物理学地发展，许多高精度地电磁铁、波导管和谐振腔应用到有关设备中，它们不仅要赋与带电粒子能量，并且要有特殊地型场去控制带电粒子地轨迹.这些都对电磁系统地设计和制造提出了新地要求，传统地分析计算方法越来越感到不足，这就促使人们发展经典地电磁场理论，促使人们用场地观点、数值计算地方法进行定量研究.资料个人收集整理，勿做商业用途 电子计算机地出现为数值计算方法地迅速发展创造了必不可少地条件.即使采用“路”地方法来计算，由于计算速度地加快和新地算法地应用，不仅使得计算精度得到了很大地提高，而且使得工程设计人员能从繁重地计算工作中解脱出来.从“场”地计算方面来看，由于很多求解偏微分方程地数值方法，诸如有限差分法、有限元法、积分方程法等等地运用，使得大量工程电磁场问题有可能利用数值计算地方法获得符合工程精度要求地解答，它使电磁系纯地设计计算地面貌焕然一新.电磁场地各种数值计算方法正是在计算机地发展、计算数学地前进和工程实际问题不断地提出地情况下取得一系列进展地.资料个人收集整理，勿做商业用途 二、电磁场数值计算方法地发展历史 电磁场数值计算已发展了许多方法，主要可分为积分法（积分方程法、边界积分法和边界元法）、微分法（有限差分法、有限元法和网络图论法等）及微分积分法地混合法.资料个人收集整理，勿做商业用途 年，利用向量位，采用有限差分法离散，求解了二维非线性磁场问题.随后和用该程序设计了同步加速器磁铁，并把它发展成为软件包.此后，采用有限差分法计算线性和非线性二维场地程序如雨后春笋般地在美国和西欧出现.有限差分法不仅能求解均匀线性媒质中地位场，还能解决非线性媒质中地场；它不仅能求解恒定场和似稳场，还能求解时变场.在边值问题地数位方法中，此法是相当简便地.在计算机存储容量许可地情况下，采取较精细地网格，使离散化模型较精确地逼近真实问题，可以获得足够精度地数值解.但是, 当场城几何特
1. 已知325413.0,325413*2*1==X X 都有6位有效数字，求绝对误差限。（4分） 解： 由已知可知6 5.0102 1 ,0,6,10325413.0016*1=?= =-=?=ε绝对误差限n k k X 2分 620*2102 1 ,6,0,10325413.0-?= -=-=?=ε绝对误差限n k k X 2分 2. 已知?? ???=0 01 A 220- ?????440求21,,A A A ∞ （6分） 解： {}, 88,4,1max 1==A 1分 {}, 66,6,1max ==∞A 1分 () A A A T max 2λ= 1分 ?????=0 1 A A T 4 2 ???? ? -420?????0 01 2 20 - ???? ?440= ?????0 01 80 ???? ?3200 2分 {}32 32,8,1max )(max ==A A T λ
1分 24322==A 3. 设32)()(a x x f -= （6分） ① 写出f(x)=0解的迭代格式 ② 当a 为何值时，)(1k k x x ?=+ （0,1……）产生的序列{}k x 收敛于 2 解： ①迭代格式为： x a x x x a x a x x a x x x f x f x x k k k k k k k k k k 665)(665)(6)()(')(2 2 32 1 += +=---=-=+? 3 分 ②时迭代收敛即当222,112 10)2(',665)('2<<-<-=-= a a x a x ?? 3分 4. 给定线性方程组，其中：?? ?=13A ?? ?2 2，?? ? ???-=13b 用迭代公式 )()()()1(k k k Ax b x x -+=+α（0,1……）求解，问取什么实数α ，可使 迭代收敛 （8分） 解： 所给迭代公式的迭代矩阵为?? ? --???--=-=ααααα21231A I B 2分
例1、 已知函数表 求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。 解： （1） 故所求二次拉格朗日插值多项式为 （2）一阶均差、二阶均差分别为 例2、 设2 ()32f x x x =++，[0,1]x ∈，试求()f x 在[0,1]上关于()1x ρ=，{}span 1,x Φ=的最佳平 方逼近多项式。 解： 若{}span 1,x Φ=，则0()1x ?=，1()x x ?=，且()1x ρ=，这样，有 所以，法方程为
011231261192 34a a ??????????=?????????? ?????????? ，经过消元得012311 62110123a a ??? ???????=???????????????????? 再回代解该方程，得到14a =，011 6 a = 故，所求最佳平方逼近多项式为* 111()46S x x =+ 例3、 设()x f x e =，[0,1]x ∈，试求()f x 在[0,1]上关于()1x ρ=，{}span 1,x Φ=的最佳平方逼近 多项式。 解： 若{}span 1,x Φ=，则0()1x ?=，1()x x ?=，这样，有 所以，法方程为 解法方程，得到00.8732a =，1 1.6902a =， 故，所求最佳平方逼近多项式为 例4、 用4n = 的复合梯形和复合辛普森公式计算积分1 ? 。 解： （1）用4n =的复合梯形公式 由于2h =，( )f x =()121,2,3k x k k =+=，所以，有 （2）用4n =的复合辛普森公式 由于2h =，( )f x =()121,2,3k x k k =+=，()12 220,1,2,3k x k k + =+=，所以，有 例5、 用列主元消去法求解下列线性方程组的解。 解：先消元 再回代，得到33x =，22x =，11x = 所以，线性方程组的解为11x =，22x =，33x = 例6、 用直接三角分解法求下列线性方程组的解。 解： 设 则由A LU =的对应元素相等，有 1114u = ，1215u =，1316u =， 2111211433l u l =?=，3111311 22 l u l =?=， 2112222211460l u u u +=?=-，2113232311 545l u u u +=?=-，
电磁场数值分析 电和磁现象在自然界普遍存在，两者相互依存形成一个不看分割的整体。电能产生磁，磁能生电。很早以前人们就注意到电现象和磁现象，但是两者之间的这种相互联系在很长的一段时间内都没有被人们认识。直到奥斯特首先发现了通电直导线周围存在磁场这一现象人们才开始把电和磁放在一起来研究。然而这个时候人们依然没有办法揭示电和磁中间的秘密，只是停留在实验研究阶段，没有形成科学的理论。1831年法拉第发现了电磁感应定律，从此电和磁的计算可以量化了，人类历史也开启了一个新的时代—电气时代。由于法拉第的杰出工作，电和磁不再是不可触摸的了，人们已经掌握了运用它的钥匙。在法拉第之后，另一位杰出的科学家麦克斯韦则更进一步，建立了麦克斯韦方程组，电和磁的理论已经到了相当完美的程度。 现代电机，不管结构多么复杂，都是基于法拉第电磁感应定律和麦克斯韦方程组的原理来运行的，其电和磁的相关量都可以利用这两个定律来进行精确地分析，在设计电机时，我们也是基于这两个定律对电机的电磁过程来进行精确的设计，从而设计出理想的电机。 学会电磁场分析，主要是基于麦克斯韦方程组的相关计算，对电机的学习非常重要。它为我们今后的学习打下基础。在学习过程中，主要要把握以下几个度之间的关系：梯度、旋度、散度，这三者的变换正体现了电和磁之间的转换。 一基本原理 电磁场的内在规律由电磁场基本方程组—麦克斯韦（Maxwell ）方程组表达。这些方程是由麦克斯韦对大量实验结果及基本概念进行了数学加工和推广归纳而成的。麦克斯韦方程组是分析和计算电磁场问题的出发点，它既可写成微分形式，又可写成积分形式。 微分形式的麦克斯韦方程组为 t D J H ??+=?? （1） t B E ??-=?? （2） 0=??B （3） ρ=??D （4）
电磁学基础知识 电场 一、场强E （矢量，与q 无关） 1．定义：E ＝ 单位：N/C 或V/m 方向：与＋q 所受电场力方向 电场线表示E 的大小和方向 2．点电荷电场：E ＝ 静电力恒量 k ＝ Nm 2/C 2 匀强电场：E ＝ d 为两点在电场线方向上的距离 3．E 的叠加——平行四边形定则 4．电场力（与q 有关） F ＝ 库仑定律：F ＝ （适用条件：真空、点电荷） 5．电荷守恒定律（注意：两个相同带电小球接触后，q 相等） 二、电势φ（标量，与q 无关） 1．定义：φA ＝ ＝ ＝ 单位：V 说明：φ＝单位正电荷由某点移到φ＝0处的W ⑴沿电场线，电势降低 ⑵等势面⊥电场线；等势面的疏密反映E 的强弱 2．电势叠加——代数和 3．电势差：U AB ＝ ＝ 4．电场力做功：W AB ＝ 与路径无关 5．电势能的变化：Δε＝W 电场力做正功，电势能 ；电场力做负功，电势能 需要解决的问题： ①如何判电势的高低以及正负（由电场线判断） ②如何判电场力做功的正负（由F 、v 方向判） ③如何判电势能的变化（由W 的正负判） 三、电场中的导体 1．静电平衡：远端同号，近端异号 2．静电平衡特点 ⑴E 内＝0；⑵E 表面 ⊥表面；⑶等势体（内部及表面电势相等）；⑷净电荷分布在外表面 四、电容器 1．定义：C ＝ （C 与Q 、U 无关） 单位：1 F ＝106 μF ＝1012 pF 2．平行板电容器： C ＝ 3．两类问题：①充电后与电源断开， 不变；②始终与电源相连， 不变 五、带电粒子在电场中的运动 1．加速：qU ＝ 2．偏转：v ⊥E 时，做类平抛运动 位移：L ＝ ； y ＝ ＝ ＝ 速度：v y ＝ ＝ ； v ＝ ； tan θ＝ 六、实验：描绘等势线 1．器材： 2．纸顺序：从上向下
数值分析试题及答案 一、单项选择题（每小题3分，共15分） 1. 3.142和3.141分别作为的近似数具有（）和（）位有效数字. A．4和3 B．3和2 C．3和4 D．4和4 2. 已知求积公式，则＝（） A． B．C．D． 3. 通过点的拉格朗日插值基函数满足（） A．＝0，B．＝0， C．＝1，D．＝1， 4. 设求方程的根的牛顿法收敛，则它具有（）敛速。 A．超线性B．平方C．线性D．三次 5. 用列主元消元法解线性方程组作第一次消元后得到的第3个方程（）. A．B． C．D． 单项选择题答案 1.A 2.D 3.D 4.C 5.B 得分评卷 人 二、填空题（每小题3分，共15分） 1. 设, 则， . 2. 一阶均差 3. 已知时，科茨系数，那么 4. 因为方程在区间上满足，所以在区间内有根。 5. 取步长，用欧拉法解初值问题的计算公式.填空题答案
1. 9和 2. 3. 4. 5. 得分评卷 人 三、计算题（每题15分，共60分） 1. 已知函数的一组数据：求分段线性插值函数，并计算的近似值. 计算题1.答案 1. 解， ， 所以分段线性插值函数为 2. 已知线性方程组 （1）写出雅可比迭代公式、高斯－塞德尔迭代公式； （2）对于初始值，应用雅可比迭代公式、高斯－塞德尔迭代公式分别计算（保留小数点后五位数字）. 计算题2.答案 1.解原方程组同解变形为 雅可比迭代公式为 高斯－塞德尔迭代法公式 用雅可比迭代公式得 用高斯－塞德尔迭代公式得 3. 用牛顿法求方程在之间的近似根 （1）请指出为什么初值应取2？ （2）请用牛顿法求出近似根，精确到0.0001. 计算题3.答案
第一章 绪论 姓名 学号 班级 习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。 1 若误差限为5105.0-?,那么近似数有几位有效数字（有效数字的计算） 解：2*103400.0-?=x ，325*102 1102 1---?=?≤-x x 故具有3位有效数字。 2 14159.3=π具有4位有效数字的近似值是多少（有效数字的计算） 解：10314159.0?= π，欲使其近似值*π具有4位有效数字，必需 41*102 1 -?≤-ππ，3*3102 1102 1--?+≤≤?-πππ，即14209.314109.3*≤≤π 3 已知2031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +， b a ?有几位有效数字（有效数字的计算） 解：3*1021 -?≤-a a ，2*102 1-?≤-b b ，而1811.2=+b a ，1766.1=?b a 2123****102 1 10211021)()(---?≤?+?≤ -+-≤+-+b b a a b a b a 故b a +至少具有2位有效数字。 2123*****102 1 0065.01022031.1102978.0)()(---?≤=?+?≤ -+-≤-b b a a a b b a ab
故b a ?至少具有2位有效数字。 4 设0>x ，x 的相对误差为δ，求x ln 的误差和相对误差（误差的计算） 解：已知δ=-* *x x x ，则误差为 δ=-= -* **ln ln x x x x x 则相对误差为 * * ** * * ln ln 1ln ln ln x x x x x x x x δ = -= - 5测得某圆柱体高度h 的值为cm h 20*=，底面半径r 的值为cm r 5*=， 已知cm h h 2.0||*≤-，cm r r 1.0||*≤-，求圆柱体体积h r v 2π=的绝对误差 限与相对误差限。（误差限的计算） 解：*2******2),(),(h h r r r h r r h v r h v -+-≤-ππ 绝对误差限为 πππ252.051.02052)5,20(),(2=??+????≤-v r h v 相对误差限为 %420 1 20525) 5,20() 5,20(),(2 ==??≤ -ππv v r h v 6 设x 的相对误差为%a ,求n x y =的相对误差。（函数误差的计算） 解：%* *a x x x =-， )%(* **** *na x x x n x x x y y y n n n =-≤-= - 7计算球的体积，为了使体积的相对误差限为%1，问度量半径r 时允许的相对误差限为多大（函数误差的计算）
第一章 电机中的电磁学基本知识 1.1 磁路的基本知识 1.1.1 电路与磁路 对于电路系统来说，在电动势E 的作用下电流I 从E 的正极通过导体流向负极。构成一个完整的电路系统需要电动势、电导体，并可以形成电流。 在磁路系统中，也有一个磁动势F （类似于电路中的电势），在F 的作用下产生一个 Φ（类似于电路中的电流），磁通Φ从磁动势的N 极通过一个通路（类似于电路中的导 体）到S 极，这个通路就是磁路。由于铁磁材料磁导率比空气大几千倍，即空气磁阻比铁磁材料大几千倍，所以构成磁路的材料均使用导磁率高的铁磁材料。然而非铁磁物质，如空气也能通过磁通，这就造成铁磁材料构成磁路的周围空气中也必然会有磁通σΦ（，由于空气磁阻比铁磁材料大几千倍，因而σΦ比Φ小的多，σΦ常常被称为漏磁通，Φ称为主磁通。因此磁路问题比电路问题要复杂的多。 1.1.2 电机电器中的磁路 磁路系统广泛应用在电器设备之中，如变压器、电机、继电器等。并且在电机和某些电器的磁路中，一般还需要一段空气隙，或者说空气隙也是磁路的组成部分。 图1—1是电机电器的几种常用磁路结构。图(a)是普通变压器的磁路，它全部由铁磁材料组成；图(b)是电磁继电器磁路，它除了铁磁材料外，还有一段空气隙。 图(c)表示电机的磁路，也是由铁磁材料和空气隙组成；图(b)是无分支的串联磁路，空气隙段和铁磁材料串联组成；图(a)是有分支的并联磁路。图中实（或虚）线表示磁通的路径。 (a) (b) (c) 图1—1 几种常用电器的典型磁路 (a) 普通变压器铁芯； (b) 电磁继电器常用铁芯； (c) 电机磁路 1.1.3 电气设备中磁动势的产生 为了产生较强的磁场，在一般电气设备中都使用电流产生磁场。电流产生磁场的方法是：把绕制好的N 匝线圈套装在铁心上，并在线圈内通入电流i ，这样在铁心和线圈周围的空间中就会形成磁场，其中大多数磁通通过铁心，称为主磁通Φ；小部分围绕线圈，称为漏磁通σΦ，如图1—2所示。套装在铁心上用于产生磁通的N 匝线圈称为励磁线圈，励磁
电磁场的数值计算方法 摘要：数值计算方法是一种研究并解决数学问题数值近似解的方法，广泛运用于电气、军事、经济、生态、医疗、天文、地质等众多领域。本文综述了电磁场数值计算方法的发展历史、分类，详细介绍了三种典型的数值计算方法—有限差分法、有限元法、矩量法, 对每种方法的解题思路、原理、步骤、特点、应用进行了详细阐述, 并就不同方法的区别进行了深入分析, 最后对电磁场数值计算方法的应用前景作了初步探讨。 关键词：电磁场；数值计算；有限差分法；有限元法；矩量法 引言 自从1864年Maxwell建立了统一的电磁场理论，并得出著名的Maxwell方程以来，经典的数学分析方法是一百多年来电磁学学科发展中一个极为重要的手段, 围绕电磁分布边值问题的求解国内外专家学者做了大量的工作。在数值计算方法之前, 电磁分布的边值问题的研究方法主要是解析法，但其推导过程相当繁琐和困难,缺乏通用性,可求解的问题非常有限。上个世纪六十年代以来,伴随着电子计算机技术的飞速发展,多种电磁场数值计算方法不断涌现,并得到广泛地应用,相对于解析法而言,数值计算方法受边界形状的约束大为减少,可以解决各种类型的复杂问题。但各种数值计算方法都有一定的局限性,一个复杂的问题往往难以依靠一种单一方法解决，因此如何充分发挥各种方法的优势,取长补短,将多种方法结合起来解决实际问题,即混合法的研究和应用已日益受到人们的关注。本文综述电磁场的数值计算方法，对三种常用的电磁场数值计算方法进行分类和比较。 1电磁场数值计算方法的发展历史 在上世纪四十年代，就有人试探用数值计算的方法来求解具有简单边界的电磁场问题，如采用Ritz法[1]，以多项式在整个求解场域范围内整体逼近二阶偏微分方程在求解域中的解。五十年代，采用差分方程近似二阶偏微分方程，诞生了有限差分数值计算方法，开始是人工计算，后来采用机械式的手摇计算机计算，使简单、直观的有限差分法得到应用和发展，该方法曾在欧、美风行一时。1964年美国加州大学学者Winslow以矢量位为求解变量，用有限差分法在计算机上成
例1、已知函数表 求() f x的Lagrange二次插值多项式和Newton二次插值多项式。 解： （1）由题可知 插值基函数分别为 故所求二次拉格朗日插值多项式为 （2）一阶均差、二阶均差分别为 均差表为
故所求Newton 二次插值多项式为 例2、 设2()32f x x x =++，[0,1]x ∈，试求()f x 在[0, 1]上关于()1x ρ=，{}span 1,x Φ=的 最佳平方逼近多项式。 解： 若{}span 1,x Φ=，则0()1x ?=，1()x x ?=，且()1x ρ=，这样，有 所以，法方程为 011231261192 34a a ??????????=?????????? ?????????? ，经过消元得012311 62110123a a ??? ???????=???????????????????? 再回代解该方程，得到14a =，011 6 a = 故，所求最佳平方逼近多项式为* 111 ()46 S x x = +
例3、 设()x f x e =，[0,1]x ∈，试求()f x 在[0, 1]上关于()1x ρ=，{}span 1,x Φ=的最佳平 方逼近多项式。 解： 若{}span 1,x Φ=，则0()1x ?=，1()x x ?=，这样，有 所以，法方程为 解法方程，得到00.8732a =，1 1.6902a =， 故，所求最佳平方逼近多项式为 例4、 用4n =的复合梯形和复合辛普森公式计算积分1?。 解： （1）用4n =的复合梯形公式 由于 2h =，()f x =，()121,2,3k x k k =+=，所以，有 （2）用4n =的复合辛普森公式 由于2h =，()f x =，()121,2,3k x k k =+=，()1 2 220,1,2,3k x k k +=+=，所以，有 例5、 用列主元消去法求解下列线性方程组的解。 解：先消元 再回代，得到33x =，22x =，11x =
数值分析试题及答案汇 总 TTA standardization office【TTA 5AB- TTAK 08- TTA 2C】
数值分析试题 一、填空题（2 0×2′） 1. ?? ????-=? ?????-=32,1223X A 设x =是精确值x *=的近似值，则x 有 2 位有效数字。 2. 若f (x )=x 7－x 3＋1，则f [20,21,22,23,24,25,26,27]= 1 ， f [20,21,22,23,24,25,26,27,28]= 0 。 3. 设，‖A ‖∞＝___5 ____，‖X ‖∞＝__ 3_____， ‖AX ‖∞≤_15_ __。 4. 非线性方程f (x )=0的迭代函数x =(x )在有解区间满足 |’(x )| <1 ，则使用该迭代函数 的迭代解法一定是局部收敛的。 5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。 6. 当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差 商公式的 前插公式 ，若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的 后插公式 ；如果要估计结果的舍入误差，应该选用插值公式中的 拉格朗日插值公式 。 7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是：=∑=n i i x a 0)( 1 ；所以当系数 a i (x )满足 a i (x )>1 ，计算时不会放大f (x i )的误差。 8. 要使 20的近似值的相对误差小于%，至少要取 4 位有效数字。 9. 对任意初始向量X (0)及任意向量g ，线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…) 收敛于方程组的精确解x *的充分必要条件是 (B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。
电磁场与电磁波课程知识点总结与主要公式 1 麦克斯韦方程组的理解和掌握 （1）麦克斯韦方程组 ??????=?=??=?=?????-=???- =?????+=???+ =??s s l s l s s d B B Q s d D D s d t B l d E t B E s d t D J l d H t D J H 0 )( ρ 本构关系： E J H B E D σμε=== （2）静态场时的麦克斯韦方程组（场与时间t 无关） ????=?=??=?=??=?=??=?=??s s l l s d B B Q s d D D l d E E I l d H J H 0 000 ρ 2 边界条件 （1）一般情况的边界条件 n n n sT t t s n s n n s n t t n B B B B a J H H J H H a D D D D a E E E E a 21212121212121210 )())(0 )==-?=-=-?=-=-?==-? （（ρρ （2）介质界面边界条件（ρs = 0 J s = 0） n n n t t n n n n t t n B B B B a H H H H a D D D D a E E E E a 21212121212121210 )(0 )0 )(0 )==-?==-?==-?==-? （（
（1）基本方程 00 2 2 =?==?- =?=?=??=?=??? ??A A p s l l d E Q s d D D l d E E ???ε ρ ?ρ 本构关系： E D ε= （2）解题思路 ● 对称问题（球对称、轴对称、面对称）使用高斯定理或解电位方程（注 意边界条件的使用）。 ● 假设电荷Q ——> 计算电场强度E ——> 计算电位φ ——> 计算能 量ωe =εE 2/2或者电容（C=Q/φ）。 （3）典型问题 ● 导体球（包括实心球、空心球、多层介质）的电场、电位计算； ● 长直导体柱的电场、电位计算； ● 平行导体板（包括双导体板、单导体板）的电场、电位计算； ● 电荷导线环的电场、电位计算； ● 电容和能量的计算。 例 ： ρ s 球对称 轴对称 面对称