python线性拟合曲线代码

python线性拟合曲线代码

import numpy as np
import matplotlib.pyplot as plt
from scipy import optimize
plt.rcParams['font.sans-serif'] = [u'SimHei']
plt.rcParams['axes.unicode_minus'] = False             ##pycharm中无法显示中文字体


def f_1(x, a, b):
    return a * x + b


plt.figure()
# 拟合点
x0 = [0 ,20,40,60,80,100 ,120,140 ,160 ,180]
y0 = [0,-28,-58,-87,-115,-144,-173,-202,-231,-260]

# 绘制散点
plt.scatter(x0[:], y0[:], 3, "red")
# 直线拟合与绘制
A1, B1 = optimize.curve_fit(f_1, x0, y0)[0]
x1 = np.arange(0,185,0.1)     #30和75要对应x0的两个端点，0.01为步长
y1 = A1 * x1 + B1
plt.plot(x1, y1, "blue")
print(A1)
print(B1)

plt.title("应变片全桥曲线")
plt.xlabel("重量（g)")
plt.ylabel("电压（mv）单臂")
plt.show()