# 利用最小二乘法
Xi, Yi, XX, YY, XY = 0, 0, 0, 0, 0
for i in range(0, int(N)):
Xi += x[i]
Yi += y[i]
XX += x[i]*x[i]
YY += y[i]*y[i]
XY += x[i]*y[i]
a = (Yi*Xi - N*XY)/(Xi*Xi - N*XX)
b = (Yi - a*Xi)/N
r = abs(Yi*Xi/N-XY)/math.sqrt((XX-Xi*Xi/N)*(YY-Yi*Yi/N))
# 进行绘制图像
mpl.rcParams['font.sans-serif'] = ['SimHei']
plt.plot(X, Y, 'ro-', color='blue', alpha=0.5, linewidth=1, label="线性拟合结果:Y={}X+{}".format(a, b))
plt.legend(loc="upper right")
plt.xlabel('截止频率V')
plt.ylabel('遏制电压U')
实验得到遏止电压后
分别输入运行:
结果得到为