fortran常用算法的设计举例导论.ppt

* * 第八章 常用算法的程序设计举例 第一章 算法 第二章 计算机和计算机程序 第四章 逻辑运算和选择结构 第五章 循环结构的实现 第六章 Fortran的数据结构 第七章 数据的输入、输出 第三章 Fortran语言程序设计初步 一、数值积分 几何意义： 近似求小曲边梯形面积的方法： （1）用小矩形代替小曲边梯形； （2）用小梯形代替小曲边梯形； （3）在小区间范围内，用一条抛物线代替该区间内的f(x)。 read(*,*) a,b,n x=a h=(b-a)/n f0=exp(x) s=0.0 do 10 i=1,n si=f0*h s=s+si x=x+h f0=exp(x) 10 continue write(*,100) a,b,n write(*,200) s 100 format(1x,a=,f10.3,3x,b=, $ f10.3,3x,n=,i4) 200 format(1x,s=,f15.8) end 1. 矩形法 read(*,*) a,b,n x=a h=(b-a)/n s=0.0 do 10 i=1,n si=(sin(x+(i-1)*h)+ $ sin(x+i*h))*h/2.0 s=s+si 10 continue write(*,100) a,b,n write(*,200) s 100 format(1x,a=,f10.3,3x, $ b=,f10.3,3x,n=,i4) 200 format(1x,s=,f15.8) end 2. 梯形法 其他几种 程序变形 . . . f1=sin(a) . . . do 10 i=1,n f2=sin(a+i*h) si=(f1+f2)*h/2.0 s=s+si f1=f2 10 continue . . . . . . x2=a . . . do 10 i=1,n x1=x2 x2=x2+h si=(sin(x1)+ $ sin(x2))*h/2.0 s=s+si 10 continue . . . . . . f0=sin(a) h=(b-a)/n s=f0 do 10 i=1,n f=sin(a+i*h) s=s+2.0*f 10 continue s=(s-sin(b))*h/2.0 . . . 3. Sinpson法 取a，b中点c—((a+b)/2，0)，通过f(a)，f(b)，f(c)三点可作唯一一条抛物线f1(x)。 根据抛物线定积分求值公式，有： 如果将(a，b)分成两个小区间(a，c) 和(c，b): 如果将(a，b)分成四个小区间: 如果将(a，b)分成n个小区间: read(*,*) a,b,n h=(b-a)/(2.0*n) s=0.0 fa=1.0/(1.0+a) fb=1.0/(1.0+b) x=a+h f2=0.0 f4=1.0/(1.0+x) do 10 i=1,n-1 x=x+h f2=f2+1.0/(1.0+x) x=x+h f4=f4+1.0/(1.0+x) 10 continue s=h/3.0*(fa+fb+4.0*f4+2.0*f2) write(*,100) a,b,n write(*,150) s 100 format(1x,a=,f8.2,2x,b=,f8.2, $ 2x,n=,i4) 150 format(1x,s=,f16.7) end 三种求定积分的方法中，矩形法的误差较大，梯形法次之，辛普生法最好。 二、解一元方程（解非线性函数） 1. 直接迭代法 read(*,*) x,m do 10 i=1,m x1=(-x**3-2.0*x*x-2.0)/2.0 write(*,100) i,x1 if(abs(x-x1).gt.1e-6) then x=x1 else stop end if 10 continue write(*,200) m 100 format(1x,i=,i3,5x,x1=,f15.7) 200 format(1x,computati