计算数值实验报告(太原理工大学).doc

《计算数值实验报告(太原理工大学).doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、本科实验报告课程名称：计算数值方法实验地点：综合楼五层506室专业班级：计科1002学号：学生姓名：xxx指导教师：王峥2012年6月20太原理工大学学生实验报告学院名称计算机科学与技术专业班级计科1002学号学生姓名xxx实验日期2012.6.5成绩课程名称计算数值方法实验题目实验一方程求根一、课题名称方程求根：熟悉使用、迭代法、牛顿法、割线法等方法对给定的方程进行根的求解。选择上述方法中的两种方法求方程：二分法f(x)=x3+4x2-10=0在[1,2]内的一个实根，且要求满足精度

2、x*-xn

3、<0.5×10-5迭代法：用迭代公式x=f(x)进行迭代计算，直到满足

4、x*-xn

5、

6、<0.5×10-5为止。二分法：设f（x）在[a,b]上连续，且f（a1）*f（x1）<0,记（a2，b2）=（x1，b1）带入计算式进行计算直到

7、x*-xn

8、<0.5×10-5为止。二、目的和意义（1）了解非线性方程求根的常见方法，如二分法、迭代法、牛顿法、割线法。（2）加深对方程求根方法的认识，掌握算法。会进行误差分析，并能对不同方法进行比较。三、计算公式　（1）迭代法1）.首先对给定的计算公式进行变形使其能够迭代或者找出相应迭代速度较快的式子。2）.带入求好的式子到循环中去比如：（2）二分法：f(x)在区间（x，y）上连续　　1）.先找到a、b属于区间（x，y），使f(a

9、)，f(b)异号，说明在区间(a,b)内一定有零点，然后求f[(a+b)/2],　　2）.如果f[(a+b)/2]=0，该点就是零点，　　如果f[(a+b)/2]<0,则在区间（(a+b)/2，b)内有零点，反之在（a，(a+b)/2）内有零点带入1）中继续。　　　四、主要仪器设备Vc++9.0C-freeCodeBlocks五、结构程序设计迭代法:#include#includemain(){inti;doublexn[15],y,x1,x2,m;printf("请输入x1,x2的值：");scanf("%lf%lf",&x1,&x2);p

10、rintf("请输入精度要求：");scanf("%lf",&m);printf("nxn");i=0;do{xn[0]=(x1+x2)/2;xn[i+1]=sqrt(10/(4+xn[i]));//迭代printf("%5d%5lf",i,xn[i]);y=fabs(xn[i+1]-xn[i]);i++;if(y#includemain(){intm,n,o,p;doublea,b,l;printf("请输入x^3,x^2,x的系数和常数p：");scanf("%d

11、%d%d%d",&m,&n,&o,&p);//140-10printf("请输入x1,x2:");scanf("%lf%lf",&a,&b);//12printf("请输入精度要求:");scanf("%lf",&l);//0.5x10^5printf("nanbnxnf(xn)");doublex,fx;inti=1;do{x=(b+a)/2;fx=m*x*x*x+n*x*x+o*x+p;printf("%5d%5f%5f%5f%5f",i,a,b,x,fx);i++;if(fx==0)break;if(fx>0)b=x;elseif(fx<0)a=x;if((

12、b-a)

13、s消元法、LU分解法、追赶法求解下列方程组：①②③④（n=5,10,100…）二、目的和意义（1）了解线性方程组常见的直接解法，如Guass消元法、LU分解法、追赶法。（2）加深对线性方程组求解方法的认识，掌握算法。（3）会进行误差分析，并能对不同方法进行比较。三、计算公式高斯分解法：⑴将原方程组化为三角形方阵的方程组：lik=aik/akkaij=aij-lik*akjk=1,2,…,n-1i=k+1,k+2,…,nj=k+1,k+2,…,n+1⑵由回代过程求得原方程组的解：x