例题7.4

合集 - 第七章例题加作业(9)

1.例题7.12024-11-062.例题7.22024-11-063.例题7.32024-11-06

4.例题7.42024-11-06

5.例题7.72024-11-066.例题7.102024-11-067.习题7.12024-11-068.习题7.32024-11-069.习题7.42024-11-06

收起

import numpy as np
import pylab as plt
from scipy.interpolate import interp1d
from scipy.interpolate import lagrange

a = np.loadtxt('F:\python数学建模与算法\源程序\《Python数学建模算法与应用》程序和数据\第7章  插值与拟合\data7_4.txt')
x0 = a[0]; y0 = a[1]
x = np.linspace(0, 15, 151)

# 分段线性插值
yx1 = interp1d(x0, y0)
y1 = yx1(x)

#lagrange插值
p2 = lagrange(x0, y0)
y2 = np.polyval(p2, x)
yx3 = interp1d(x0, y0, 'cubic')
y3 = yx3(x)
dx = np.diff(x); dy = np.diff(y3)
dyx = dy / dx; dyx0 = dyx[0]
xt = x[130:]; yt = y3[130:]
ymin = min(yt)
xmin = [xt[ind] for ind, v in enumerate(yt) if v == ymin]
print('x=0处斜率的数值解为:', dyx0)
print('xmin=', xmin); print('ymin=', ymin)
plt.rc('font', family='SimHei')
plt.rc('axes', unicode_minus=False)
plt.subplots_adjust(wspace=0.5)
plt.subplot(131); plt.plot(x, y1)
plt.title('分段线性插值')
plt.subplot(132); plt.plot(x, y2)
plt.title('拉格朗日插值')
plt.subplot(133); plt.plot(x, y3)
plt.title('三次样条插值')
plt.show()