    # 利用最小二乘法Xi, Yi, XX, YY, XY = 0, 0, 0, 0, 0for 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)/Nr = 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')

实验得到遏止电压后

分别输入运行：

结果得到为

完整代码：

import math
from pylab import *
def 最小二乘法拟合(x, y):N = float(len(x))Xi, Yi, XX, YY, XY = 0, 0, 0, 0, 0for 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)/Nr = abs(Yi*Xi/N-XY)/math.sqrt((XX-Xi*Xi/N)*(YY-Yi*Yi/N))return a, b, rif __name__ == '__main__':Y = eval(input("请分别输入不同波长对应不同遏制电压:"))X = [5.198e14, 5.492e14, 6.882e14, 7.410e14, 8.216e14]Z = ["577.0nm 截止频率：", "546.1nm 截止频率：", "435.8nm 截止频率：", "404.7nm 截止频率：", "365.0nm 截止频率："]Z2 = ["577.0nm 截止频率：", "546.1nm 截止频率：", "435.8nm 截止频率：", "404.7nm 截止频率：", "365.0nm 截止频率："]a, b, r = 最小二乘法拟合(X, Y)H = 6.6260693e-34K = -a * 1.6e-19for i, j in enumerate(X):z = Z.pop(0)print("{} v{}: {}e14 Hz".format(z, i, j/1e14))for i, j in enumerate(Y):z2 = Z2.pop(0)print("{} U{}: {} V".format(z2, i, j))print("线性拟合结果:Y = {} X +{}, 相关系数={}".format(a, b, r))print("故求得普朗克常量为(h):|h| = k*e = {}J·s".format(K))print("公认普朗克值(H):{}J·s".format(H))print("相对误差: (h-H)/H = {}%".format(((K-H)/H)*100))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')plt.show()
# plt.savefig('demo.jpg')  # 保存该图片
plt.plot()

光电效应测普朗克常数实验，科大奥瑞平台计算Python程序相关推荐

科大奥锐思考题_科大奥锐光电效应测普朗克常数
姓名:学号:班级:教师:信箱号: 预约时间:第 _____ 周.星期 _____ .第 _____~_____ 节:座位号: _______ 科大奥锐<大学物理实验>报告一. 实验名称: ...
用matlab来画表格（实例：处理光电效应及普朗克常数的实验报告）
在处理光电效应及普朗克常数的数据时, 需要将对应数据绘制成表格,并画图.画图较为简单,在这就不提啦.表格的话,可以用excel这也可以.那么matlab该怎么画出那样的表格呢? 在matlab中有现成 ...
光电效应matlab处理图像,用matlab来画表格（实例：处理光电效应及普朗克常数的实验报告）...
在处理光电效应及普朗克常数的数据时, 需要将对应数据绘制成表格,并画图.画图较为简单,在这就不提啦.表格的话,可以用excel这也可以.那么matlab该怎么画出那样的表格呢? 在matlab中有现成 ...
密立根油滴实验科大奥瑞计算代码
静态法: #include <iostream> #include <cmath> using namespace std; int main(){ double q, ...
python程序项目报告_python程序设计实验报告
安徽理工大学Python编程实验报告类物流192名陆代学生号码3190505207年级日期教师秀于实验名称实验一个熟悉闲置和在线编程平台实验目的1.掌握python闲置的安装和使用集成开发环境2.熟悉 ...
科大奥瑞物理实验——光电效应和普朗克常量的测定
实验名称:光电效应和普朗克常量的测定 1. 实验目的: 通过实验了解光电效应的基本规律,并用光电效应法测量普朗克常量. 2. 实验器材: 光电管检流计(或微电流计) 直流电压计直流电源滤波片组 ...
科大奥瑞物理实验——光强调制法测光速
实验名称:光强调制法测光速 1. 实验目的: 掌握光强调制法测光速原理. 了解光速测定仪结构,并理解李萨如图像的成像原理. 2. 实验器材: 光速测定仪 1m长的水管直角反光镜透镜2个信号发生器 ...
matlab 输入普朗克常量,利用matlab和excel进行光电效应测普朗克常量实验中的数据处理...
利用matlab和excel进行光电效应测普朗克常量实验中的数据处理利用MATLAB和EXCEL 进行光电效应测普朗克常量实验中的数据处理 (东南大学学院南京211189) 摘要:本文回顾了光电 ...
科大奥瑞物理实验——迈克尔逊干涉仪实验
实验名称:迈克尔逊干涉仪实验 1. 实验目的: 掌握迈克尔逊干涉仪的干涉原理: 了解非定域干涉和时间相干性: 学会测量激光波长和介质的折射率. 学习一种测定光波长的方法,加强对等倾干涉的理解. 2. ...

光电效应测普朗克常数实验，科大奥瑞平台计算Python程序

实验得到遏止电压后

完整代码：

光电效应测普朗克常数实验，科大奥瑞平台计算Python程序相关推荐

最新文章

热门文章