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

计算机数值方法计算机数值方法 实验指导书 (2017 版) 太原理工大学 计算机科学与技术学院 二一七年三月 0 目目 录录 实验一实验一 方程求根方程求根2 2 实验二实验二 线性方程组的直接解法线性方程组的直接解法4 4 实验三实验三 线性方程组的迭代解法线性方程组的迭代解法9 9 实验四实验四 代数插值代数插值1010 实验五实验五 最小二乘法拟合多项式最小二乘法拟合多项式1111 计算数值方法实验指导书 1 实验一 方程求根 【目的与要求目的与要求】 1了解方程求根的基本方法、基本原理、误差估计； 2能够按照工程实际要求，选择适当的算法； 3通过编写程序，进行算法设计和数值求解。 【实验内容实验内容】 1. 熟悉使用二分法、迭代法、牛顿法、割线法等方程求根的基本方法、基 本原理、误差估计的相关理论。 2. 选择方程求解方法中的两种方法求方程：f(x)=x3+4x2-10=0 在1,2内的一 个实根，且要求满足精度|x*-xn|0，说明xn+1, xn内无根； (3) f(xn+1)* f(xn) double fun(double x) double y=x*x; return y*x+4*y-10; void main(void) 计算数值方法实验指导书 2 double x0=-4,x; double h=1; double step=0; x=x0; while(h0.000001) cout0) x+=h; continue; else if(fun(x)*fun(x+h)=1;k-) tmp=0; for(j=k+1;j(%i,%i)“,k,k,maxL,maxC); /进行行变换 if(maxL!=k) for(i=1;i=1;k-) sum=0; for(j=k+1;j=0?“+%lfx%i“:“%lfx%i“,xi,i-1); printf(“n“); 程序运行结果如下： 已知点的个数 n=9 No. x y 1 1.000000 10.000000 2 3.000000 5.000000 3 4.000000 4.000000 4 5.000000 2.000000 5 6.000000 1.000000 6 7.000000 1.000000 7 8.000000 2.000000 8 9.000000 3.000000 9 10.000000 4.000000 sumX1=53.000000 sumY1=32.000000 拟和次数=2 sumX2=381.000000 sumX3=3017.000000 sumX4=25317.000000 sumY2=147.000000 sumY3=1025.000000 3 9.000000 53.000000 381.000000 32.000000 53.000000 381.000000 3017.000000 147.000000 381.000000 3017.000000 25317.000000 1025.000000 9.000000 53.000000 381.000000 32.000000 53.000000 381.000000 3017.000000 147.000000 381.000000 3017.000000 25317.000000 1025.000000 1.000000 2.000000 3.000000 4.000000 (1,1)(3,3) 计算数值方法实验指导书 15 25317.000000 3017.000000 381.000000 1025.000000 3017.000000 381.000000 53.000000 147.000000 381.000000 53.000000 9.000000 32.000000 3.000000 2.000000 1.000000 4.000000 l21=0.119169 25317.000000 3017.000000 381.000000 1025.000000 0.000000 21.467314 7.596635 24.851839 381.000000 53.000000 9.000000 32.000000 3.000000 2.000000 1.000000 4.000000 l31=0.015049 25317.000000 3017.000000 381.000000 1025.000000 0.000000 21.467314 7.596635 24.851839 0.000000 7.596635 3.266264 16.574594 3.000000 2.000000 1.000000 4.000000 (2,2)(2,2) 25317.000000 3017.000000 381.000000 1025.000000 0.000000 21.467314 7.596635 24.851839 0.000000 7.596635 3.266264 16.574594 3.000000 2.000000 1.000000 4.000000 l32=0.353870 25317.000000 3017.000000 381.000000 1025.000000 0.000000 21.467314 7.596635 24.851839 0.000000 0.000000 0.578044 7.780278 3.000000 2.000000 1.000000 4.000000 x1=13.459664 x2=-3.60530