数值分析（一） 牛顿插值法及matlab代码

数学： 数值分析

  刚上完数值分析课在其中学习了不少的知识，课后还做了一些课程实验主要都是利用matlab编程来解决问题，接下先讲插值法中的牛顿插值法

一、牛顿插值法原理

1.牛顿插值多项式

  定义牛顿插值多项式为：
$$N_n\left(x\right)=a_0+a_1\left(x-x_0\right)+a_2\left(x-x_0\right)\left(x-x_1\right)+\cdots +a_n\left(x-x_0\right)\left(x-x_1\right)\cdots \left(x-x_{n-1}\right)$$  其中$a_k\left(k=0,1,2,\cdots ,n\right)$为待定系数

  可见，牛顿插值多项式$N\left(x\right)$是插值多项式$P\left(x\right)$的另一种表示形式, 与Lagrange多项式相比它不仅克服了“增加一个节点时整个计算工作重新开始”的缺点, 且可以节省乘除法运算次数, 同时在Newton插值多项式中用到差分与差商等概念，又与数值计算的其他方面有密切的关系.

2.差商

2.1 定义

  自变量之差与因变量之差之比叫差商
  定义: 函数$y=f\left(x\right)$在区间$\left[x_i,x_{i+1}\right]$上的平均变化率$$f\left[x_i,x_{i+1}\right]=\frac{f\left(x_{i+1}\right)-f\left(x_i\right)}{x_{i+1}-x_i}$$  称为$f\left(x\right)$关于$x_i,x_{i+1}$的一阶差商，并记为$f\left[x_i,x_{i+1}\right]$
  二阶差商:$$f\left[x_i,x_{i+1},x_{i+2}\right]=\frac{f\left[x_{i+1},x_{i+2}\right]-f\left[x_i,x_{i+1}\right]}{x_{i+2}-x_i}$$  m阶差商:$$f\left[x_0,x_1,\cdots ,x_m\right]=$$