《计算机数值方法B》实验指导书.doc

《计算机数值方法》实验指导书计算机科学与技术与学院计算机科学与技术系目录实验一方程求根 1实验二线性方程组的直接解法 3实验三线性方程组的迭代解法 8实验四代数插值 9实验一方程求根【目的与要求】、基本原理、误差估计;,选择适当的算法;,进行算法设计和数值求解。【实验内容】、迭代法、牛顿法、割线法等方程求根的基本方法、基本原理、误差估计的相关理论。:f(x)=x3+4x2-10=0在[1,2]内的一个实根,且要求满足精度|x*-xn|<×10-5。【示例程序】增值寻根法又叫做逐步搜索法,是求解方程根的一种方法,也是初步确定方程隔根区间的一种方法。增值寻根法的基本思想是:从初值x0开始,按规定的一个初始步长h来增值:xn+1=xn+hn=0,1,2,3……同时计算f(xn+1)可能遇到三种情形:f(xn+1)=0,此时,xn+1即为方程的根x*;f(xn+1)*f(xn)>0,说明[xn+1,xn]内无根;f(xn+1)*f(xn)<0,说明[xn+1,xn]内有根;搜索过程,可从a开始,也可从b开始,这时应取更小的步长h,直到有根区间的长度|xn+1-xn|<e。(e为所要求的精度),此时f(xn)或f(xn+1)就可近似认为是零,xn或xn+1就是满足精度的方程的近似根。基于增值寻根法,设计程序如下:#include<>doublefun(doublex){ doubley=x*x; returny*x+4*y-10; }voidmain(void){ doublex0=-4,x; doubleh=1; doublestep=0; x=x0; while(h>0.){ cout<<"-------h="<<h<<"------"<<endl; while(true){ step++; if(fun(x)==0){ cout<<"x="<<x<<endl; h=0; break; } elseif(fun(x)*fun(x+h)>0){ x+=h; continue; }elseif(fun(x)*fun(x+h)<0){ cout<<"x1="<<x<<endl; cout<<"x2="<<(x+h)<<endl; break; } } h=h/10; } cout<<"step="<<step<<endl; }实验二线性方程组的直接解法【目的与要求】、LU分解法、追赶法等线性方程组直接求解的基本方法、基本原理;,选择适当的算法;,进行算法设计和数值求解。【实验内容】合理利用主元素消元法、LU分解法、追赶法求解下列方程组:①②③④(n=5,10,100,…)【示例程序】高斯消元法的基本思想是将原有线性方程组化为与之等价的三角形方程组,这个过程也称为消元过程,再通回代过程进行求解。消元过程使用公式为:回代过程使用公式如下:首先,,如下:322334771-245-7其中的3表示三阶线性方程组,下面是三行四列的系数矩阵。编写函数displayA()以显示系数矩阵及其的运算中的变化如下:voiddisplayA(){ printf("\n"); for(intj=1;j<=n;j++) { for(inti=1;i<=n+1;i++) printf("%lf",a[j][i]); printf("\n"); }}定义全局变量如下:doublea[15][16],a0[15][16];//a和a0用以记录方程组系数,//其中:a为本次运算的结果,a0为上次运算的结果doublel[15],tmp;//l用做消元过程使用系数l的计算和存储,tmp为零时变量intn;//n为线性方程组阶数完整的程序如下:#include""#include""doublea[15][16],a0[15][16];doublel[15],tmp;intn;voiddisplayA(){ printf("\n"); for(intj=1;j<=n;j++) { for(inti=1;i<=n+1;i++) printf("%lf",a[j][i]); printf("\n"); }}voidmain(){ FILE*f; inti,j,k; f=fopen("","r"); fscanf(f,"%d",&n); for(i=1;i<=n;i++) { for(j=1;j<=n+1;j++) { fscanf(f,"%lf",