Python实现牛顿插值法(差商表)
admin
2023-07-18 04:42:35
0

def func(x,y,X,infor=True):
    list2=[y[0]]       #  差商表的对角线的第一个元素始终是y0
    count=1
    while(True):
        if len(y)>1:
            list=[]                 # 空列表用来保存,每次计算后差商表的行
            for i in range(len(y)-1):
                n=x[i+count]-x[i]
                m=y[i+1]-y[i]
                l=m/n
                list.append(l)
            list2.append(list[0])        # list2用来记录差商表的对角线元素,每计算一次,取行的第一个元素
            count += 1
            y = list
        else:
            break
    if infor:                           # 判断是否要继续计算,结果
        W=0
        for i in range(len(list2)):
            if i==0:
                w=list2[i]
            else:
                w = list2[i]
                for j in range(i):
                    w*=(X-x[j])
            W+=w
        print('牛顿插值:', W)
    return '牛顿差商表对角线列:%s' %list2

ret=func([0.32, 0.34, 0.36],[0.314567, 0.333487, 0.352274],'',infor=False)
print(ret)
ret=func([0.32, 0.34, 0.36],[0.314567, 0.333487, 0.352274],0.3367)
print(ret)

运行结果:
Python实现牛顿插值法(差商表)

相关内容

热门资讯

中国气象局升级暴雨应急响应为三... 预计6月7日至8日,西南地区东部、江南、华南等地有大到暴雨,其中,四川东南部、重庆中部、贵州东部和南...
办理临时身份证、免费乘车……各... 高考在即,为了给高考考生营造良好的考试环境,各地在交通、供电、食品安全、天气等多方面推出暖心服务。在...
菲防长绕后厨躲中国记者幕后:代... 据央视·玉渊谭天报道,菲律宾防长更换路线,绕道后厨,避开中国媒体,这几天相信很多人都看到了这条视频。...
荔枝和龙眼杂交成功了,预计五六... 荔枝和龙眼有“宝宝”了!近日,在广州天河区的国家荔枝种质资源圃内,广东省农业科学院果树研究所研究人员...
特朗普女儿女婿,在欧洲被“火烈... 坐拥权势圈层的顶级荣耀与光环,有时也必然要承受盛名附带的舆论反噬与风险考验。特朗普的女儿女婿估计也没...
两次高考两度考研,三甲医院科室... 澎湃新闻记者 陈斯斯 张呈君【编者按】又到一年高考时。对很多人来说,这是一次重要经历,承载着独一无二...
2026年高考今天开考:报名总... 2026年全国统一高考今天开考,今年高考报名总人数为1290万人,全国共设7981个考点,34.8万...
中企高管被捕后,下步怎么办? 截至今天,中国集装箱制造企业胜狮货柜高管马南庆在巴黎戴高乐机场,被法国警方以美国引渡申请为由抓捕,已...
人民日报钟声:日菲勾连,祸乱亚... 日菲勾连,祸乱亚太的危险一步(钟声)《人民日报》(2026年06月07日 第 03 版)亚太地区是和...
原创 华... 华为在影像方面,发展越来越成熟,而且推出红枫原色镜头、2亿像素、AI智拍、达芬奇人像引擎、实况功能等...