实验一误差的基本性质与处理
(一)问题与解题思路:假定该测量列不存在固定的系统误差,则可按下列步骤求测量结
果
1、算术平均值
2、求残余误差
3、校核算术平均值及其残余误差
4、判断系统误差
5、求测量列单次测量的标准差
6、判别粗大误差
7、求算术平均值的标准差
8、求算术平均值的极限误差
9、写出最后测量结果
(二)在matlab中求解过程:
a = [24.674,24.675,24.673,24.676,24.671,24.678,24.672,24.674] ;%试验测得数据
x1 = mean(a) %算术平均值
b = a -x1 %残差
c = sum(b) %残差和
c1 = abs(c) %残差和的绝对值
bd = (8/2) *0.0001 %校核算术平均值及其误差,利用c1(残差和的绝对值)<=(n/2)*A 时,以上计算正确
% 3.5527e-015(c1) < 4.0000e-004(bd),以上计算正确
xt = sum(b(1:4)) - sum(b(5:8)) %判断系统误差,算的xt= 0.0030.由于xt较小,不存在系统误差
dc = sqrt(sum(b.^2)/(8-1)) %求测量列单次的标准差dc = 0.0022
sx = sort(a) %根据格罗布斯判断准则,先将测得数据按大小排序,进而判断粗大误差。
g0 = 2.03 %查表g(8,0.05)的值
g1 = (x1 - sx(1))/dc %解得g1 = 1.4000
g8 = (sx(8) - x1)/dc %解得g8 = 1.7361 由于g1和g8都小于g0,故判断暂不存在粗大误差sc = dc/sqrt(8) %算术平均值得标准差sc = 7.8916e-004
t=2.36; %查表t(7,0.05)值
jx = t*sc %算术平均值的极限误差jx = 0.0019
l1 = x1 - jx %测量的极限误差l1 = 24.6723
l2 = x1 + jx %测量的极限误差l2 = 24.6760
(三)在matlab中的运行结果