FORTRAN 课程设计zwd【DOC精选】.docVIP

一、求一元方程的根 1、采用函数子程序定义一元方程； 2、程序选择以下三种方法求该方程的根； METHOD = 1 牛顿迭代法 METHOD = 2 二分法 METHOD = 3 弦截法 3、对于不同的近似算法分别编写子程序，精度要求10-6。 一.用二分法、弦解法和牛顿迭代法求x2-6x-7=0的根。 二分法的基本思路： 现任取两个值x1和x2，使得f（x1）* f（x2）0，也就是f（x1）和f（x2）必须异号。这才能保证在［x1，x2］区间有解，即存在一个x使得f（x）=0。 令x=（x1+x2）/2，如果f（x）=0，就找到了这个解，计算完成。由于f（x）是一个实型数据，所以在判断f（x）是否等于0时，是通过判断｜f（x）｜是否小于一个很小的数ε，如果是就认为f（x）=0。 若f（x）不等于0，判断如果f（x1）和f（x）异号，就说明解在［x1，x］区间，就以x1，x为新的取值重复步骤（2），这时用x代替否则x2，否则反之，直到找到满足条件的解为止。 程序编写如下： real x1,x2,x real bisect,func do print*,输入x1,x2的值 read*,x1,x2 if(func(x1)*func(x2)0.0)exit print*,不正确的输入 enddo x=bisect(x1,x2) print 10,x=,x 10 format(a,f15.7) end real function bisect(x1,x2) real x1,x2,x,f1,f2,fx x=(x1+x2)/2.0 fx=func(x) do while(abs(fx)1e-6) f1=func(x1) if(f1*fx0)then x2=x else x1=x endif x=(x1+x2)/2.0 fx=func(x) enddo bisect=x end function func(x) real x func=x**2-6*x-7 end 二分法运行结果 弦解法的基本思路： 现任取两个值x1和x2，使得f（x1）* f（x2）0。 做一条通过（x1，f（x1））和（x2，f（x2））两点的直线，这条直线与x轴的交点为x。可用以下公式求出 X=x2-（x2- x1）* f（x2）/（f（x1）- f（x2））， 代入函数求得f（x），判断｜f（x）｜是否小于一个很小的数ε，如果是就认为f（x）=0。 （3）否则，判断如果f（x1）和f（x）异号，就说明解在［x1，x］区间，就以x1，x为新的取值重复步骤（2），否则反之，然后以同样的办法再进一步缩小范围，直到｜f（x）｜ε。 程序编写如下： real x1,x2,x real secant,func do print*,输入x1,x2的值 read*,x1,x2 if(func(x1)*func(x2)0)exit print*,不正确的取值 enddo x=secant(x1,x2) print 10,x=,x 10 format(a,f15.7) end real function secant(x1,x2) implicit none real x1,x2,x,f1,f2,fx real func x=x2-(x2-x1)/(func(x2)-func(x1))*func(x2) fx=func(x) do while(abs(fx)1e-6) f1=func(x1) if(f1*fx0)then x2=x else x1=x endif x=x2-(x2-x1)/(func(x2)-func(x1))*func(x2) fx=func(x) enddo secant=x end real function func(x) real x func=x**2-6*x-7 end 弦解法运行结果 牛顿迭代法基本思路： （1）现任取一个值x1 （2）做一条通过（x1，f（x1））的切线，即以f＇（x1）为斜率作直线，直线与x轴的交点为x2， 因为f＇（x1）= f（x1）/(x1- x2） x2=x1- f（x1）/ f＇（x1） 判断｜f（x2）｜ε是否成立，如 果是就找到了这个解，计算完成。 （3）否则，重复步骤（2），以f＇（x1）为斜率做一条通过（x2，f（x2））的切线，直线与x轴的交点为x3，······，直到｜f（xn）｜ε，即xn为所得解。 程序编写如下： real x integer m print*,输入初值 read*,x call newton(x) end subrout