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

数值计算方法 实验报告 学院：计算机科学与技术学院 专业：计算机科学与技术 班级：0802班 学号：2008001424 姓名：王春生 太原理工大学学生实验报告 学院名称 计算机科学与技术学院 专业班级 计算机0802 学号 2008001424 学生姓名 王春生 实验日期 成绩 课程名称 数值计算方法 实验题目 实验一 方程求根 一、实验目的及意义： （1）了解非线性方程求根的常见方法，如二分法、牛顿法、割线法。 （2）加深对方程求根方法的认识，掌握算法。 （3）会进行误差分析，并能对不同方法进行比较。 二、实验任务： 选择上述方法中的两种方法求方程f(x)=x3+4x2-10=0在[1,2]内的一个实根，且要求满足精度|x*-xn|0.5×10-5 。 三、计算公式： （1）二分法公式 f（x）在区间[a,b]上连续，且f(a)·f(b)0,令x1=（a+b）/2，计算f(x1)。 若f(x1)=0,则x﹡=x1； 若f(a) ·f(x1)0,则记（a,b）=（a,x1）； 若f(a) ·f(x1)0,则记（a,b）=（x1,b）； 重复此过程，直至满足精度要求，（b-a）/2ε。 （2）牛顿法 迭代公式： 四、结构程序设计 1）二分法源代码： //二分法.cpp #include iostream #include iomanip using namespace std; double fx(double x) { return ((x*x*x)+(4*x*x)-10); } void main() { double A,B,mid,result,accuracy=0.000005; double a0=1.0,b0=2.0; coutf(x)=x3+4x2-10=0 [1,2] |x*-xn|0.5*10-5endl; cout************************************************endl; do { A=a0,B=b0; mid=(A+B)/2; if( fx(mid)==0) { result=mid; a0=A,b0=B; } else { if (fx(A)*fx(mid)0) { a0=A;b0=mid; coutA=setw(10)setprecision(10)a0 B=setw(10)setprecision(10)b0endl; } else if(fx(mid)*fx(B)0) { a0=mid;b0=B; coutA=setw(10)setprecision(10)a0 B=setw(10)setprecision(10)b0endl; } result=(a0+b0)/2; } } while((b0-a0)=accuracy); cout************************************************endl; coutthe result is X=setw(10)setprecision(10)resultendl; coutA=setw(10)setprecision(10)a0endl; coutB=setw(10)setprecision(10)b0endl; } 2）牛顿法源代码： //牛顿法.cpp //f(x)=x^3+4x^2-10=0 |x-xn|0.5*10^(-5) #include iostream #include iomanip using namespace std; double f(double x) { return (x*x*x+4*x*x-10); } double fo(double x) { return (3*x*x+8*x); } int main() { double x0=1,y; double temp,x1,xn; x1=x0; do { xn=x1-f(x1)/fo(x1); temp=x1; x1=xn; y=xn-temp; if (y0) y=(-y); coutxn=setprecision(6)xnendl; }while (y=0.000005); cout**************************************endl; coutresult is xnendl; return 0;