「学习记录」《数值分析》第二章计算实习题(Python语言)

在假期利用Python完成了《数值分析》第二章的计算实习题,主要实现了牛顿插值法和三次样条插值,给出了自己的实现与调用Python包的实现——现在能搜到的基本上都是MATLAB版,或者是各种零碎的版本。
代码如下:
(第一题使用的自己的程序,第二第三题使用的Python自带库)

import math

import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
from numpy.linalg import solve
from scipy import interpolate
from scipy.interpolate import lagrange

plt.rc('figure',figsize=(20,15))

print("Problem I:")

given_x=[0.2,0.4,0.6,0.8,1.0]
given_y=[0.98,0.92,0.81,0.64,0.38]
given_times=4
x_range=(0,1.1,0.02)

#@brief: Convert(begin,end,interval) to a list, but interval can be float numbers.
def process_xpara(xpara):
    max_times=0
    if 0<xpara[2]<1:
        tmp_xpara_interval=xpara[2]
        while tmp_xpara_interval-int(tmp_xpara_interval)!=0:
            max_times=max_times+1
            tmp_xpara_interval=tmp_xpara_interval*10
    max_times=10**max_times
    return [i/max_times for i in range(int(xpara[0]*max_times),int(xpara[1]*max_times),int(xpara[2]*max_times))]


def divide_difference(x,y,times):
    now=[(x[i],y[i]) 
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值