数值计算实验精读.docx

线性方程的数值解法一、课题名称：二分法求解非线性方程二、引言：目的：求方程f(x)=0的根二分法所使用的方法主要是介值定理，它的基本思想是逐步缩小根存在的区间，搜索根所在的位置许多实际算法表现为某种无穷递推过程的阶段，实现这类算法，不但需要建立计算公式，还需要解决精度控制问题，而我们下面介绍的二分法就是用来解决这一类问题，它是电子计算机上一种常用的算法。三、算法：首先，满足f(a)f(b)0,给定ε(根的容许误差)或δ(|f(x)|的容许误差)。 1、对分【a,b】，取x=(a+b)/2,计算f()； 2、若|b-a|ε,或|f()|δ,则取ζ=，停机； 3、判断根所在新区间 若f(a)f()≤0，则【】=【a,】 若f(a)f()0 ,则【】=【,b】=(b-a)/2 4、转向1，重复上述步骤。四、结构程序设计（求的根）program mainimplicit nonereal :: a=1.0000real :: b=1.5000real,parameter :: c=0.0050real x,y0,ydo while( b-a=c)y0=a**3-a-1x=(a+b)/2y=x**3-x-1if( y*y00) thena=xelseb=xend ifwrite (*,*)x,yend dostopend五、结果及讨论分析1、结果：2、误差分析：第k次二分后，则在二分过程中可以获得一个近似根的序列，该序列以根为极限 在实际计算时，人们不可能也没有必要完成这种无穷过程，因为计算结果允许带有一点的误差。由于只要二分足够多次（即k充分大），便有这里ε为预定精度，数值计算中的近似是正常的，计算误差是不可避免的。3、收敛性分析：对于二分法，只要能够保证在给定的区间内有根，使能够收敛的，当时收敛的速度和给定的区间有关，而且总体上来说速度比较慢。牛顿法(程序见附)，收敛速度要比二分法快，但是最终其收敛的结果与初值的选取有关，初值不同，收敛的结果也可能不一样，也就是结果可能不时预期需要得结果。4、算法评价：就二分法来说，它的优点有程序简单、方法可靠、收敛；缺点有无法解决多根、复根、一个切点的问题。附：1、简单迭代结构程序设计(求的根)：program mainimplicit nonereal::x0integer::kreal::x1integer,parameter::N=9real,external::fk=0x0=1.5do while(k=N)write(*,*)k,x0x0=f(x0)k=k+1end dofunction f(x)real::x,ff=(x+1)**(real(1)/real(3))returnend2、牛顿迭代结构程序设计（求）：program mainimplicit nonereal::x0,x1,x2real,parameter ::a=0.000005integer kinteger,parameter :: n=20read(*,*)x0x1=2*x0if(x1/=0) thendo k=0,nwrite (*,*)k,x0x2=x0-(x0**2-115)/x1if(abs(x2-x0)a) thenexitelse if(k/=n) then x0=x2end ifend doelsewrite(*,*)cuowuend ifstopend3、埃特金算法结构流程图（求x=在x=0.5附近的根）：program mainimplicit nonereal x0,x1,x2real,parameter ::a=0.0000005integer kinteger,parameter :: n=20read(*,*)x0do k=0,nx1=exp(-x0)x2=exp(-x1)x2=x2-(x2-x1)**2/(x2-2*x1+x0)write (*,*)k,x1,x2,x0if(abs(x2-x0)a) thenexitelse if(k/=n) thenx0=x2end ifend dostopend线性方程组的迭代法一、课题名称：雅可比迭代法求解线性方程组二、引言：雅可比迭代的目的是利用此方法计算出线性方程组符合一定精度的解，它的基本思想是将联立方程组的求解，归结为重复计算一组彼此独立的线性表达式，雅可比迭代的意义：在科学计算与工程计算中，线性方程组会经常遇到，因此，线性方程组的解法在数值分析中占有极其重要的地位。三、算法：考察一般形式的线性方程组，i=1,2,…,n (1)设从式（1）中分离出变量，将它改写成 ，i=1,2,┈，n (2)据此建立的迭代公式