计算数值方法实验报告.docVIP

课程名称： 计算机数值方法 实验项目：方程求根，线性方程组的直接解法与迭代 解法，代数插值，最小二乘法拟合多项式 实验地点： 逸夫楼402 专业班级： 学号： 学生姓名： 指导教师： 于亚男 2012年4月26日 太原理工大学学生实验报告 学院名称 软件学院 专业班级 学号 学生姓名 实验日期 4月26日 成绩 课程名称 计算机数值方法 实验题目 方程求根 实验目的和要求 （1）了解非线性方程求根的常见方法，如二分法、牛顿法、割线法。 （2）加深对方程求根方法的认识，掌握算法。 二、实验内容和原理 熟悉使用二分法、迭代法、牛顿法、割线法等方法对给定的方程进行根的求解。选择上述方法中的两种方法求方程：f(x)=x3+4x2-10=0在[1,2]内的一个实根，且要求满足精度|x*-xn|0.5×10-5 主要仪器设备 HP笔记本，VC++6.0 操作方法与实验步骤 1.二分法：#include stdio.h #include math.h double f(double x) { double y; y=pow(x,3)+4*pow(x,2)-10; return (y); } main() { double a,b,yl,yr,ym,x1,xr,xm; a=1; b=2; yl=f(a); yr=f(b); while(yl*yr0) { yl=f(a); yr=f(b); } xr=a; x1=b; while(fabs(xr-x1)0.000001) { xm=(x1+xr)/2; ym=f(xm); yl=f(x1); yr=f(xr); if(yl*ym0yr*ym0) xr=xm; else x1=xm; printf(%f\t%f\t%f\n,x1,xr,xm); } printf(\n所求根为x=%f\n,xm); return 0; } 2.牛顿法： #include stdio.h #include math.h double f(double x) { return (pow(x,3)+4*pow(x,2)-10); } double f1(double x) { return (3*pow(x,2)+10*x); } int main() { double x,x1,y1,y2; printf(请输入一个任意实数:X=); scanf(%lf,x); printf(方程的解为：\n); do{ x1=x; y1=f(x); y2=f1(x1); x=x1-y1/y2; } while (fabs(x-x1)=5e-6); printf( %lf\n,x1); } 3.追赶法 实验数据记录和处理 二分法 .牛顿法 实验结果与分析 通过这个两个程序可看出，二分法的计算量更大一些。 讨论、心得 通过这个实验，我了解了线性方程的一些求根方法，对于方程近似值的求解有了更多的理解。 实验地点 指导教师 于亚男 太原理工大学学生实验报告 学院名称 专业班级 学号 学生姓名 实验日期 4月26日 成绩 课程名称 计算机数值方法 实验题目 求解线性方程组 一、实验目的和要求 了解直接法解线性方程组 对Gauss消元法、LU分解法、追赶法有更深刻的理解 二、实验内容和原理 合理利用Gauss消元法、LU分解法、追赶法求解下列方程组： ① ② ③ ④（n=5,10,100,…） 三．主要仪器设备 HP笔记本，VC++6.0 四．操作方法与实验步骤 1.高斯： #includestdio.h void main() { int i,j,f,n; float a[5][5],b[5],x[5],l,k; printf(Input n:); scanf(%d,n); printf(Input a[][]:); for(i=0;in;i++) for(f=0;