分类：Python牛顿插值法

2021
06-23

详解Python牛顿插值法

一、牛顿多项式拉格朗日多项式的公式不具备递推性，每个多项式需要单独构造。但很多时候我们需要从若干个逼近多项式选择一个。这个时候我们就需要一个具有递推关系的方法来构造——牛顿多项式这里的的a0，a1…等可以通过逐一带入点的值来求得。但是当项数多起来时，会发现式子变得很大，这个时候我们便要引入差商的概念（利用差分思想）具体见式子能更好理解这里在编程实现中我们可以推出相应的差商推导方程d(k,0)=y(k)d(k,j)=(d(k... 继续阅读 >

2021
05-02

用Python实现Newton插值法

1.n阶差商实现defdiff(xi,yi,n):"""paramxi:插值节点xiparamyi:插值节点yiparamn:求几阶差商return:n阶差商"""iflen(xi)!=len(yi):#xi和yi必须保证长度一致returnelse:diff_quot=[[]foriinrange(n)]forjinrange(1,n+1):ifj==1:foriinrange(n+1-j):diff_quot[j-1].append((yi[i]-yi[i+1])/... 继续阅读 >