还剩10页未读,继续阅读
本资源只提供10页预览,全部文档请下载后查看!喜欢就下载吧,查找使用更方便
文本内容:
第三章测量数据处理第一节测量误差的处理
一、系统误差的发现和减小系统误差的方法1,减小系统误差的方法
(1)采用修正的方法
(2)在实验过程中尽可能减少或消除一切产生系统误差的因素
(3)选择使系统误差抵消而不致带入测得值的测量方法2,试验和测量中常用的几种减小系统误差的测量方法
(1)恒定系统误差消除法
①异号法
②交换法
③替代法
(2)可变系统误差消除法
①对称测量法消除线性系统误差替代方案采用按“标准被校被校标准”顺序进行〜〜〜
②半周期偶数测量法消除周期性系统误差一一这种方法广泛用于测角仪上3,修正系统误差的方法
(1)在测得值上加修正值
(2)对测得值乘修正因子
(3)画修正曲线;
(4)制定修正值表4,获得修正值或修正因子的注意事项⑴修正值或修正因子的获得,最常用的方法是将测得值与计量标准的标准值比较得到,也就是通过校准得到修正曲线往往还需要采用实验方法获得
(2)修正值和修正因子都是有不确定度的在获得修正值或修正因子时,需要评定这些值的不确定度
(3)使用已修正测得值时,该测得值的不确定度中应该考虑由于修正不完善引入的不确定度分量
二、实验标准偏差的估计方法1,几种常用的实验标准偏差的估计方法
(1)贝塞尔公式法——适合于测量次数较多的情况⑶二任尹(3-6)计算步骤如下1)计算算术平均值2)计算10个残差3)计算残差平方和4)计算实验标准偏差
(2)极差法一般在测量次数较小时采用该法S(X)=(Xmax—Xmm)/C(3-8)
(3)较差法——适用于频率稳定度测量或天文观测等领域2,各种实验标准偏差估计方法的比较贝塞尔公式法是一种基本的方法,但n很小时其估计的不确定度较大,例如n=9时,由这种方法获得的标准偏差估计值的标准不确定度为25%,而n=3时标准偏差估计值的标准不确定度达50%,因此它适合于测量次数较多的情况极差法使用起来比较简便,但当数据的概率分布偏离正态分布较大时,应当以贝塞尔公式法的结果为准在测量次数较少时常采用极差法较差法更适用于随机过程的方差分析,如适用于频率稳定度测量或天文观测等领域
三、算术平均值及其实验标准偏差的计算、异常值的判别和剔除I1,什么是异常值异常值又称离群值,指在对一个被测量重复观测所获的若干观测结果中,出现了与其他值偏离较远且不符合统计规律的个别值,他们可能属于来自不同的总体,或属于意外的、偶然的测量错误2,用合成标准不确定度报告测量结果1基础计量学研究;2基本物理常量测量;3复现国际单位制单位的国际比对合成标准不确定度可以表示测量结果的分散性大小,便于测量结果间的比较3,带有扩展不确定度的测量结果报告的表示⑵测量结果及其扩展不确定度的报告形式扩展不确定度的报告有〃或两种
①历与的报告例如标准祛码的质量为m,,测量结果为
100.02147g,合成标准不确定度ucmJ为
0.35mg,取包含因子k=2,U=kuc仞=2X
0.35mg=
0.70mg,,一般,〃可用以下两种形式之一报告a.ms=
100.02147g;6^
0.70mg,62b.m=
100.02147±
0.00070g,62s
②如士儿仞的报告U.可用以下四种形式之一报告a.m=
100.02147;1195=
0.79mg,vf=
9.sc.ms=
100.02147±
0.00079g,p修9,括号内第二项为飒的值d.m=
100.0214779g,厅9,括号内为出的值,其末位与前面结果末位数对齐e.ms=
100.
021470.00079g,k9,括号内为45的值,与前面结果有相同的计量单位另外,给出扩展不确定度〃时,为了明确起见,推荐以下说明方式m=
100.02147土
0.00079g式中,s正负号后的值为扩展不确定度U95=k u,而合成标准不确定度u m=
0.35mg,自由度v f=9,包含因子95c cs efk5=t9=
2.26,从而具有约为95%概率的包含区间9954,相对扩展不确定度的表示
①相对扩展不确定度%el=u/y
②相对不确定度的报告形式举例a.m=
100.02147g;el=
0.70X106,4=2sb.m=
100.02147g;%⑹=
0.79X106osc.m=
100.021471+
0.79X
10.6应;片95%,々尸9,括号内第二项为相对扩展不确定度^rels5o也称为存在着“粗大误差”所以必须正确地判别和剔除异常值在测量过程中,记错、读错、仪器突然跳动、突然震动等异常情况引起的已知原因的异常值,应该随时发现,随时剔除,这就是物理判别法有时,仅仅是怀疑某个值,对于不能确定哪个是异常值时,可采用统计判别法进行判别2,判别异常值常用的统计方法⑴拉依达准则|Xd-x|三3s2格拉布斯准则:|Xd-x|/sNG a,n⑶狄克逊准则考前加强,出的可能性不大3,三种异常值判别准则的比较1当n〉50的情况下,3o准则较简便;3n50的情况下,格拉布斯准则效果较好,适用于单个异常值;有多于一个异常值时狄克逊准则较好2实际工作中,有较高要求的情况下,可选用多种准则同时进行,若结论相同,可以放心当结论出现矛盾,则应慎重,此时通常需选a=
0.01当出现既可能是异常值,又可能不是异常值的情况时,一般以不是异常值处理较好
五、测量重复性和测量复现性的评定
六、加权算术平均值及其实验标准偏差的计算方法
七、计量器具误差的表示与评定1,最大允许误差可以用绝对误差,相对误差,引用误差或它们的组合形式表示2,计量器具的示值误差是指计量器具即测量仪器的示值与相应测量标准提供的量值之差在计量检定时,用高一级计量标准所提供的量值作为约定值,称为标准值;被检仪器的指示值或标称值统称为示值则示值误差可以用下式表示示值误差二示值一标准值根据被检仪器的情况不同,示值误差的评定方法有比较法、分部法和组合法几种3,测量仪器示值误差符合性评定的基本要求评定示值误差的测量不确定度如或k=2时的与被评定测量仪器的最大允许误差的绝对值MPEV之比小于或等于13,即满足M1/3MPEV,时示值误差评定的测量不确定度对符合性评定的影响可忽略不计也就是合格评定误判概率很小合格14IWMPEV;不合格[△|〉MPEV对于型式评价和仲裁鉴定,必要时处与MPEV之比也可取小于或等于15o4,考虑示值误差的测量不确定度后的符合性评定当示值误差的测量不确定度历5或是k=2时的9与被评定测量仪器的最大允许误差的绝对值MPEV之比不满足小于或等于13的要求时,必须要考虑示值误差的测量不确定度对符合性评定的影响D合格判据当被评定的测量仪器的示值误差△的绝对值小于或等于其最大允许误差的绝对值MPEV与示值误差的扩展不确定度处之差时可判为合格,即[AlWMPEV-^5判为合格2不合格判据当被评定的测量仪器的示值误差△的绝对值大于或等于其最大允许误差的绝对值MPEV与示值误差的扩展不确定度则之和时可判不合格,即[△INMPEV+%判为不合格⑶待定区当被评定的测量仪器的示值误差既不符合合格判据又不符合不合格判据时,为处于待定区这时不能下合格或不合格的结论,即MPEV-65|△|〈MPEV+及5判为待定区当测量仪器示值误差的评定处于不能做出符合性判定时,可以通过采用准确度更高的计量标准、改善环境条件、增加测量次数和改善测量方法等措施,以降低示值误差评定的测量不确定度处后再进行合格评定对于只具有不对称或单侧允许误差限的被评定测量仪器,仍可按照上述原则进行符合性评定5,计量器具其他一些计量特性的评定
(一)准确度等级测量仪器的准确度等级应根据检定规程的规定进行评定有以下几种情况
(1)按最大允许误差评定准确度等级
(2)按示值的标准值的测量不确定度评定准确度等级
(3)测量仪器多个测量范围成多个参数时准确度等级的评定当被评定的测量仪器包含两个或两个以上的测量范围,并对应不同的准确度等级时,应分别评定各个测量范围的准确度等级对多参数的测量仪器,应分别评定各测量参数的准确度等级
(二)分辨力对测量仪器分辨力的评定,可以通过测量仪器的显示装置或读数装置能有效辨别的最小示值来确定
(1)带数字显示装置的测量仪器的分辨力为最低位数字显示变化一个步进量时的示值差
(2)用标尺读数装置(包括带有光学机构的读数装置)的测量仪器的分辨力为标尺上任意两个相邻标记之间最小分度值的一半
(三)灵敏度对被评定测量一起,在规定的某激励值上通过一个小的激励变化以,得到相应的响应变化△¥,,则比值s=4y厂x,即为该激励值时的灵敏度对线性测量仪器来说,灵敏度是一个常数
(四)鉴别阈对被评定测量仪器,在一定的激励和输出响应下,通过缓慢单方向地逐步改变激励输入,观察其输出响应使测量仪器产生恰能察觉有响应变化时的激励变化,就是该测量仪器的鉴别阈
(五)稳定性
(1)方法一通过测量标准观测被评定测量仪器计量特性的变化,当变化达到某规定值时,其变化量与所经过的时间间隔之比即为被评定测量仪器的稳定性
(2)方法二通过测量标准定期观测被评定测量仪器计量特性随时间的变化,用所记录的被评定测量仪器计量特性在观测期间的变化幅度除以其变化所经过的时间间隔,即为被评定测量仪器的稳定性
(3)方法三频率源的频率稳定性用阿伦方差的正平方根值评定,称频率稳定度当稳定性不是对时间而言时,应根据检定规程、技术规范或仪器说明书等有关技术文件规定的方法评定
(六)漂移根据技术规范要求,用测量标准在一定时间内观测被评定测量仪器计量特性随时间的慢变化,记录前后的变化值或画出观测值随时间变化的漂移曲线当测量仪器计量特性随时间呈线性变化时,漂移曲线为直线,该直线的斜率即漂移率在测得随时间变化的一系列观测值后,可以用最小二乘法拟合得到最佳直线,并根据直线的斜率计算出漂移率
(七)响应特性在确定条件下,激励与对应响应之间的关系称为测量仪器的响应特性评定方法是在确定条件下,对被评定测量仪器的测量范围内不同测量点输入信号,并测量输出信号当输入信号和输出信号不随时间变化时,记下被评定测量仪器的不同激励输入时的输出值,列成表格、画出曲线或得出输入输出量的函数关系式,即为测量仪器静态测量情况下的响应特性第二节测量不确定度的评定与表示
一、统计技术应用1(-)期望期望是概率分布曲线与横坐标轴所构成面积的重心所在的横坐标,所以期望是决定概率分布曲线位置的量
(二)方差方差就是随机误差平方的期望值方差说明了随机误差的大小和测量值的分散程度2,标准偏差是表明测量值分散性的参数,小表明测量值比较集中,大表明测量值比较分散用期望与标准偏差表征概率分布期望和方差是表征概率分布的两个特征参数由于方差不便使用,通常用期望和标准偏差来表征一个概率分布——u影响概率分布曲线的位置;对于单峰、对称的概率分布来说,期望值在分布曲线峰顶对应的横坐标处——o影响概率分布曲线的形状,表明测量值的分散性(小表明测量值比较集中,大表明测量值比较分散)期望与标准偏差都是以无穷多次测量的理想情况定义的,无法由测量得到U和2,因此都是概念性的术语3,实验标准偏差用有限次测量的数据得到的标准偏差的估计值称为实验标准偏差,用符号s表示在给出标准偏差的估计值时,自由度越大,表明估计值的可信度越高[n—1越大,1/n—1值越小,则其sx值也越小]4,正态分布曲线正态分布图,具有如下特征
①单峰概率分布曲线在均值u处具有一个极大值;
②对称分布正态分布以x=—U为其对称轴,分布曲线在均值u的两侧是对称的;
③当x—8时,概率分布曲线以x轴为渐近线;
④概率分布曲线在离均值等距离即x二u±o处两边各有一个拐点;
⑤分布曲线与x轴所围面积为1,即各样本值出现概率的总和为1;
⑥u为位置参数,为形状参数由于u,能完全表达正态分布的形态,所以常用简略符号x-nu表示正态分布当u=0,o=1时表示为X-n0,1,称为标准正态分布5,相关系数与协方差的关系1相关系数是一个纯数字,相关系数的值在一1到+1之间,它表示两个量的相关程度,通常比协方差更直观2协方差估计值sx,y与相关系数估计值rx,y的关系s[x9y=r[^ys[xsy^式中,sx和sy分别为X和y的实验标准偏差
二、法评定测量不确定度的步骤和方法GUM1,GUM法评定测量不确定度的步骤1明确被测量,必要时给出被测量的定义及测量过程的简单描述;2分析不确定度来源并写出测量模型;⑶评定测量模型中各输入量的标准不确定度ux.九计算灵敏系数ci从而给出与各输入量对应的输出量y的不确定度分量〃仞=1ci|ux\4计算合成标准不确定度〃仞,计算时应考虑各输入量之间是否存在值得考虑的相关性,对于非线性数学模型则应考虑是否存在值得考虑的高阶项;5列出不确定度分量的汇总表,表中应给出每一个不确定度分量的详细信息;6对被测量的概率分布进行估计,并根据概率分布和所要求的置信水平夕确定包含因子心⑺在无法确定被测量y的概率分布时,或该测量领域有规定时,也可以直接取包含因子k=2;8由合成标准不确定度〃0口和包含因子女或人的乘积,分别得到扩展不确定度〃或%;9给出测量不确定度的最后陈述,其中应给出关于扩展不确定度的足够信息利用这些信息,至少应该使用户能从所给的扩展不确定度进而评定其测量结果的合成标准不确定度简化步骤分析不确定度来源和建立测量模型一评定标准不确定度分量一计算合成标准不确定度火一确定扩展不确定度〃或%—报告测量结果’’2,通常测量不确定度来源从以下方面考虑1被测量的定义不完整⑵复现被测量的测量方法不理想⑶取样的代表性不够,即被测样本不能代表所定义的被测量⑷对测量过程受环境影响的认识不恰如其分或对环境的测量与控制不完善⑸对模拟式仪器的读数存在人为偏移⑹测量仪器的计量性能的局限性7测量标准或标准物质提供的量值的不准确⑻引用的数据或其他参量值的不准确⑼测量方法和测量程序的近似和假设10在相同条件下被测量在重复观测中的变化通常,在分析测量结果的不确定度来源时,可以从测量仪器、测量环境、测量方法、被测量等方面全面考虑,应尽可能做到不遗漏、不重复特别应考虑对测量结果影响较大的不确定度来源3,规范化常规测量时a类标准不确定度评定在规范化的常规测量中,测量m个同类被测量,得到m组数据,每组测量n次,第j组的平均值为x.,则合并样本标准偏差s」P〃一£,£Z Xq=4____________‘V/77H-13-58对每个量的测量结果的a类标准不确定度〃■1J=、/359359自由度为v=m n—1o3-60若对每个被测件的测量次数
4.不同,即各组的自由度工不等,各组的实验标准偏差为号则对于常规的计量检定或校准,当无法满足nN10时,为使得到的实验标准差更可靠,如果有可能,建议采用合并样本标准差Sp作为由重复性引入的标准不确定度分量4,标准不确定度的B类评定是借助于一切可利用的有关信息进行科学判断,得到估计的标准偏差
①根据有关信息或经验,判断被测量的可能值区间一a,a;
②假设测得值在区间内的概率分布;
③根据概率分布和要求的包含概论夕估计包含因子匕则B类评定的标准不确定度u为式中,v.=n.11ou=a/k3-62式中》为被测量可能值区间的半宽度;4为包含因子5,区间半宽度a值是根据有关信息确定的,一般情况下,可利用的信息包括:
①以前的观测数据;
②对有关技术资料和测量仪器特性的了解和经验;
③生产部门提供的技术说明文件制造厂的技术说明书;
④校准证书、检定证书、测试报告或其他提供的数据、准确度等级等;
⑤手册或某些资料给出的参考数据及其不确定度;
⑥规定测量方法的校准规范、检定规程或测试标准中给出的数据;
⑦其他有用信息6,B类评定时如何建设可能值的概率分布和确定k值
①概率分布的假设a.若被测量受许多相互独立的随机影响量的影响,这些影响量变化的概率分布各不相同,但各个变量的影响均很小时,被测量的随机变化服从正态分布b.如果有证书或报告给出的扩展不确定度是优、U绫或U.除非另有说明,可以按正态分布来评定B类标准不确定度c.一些情况下,只能估计被测量的可能值区间的上限和下限,测量值落在区间外的概率几乎为零若测量值落在该区间内的任意值的可能性相同,则可假设为均匀分布d.若落在该区间中心的可能性最大,则假设为三角分布e.若落在该区间中心的可能性最小,而落在该区间上限和下限处的可能性最大,则假设为反正弦分布£对被测量的可能值落在区间内的情况缺乏了解时,一般假设为均匀分布实际工作中,可依据同行专家的研究和经验来假设概率分布例如无线电计量中失配引起的不确定度为反正弦分布;几何量计量中度盘偏心引起的测角不确定度为反正弦分布;测量仪器最大允许误差、分辨力、数据修约、度盘或齿轮回差、平衡指示器调零不准等导致的不确定度按均匀分布考虑;两个独立量值之和或之差的概率分布为三角分布;按级使用量块时,中心长度偏差导致的概率分布为两点分布
②k值的确定a.已知扩展不确定度是合成标准不确定度的若干倍时,则该倍数包含因子就是4值b.假设概率分布后,根据要求的置信概率查表得到置信因子A值
③常用的概率分布与置信因子的关系见表3-12和表3-13o7,输入量间相关时的处理方法1在以下情况时可取协方差为零或忽略不计
①在不同实验室用不同测量设备、在不同时间测得的量值
②Xi与x.中任意一个量可作为常数处理
③独立测量的不同量的测量结果⑵用同时观测两个量的方法确定协方差估计值对两个输入量Xi及X.进行同时重复观测,设Xik,Xjk分别是输入量Xi及X.的观测值k为测量次数k=l,2,…,n「乙分别为第i个输入量和第j个输入量的k次测量的算术平均值;xi与x.的协方差估计值可o7------X3-74⑶用同时观测两个量的方法确定相关系数的估计值;3-75a⑷用经验公式估计相关系数如果两个输入量Xi,x.相关,土变化81会使x.相应变化变化3j,则xi和x.的相关系数可用经验公式3-76估计rx X硬®〜r:3-75b式中,以音和ux.分别Xi和X.的标准不确定度⑸当两个量均因与同一个量有关而相关时,协方差的估计方法设xi=fq,xj=gq,q是为使Xi与x.相关的变量q的估计值,f,g分别表示两个量与q的测量函数则xi与xj的协方差按式3—76a计算oq oq3—76a如果多个变量使Xi与X.相关,当、dFdGi^二尸肛耳吃二X,
2....%Gqq.©/时,则协方差按式3—76b计算/、+HF dG、3丁丁3-76bZ2/16采用适当方法去除相关性
①将引起相关的量作为独立的附加输入量进入数学模型
②采取有效措施变换输入量8,合成标准不确定度的有效自由度的计算合成标准不确定度,
①的自由度称为有效自由度,用符号表示在以下情况时需要计算有效自由度-纷’1当需要评定Up时为求得kp而必须计算,
①的自由度勺2当用户为了解所评定的不确定度的可靠程度而提出要求时有效自由度的计算公式当各分量间相互独立且输出量解决正态分布或t分布时,合成标准不确定度的有效自由度通常可按式3—77计算得到3-77Y二八乂V-
11.文八当测量模型为〜时,有效自由度可用相对标准不确定度的形式计算,式3—3-78实际计算中,得到的有效自由度Veff不一定是一个整数如果不是整数,可以采用将外ff数字舍位到最接近的一个较低的整数例如计算得到Ve产
12.65,则取V f=12oef9,明确规定包含概率时扩展不确定度如的评定方法当要求扩展不确定度所确定的区间具有接近于规定的包含概率夕时,扩展不确定度用符号即表示3—80Up^kp Ue片是包含概率为0时的包含因子⑴接近正态分布时kp的确定根据中心极限定理,当不确定度分量很多,且每个分量对不确定度的影响都不大时,其合成分布接近正态分布,此时若以算术平均值作为测量结果y,通常可假设概率分布为t分布,可以取人值为t值即3~81kp=tpMf根据合成标准不确定度wy的有效自由度r附和需要的置信水平口查表得到的t值即置信水平为夕的包含因子儿扩展不确定度〃二左ucy提供了一个具有包含概率置信水平为夕的区间y士Up.获得kp的计算步骤为
①求得测量结果的估计值y及其合成标准不确定度4仞
②按式3—82计算“仞的有效自由度Veff〃4q X/3-82式中,C为灵敏系数,〃氏/为输入量%的标准不确定度,力为〃出的自由度当〃为A类标准不确定度时是由〃次观测得到的S
⑨或SZ其自由度为片77-/;当〃6为B类标准不确定度时,用式3—83估计自由度内‘△〃・1£3-83式中,§u2)//的)是标准不确定度〃出)的相对不确定度,是所评定的〃出)的不可靠程度在实际工作中,B,类标准不确定度通常根据区间[—a,a]的信息来评定若可假设被测量值落在区间外的率极小,则可认为的评定是很可靠的,即§〃
⑨趋于0,此时,可假设忧幻的自由度匕-8
③根据要求的置信水平夕和计算得到的有效自由度欣查t分布的t值表得到卬/值
④取kp二讣Vejf),并计算Up=k Ucop
三、蒙特卡洛法评定测量不确定度的步骤和方法1,蒙特卡洛法简称MCM,是用概率分布传播的方法来评定测量不确定度蒙特卡洛法评定测量不确定度的方法
(1)建立测量模型;
(2)对每个输入量设定概率密度函数(PDF);
(3)选定蒙特卡洛试验数;⑷输入量概率分布的抽样及模型值计算;
(5)输出量分布函数的离散表示;
(6)输出量的估计值及其标准不确定度和包含区间;
(7)报告评定结果;
四、法与蒙特卡洛法的比较GUM1,GUM法与蒙特卡洛法
(1)GUM法通过不确定传播率计算合成标准不确定度,从而得到被测量估计值的测量不确定度的方法成为GUM法,即不确定度指南的方法GUM法的使用详见JJF
1059.1—2012《测量不确定度表示与评定》GUM法主要适用条件
1.可以假设输入量的概率分布呈对称分布;
2.可以假设输入量的概率分布近似为正态分布或t分布;
3.测量模型为线性模型、可以转换为线性的模型或可用线性模型近似的模型⑵蒙特卡洛法适用GUM法的条件MCM也都适用,除此之外,MCM对以下情况尤为有利1,测量模型明显呈非线性2,输入量的概率分布明显非对称3,输出量的概率分布较大程度地偏离正态分布或t分布,尤其是明显非对称分布第三节测量结果的处理和报告
一、最终报告时测量不确定度的有效位数1,在报告测量结果时,不确定度以〃或见0)都只能是『2位有效数字也就是说,报告的测量不确定度最多为2位有效数字建议当第一位有效数字是1或2时,应保留2位有效数字除此之外,对测量要求不高的情况可以保留1位有效数字测量要求较高时,一般取二位有效数字“近似值修约误差限的绝对值不超过末位的单位量值的一半”
二、报告测量结果的最佳估计值的有效位数的确定1,测量结果(即被测量的最佳估计值)的末位一般应修约到与其测量不确定度的末位对齐即同样单位情况下,如果有小数点,则小数点后的位数一样;如果是整数,则末位一致
三、测量结果的表示和报告1,完整的测量结果应包含
(1)被测量的最佳估计值⑵测量不确定度。