《数值积分实验报告(共5页).doc》由会员分享，可在线阅读，更多相关《数值积分实验报告(共5页).doc（6页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上计算方法实验报告实验四 数值积分二级学院： 计算机学院 专 业： 计算机科学与技术 指导教师： 爨莹 班级学号： 姓 名： 实验一 非线性方程的迭代数值解法1、 实验目的： 通过编程和插值与拟合中的某种具体算法解决具体问题，更深一步的体会计算方法这门课的重要性，同时加深对插值与拟合公式某种具体算法的理解。 熟悉编程环境。2、实验要求: 实现插值与拟合中的某种具体算法编写并执行3、实验内容：1)用牛顿法求解的根，取初始值为10。 2) 用弦截法求解数学方程。 4、题目： 非线性方程的迭代数值解法5、原理：1) 用牛顿法求解的根的原理：牛顿迭代法是以微分为基础，用直线来代
2、替曲线，由于曲线不规则，那么研究直线代替曲线后，剩下的差值是不是高阶无穷小，如果是高阶无穷小，只用直线就可以了.。牛顿迭代法是取初始值x0之后，过曲线y= f(x)上的点（x0, f(x0)）做切线，切线方程为:y=f(x0)+f(x0)(x-x0), 它与x轴交点横坐标x1 为 x1=x0-f(x0)/f(x0)再过曲线y= f(x)上的点（x1, f(x1)）做切线，切线方程为y=f(x1)+f(x1)(x-x1) 它与x轴交点横坐标x2 为: x2=x1-f(x1)/f(x1)如此做下去，第n+1条切线方程为：y=f(xn)+f(xn)(x-xn) 它与x轴交点横坐标xn 为 Xn+1=
3、Xn-f(Xn)/f(Xn)在这个基础上，从而找到更接近方程根的近似跟。2) 用弦截法求解数学方程的原理：设 xk、xk+1是f(x)=0的近似根,我们利用f(xk),f(xk+1)构造一次插值多项式p1(x), 并用p1(x)=0的根作为f(x)=0的新的近似根xk+1,由于 p1(x)=f(xk)+f(xk)-f(xk-1)xk-xk-1(x-xk)(1)因此有 xk+1=xk-f(xk)f(xk)-f(xk-1)(xk-xk-1)(2)所以弦截法的几何意义为:依次用弦线代替曲线，用线性函数的零点作为函数零点的近似值6、设计思想： 1)用牛顿法求解的根的基本思想是：将非线性方程f(x)=0
4、逐步转化为线性方程来求解，即依次用切线代替曲线，用线性函数的零点作为函数f(x)=0的近似值。2)用弦截法求解根的基本思想是：依次用弦线代替曲线，用线性函数的零点作为函数零点的近似值。7、对应程序：用牛顿法求解#include #include double F1(double); double F2(double); double Newton(double,double);int main(int argc, int *argv)double x0 = 10.0; double e = pow(10,-5); printf(the result is %fn,Newton(x0,e); p
5、rintf(pause);double F1(double x)return x*x*x-155 ;double F2(double x) return 3*x*x ;double Newton(double x0, double e) double x1; do x1 = x0; x0 = x1 - F1(x1) / F2(x1); while (fabs(x1 - x0) e); return x0;用弦截法求解数学方程：#include#include#includeusing namespace std;double f(double); double xpoint(double,do
6、uble); double root(double,double);int main() double x1,x2,f1,f2,x; do cout x1 x2; f1=f(x1); f2=f(x2); while(f1*f2 = 0); x = root(x1,x2); cout setiosflags(ios:fixed) setprecision(7); cout A root of equation is x 0) y1 = y; x1 = x; else x2 = x;while(fabs(y)= 0.00001); return x;8、实验结果： 的根是 X= 5. 的根是 X= 37. 9、图形（如果可视化）10、实验体会： 插值与拟合中的牛顿法和弦截法都是通过对图形的观察和实际计算方法的运用，从而解决一些复杂的数值问题，我们应熟练掌握它，将其运用于我们的实际生活中。专心-专注-专业