5.3 函数逼近与拟合

import numpy as np
import matplotlib.pyplot as plt
from scipy import optimize

plt.rcParams['font.sans-serif']=['SimHei']    # 用来正常显示中文
plt.rcParams['axes.unicode_minus']=False      # 用来正常显示负号

#代码5-9
# 拉伸倍数
x = np.array([1.9,2,2.1,2.5,2.7,2.7,3.5,3.5,4,4,4.5,4.6,
              5,5.2,6,6.3,6.5,7.1,8,8,8.9,9,9.5,10])
# 强度
y = np.array([1.4,1.3,1.8,2.5,2.8,2.5,3,2.7,4,3.5,4.2,3.5,
              5.5,5,5.5,6.4,6,5.3,6.5,7,8.5,8,8.1,8.1])
x1 = np.arange(-1,10.0,0.1)
plt.figure(figsize=(8,6))
plt.scatter(x,y,label= '真实值')             # 绘制原来的点
plt.plot(x1,0.86*x1,'r',label = '近似直线')    # 近似直线
plt.legend()
#plt.xlim(0,10)
#plt.ylim(0,9)
plt.title('纤维线性拟合')
plt.xlabel('拉伸倍数')
plt.ylabel('强度')
#plt.savefig('../tmp/纤维线性拟合1.png')
plt.show()