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爱片免费69 精品国产鲁一鲁一区二区张丽 | 精品产国自在拍 | 国产福利免费视频 | 国产成人精品一区二区三区视频 | 国产成人综合在线 | 亚洲午夜网 | 99国产精品99久久久久久 | 国产一区二区三区www | 日本视频在线播放 | 久久精品天堂 | 国产精品久久国产精品 | 日韩中文字幕在线观看 | 亚洲色图50p| 亚洲成人精品国产 | 四虎影院欧美 | 午夜专区 | 99精品久久久久 | 久草在线 | 亚洲国产成人av | 欧美三级电影在线播放 | 97精品国产一区二区三区 | 91精品国产综合久久久久久丝袜 | 91视频日本 | www.99热.com| 久久国产欧美日韩精品 | 亚洲综合视频一区 | 妞干网av | av网站免费观看 | 精品久久国产 | www.久草.com | 黄色一级大片在线免费看产 | 国产精品久久久久久吹潮 | 毛片免费观看 | 天天干成人网 | 香蕉二区| 中文字幕1区2区3区 日韩在线视频免费观看 | 国产精品久久久久久福利一牛影视 | 中文字幕免费视频 | 麻豆天堂 | 国内精品视频一区二区三区 | 欧美成人一区二区 |