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) # 多项式