第四章数值积分..docVIP

第四章 数值积分 主要涉及的数值计算。 由Leibniz公式论计算。 有时被“积函数”就不是一个函数，很复杂的计算结果。 §1 机械求积分公式 求的数值左矩形公式右矩形公式 以上公式都是以直代曲，一般式可写成： 理解为去[a,b]之间某点的值 或者经常用到的还有simpson公式， 可理解为按照权重取，。 三个点作抛物线：核心是进行插值用插值函数代替原积分函数（插值求积分）。 具体积分时将积分区间分成若干小区间，在每个小区间中应用simpson公式。 程序: Fx…….WRITE*,* ‘INPUT A,B ,N’ READ*,* A, B, N HB-A/N SUM0.0 DO 100 I0,N-1 A1A+H*I SUMSUM+H/6*FA1+4*FA1+H/2+FA1+H 100CONTINUEWRITE*,*SUMSTOP END §2 变步长积分公式梯形公式辛普森公式精度得到改善，两者都是定步长公式，在实际计算中，步长的选择成为一件困难的事。本节变步长积分。 变步长梯形积分 在区间[a,b]上使用梯形公式计算， 得： 得：规律： 第k-1次，将区间[a，b]n 等分。，， 做到第k次，2n等分。步长为：， ，我们可以比较和，如果差值很小认为精度已够。再减小步长已没多大意义。 程序： Fx……. 定义函数 WRITE*,* ‘INPUT A,B ,N’ READ*,* A, B, N HB-A/N T0B-A/2.0*FA+FB DO 100 K1,NTT0.0DO 200 N1,2**K-1 200TTTT+FA+2.0*N-1.0/2.0**K*B-AT1T0/2+TT*B-A/2**KIFABST1-T0/t1.LT.0.00001 GOTO 300 T0t1 100 WRITE*,*K,T1STOP 300 WRITE*,* ‘THE RESULT IS’,T1 STOP END 数值微分 我们讨论的导数问题。实际计算及讨论函数在具体某点的导数值。 我们自然想到微积分中的导数定义： 用作为的近似。 理论上这种近似，h越小越好，但实际计算相近的数相减会损失有效数字，实际计算中要注意这些问题。 二阶导数的三点公式：精度高一个量级 往往利用插值方法先求出插值函数，在对函数求导，用之代替原函数的导数。 数值积分我们也曾这样做。