第七章 静电场与积分计算.pptVIP

静电场 线电荷电场 取: 则其等势线为: 取: 7.1 矩形梯形和抛物线形积分近似计算 定积分: 数值积分的原因: (1) f(x) 形式很复杂难于求其原函数F(x); 如: (2) f(x) 的原函数不能由初等函数表达; 如: (3) f(x) 没有给出表达式,只给出数据值. 数值积分原理: 把区间等分,求所有小条形面积的和,当?x趋于0时就是函数f(x)在整个区间上的积分值. 矩形公式: 对于积分: 把区间[a,b]等分成n-1份,等分点记为x1, x2, …, xn 则等分间隔为: 矩形公式就是用矩形面积近似小条形的面积 每个小矩形面积: 矩形法计算公式: 例: 求sin(x) 在[0, ?]的积分. function v = jx1(n) % 矩形法求sin(x)在[0,pi]上的积分 dx = pi/(n-1); x = 0:dx:pi; v = 0; for i = 1:n-1 v = v + sin(x(i))*dx; end function v = jx2(n) % 矩形法求sin(x)在[0,pi]上的积分 dx = pi/(n-1); x = 0:dx:pi; y = sin(x); v = sum(y(1:n-1))*dx; jx1(100) ans = 1.9998 jx2(100) ans = 1.8138 梯形公式: 梯形公式就是用梯形面积近似小条形的面积 每个小梯形面积: 梯形法计算公式: 或: 即: function v = tx(n) % 梯形法求sin(x)在[0,pi]上的积分 dx = pi/(n-1); x = 0:dx:pi; y = sin(x); v = (0.5*(y(1)+y(n))+sum(y(2:n-1)))*dx; 抛物线形公式(辛卜生): 抛物线公式就是用抛物线代替f(x)进行积分 梯形公式相当于用两点(xi-1,f(xi-1)), (xi,f(xi))上的线形插值函数代替f(x)进行积分. 抛物线公式就是用三点(xi-1,f(xi-1)), (xi,f(xi)), (xi+1,f(xi+1)) 二次插值函数代替f(x)进行积分. 二次插值公式: 即: 取 ，在[a,b]上用 ，得 为了便于计算，做变量代替 于是 从而有： 那么: 或 要求n为奇数, 即有(n-1)偶数个小条形. 梯形法和抛物线法都可以写成统一格式: 称为牛顿-柯特斯公式, 称为牛顿-柯特斯系数,与 f(x)无关为常数. h 为区间范围. 梯形公式(n=1): 抛物线形公式(n=2): 例: 求sin(x) 在[0, ?]的积分. function v = pwx(n) % 抛物线形法求sin(x)在[0,pi]上的积分 if mod(n,2)~=1 % 判断n是否是奇数. mod()求余数函数. error(n should be odd); end dx = pi/(n-1); x = 0:dx:pi; y = sin(x); v = (y(1)+y(n)+4*sum(y(2:2:n-1))+… 2*sum(y(3:2:n-2)))*dx/3; 7.2变步长辛卜生积分 各种算法的积分运算，要想得到与实际积分足够接近的数值结果，都必须取足够小的步长，也即足够多的区间分割。 变步长辛卜生积分就是“逐次加倍等分区间的个数，直至误差符合精度要求”的算法： 第1次, 把积分区间 2 等分, n=21+1=3, 结果记 S1 第2次, 把积分区间 4 等分, n=22+1=5, 结果记 S2 … 第i次, 把积分区间 2i 等分, n=2i+1, 结果记 Si 第i+1次, 把积分区间 2i+1 等分, n=2i+1+1, 结果 Si+1 ... 这样得到的序列 {S1, S2, S3,…, Si, Si+1,…} 肯定是收敛于实际积分值的。然而具体计算到哪一步结束，要看所需结果的精度而定。 预先设置误差精度ε, 令 如 |D| ?, 则停止, 否则 i+=1, 继续。 稍加分析，我们会发现， Si+1 与 Si 的计