python + numpy + np.polyfit()（最小二乘多项式拟合曲线）

import numpy as np

x = np.arange(1994, 2004, 1)
y = np.array([67.052, 68.008, 69.803, 72.024, 73.400, 72.063, 74.669, 74.487, 74.065, 76.777])

# 10代表拟合10次多项式，可以自由更改
z1 = np.polyfit(x, y, 10)  # 系数的集合
p1 = np.poly1d(z1)  # 多项式