第四章数值积分out.docVIP

第四章 数值积分 主要涉及的数值计算。 由Leibniz公式可理论计算。但原函数有时很难求，或可根本不存在。 例如： 有时被“积函数”就不是一个函数，很复杂的计算结果。 §1 机械求积分公式 求的数值 左矩形公式 右矩形公式 梯形公式 以上公式都是以直代曲，一般式可写成： 理解为去[a,b]之间某点的值 或者 理解为按权重取 经常用到的还有simpson公式， 可理解为按照权重取，。 也可看作过 三个点作抛物线： 核心是进行插值用插值函数代替原积分函数（插值求积分）。 具体积分时将积分区间分成若干小区间，在每个小区间中应用simpson公式。 程序: F(x)=……. 定义函数 WRITE(*,*) ‘INPUT A,B ,N’ READ(*,*) A, B, N H=(B-A)/N SUM=0.0 DO 100 I=0,N-1 A1=A+H*I SUM=SUM+H/6*(F(A1)+4*F(A1+H/2)+F(A1+H)) 100 CONTINUE WRITE(*,*)SUM STOP END §2 变步长积分公式 梯形公式辛普森公式精度得到改善，两者都是定步长公式，在实际计算中，步长的选择成为一件困难的事。本节变步长积分。 变步长梯形积分 在区间[a,b]上使用梯形公式计算， 首先：步长取 得： 第一次，将区间[a，b]二等分，即取步长为 得： 第二次，将区间[a，b]四等分，即取步长为 规律： 第k-1次，将区间[a，b]n 等分。，步长为：，点为： 做到第k次，2n等分。步长为：， ，点为 规律：双数的点是新加的点（红色）。 我们可以比较和，如果差值很小认为精度已够。再减小步长已没多大意义。 程序： F(x)=……. 定义函数 WRITE(*,*) ‘INPUT A,B ,N’ READ(*,*) A, B, N H=(B-A)/N T0=(B-A)/2.0*(F(A)+F(B)) DO 100 K=1,N TT=0.0 DO 200 N=1,2**(K-1) 200 TT=TT+F(A+(2.0*N-1.0)/2.0**K*(B-A)) T1=T0/2+TT*(B-A)/2**K IF(ABS((T1-T0)/t1).LT.0.00001) GOTO 300 T0=t1 100 WRITE(*,*)K,T1 STOP 300 WRITE(*,*) ‘THE RESULT IS’,T1 STOP END C==================================================================C C 多重积分。计算积分区域为圆心在（Y0，Y0）半径为R的圆。 C 如所解问题不同，请改红色部分。分别代表积分函数及积分区域。 C 2009.5.21. C================================================================== DOUBLE PRECISION F,X,Y,X0,Y0,HX,HY,SUM,R F(x,Y)=1 WRITE(*,*) INPUT x0,y0,R,N READ(*,*) X0,Y0,R,N HX=2.D0*R/N HY=2.D0*R/N Y1=Y0-R Y2=Y0+R SUM=0.D0 C WRITE(*,*)HX,HY,SUM DO Y=Y1,Y2,HY X1=X0-SQRT(R*R-(Y-Y0)*(Y-Y0)) X2=X0+SQRT(R*R-(Y-Y0)*(Y-Y0)) DO X=X1,X2,HX SUM=SUM+HX*HY* ! (F(X,Y)+4.D0*F(X+HX/2.,Y+HY/2.)+F(X+HX,Y+HY))/6.0 C WRITE(*,*)SUM ENDDO ENDDO WRITE(*,*)SUM END 数值微分 我们讨论的导数问题。实际计算及讨论函数在具体某点的导数值。 我们自