av一区二区在线观看_亚洲男人的天堂网站_日韩亚洲视频_在线成人免费_欧美日韩精品免费观看视频_久草视

您的位置:首頁技術(shù)文章
文章詳情頁

用Python實現(xiàn)Newton插值法

瀏覽:2日期:2022-06-22 10:19:10
1. n階差商實現(xiàn)

def diff(xi,yi,n): ''' param xi:插值節(jié)點xi param yi:插值節(jié)點yi param n: 求幾階差商 return: n階差商 ''' if len(xi) != len(yi): #xi和yi必須保證長度一致return else:diff_quot = [[] for i in range(n)]for j in range(1,n+1): if j == 1:for i in range(n+1-j): diff_quot[j-1].append((yi[i]-yi[i+1]) / (xi[i] - xi[i + 1])) else:for i in range(n+1-j): diff_quot[j-1].append((diff_quot[j-2][i]-diff_quot[j-2][i+1]) / (xi[i] - xi[i + j])) return diff_quot

測試一下:

xi = [1.615,1.634,1.702,1.828]yi = [2.41450,2.46259,2.65271,3.03035]n = 3print(diff(xi,yi,n))

返回的差商結(jié)果為:

[[2.53105263157897, 2.7958823529411716, 2.997142857142854], [3.0440197857724347, 1.0374252793901158], [-9.420631485362996]]

2. 牛頓插值實現(xiàn)

def Newton(x): f = yi[0] v = [] r = 1 for i in range(n):r *= (x - xi[i])v.append(r)f += diff_quot[i][0] * v[i] return f

測試一下:

x = 1.682print(Newton(x))

結(jié)果為:

2.5944760289639732

3.完整Python代碼

def Newton(xi,yi,n,x): ''' param xi:插值節(jié)點xi param yi:插值節(jié)點yi param n: 求幾階差商 param x: 代求近似值 return: n階差商 ''' if len(xi) != len(yi): #xi和yi必須保證長度一致return else:diff_quot = [[] for i in range(n)]for j in range(1,n+1): if j == 1:for i in range(n+1-j): diff_quot[j-1].append((yi[i]-yi[i+1]) / (xi[i] - xi[i + 1])) else:for i in range(n+1-j): diff_quot[j-1].append((diff_quot[j-2][i]-diff_quot[j-2][i+1]) / (xi[i] - xi[i + j])) print(diff_quot)f = yi[0] v = [] r = 1 for i in range(n):r *= (x - xi[i])v.append(r)f += diff_quot[i][0] * v[i] return f

到此這篇關(guān)于用Python實現(xiàn)牛頓插值法的文章就介紹到這了,更多相關(guān)python牛頓插值法內(nèi)容請搜索好吧啦網(wǎng)以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持好吧啦網(wǎng)!

標(biāo)簽: Python 編程
相關(guān)文章:
主站蜘蛛池模板: 国产高清美女一级a毛片久久w | 欧美www在线观看 | 永久精品| 亚洲人成人一区二区在线观看 | 中文在线视频 | 成人欧美一区二区三区黑人孕妇 | 日日夜夜精品视频 | 电影91久久久 | 国产精品视频导航 | 欧美精品久久久 | 亚洲欧美国产视频 | 日日夜夜操天天干 | 日韩三级 | 日本三级网站在线 | 国产精品精品视频一区二区三区 | 成人三级视频在线观看 | 日韩视频在线播放 | 久在线精品视频 | 高清一区二区三区 | 精品福利在线 | 精品亚洲一区二区 | 成人毛片一区二区三区 | 亚洲欧洲在线观看视频 | 激情视频网站 | 久久新 | 一级特黄色毛片 | 国产伦精品一区二区三区四区视频 | 超碰人人艹| 天天射网站 | 亚洲精品一区二区在线观看 | 欧美日韩福利视频 | 91视频电影 | 国产原创视频 | 欧美精品日韩 | 国产免费一区二区 | 日韩精品一区二区三区高清免费 | 手机在线观看av | 成人在线一级片 | 精品一区二区三区在线视频 | 精品国产一区二区三区免费 | 黄色在线免费观看 |