matlab误差分析报告,误差分析实验报告

实验一误差的基本性质与处理

(一)问题与解题思路：假定该测量列不存在固定的系统误差，则可按下列步骤求测量结

果

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中的运行结果