实验1插值和拟合.docVIP

《计算方法》实验报告 二级学院： 计算机学院 专 业： 计算机科学与技术 指导教师： 爨莹 班级学号： 计0901 姓 名： 张微玉  实验一 插值与拟合 实验目的： ① 通过编程和插值与拟合中的某种具体算法解决具体问题，更深一步的体会计算方法这门课的重要性，同时加深对插值与拟合公式某种具体算法的理解。 ② 熟悉编程环境。 2、实验要求: 实现插值与拟合中的某种具体算法编写并执行 3、实验内容： 1)用牛顿法求解的根，取初始值为10。 2) 用弦截法求解数学方程。 4、题目： 插值与拟合 5、原理： 1)用牛顿法求解的根的原理： 牛顿迭代法是以微分为基础的，微分就是用直线来代替曲线，由于曲线不规则，那么我们来研究直线代替曲线后，剩下的差值是不是高阶无穷小，如果是高阶无穷小，那么这个差值就可以扔到不管了，只用直线就可以了，这就是微分的意义。 牛顿法是牛顿在17世纪提出的一种求解方程f(x)=0.多数方程不存在求根公式，从而求精确根非常困难，甚至不可能，从而寻找方程的近似根就显得特别重要。 牛顿迭代法是取x0之后，在这个基础上，找到比x0更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x) = 0的根。 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逐步转化为线性方程来求解，即是依次用切线代替曲线，用线性函数的零点作为函数f(x)=0的近似值。 2)用弦截法求解根的基本思想是： 依次用弦线代替曲线，用线性函数的零点作为函数零点的近似值。 7、对应程序： 用牛顿法求解的根，取初始值为10 #include stdio.h #include math.h 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 %f\n,Newton(x0,e)); printf(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; } 用弦截法求解数学方程： #includeiostream #includeiomanip #includecmath using namespace std; double f(double); double xpoint(double,double); double root(double,double); int main() { double x1,x2,f1,f2,x; do { cout input x1,x2:; cin x1 x2; f1=f(x1);