数值实验报告.docVIP

本科实验报告 课程名称： 计算机数值方法 实验项目： 方程求根、线性方程组的直接求解、线性方程组的迭代求解、最小二乘法拟合和代数插值 实验地点： 专业班级： 学号： 学生姓名： 指导教师： 2014年 5月 日  学生姓名 实验成绩 实验名称 实验一 方程求根 实验目的和要求（必填） （1）了解非线性方程求根的常见方法，如二分法、牛顿法、割线法。? （2）加深对方程求根方法的认识，掌握算法。? （3）会进行误差分析，并能对不同方法进行比较 实验内容和原理（必填） 方程求根： 熟悉使用二分法、牛顿法、割线法等方法对给定的方程进行根的求解。选择上述方法中的两种方法求方程：f(x)=x3+4x2-10=0在[1,2]内的一个实根，且要求满足精度|x*-xn|0.5×10-5 二分法原理： f(x)在区间（x，y）上连续????先找到a、b属于区间（x，y），使f(a)，f(b)异号，说明在区间(a,b)内一定有零点，然后求f[(a+b)/2],?? ?? 现在假设f(a)0,f(b)0,ab?? 如果f[(a+b)/2]=0，该点就是零点，??如果f[(a+b)/2]0,则在区间（(a+b)/2，b)内有零点，(a+b)/2=a，从①开始继续使用 割线法原理： ? 主要仪器设备 惠普 ProBook 6470b 实验记录(写出实验内容中的程序代码和运行结果)(可分栏或加页) 二分法： #include iostream using namespace std; #includestdio.h #include math.h int i; double a[20],b[20]; double x[20]; double A,y; int main() { a[0]=1.0,b[0]=2.0; for(i=0;i11;i++) { printf( a[i] b[i] x[i]y \n); x[i]=(b[i]+a[i])/2; y=x[i]*x[i]*x[i]+4*x[i]*x[i]-10; A=a[i]*a[i]*a[i]+4*a[i]*a[i]-10; if((y*A)0) a[i+1]=a[i],b[i+1]=x[i]; else a[i+1]=x[i],b[i+1]=b[i]; printf(%f, %f, %f, %f\n,a[i],b[i],x[i],y); } } 割线法 #includestdio.h #includemath.h #includeiostream using namespace std; float main() { float c,a=1.0,b=2.0; while(1) { c=b-(b*b*b+4*b*b-10)*(b-a)/(b*b*b+4*b*b-(a*a*a+4*a*a)); if(fabs(b-c)0.5*0.000001) break; b=c; } coutc; } 实验结果和分析 心得体会（遇到的问题和解决方法） 使用不同的方法，可以不同程度的求得方程的解，不同的方法速度不同，求得的结果也稍有区别，当然和要求精度也有关系。刚开始的时候用数组对二分法进行求解，发现?循环到第二次?就无法实现值的传递，于是换了另外一种方法代替了数组。  实验名称 实验二 线性方程组的直接求解 实验目的和要求 （1）了解线性方程组常见的直接解法，如Guass消元法、LU分解法、追赶法。? （2）加深对线性方程组求解方法的认识，掌握算法。? （3）会进行误差分析，并能对不同方法进行比较。 实验内容 合理选择利用Gauss消元法、LU分解法、追赶法求解下列方程 ?高斯分解法：?? ⑴将原方程组化为三角形方阵的方程组：? lik=aik/akk? ???aij=?aij-?lik*?akj??????k=1,2,…,n-1???? ???????????????????????????i=k+1,k+2,?…,n???j=k+1,k+2,?…,n+1 ⑵由回代过程求得原方程组的解：? ???xn=?ann+1/?ann? ???xk=(?akn+1-∑akj?xj)/?akk?????(k=n-1,n-2,?…,2,1) ?LU分解