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

(Of 它3少 TAIYUAN UNIVERSITY OF TECHNOLOGY 本科实验报告 课程名称: 计算数值方法 实验地点: 综合楼五层506室 专业班级: 计科 1002 学号： 2010001414 学生姓名: XXX 指导教师: 王峥 2012 年 6 月 20 太原理工大学学生实验报告 学院名称 计算机科学与技术 专业班级 计科1002 学号 2010001414 学生姓名 xxx 实验日期 2012.6.5 成绩 课程名称 计算数值方法 实验题目 实验一方程求根 、课题名称 方程求根：熟悉使用、迭代法、 牛顿法、割线法等方法对给定的方程进行根的求解。 二分法f(x)=x 3+4x2-10=0在[1,2]内的一个实根，且要求 选择上述方法中的两种方法求方程： 满足精度 |x*-xn|0.5 W-5 迭代法：用迭代公式 x=f(x) 二分法：设 f (x)在[a,b]上连续，且 f (a1) *f (x1) 0,记(a2, b2) = (x1, b1)带 进行迭代计算，直到满足|x -xn|0.5 W5为止。 入计算式进行计算 直到|x*-xn|v0.5 W5为止。 二、目的和意义 了解非线性方程求根的常见方法，如二分法、迭代法、牛顿法、割线法。 加深对方程求根方法的认识，掌握算法。会进行误差分析，并能对不同方法进行 比较。 、计算公式 (1)迭代法 ).首先对给定的计算公式进行变形使其能够迭代或者找出相应迭代速度较快的 式子。 ).带入求好的式子到循环中去比如： Xk4l=W(Xk) (k = 0,1,2,in) (2)二分法： f(x)在区间(x, y)上连续 .先找到a、b属于区间(x, y)，使f(a) , f(b) 定有零点，然后求 f[(a+b)/2], .如果f[(a+b)/2]=0 ，该点就是零点， 如果f[(a+b)/2]0, 则在区间((a+b)/2 , b)内有零点，反之在(a, (a+b)/2 )内 有零点带入1 )中继续。 异号, 说明在区间(a,b)内一 四、主要仪器设备 Vc++ 9.0 C-free CodeBlocks 五、结构程序设计 迭代法： #in cludestdio.h #in cludemath.h mai n() { int i; double xn [15],y,x1,x2,m ; printf(” 请输入 x1,x2 的值：\n); sca nf(%lf%lf, x1, x2); printf(”请输入精度要求：\n); sca nf(%lf,m); prin tf( n xn\n); i=0; do{ xn [0]=(x1+x2)/2 ; xn [i+1]= sqrt(10/(4+x n[i])); // 迭代 prin tf(%5d %5lf\n,i,x n[i]); y= fabs(x n[i+1]-x n[i]) ; i++; if(ym)break; }while(1); } 二分法： #in cludestdio.h #in cludemath.h mai n() { int m,n ,o, p; double a,b,l; printf(”请输入x^3, x^2, x的系数和常数 scan f(%d%d%d%d,m,n,o, p); //1 printf(请输入 x1,x2:\n); scan f(%lf%lf,a,b); //1 2 printf(请输入精度要求:\n); scan f(%lf, I); //0.5X10A5 printf(” n an double x,fx; int i=1; do { x=(b+a)/2; fx=m*x*x*x+ n*x*x+o*x+ p; p： \n”)； 4 0 -10 bn xn f(xn )\n); 分析讨论： 使用不同的方法，可以不同程度的求得方程的解，不同的方法速度不同，求得的结果也 稍有区别，当然和要求精度也有关系。刚开始的时候用数组对二分法进行求解，发现 循环 到第二次就无法实现值的传递，于是换了另外一种方法代替了数组。 流程图: 实验地点 综合楼五层506室 指导教师 王峥 { { 太原理工大学学生实验报告 学院名称 计算机科学与技术 专业班级 计科1002 学号 2010001414 学生姓名 xxx 实验日期 2012.6.5 成绩 课程名称 计算数值方法 实验题目 实验二线性方程组的直接解法 、课题名称 线性方程组的直接解法合理利用Gauss消元法、LU分解法、追赶法求解下列方程组:「1①0L23[「X1[「14[2HLxs」L13JX25.291-6.130-1 线性方程组的直接解法 合理利用Gauss消元法、 LU分解法、追赶法求解下列方程组: 「1 ①0 L2 3[「X1[「1