实验报告误差(共9篇).doc

实验报告误差(共9篇) 误差分析实验报告 实验一 误差的基本性质与处理 (一) 问题与解题思路：假定该测量列不存在固定的系统误差，则可按下列步骤求测量结果 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.)/(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中的运行结果 实验二 测量不确定度 一、 测量不确定度计算步骤: 1. 分析测量不确定度的来源，列出对测量结果影响显著的不确定度分量； 2. 评定标准不确定度分量，并给出其数值 和自由度 ； 3. 分析所有不确定度分量的相关性，确定各相关系数 ； 4. 求测量结果的合成标准不确定度 及自由度 ； 5. 若需要给出伸展不确定度，则将合成标准不确定度 乘以包含因子k，得伸展不确定度 ； 二、 求解过程：用matlab编辑以下程序并运行 clc clear all close all D=[8.075 8.085 8.095 8.085 8.080 8.060]; h=[8.105 8.115 8.115 8.110 8.115 8.110]; D1=sum(D)/length(D);%直径的平均数 h1=sum(h)/length(D);%高度的平均数 V=pi*D1*h1/4; %体积 fprintf(#39;体积V的测量结果的估计值=%.1fmm#39;,V); fprintf(#39;不确定度评定： #39;); fprintf(#39;对体积V的测量不确定度影响显著的因素主要有：\n#39;); fprintf(#39;直径和高度的测量重复性引起的不确定度u1、u2，采用A类评定\n#39;);fprintf(#39;测微仪示值误差引起的不确定度u3，采用B类评定\n#39;); %%下面计算各主要因素引起的不确定度分量 fprintf(#39;直径D的测量重复性引起的标准不确定度分量u1,自由度v1\n#39;);M=std(D)2+(pi*D1/4))*(0.01/sqrt(3)) v3=round(1/(2*0.35*0.35)) fprintf(#39;不确定度合成：\n#39;); fprintf(#39;不确定度分量u1,u2,u3是相互独立的\n#39;); uc=round(sqrt(u1+u2+u3)*10)/10%标准不确定度 v=round(uc/(u1/v1+u2/v2+u3/v3))%自由度 fprintf(#39;展伸不确定度：\n#39;); fprintf(#39;取置信概率P=0.95,可查表得t=2.31，即包含因子k=2.31\n#39;);fprintf(#39;体积测量的展伸不确定度：\n#39;); P=0.95 k=2.31 U=ro