数值分析牛顿插值法实验报告mrgp4.docVIP

数值分析牛顿插值法实验报告 数值分析牛顿插值法实验报告 篇一： 数值分析课程实验报告-拉格朗日和牛顿插值法 《数值分析》课程实验报告 用拉格朗日和牛顿插值法求解函数值 算法名称 用拉格朗日和牛顿插值法求函数值 学科专业 机械工程 作者姓名 程习康 作者学号 153711006 作者班级 机电院研究生15级C2班 中南大学 二〇一五年十二月 《数值分析》课程实验报告 篇二： 数值分析实验报告： 拉格朗日插值法和牛顿插值法 实验一报告 拉格朗日插值法 一、实验目的 1、学习和掌握拉格朗日插值多项式 2、运用拉格朗日插值多项式进行计算 二、实验原理 根据x0,x1,…xn；y0,y1,…yn构造插值多项式其表达式为： 将插值点x代入上式，就可得到函数f(x)在点x处的函数值的近似值。 三、运行结果 四、代码 using System; using System.Cllectins.Generic; using System.Linq; using System.Text; namespace CnsleApplicatin3 { class Prgram { static duble lglr(duble[] x, duble[] y, duble x1, int n) { duble result = 0.0; fr (int i = 0; i i++) { duble temp = y[i]; fr (int j = 0; j j++) { if (j == i) cntinue; temp = temp * (x1 - x[j]); temp = temp / (x[i] - x[j]); } result = result + temp; } return result; } static vid Main(string[] args) { duble[] x; duble[] y; Cnsle.riteLine( 请输入插值点数： int length = Cnvert.TInt32(Cnsle.ReadLine); x = ne duble[length]; y = ne duble[length]; fr (int i = 0; i length; i++) { Cnsle.rite( 请输入第{0}个点的x值： , i + 1); x[i] = Cnvert.TDuble(Cnsle.ReadLine); Cnsle.rite( 请输入第{0}个点的y值： , i + 1); y[i] = Cnvert.TDuble(Cnsle.ReadLine); } Cnsle.riteLine( 请输入x1值： duble x1 = Cnvert.TDuble(Cnsle.ReadLine); duble result=lglr(x,y,x1,length); Cnsle.rite( 插值计算结果为： {0}： , result); Cnsle.ReadLine; } 牛顿插值法 一、实验目的 体会并了解牛顿插值法，用计算机插入x值，输出相应的y值。 二、实验原理 根据x0,x1,…xn；y0,y1,…yn构造插值多项式Nn(x)=f(x0)+f(x0,x1)(x-x0)+… +f(x0,x1,…xn)(x-x0)(x-x1)…(x-xn-1)牛顿差值公式中各项的系数就是函f(x)的各阶均差（差商）f(x0)，f(x0,x1)，f(x0,x1,…xn)，因此，在构造牛顿差值公式时，常常先把均差列成一个表，此表称为均差表。 三、运行结果 四、代码 using System; using System.Cllectins.Generic; using System.Linq; using System.Text; namespace netn { class Prgram { static duble chashang(int n, duble[] x, duble[] y) { duble f = 0; duble temp = 0; fr (int i = 0; i n + 1; i++) { temp = y[i]; fr (int j = 0; j n + 1; j++) if (i != j) temp /= (x[i] - x[j]); f += temp; } return f; } static duble niudun(duble[] x, duble[] y, duble x1) { duble r