用MATLAB计算光的等厚干涉实验中的不确定度

这学期终于有大物实验了，而且一个学期只需要做4个实验，交4份报告，考一场试，就没有别的任务了，那么周二下午的课简直形同虚设啊。快乐快乐~
又及：我们大物实验的老师是个用路飞头像的胖小哥，他总是让我们很快地把实验做掉，留下原始数据给他看，然后数据处理叫我们回去再算

0.FOREWORD

做实验自然有原始数据记录，以及数据处理，包括不确定度等的计算。正经人是不会去按计算器的（因为不会按…说真的，你如果真的懂手机计算器，按个arctan45°试试？ 因为一来操作不便，二来不利于数据的存储和整理）。

所以这里记录一些 matlab 的计算代码，虽然我也是现学现卖，但用代码写确实比按计算器来得有趣啊。

1.TASK

(直接拿课件截图多方便) 所以待计算的有：

1. 曲率半径均值
2. A类不确定度
3. B类不确定度
4. 合成不确定度
5. 相对不确定度

2.MATLAB CODE

光的等厚干涉实验的相关参数

R = [894.7 , 914.8 , 911.8 , 888.8 , 877.3];   % 曲率半径：待处理数据
len = length(R)
Dm = 6.697  % 20级暗纹直径
Dn = 5.857 %15级暗纹直径
m = 20
n = 15
lamda = 5.893*10.^(-4) % 钠灯光源波长
Ud = 0.005;    %实验仪器误差

均值计算

R_mean = mean(R);

A类不确定度

Ua = sqrt(sum((R(1,:)-R_mean).^2)/(len*(len-1)));

B类不确定度

Ub = Ud * sqrt(4 * Dm.^2 + 4 * Dn.^2)/(4*(m - n)*lamda);

合成不确定度

Ur = sqrt(Ua.^2 + Ub.^2); 

相对不确定度

Er = Ur/R_mean * 100; 

3.RESULTS

所有变量的值均可以在右边的工作区中找到。
以上。