工程电磁场数值方法编程实验3－数值积分方法.pptVIP

* * * * 电压源激励的磁场 螺线管电阻 cem@uestc.edu.cn ρ电阻率 电压源激励的磁场 螺线管磁场 cem@uestc.edu.cn 编程实践四 编写轴对称线圈的矢量位计算计算函数 cem@uestc.edu.cn * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 第3章 数值积分法 工程电磁场数值方法编程实验－数值积分法 电子科技大学 赖生建 cem@uestc.edu.cn 主要内容 概述 数值积分 基于场量积分公式的数值积分法 基于场源离散化的数值积分法 编程实践 cem@uestc.edu.cn 3.1 概述 数值积分法是数值计算方法应用中的基本内容之一。 在电磁场分析计算中，对于无限大、均匀各向同性的媒质，当已知场源分布求场分布时，基于库仑定律或比奥－沙伐定律均可以导出关于场量、位函数的积分表达式。可以计算场分布以及有关电磁参数、能量和力等积分量。 它是多种电磁场数值计算方法(等参数有限元法、边界元法和模拟电荷法)数值解的必要基础。 cem@uestc.edu.cn 数值积分实质 数值积分实质上是一种近似的求积方法，即通过构造被积函数的某种线性组合的逼近函数来近似求其积分值。 如函数f(x)的定积分 ，当被积函数f(x)的原函数F(x)无法用初等函数表达，则可以用另一个具有足够逼近精度的简单函数g(x)来近似替代函数f(x)，函数g(x)的构造通常取f(x)的代数插值函数或样条函数。 cem@uestc.edu.cn 数值积分实质 函数g(x)不同的构造函数，产生不同的数值求积公式。 为提高数值积分精度，常将积分空间分成n等份，在每个小区间上采用相应的求积公式计算，称为复合的数值求积公式。 cem@uestc.edu.cn 3.2 梯形求积公式 梯形求积是一种直观的近似方法，由边界上近似直线围成的梯形面积近似代替原来曲边梯形面积。 上面误差较大，为得到求积的精度，可将积分区间细分k次，得到分段数 ，其梯形积分公式为 cem@uestc.edu.cn 梯形求积公式 如果精度不够，再各子区间平分，得2n=2k+1 当满足 cem@uestc.edu.cn 辛普生求积公式 如用二次插值多项式－－抛物线g(x)所围成的曲边梯形面积近似代替y=f(x)围成的曲边梯形面积 为得到求积的精度，将积分区间细分k次，得到分段数 ，其复合辛普生求积公式为 cem@uestc.edu.cn 辛普生求积公式 计算二重积分时，数值积分的处理是将二重积分分解为两个单积分，每个积分使用辛普生求积公式，即在第一重积分内采用辛普生求积公式，公式中每产生一个固定某变量值x，在另一重积分也用辛普生求积公式计算。 分解为两个单积分 cem@uestc.edu.cn 高斯求积公式 高斯求积法，在积分区间，选择某些积分点（设为n个积分点），计算出函数在这些积分点上的数值，然后用相应的权系数乘以这些函数值，并求和。 在相同数量积分点的选取条件下，高斯积分的精度最佳。 求积节点 为高斯积分点。 称为积分区间[a,b]上的权函数，不同的权函数选取，不同的高斯型求积公式，其中 是最常用的权函数。 cem@uestc.edu.cn 高斯求积公式 当给定权函数 ，可假定积分区间为[-1,1],利用正交多项式来确定求积点（高斯积分点）时，该正交多项式为勒让德多项式，这样积分点为 n次勒让德多项式的n个零点，由此得到相应的求积公式的权系数 上式为高斯－勒让德求积公式。高斯积分点与权系数的值参考表。 cem@uestc.edu.cn 高斯求积点 cem@uestc.edu.cn 参考高斯积分表 高斯求积公式 对于一般区间[a,b]，利用积分变量代换，转换为[-1,1]积分区间。 其中 cem@uestc.edu.cn 高斯求积公式 对于多重积分，可以化重积分为多次积分的方法，与前面辛普生的多重积分相似，每一重积分采用相同或不同的高斯积分点 cem@uestc.edu.cn 编程实践一 编写辛普生积分计算函数 cem@uestc.edu.cn 编程实践二 编写高斯积分计算函数 cem@uestc.edu.cn 椭圆积分的数值计算 在电磁场的数值计算中，椭圆积分的应用广泛。应用较多的是第一与第二类椭圆积分 第一类完全椭圆积分 第二类完全椭圆积分 称为积分模数 cem@uestc.edu.cn 椭圆积分的数值计算 级数展开