（十一）牛顿差值高阶拟合

最新推荐文章于 2024-10-28 17:33:40 发布

weixin_30487701

最新推荐文章于 2024-10-28 17:33:40 发布

阅读量225

点赞数

CC 4.0 BY-SA版权

文章标签： python

原文链接：http://www.cnblogs.com/the-wang/p/8021536.html

 1 #coding=utf-8
 2 from numpy import  *
 3 import  numpy as  np
 4 
 5 x = [0.40,0.55,0.65,0.80,0.90,1.05];
 6 y = [[0 for a in range(6)] for b in range(10)];
 7 y[0][0] = 0.41075;
 8 y[1][0] = 0.57815;
 9 y[2][0] = 0.69675;
10 y[3][0] = 0.88811;
11 y[4][0] = 1.02653;
12 y[5][0] = 1.25382;
13 
14 # 计算 一列一列的计算 i为列 j为行
15 for i in range(1,len(x)):
16     for j in range(0,len(x)-i):
17         x2 = j+1;
18         y2 = i-1;
19         x1 = j;
20         y1 = i-1;
21         zi = (y[x2][y2]-y[x1][y1])/(x[j+i]-x[j]);
22         y[j][i] = zi;
23 # #输出  第一行为差商
24 for i in range(1,len(x)):
25     show="第 "+str(i)+" 阶差商";
26     print show;
27     print y[0][i];