第二十一讲广义特征值与极小极大原理

特征值解法

《结构动力学》大作业 结构大型特征值问题的求解 0810020035 吴亮秦 1振动系统的特征值问题 1.1实特征值问题 n 自由度无阻尼线性振动系统的运动微分方程可表示为： []{}[]{}()M u K u F t += （1.1） 其中，{}u 是位移向量，[]M 和[]K 分别是系统的质量矩阵和刚度矩阵，都是n 阶正定矩阵，()F t 是激励向量。 此系统的自由振动微分方程为 []{}[]{}0M u K u += （1.2） 设其主振型为： {}{}sin()u v t ω?=+ （1.3） 其中，{}v 为振幅向量，ω为圆频率，?为初相位。将(1.3)代入自由振动微分方程(1.2)， 得： []{}[]{}K v M v λ= （1.4） 其中2 λω=，(1.4)具有非零解的条件是 ()[][]det 0M K λ-= (1.5) 式(1.4)称为系统的特征方程，由此可以确定方程的n 个正实根1{}n i i λ=，称为系统的特征值，1{}n i i ω=称为系统的固有频率，{}i v （i=1,2,…..n ）为对应于特征值的特征向量或称为系统的振型或模态。 因为[]M 矩阵正定，则[]M 有Cholesky 分解： [][][]T M L L = (1.6) 其中，[]L 是下三角矩阵。引入向量{}x 满足：{}[]{}T x L v =，则： 1 {}([]){}T v L x -= (1.7) 代入(1.4)，得： ([][]){}0I P x λ-= (1.8) 其中，( ) 1 1 [][][][] T P L K L --=，式(1.8)称为标准实特征值问题。 1.2复特征值问题 多自由度阻尼自由振动系统的运动方程为如下二阶常系数微分方程组： []{()}[]{()}[]{()}0 M x t C x t K x t ++= (1.9) 其中 []M ，[]C ，[]K 分别是n 阶的质量、阻尼和刚度矩阵，{()}q t 是n 维可微向量函数。用分离变量法，设{()}{}t x t e λφ=，其中{}φ是与时间t 无关的常向量，λ为待定参数。将

计算电磁学作业_二)

计算电磁学课程作业(二) 1. 电磁场的线性系统（满足标量亥姆霍兹方程的系统）与一般电 子线性系统有何异同点？ 2. 试阐述格林函数对工程电磁场计算和求解的意义。 3. 任何源函数都可很方便地表示为基本函数（一般为函数）的线 性组合。任何波函数都可很方便地表示为基本函数（各种谐函 数）的线性组合。利用电磁场线性系统的函数和格林函数， 对于矢量磁位的亥姆霍兹方程： ，其在自由空间的解为 试写出两个有关矢量磁位的结论。 4. 对于无源区，电场、磁场、矢量磁位、标量电位、矢量电 位、标量磁位以及德拜位、赫兹矢量位等波函数，在时 域均可以写成矢量达朗伯方程的形式： 或标量达朗伯方程的形式。 对于矢量达朗伯方程，也常常只对标量达朗伯方程进行讨论和求解。这是因为：一方面矢量方程可以通过分离变量法后看做各个坐标分量标量方程的叠加；另一方面不同的波函数（平面波、柱面波、球面波）之间可以相互转换表达或相互展开表示（通过广义傅里叶变换）。 试写出无源区标量达朗伯方程的一个通解形式及其推导过程，并阐述通解的物理含义。 5. 类似地，在无源区，频域中波函数的波动方程可以表达为标量 亥姆霍兹方程（谐方程）： （） 其解在为谐函数（正弦函数、余弦函数、指数函数或柱谐函数、 球谐函数）。 电磁波在无限空间传播与存在的是连续谱；而电磁波在有限空 间传播与存在的是分立谱。试分别写出无源区的标量亥姆霍兹方程在直

角坐标、柱坐标和球坐标下的的一般解（通解）形式。 以下题目需提交作业： 6. 当矢量位为 （1），； （2），； 时，分别推导由矢量位计算电磁场各直角坐标和圆柱坐标分量的关系式，并且讨论其电磁场特点。 7. 对于TEM 波（横电磁波），标量电位函数满足拉普拉斯方 程：，即在横街面上具有静电场的行为特征，这种特征给电磁场 的数值计算带来很大的方便，试证明之。 电场E和磁场H满足此关系吗？ TE波（横电波）和TM 波（横磁波）的情况如何呢？ 8. 电磁场中的标量格林函数满足亥姆霍兹方程： 对于无界空间，标量格林函数是关于源点球对称的，标量格林函数对应的亥姆霍兹方程可以变化为： 其中。其通解为：，试将通解代入上式求出。注意到一般边值问题的特解是将通解代入到边界条件（时域还需知道初始条件）中得到的，此问题的另外一个边界在无限远。能不能利用索莫菲辐射条件求出？为什么？ 下题选做： 9. 试说明准静态场的概念，并分别推导磁准静态场和电准静态场的场波动方程及其通过矢量磁位求解的过程。

几种数学计算方法的比较

有限元法，有限差分法和有限体积法的区别 有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分为中心格式和逆风格式。考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法，从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计算域内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程，即在配置点上令方程余量为0。插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘积表示，但最常用的多项式插值函数。有限元插值函数分为两大类，一类只要求插值多项式本身在插值点取已知值，称为拉格朗日(Lagrange)多项式插值；另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值，称为哈密特(Hermite)多项式插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状，一维看作长度比，二维看作面积比，三维看作体积比。在二维有限元中，三角形单元应用的最早，近来四边形等参元的应用也越来越广。对于二维三角形和四边形电源单元，常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。 对于有限元方法，其基本思路和解题步骤可归纳为 (1)建立积分方程，根据变分原理或方程余量与权函数正交化原理，建立与微分方程初边值

12 特征值估计、广义特征值与极大极小原理

第十二讲 矩阵特征值估计 特征值计算较困难，希望找到简便的特征值界限或分布范围的估计方法。 一、 特征值界的估计 定理1. 设n n A R ?∈，λ为A 的任意特征值，则有 () Im M λ≤其中，ij ji 1i ,j n a a M m a x 2 ≤≤-= 证明：设x 为A 的属于特征值λ的单位特征向量，即A x x =λ， H x x 1=， 则 H x A x λ= → ( ) () H H H H H x A x x A x x A x λ== = () ()()H H H T 2jIm x A A x x A A x λ-λ=λ=-=- 将x 写成[] T 12n x ,,,=ξξξ ()()n n H T i ij ji j i 1 j 1 x A A x a a ==-=ξ-ξ∑∑ () ()()n n i ij ji j i 1j 1 n n i ij ji j i 1 j 1 2I m a a a a ====λ= ξ-ξ≤ ξ-ξ∑∑ ∑∑ n ' i j ij ji i ,j 1 a a == ξξ-∑ （'∑表示不含i ＝j ） n ' i j i ,j 1 2M =≤ξξ∑ () 2 n 2 2 ' i j i ,j 1 I m M =? ?λ≤ξξ ? ? ? ∑

() n 2 2 ' i j i ,j 1M n n 1=≤-ξξ∑ () n 2 2 2 ' i j i ,j 1M n n 1==-ξξ∑ n n n n n 2 2 2 2 4 2 4 ' i j i j i i i i ,j 1 i ,j 1 i 1 i 1 i 1 =====ξξ= ξξ- ξ≤ ξ- ξ∑ ∑ ∑ ∑ ∑ ( )n 2 2 i i i 11== ξ-ξ∑ 不妨写为： ( ) ( ) ( )n 2 222 2 2 1 1 2 2 i i i 3 111==ξ-ξ +ξ -ξ + ξ -ξ∑ ( )( )( )2 2 2 2 2 2 n 11 22 2 2 i i i 3 1112 2 =????ξ +-ξξ +-ξ ? ? ≤++ ξ-ξ ? ? ? ???? ? ∑ 12 ≤ 取等号的条件为2 2 1 2 12 ξ=ξ= ，但 2 x 1 =，所以其它2 i ξ= ∴ () Im M λ≤定理2. 设n n A R ?∈，λ为A 的任意特征值，则有 n λ≤ρ ()R e n λ≤τ () I m n s λ≤ 其中，ij 1i,j n m a x a ≤≤ρ =，ij ji 1i,j n m a x a a ≤≤τ =+，ij ji 1i,j n s m a x a a ≤≤=- 二、 盖尔圆法 定义：设() n n ij n n A a C ??= ∈,由方程 n ii i ij j 1 i j z a R a =≠-≤= ∑ 所确定的圆称 为A 的第i 个盖尔圆，i R 称为盖尔圆的半径。

电磁学主要公式、定理、定律

电磁学主要公式、定理、定律 一. 电场 1.库仑定律：212 q q F K r = 2.电场强度定义式：F E q = 3.点电荷电场强度决定式：2 Q E K r = 4.电势定义式：P E q ?= 5.两点间电势差：AB A B U ??=- 6.场强与电势差的关系式：AB U Ed = （只适用于匀强电场） 7.电场力移动电荷做功：AB W U q =? 8平行板电容器电容定义式：Q C U = (U 就是电势差AB U ) 9.平行板电容器电容决定式：4S C Kd επ= （ 式中，ε为介质的介电常数，S 为两板正对面积， K 为静电力恒量，d 为板间距离） 10.带电粒子在匀强电场中被加速：21 2mv qU = 11.带电粒子在匀强电场中偏转：2 2 02qL U y mv d = （U 为两板间电压） 二.恒定电流 1.电流强度定义式：q I t = 2.电流微观表达式：I nqSv = （其中n 为单位 体积内 的自由 电荷数，q 为每个电荷的电量值，S 为导体的横截面积，v 为 自由电荷定向移动速率。） 3.电动势定义式：W E q = （W 为非静电力移送电荷做的功，q 为被移送的电荷量） 4.导线电阻决定式：L R S ρ = ( 式中ρ为电阻率，由导线材料、温度决定，L 为导线长，S

为导线横截面积。) 5.欧姆定律：U I R = (只适用于金属导电和电解液导电的纯电阻电路，对含电动机、电解槽 的非纯电阻电路，气体导电和半导体导电不适用） 6.串联电路： （1） 总电阻 12......R R R =++总 （2） 电流关系 123.....I I I I === （3） 电压关系 123......U U U U =++总 7.并联电路： （1）总电阻 123 1111 ......R R R R =+++总 ①只有两个电阻并联时用 12 12 R R R R R = +总 更方便快捷； ②若是n 个相同的电阻并联。可用1= R R n 总 （2） 电流关系 123=......I I I I +++总 (3) 电压关系 123=......U U U U ===总 8.电功的定义式：W qU UIt == （ 在纯电阻电路中 ，2 2 U W UIt I Rt t R ===） 9.电功率定义式：W P UI t == （ 在纯电阻电路中 , 22 U P I R R ==） 10.焦耳定律(电热计算式)：2Q I Rt = 11.电热与电功的关系 ： （1）在纯电电路中，W Q = （2）在非纯电阻电路中 W qU UIt == ＞Q 2I Rt = 12.电功率定义式：W P t = 13.电功率通用式：W P t = 和 P UI = （对纯电阻电路，22 W U P UI I R t R ====） 14.闭合电路欧姆定律：E I R r =+ （变形：E U U =+外内 ；E IR Ir =+； E U Ir =+外） 三. 磁场

第8章 特征值和特征向量

第8章特征值和特征向量 M A T L A B中的命令计算特征值和特征向量很方便，可以得到不同的子结果和分解，这在线性代数教学时很有用。注意，本章中的命令只能对二维矩阵操作。 8.1 特征值和特征向量的计算 假设A是一个m×n的矩阵，A的特征值问题就是找到方程组的解： 其中λ是一个标量，x是一个长度为n的列向量。标量λ是A的特征值，x是相对应的特征向量。对于实数矩阵A来说，特征值和特征向量可能是复数。一个n×n的矩阵有n个特征值，表，λ2，. ..，λn。 示为λ 1 M A T L A B中用命令e i g来确定矩阵A的特征值和特征向量。特征向量的规格化，就是每个特征向量的欧几里得范数为1；参见7 .6节。 命令e i g自动完成对矩阵A的平衡化。这就要求M A T L A B找出一个相似变换矩阵Q，满足 条件。求的特征值比求A的特征值条件更好些。万一A有一个和机器错误大小一样的元素，平衡化对于计算过程是没有好处的。带有参数n o b a l a n c e的命令e i g可用来计算没有这个变换矩阵的特征值和特征向量。 命令集7 9特征值和特征向量 e i g(A)求包含矩阵A的特征值的向量。 [ X,D]=e i g(A)产生一个矩阵A的特征值在对角线上的对角矩阵D和矩阵 X，它们的列是相应的特征向量，满足A X=X D。为了得到 有更好条件特征值的矩阵要进行相似变换。 [ X,D]=不经过平衡处理求得矩阵A的特征值和特征向量，也就是 e i g(A,’n o b a l a n c e’)不进行平衡相似变换。 b a l a n c e(A)求平衡矩阵。 [ T,B]=b a l a n c e(A)找到一个相似变换矩阵T和矩阵B，使得它们满足B=T－1AT。 B是用命令b a l a n c e求得的平衡矩阵。 e i g s(A)返回一个由矩阵A的部分特征值组成的向量，和命令e i g 一样，但是不返回全部的特征值。如果不带有参量，则计 算出最大的特征值。当计算所有特征值时，如果矩阵A的 秩不小于6，则计算出6个特征值来。 e i g s(f,n)求出矩阵A的部分特征值。在使用一个矩阵列的线性运算 符时，字符串f中包含的是M文件的文件名，n指定问题的 阶次。用这种方法来求特征值比开始就用运算符来求要快。

电磁学作业及解答

电磁学习题 1 (1)在没有电流的空间区域里，如果磁感应线是平行直线，磁感应强度B 的大 小在沿磁感应线和垂直它的方向上是否可能变化(即磁场是否一定是均匀的)? (2)若存在电流，上述结论是否还对? 2 如题图所示，AB 、CD 为长直导线，C B 为圆心在O 点的一段圆弧形导线， 其半径为R ．若通以电流I ，求O 点的磁感应强度． 图 3 在半径为R 的长直圆柱形导体内部，与轴线平行地挖成一半径为r 的长直圆柱形空腔，两轴间距离为a ,且a ＞r ，横截面如题9-17图所示．现在电流I 沿导体管流动，电流均匀分布在管的横截面上，而电流方向与管的轴线平行．求： (1)圆柱轴线上的磁感应强度的大小； (2)空心部分轴线上的磁感应强度的大小． 4 如图所示，长直电流1I 附近有一等腰直角三角形线框，通以电流2I ，二者 共面．求△ABC 的各边所受的磁力． 图 5 一正方形线圈，由细导线做成，边长为a ，共有N 匝，可以绕通过其相对两边中点的一个竖直轴自由转动．现在线圈中通有电流I ，并把线圈放在均匀的水平

外磁场B 中，线圈对其转轴的转动惯量为J .求线圈绕其平衡位置作微小振动时 的振动周期T . 6 电子在B =70×10-4 T 的匀强磁场中作圆周运动，圆周半径r =3.0cm ．已知B 垂直于纸面向外，某时刻电子在A 点，速度v 向上，如图． (1) 试画出这电子运动的轨道； (2) 求这电子速度v 的大小； (3)求这电子的动能k E ． 图 7 在霍耳效应实验中，一宽1.0cm ，长4.0cm ，厚1.0×10-3cm 的导体，沿长度 方向载有3.0A 的电流，当磁感应强度大小为B =1.5T 的磁场垂直地通过该导体时，产生1.0×10-5V 的横向电压．试求： (1) 载流子的漂移速度； (2) 每立方米的载流子数目． 8 如图所示，载有电流I 的长直导线附近，放一导体半圆环MeN 与长直导线共面，且端点MN 的连线与长直导线垂直．半圆环的半径为b ，环心O 与导线相距a ．设半圆环以速度v 平行导线平移．求半圆环内感应电动势的大小和方向及MN 两端的电压 N M U U ． 图 9 如图所示，用一根硬导线弯成半径为r 的一个半圆．令这半圆形导线在磁场

物理电磁学论文

物理电磁学论文 现代人的生活已经离不开电，与此同时，电磁也充斥着我们生活中的每一个角落。随着电磁学，电磁技术的发展，我们已经离不开它了，在越来越多的领域，越来越多的角落，电磁学都在发挥着它的作用。1电磁对家庭输电的影响 现在人们越来越关注周围的生活环境了，所谓的污染已经不再是我们的眼睛所能看到的垃圾，耳朵听到的噪声，鼻子闻到的恶臭，还有我们看不见，摸不着的电磁辐射。随着科学技术的发展和信息社会的到来，我们的居室内不仅有冰箱，彩色电视机，洗衣机，微波炉和空调机等家用电器，而且不少家庭中还有计算机，传真机等多种信息交流的工具，相应地，进入每个家庭的输电线强磁场对人体也特别有害处。 摘要：介绍了电磁学计算方法的研究进展和状态，对几种富有代表性的算法做了介绍，并比较了各自的优势和不足，包括矩量法、有限元法、时域有限差分方法以及复射线方法等。 关键词：矩量法；有限元法；时域有限差分方法；复射线方法 1 引言 1864年Maxwell在前人的理论（高斯定律、安培定律、法拉第定律和自由磁极不存在）和实验的基础上建立了统一的电磁场理论，并用数学模型揭示了自然界一切宏观电磁现象所遵循的普遍规律，这就是著名的Maxwell方程。在11种可分离变量坐标系求解Maxwell方程组或者其退化形式，最后得到解析解。这种方法可以得到问题的准确解，而且效率也比较高，但是适用范围太窄，只能求解具有规则边界的简单问题。对于不规则形状或者任意形状边界则需要比较高的数学技巧，甚至无法求得解析解。20世纪60年代以来，随着电子计算机技术的发展，一些电磁场的数值计算方法发展起来，并得到广泛地应用，相对于经典电磁理论而言，数值方法受边界形状的约束大为减少，可以解决各种类型的复杂问题。但各种数值计算方法都有优缺点，一个复杂的问题往往难以依靠一种单一方法解决，常需要将多种方法结合起来，互相取长补短，因此混和方法日益受到人们的重视。 2 电磁场数值方法的分类 电磁学问题的数值求解方法可分为时域和频域2大类。频域技术主要有矩量法、有限差分方法等，频域技术发展得比较早，也比较成熟。时域法主要有时域差分技术。时域法的引入是基于计算效率的考虑，某些问题在时域中讨论起来计算量要小。例如求解目标对冲激脉冲的早期响应时，频域法必须在很大的带宽内进行多次采样计算，然后做傅里叶反变换才能求得解答，计算精度受到采样点的影响。若有非线性部分随时间变化，采用时域法更加直接。另外还有一些高频方法，如GTD，UTD和射线理论。 从求解方程的形式看，可以分为积分方程法（IE）和微分方程法（DE）。IE和DE相比，有如下特点：IE法的求解区域维数比DE法少一维，误差限于求解区域的边界，故精度高；IE法适合求无限域问题，DE法此时会遇到网格截断问题；IE法产生的矩阵是满的，阶数小，DE法所产生的是稀疏矩阵，但阶数大；IE法难以处理非均匀、非线性和时变媒质问题，DE 法可直接用于这类问题〔1〕。 3 几种典型方法的介绍 有限元方法是在20世纪40年代被提出，在50年代用于飞机设计。后来这种方法得到发展并被非常广泛地应用于结构分析问题中。目前，作为广泛应用于工程和数学问题的一种通用方法，有限元法已非常著名。

电磁学计算题题库(附答案)

《电磁学》练习题（附答案） 1. 如图所示，两个点电荷＋q 和－3q ，相距为d . 试求： (1) 在它们的连线上电场强度0=E ? 的点与电荷为＋q 的点电荷相距多远？ (2) 若选无穷远处电势为零，两点电荷之间电势U =0的点与电荷为＋q 的点电荷相距多远？ d 2. 一带有电荷q ＝3×10-9 C 的粒子，位于均匀电场中，电场方向如图所示．当该粒子沿水平方向向右方运动5 cm 时，外力作功6×10-5 J ，粒子动能的增量为4.5×10-5 J ．求：(1) 粒子运动过程中电场力作功多少？(2) 该电场的场强多大？ 3. 如图所示，真空中一长为L 的均匀带电细直杆，总电荷为q ，试求在直杆延长线上距杆的一端距离为d 的P 点的电场强度． 4. 一半径为R 的带电球体，其电荷体密度分布为 ρ =Ar (r ≤R ) ， ρ =0 (r ＞R ) A 为一常量．试求球体内外的场强分布． 5. 若电荷以相同的面密度σ均匀分布在半径分别为r 1＝10 cm 和r 2＝20 cm 的两个同心球面上，设无穷远处电势为零，已知球心电势为300 V ，试求两球面的电荷面密度σ的值． (ε0＝8.85×10-12 C 2 / N ·m 2 ) 6. 真空中一立方体形的高斯面,边长a ＝0.1 m ，位于图中所示位置．已知空间的场强分布为： E x =bx , E y =0 , E z =0． 常量b ＝1000 N/(C ·m)．试求通过该高斯面的电通量． 7. 一电偶极子由电荷q ＝1.0×10-6 C 的两个异号点电荷组成，两电荷 相距l ＝2.0 cm ．把这电偶极子放在场强大小为E ＝1.0×105 N/C 的均匀电场中．试求： (1) 电场作用于电偶极子的最大力矩． (2) 电偶极子从受最大力矩的位置转到平衡位置过程中，电场力作的功． 8. 电荷为q 1＝8.0×10-6 C 和q 2＝－16.0×10-6 C 的两个点电荷相距20 cm ，求离它们都是20 cm 处的电场强度． (真空介电常量ε0＝8.85×10 -12 C 2N -1m -2 ) 9. 边长为b 的立方盒子的六个面，分别平行于xOy 、yOz 和xOz 平面．盒子的一角在坐标原点处．在此区域有 一静电场，场强为j i E ? ??300200+= .试求穿过各面的电通 量． 10. 图中虚线所示为一立方形的高斯面，已知空间的场强分布为： E x ＝bx ， E y ＝0， E z ＝0．高斯面边长a ＝0.1 m ，常量b ＝1000 N/(C ·m)．试求该闭合面中包含的净电荷．(真空介电常数ε0＝8.85×10-12 C 2 ·N -1 ·m -2 ) 11. 有一电荷面密度为σ的“无限大”均匀带电平面．若以该平面处为电势零点，试求带电平面周围空间的电势分 布． 12. 如图所示，在电矩为p ? 的电偶极子的电场中，将一电荷为q 的点电荷从A 点沿半径为R 的圆弧(圆心与电偶极子中心重合，R >>电偶极子正负电荷之 间距离)移到B 点，求此过程中电场力所作的功． 13. 一均匀电场，场强大小为E ＝5×104 N/C ，方向竖直朝上，把一电荷为q ＝ 2.5×10-8 C 的点电荷，置于此电场中的a 点，如图所示．求此点电荷在下列过程中电场力作的功． (1) 沿半圆路径Ⅰ移到右方同高度的b 点，ab ＝45 cm ； (2) 沿直线路径Ⅱ向下移到c 点，ac ＝80 cm ； (3) 沿曲线路径Ⅲ朝右斜上方向移到d 点，ad ＝260 cm(与水平方向成45°角)． 14. 两个点电荷分别为q 1＝＋2×10-7 C 和q 2＝－2×10-7 C ，相距0.3 m ．求距q 1为0.4 m 、距q 2为0.5 m 处P 点的电场强度． ( 41επ=9.00×109 Nm 2 /C 2 ) 15. 图中所示， A 、B 为真空中两个平行的“无限大”均匀带电平面，A 面上电荷面密度σA ＝－17.7×10-8 C ·m -2 ，B 面的电荷面密度σB ＝35.4 ×10-8 C ·m -2 ．试计算两平面之间和两平面外的电场强度．(真空介电常量ε0＝8.85×10-12 C 2 ·N -1 ·m -2 ) 16. 一段半径为a 的细圆弧，对圆心的张角为θ0，其上均匀分布有正电荷q ，如图所示．试以a ，q ，θ0表示出圆心O 处的电场强度． 17. 电荷线密度为λ的“无限长”均匀带电细线，弯成图示形状．若半圆弧AB R ，试求圆心O 点的场强． E ? q L d q O x z y a a a a A B R ? Ⅰ Ⅱ Ⅲ d b a 45?c E ? σA σB A B O a θ0 q A R ∞ ∞ O

各种计算电磁学方法比较和仿真软件

各种计算电磁学方法比较和仿真软件 各种计算电磁学方法比较和仿真软件微波EDA 仿真软件与电磁场的数值算法密切相关，在介绍微波EDA 软件之前先简要的介绍一下微波电磁场理论的数值算法。所有的数值算法都是建立在Maxwell 方程组之上的，了解Maxwell 方程是学习电磁场数值算法的基础。计算电磁学中有众多不同的算法，如时域有限差分法（FDTD ）、时域有限积分法（FITD ）、有限元法（FE）、矩量法（MoM ）、边界元法（BEM ）、谱域法（SM）、传输线法（TLM ）、模式匹配法（MM ）、横向谐振法（TRM ）、线方法（ML ）和解析法等等。在频域，数值算法有：有限元法（ FEM -- Finite Element Method）、矩量法（MoM -- Method of Moments ），差分法（ FDM -- Finite Difference Methods ），边界元法（ BEM --Boundary Element Method ），和传输线法 （ TLM -Transmission-Line-matrix Method ）。在时域，数值算法有：时域有限差分法（ FDTD - Finite Difference Time Domain ），和有限积分法（ FIT - Finite Integration Technology ）。这些方法中有解析法、半解析法和数值方法。数值方法中又分零阶、一阶、二阶和高阶方法。依照解析程度由低到高排列，依次是：时域有限差分法（FDTD ）、传输线法（TLM ）、时域有限积分法（FITD ）、有限元法（FEM ）、矩量法（MoM ）、线方法（ML ）、边界元法（BEM ）、谱域法（SM ）、模式匹配法

eig求所有特征值和特征向量

最近看了看matlab求特征值的函数，记下来备用。 eig求所有特征值和特征向量。 d = eigs(A) %求稀疏矩阵A的6个绝对值最大特征值d，d以向量形式存放。 d = eigs(A,B) %求稀疏矩阵的广义特征值问题。满足A V=BVD，其中D为特征值对角阵，V为特征向量矩阵，B必须是对称正定阵或Hermitian正定阵。 d = eigs(A,k) %返回k个最大特征值 d = eigs(A,B,k) %返回k个最大特征值 d = eigs(A,k,sigma) %sigma取值：'lm'表示绝对值最大的特征值；'sm'绝对值最小特征值；对实对称问题：'la'表示最大特征值；'sa'为最小特征值；对非对称和复数问题：'lr'表示最大实部；'sr'表示最小实部；'li'表示最大虚部；'si'表示最小虚部 d = eigs(A,B,k,sigma) %同上 d = eigs(A,k,sigma,opts) % opts为指定参数：参见eigs帮助文件。opts为一个向量 参数描述value opts.issym =1:如果A对称 =0：A不对称 {0|1} opts.isreal =1:A为实数 =0：otherwise {0|1} opts.tol 收敛？？？（没看懂）**估计

d = eigs(A,B,k,sigma,options) %同上。以下的参数k、sigma、options相同。 d = eigs(Afun,n) %用函数Afun代替A，n为A的阶数，D为特征值。 d = eigs(Afun,n,B) d = eigs(Afun,n,k) d = eigs(Afun,n,B,k) d = eigs(Afun,n,k,sigma) d = eigs(Afun,n,B,k,sigma) d = eigs(Afun,n,k,sigma,options) d = eigs(Afun,n,B,k,sigma,options) [V,D] = eigs(A,…) %D为6个最大特征值对角阵，V的列向量为对应特征向量。 [V,D] = eigs(Afun,n,…) [V,D,flag] = eigs(A,…) %flag表示特征值的收敛性，若flag=0，则所有特征值都收敛，否则，不是所有都收敛。 [V,D,flag] = eigs(Afun,n,…)

电磁学复习计算题(附答案)

《电磁学》计算题（附答案） 1. 如图所示，两个点电荷＋q 和－3q ，相距为d . 试求： (1) 在它们的连线上电场强度0=E ? 的点与电荷为＋q 的点电荷相距多远？ (2) 若选无穷远处电势为零，两点电荷之间电势U =0的点与电荷为＋q 的点电荷相距多远？ d +q 2. 一带有电荷q ＝3×10-9 C 的粒子，位于均匀电场中，电场方向如图所示．当该粒子沿水平方向向右方运动5 cm 时，外力作功6×10-5 J ，粒子动能的增量为4.5×10-5 J ．求：(1) 粒子运动过程中电场力作功多少？(2) 该电场的场强多大？ 3. 如图所示，真空中一长为L 的均匀带电细直杆，总电荷为q ，试求在直杆延长线上距杆的一端距离为d 的P 点的电场强度． 4. 一半径为R 的带电球体，其电荷体密度分布为 =Ar (r ≤R ) ， =0 (r ＞R ) A 为一常量．试求球体内外的场强分布． 5. 若电荷以相同的面密度均匀分布在半径分别为r 1＝10 cm 和r 2＝20 cm 的两个同心球面上，设无穷远处电势为零，已知球心电势为300 V ，试求两球面的电荷面密度的值． (0 ＝8.85× 10-12C 2 / N ·m 2 ) 6. 真空中一立方体形的高斯面,边长a ＝0.1 m ，位于图中所示位 置．已知空间的场强分布为： E x =bx , E y =0 , E z =0． 常量b ＝1000 N/(C ·m)．试求通过该高斯面的电通量． 7. 一电偶极子由电荷q ＝1.0×10-6 C 的两个异号点电荷组成，两电荷相距l ＝2.0 cm ．把这电偶极子放在场强大小为E ＝1.0×105 N/C 的均匀电场中．试求： (1) 电场作用于电偶极子的最大力矩． (2) 电偶极子从受最大力矩的位置转到平衡位置过程中，电场力作的功． 8. 电荷为q 1＝8.0×10-6 C 和q 2＝－16.0×10-6 C 的两个点电荷相距20 cm ，求离它们都是20 cm 处的电场强度． (真空介电常量 ＝8.85×10 -12 C 2N -1m -2 ) 9. 边长为b 的立方盒子的六个面，分别平行于xOy 、yOz 和xOz 平面．盒子的一角在坐标原点处．在 此区域有一静电场，场强为j i E ? ??300200+= .试求穿过各面的电通量． E ? q L d q O x z y a a a a

计算电磁学入门基础介绍

计算电磁学入门基础介绍 一. 计算电磁学的重要性 在现代科学研究中，“科学试验，理论分析，高性能计算”已经成为三种重要的研究手段。在电磁学领域中，经典电磁理论只能在11 种可分离变量坐标系中求解麦克斯韦方程组或者其退化形式，最后得到解析解。解析解的优点在于： ①可将解答表示为己知函数的显式，从而可计算出精确的数值结果; ②可以作为近似解和数值解的检验标准; ③在解析过程中和在解的显式中可以观察到问题的内在联系和各个参数对数值结果所起的作用。 这种方法可以得到问题的准确解，而且效率也比较高，但是适用范围太窄，只能求解具有规则边界的简单问题。当遇到不规则形状或者任意形状边界问题时，则需要比较复杂的数学技巧，甚至无法求得解析解。20 世纪60 年代以来，随着电子计算机技术的发展，一些电磁场的数值计算方法也迅速发展起来，并在实际工程问题中得到了广泛地应用，形成了计算电磁学研究领域，已经成为现代电磁理论研究的主流。简而言之，计算电磁学是在电磁场与微波技术学科中发展起来的，建立在电磁场理论基础上，以高性能计算机技术为工具，运用计算数学方法，专门解决复杂电磁场与微波工程问题的应用科学。相对于经典电磁理论分析而言，应用计算电磁学来解决电磁学问题时受边界约束大为减少，可以解决各种类型的复杂问题。原则上来讲，从直流到光的宽广频率范围都属于该学科的研究范围。近几年来,电磁场工程在以电磁能量或信息的传输、转换过程为核心的强电与弱电领域中显示了重要作用。 二. 电磁问题的分析过程 电磁工程问题分析时所经历的一般过程为： 三. 计算电磁学的分类 (1) 时域方法与谱域方法 电磁学的数值计算方法可以分为时域方法(Time Domain或TD)和频域方法(Frequeney Domain或FD)两大类。 时域方法对Maxwell方程按时间步进后求解有关场量。最著名的时域方法是时域有限差分法(Finite Difference Time Domain或FDTD)。这种方法通常适用于求解在外界激励下场

电磁学第二版答案(DOC)

第一章静电场 §1.1 静电的基本现象和基本规律 思考题： 1、给你两个金属球，装在可以搬动的绝缘支架上，试指出使这两个球带等量异号电荷的方向。你可以用丝绸摩擦过的玻璃棒，但不使它和两球接触。你所用的方法是否要求两球大小相等？ 答：先使两球接地使它们不带电，再绝缘后让两球接触，将用丝绸摩擦后带正电的玻璃棒靠近金属球一侧时，由于静电感应，靠近玻璃棒的球感应负电荷，较远的球感应等量的正电荷。然后两球分开，再移去玻璃棒，两金属球分别带等量异号电荷。本方法不要求两球大小相等。因为它们本来不带电，根据电荷守恒定律，由于静电感应而带电时，无论两球大小是否相等，其总电荷仍应为零，故所带电量必定等量异号。 2、带电棒吸引干燥软木屑，木屑接触到棒以后，往往又剧烈地跳离此棒。试解释之。答：在带电棒的非均匀电场中，木屑中的电偶极子极化出现束缚电荷，故受带电棒吸引。但接触棒后往往带上同种电荷而相互排斥。 3、用手握铜棒与丝绸摩擦，铜棒不能带电。戴上橡皮手套，握着铜棒和丝绸摩擦，铜棒就会带电。为什么两种情况有不同结果？ 答：人体是导体。当手直接握铜棒时，摩擦过程中产生的电荷通过人体流入大地，不能保持电荷。戴上橡皮手套，铜棒与人手绝缘，电荷不会流走，所以铜棒带电。 7、两个点电荷带电2q 和q，相距l，第三个点电荷放在何处所受的合力为零？ 解：设所放的点电荷电量为Q。若Q与q同号，则三者互相排斥，不可能达到平衡；故Q 只能与q异号。当Q在2q和q联线之外的任何地方，也不可能达到平衡。由此可知，只有Q与q异号，且处于两点荷之间的联线上，才有可能达到平衡。设Q到q的距离为x. 8、三个相同的点电荷放置在等边三角形的各顶点上。在此三角形的中心应放置怎样的电荷，才能使作用在每一点电荷上的合力为零？ 解：设所放电荷为Q，Q应与顶点上电荷q异号。中心Q所受合力总是为零，只需考虑q 受力平衡。 平衡与三角形边长无关，是不稳定平衡。 9、电量都是Q的两个点电荷相距为l，联线中点为O；有另一点电荷q，在联线的中垂面上距O为r处。（1）求q所受的力；（2）若q开始时是静止的，然后让它自己运动，它将如何运动？分别就q与Q同号和异号两种情况加以讨论。 解： (1) (2)q与Q同号时，F背离O点，q将沿两Q的中垂线加速地趋向无穷远处。 q与Q异号时，F指向O点，q将以O为中心作周期性振动，振幅为r . <讨论>：设q 是质量为m的粒子，粒子的加速度为 因此，在r<

电磁学练习题积累(含部分答案)

一.选择题（本大题15小题，每题2分） 第一章、第二章 1.在静电场中，下列说法中哪一个是正确的 [ ] (A)带正电荷的导体，其电位一定是正值 (B)等位面上各点的场强一定相等 (C)场强为零处，电位也一定为零 (D)场强相等处，电位梯度矢量一定相等 2.在真空中的静电场中，作一封闭的曲面，则下列结论中正确的是［］ (A)通过封闭曲面的电通量仅是面内电荷提供的 (B) 封闭曲面上各点的场强是面内电荷激发的 (C) 应用高斯定理求得的场强仅是由面内电荷所激发的 (D) 应用高斯定理求得的场强仅是由面外电荷所激发的 3.关于静电场下列说法中正确的是 [ ] (A)电场和试探电荷同时存在和消失 (B)由E＝F/q知道，电场强度与试探电荷成反比 (C)电场强度的存在与试探电荷无关 (D)电场是试探电荷和场源电荷共同产生的 4.下列几个说法中正确的是： [ ] (A)电场中某点场强的方向，就是将点电荷放在该点所受电场力的方向 (B)在以点电荷为中心的球面上，由该点电荷所产生的场强处处相同 (C)场强方向可由E=F/q定出，其中q为试验电荷的电量，q可正、可负， F为试验电荷所受的电场力 (D)以上说法全不对。 5.一平行板电容器中充满相对介电常数为的各向同性均匀电介质。已知介 质两表面上极化电荷面密度为，则极化电荷在电容器中产生的电 场强度的大小为 [ ]

(A) 0εσ' (B) 02εσ' (C) 0εεσ' (D) ε σ' 6. 在平板电容器中充满各向同性的均匀电介质，当电容器充电后，介质中 D 、 E 、P 三矢量的方向将是 [ ] (A) D 与E 方向一致，与P 方向相反 (B) D 与E 方向相反，与P 方向一致 (C) D 、E 、P 三者方向相同 (D) E 与P 方向一致，与D 方向相反 7. 在一不带电荷的导体球壳的球心处放一点电荷，并测量球壳内外的场强分 布，如果将此点电荷从球心移到球壳内其它位置，重新测量球壳内外的场强分布，则将发现： [ ] (A) 球壳内、外场强分布均无变化 (B) 球壳内场强分布改变，球壳外的不变 (C) 球壳外场强分布改变，球壳内的不变 (D) 球壳内、外场强分布均改变 8. 一电场强度为E 的均匀电场，E 的方向与x 轴正向平行，如图所示，则通过 图中一半径为R 的半球面的电场强度通量为 [ ] (A) 2R E π；(B) 21 2 R E π； (C) 22R E π；(D ) 0。 9. 在静电场中，电力线为均匀分布的平行 直线的区域内，在电力线方向上任意两点的电场强度E 和电势U 相比较 [ ] (A) E 相同，U 不同 (B) E 不同，U 相同 (C) E 不同，U 不同 (D) E 相同，U 相同

计算电磁学结课论文

《计算电磁学》学习心得 姓名：桑dog 学号： 班级： 联系方式：

前言 计算电磁学是科技的重要领域它的研究涉及到应用计算机求解电磁方程它的重要性基于麦克斯韦方程——唯一的可以描述小到亚原子大到天体尺度的所有物理现象的方程, 。而且, 麦克斯韦方程式对于结果拥有很强的预测能力: 对于一个复杂问题的麦克斯韦方程的解通常可以准确的预知实验结果。因此, 麦克斯韦方程的解对于提高我们对复杂系统之物理现象的洞察力和设计复杂系统的能力均有极大帮助所以, 成功求解麦克斯韦方程式拥有广泛的应用前景: 例如纳米技术, 电脑微电子电路, 电脑芯片设计, 光学, 纳米光学, 微波工程, 遥感, 射电天文学, 生物医学工程, 逆散射和成象等等。 这篇文章的安排如下：第一章介绍了计算电磁学的重要意义以及发展状况。第二章介绍了计算电磁学中解决问题的方法分类。第三章对主要的数值方法进行了简介。第四章展望了计算电磁学的发展趋势。

第1章计算电磁学的重要性 在现代科学研究中，“科学试验，理论分析，高性能计算”已经成为三种重要的研究手段[1]。在电磁学领域中，经典电磁理论只能在11 种可分离变量坐标系中求解麦克斯韦方程组或者其退化形式，最后得到解析解。解析解的优点在于： ●可将解答表示为己知函数的显式，从而可计算出精确的数值结果; ●可以作为近似解和数值解的检验标准; ●在解析过程中和在解的显式中可以观察到问题的内在联系和各个参数对数值 结果所起的作用。 这种方法可以得到问题的准确解，而且效率也比较高，但是适用范围太窄，只能求解具有规则边界的简单问题[2]。当遇到不规则形状或者任意形状边界问题时，则需要比较复杂的数学技巧，甚至无法求得解析解。20 世纪60 年代以来，随着电子计算机技术的发展，一些电磁场的数值计算方法也迅速发展起来，并在实际工程问题中得到了广泛地应用，形成了计算电磁学研究领域，已经成为现代电磁理论研究的主流。简而言之，计算电磁学是在电磁场与微波技术学科中发展起来的，建立在电磁场理论基础上，以高性能计算机技术为工具，运用计算数学方法，专门解决复杂电磁场与微波工程问题的应用科学。相对于经典电磁理论分析而言，应用计算电磁学来解决电磁学问题时受边界约束大为减少，可以解决各种类型的复杂问题。原则上来讲，从直流到光的宽广频率范围都属于该学科的研究范围。近几年来,电磁场工程在以电磁能量或信息的传输、转换过程为核心的强电与弱电领域中显示了重要作用。[3]