测量不确定度评定指导书

测量不确定度评定指导书一、测量不确定度的基本概念（一）定义测量不确定度是表征合理地赋予被测量之值的分散性，与测量结果相联系的参数。它反映了测量结果不能肯定的程度，是对测量结果质量的定量表征。例如，在长度测量中，若得到测量结果为$(10.00\pm0.02)mm$，其中$\pm0.02mm$就是测量不确定度，表明被测量的真值以较高的概率落在$9.98mm$到$10.02mm$之间。（二）与测量误差的区别测量误差是测量结果与被测量真值之间的差值，是一个确定的值，但由于真值通常无法准确得知，所以测量误差也无法准确获取。而测量不确定度是基于对测量过程的分析和评定，给出的一个区间，用于合理地表征被测量值的分散性。例如，测量一个标准砝码的质量，假设真值为$100g$，测量结果为$100.01g$，测量误差为$0.01g$；但由于测量过程中存在各种不确定因素，评定得到的测量不确定度可能为$\pm0.005g$，这意味着被测量的真值有很大可能落在$99.995g$到$100.005g$之间。（三）重要性在科学研究、工业生产、贸易结算等众多领域，测量不确定度都具有重要意义。它是判断测量结果可靠性的重要依据，也是进行测量结果比对、数据融合和质量控制的基础。比如在国际贸易中，准确评定测量不确定度可以避免因测量结果的差异而引发的贸易纠纷；在工业生产中，通过评定测量不确定度可以优化生产工艺，提高产品质量。二、测量不确定度评定的基本程序（一）明确测量要求在进行测量不确定度评定之前，必须明确测量的目的、被测量的定义、测量范围、允许误差等要求。例如，在测量某机械零件的尺寸时，需要明确是测量长度、宽度还是高度，测量的精度要求是多少，以及测量结果的用途是用于产品验收还是质量控制等。（二）建立数学模型数学模型是描述被测量与各输入量之间关系的数学表达式。设被测量为$Y$，各输入量为$X_1,X_2,\cdots,X_n$，则数学模型可以表示为$Y=f(X_1,X_2,\cdots,X_n)$。例如，在测量长方形的面积$A$时，数学模型为$A=l\timesw$，其中$l$为长方形的长度，$w$为长方形的宽度。在建立数学模型时，需要考虑所有对测量结果有影响的输入量，包括直接测量的量和间接影响的量。同时，要对数学模型进行合理性验证，确保其能够准确反映测量过程。（三）识别不确定度来源测量不确定度的来源是多方面的，主要包括以下几个方面：测量设备：测量仪器的精度、稳定性、重复性等都会引入不确定度。例如，游标卡尺的分度值、示值误差，电子天平的灵敏度漂移等。测量环境：温度、湿度、气压、振动等环境因素的变化会对测量结果产生影响。比如在测量金属材料的长度时，温度的变化会导致材料热胀冷缩，从而引入不确定度。测量人员：测量人员的操作技能、读数习惯、视觉误差等也会带来不确定度。例如，在读取指针式仪表的示值时，不同的人员可能会有不同的读数结果。测量方法：测量方法的不完善、近似性等会引入不确定度。比如在采用近似公式进行计算时，会产生方法误差。被测量本身：被测量的不均匀性、稳定性等也会导致测量不确定度。例如，测量某合金材料的成分含量，由于材料成分分布不均匀，会引入不确定度。（四）评定标准不确定度A类评定：通过对观测列进行统计分析来评定标准不确定度，通常采用贝塞尔公式计算。对于一组重复测量的数据$x_1,x_2,\cdots,x_n$，样本标准差$s$的计算公式为：$$s=\sqrt{\frac{\sum_{i=1}^{n}(x_i-\bar{x})^2}{n-1}}$$其中$\bar{x}$为样本均值，$n$为测量次数。标准不确定度$u(x_i)$等于样本标准差$s$。例如，对某一长度进行10次重复测量，得到的数据分别为$10.01mm,10.02mm,\cdots,10.10mm$，计算得到样本标准差为$0.03mm$，则A类标准不确定度为$0.03mm$。B类评定：基于经验或其他信息的概率分布来评定标准不确定度。当无法通过重复测量得到观测列时，需要采用B类评定方法。例如，根据测量仪器的说明书，已知其最大允许误差为$\Delta$，假设其服从均匀分布，则标准不确定度$u(x)$为$\Delta/\sqrt{3}$。又如，根据以往的测量经验，某输入量的估计值的不确定度为$a$，且其服从正态分布，置信水平为$95%$，则标准不确定度$u(x)$为$a/1.96$。（五）计算合成标准不确定度当被测量$Y$由多个输入量$X_1,X_2,\cdots,X_n$通过数学模型$Y=f(X_1,X_2,\cdots,X_n)$计算得到时，合成标准不确定度$u_c(y)$可以通过以下公式计算：$$u_c(y)=\sqrt{\sum_{i=1}^{n}\left(\frac{\partialf}{\partialX_i}\right)^2u^2(X_i)+2\sum_{1\leqi<j\leqn}\frac{\partialf}{\partialX_i}\frac{\partialf}{\partialX_j}u(X_i,X_j)}$$其中$\frac{\partialf}{\partialX_i}$为灵敏系数，表示被测量$Y$对输入量$X_i$的变化的敏感程度；$u(X_i)$为输入量$X_i$的标准不确定度；$u(X_i,X_j)$为输入量$X_i$和$X_j$之间的协方差。在实际计算中，如果各输入量之间相互独立，协方差$u(X_i,X_j)=0$，则合成标准不确定度的计算公式可以简化为：$$u_c(y)=\sqrt{\sum_{i=1}^{n}\left(\frac{\partialf}{\partialX_i}\right)^2u^2(X_i)}$$例如，在测量长方形面积$A=l\timesw$的例子中，假设长度$l$的标准不确定度为$u(l)=0.01mm$，宽度$w$的标准不确定度为$u(w)=0.02mm$，灵敏系数$\frac{\partialA}{\partiall}=w$，$\frac{\partialA}{\partialw}=l$，若$l=100mm$，$w=50mm$，则合成标准不确定度为：$$u_c(A)=\sqrt{(50\times0.01)^2+(100\times0.02)^2}=\sqrt{0.25+4}=\sqrt{4.25}\approx2.06mm^2$$（六）扩展不确定度的评定扩展不确定度$U$是将合成标准不确定度$u_c(y)$乘以包含因子$k$得到的，即$U=k\timesu_c(y)$。包含因子$k$的取值取决于被测量的概率分布和所要求的置信水平。正态分布：当被测量服从正态分布时，若要求置信水平为$95%$，包含因子$k=1.96$；若要求置信水平为$99%$，包含因子$k=2.58$。例如，合成标准不确定度为$0.05mm$，置信水平为$95%$，则扩展不确定度$U=1.96\times0.05=0.098mm$。均匀分布：当被测量服从均匀分布时，若要求置信水平为$100%$，包含因子$k=\sqrt{3}\approx1.73$。（七）报告测量不确定度测量不确定度的报告应包括以下内容：被测量的定义和测量结果。合成标准不确定度和扩展不确定度的数值。包含因子$k$的取值和置信水平。测量不确定度的评定方法和依据。对测量结果有影响的主要不确定度来源。例如，测量某零件的长度，测量结果为$50.00mm$，合成标准不确定度为$0.02mm$，扩展不确定度为$0.04mm$（$k=2$，置信水平为$95%$），则测量不确定度报告可以表述为：“被测量为零件的长度，测量结果为$50.00mm$，合成标准不确定度$u_c=0.02mm$，扩展不确定度$U=0.04mm$（$k=2$，置信水平$95%$）。测量不确定度主要来源于测量仪器的示值误差和测量环境的温度变化。”三、测量不确定度评定的实例分析（一）实例一：长度测量测量要求：测量某机械零件的长度，要求测量精度达到$\pm0.01mm$，测量结果用于产品验收。数学模型：$L=x$，其中$L$为被测量的长度，$x$为测量仪器的示值。不确定度来源识别：测量设备：使用的游标卡尺的分度值为$0.02mm$，示值误差为$\pm0.01mm$。测量环境：测量环境温度为$(20\pm1)^{\circ}C$，零件材料的线膨胀系数为$12\times10^{-6}/^{\circ}C$，由于温度变化会导致零件长度发生变化，从而引入不确定度。测量人员：测量人员的读数误差估计为$\pm0.005mm$。标准不确定度评定：A类评定：对零件长度进行10次重复测量，得到的数据分别为$20.01mm,20.02mm,\cdots,20.10mm$，计算样本均值$\bar{x}=20.055mm$，样本标准差$s=0.03mm$，则A类标准不确定度$u_A=0.03mm$。B类评定：游标卡尺的示值误差服从均匀分布，标准不确定度$u_1=0.01/\sqrt{3}\approx0.0058mm$。温度变化引入的不确定度：温度变化量$\DeltaT=1^{\circ}C$，零件长度的变化量$\DeltaL=L\times\alpha\times\DeltaT=20\times12\times10^{-6}\times1=2.4\times10^{-4}mm$，假设温度变化服从均匀分布，标准不确定度$u_2=2.4\times10^{-4}/\sqrt{3}\approx1.39\times10^{-4}mm$。测量人员的读数误差服从均匀分布，标准不确定度$u_3=0.005/\sqrt{3}\approx0.0029mm$。合成标准不确定度计算：由于各输入量之间相互独立，合成标准不确定度$u_c=\sqrt{u_A^2+u_1^2+u_2^2+u_3^2}=\sqrt{0.03^2+0.0058^2+(1.39\times10^{-4})^2+0.0029^2}\approx0.0307mm$。扩展不确定度评定：假设被测量服从正态分布，置信水平为$95%$，包含因子$k=1.96$，则扩展不确定度$U=1.96\times0.0307\approx0.06mm$。测量不确定度报告：“被测量为机械零件的长度，测量结果为$20.06mm$，合成标准不确定度$u_c=0.03mm$，扩展不确定度$U=0.06mm$（$k=1.96$，置信水平$95%$）。测量不确定度主要来源于测量仪器的示值误差、测量环境的温度变化和测量人员的读数误差。”（二）实例二：电学量测量测量要求：测量某电阻的阻值，要求测量精度达到$\pm0.1%$，测量结果用于电路设计。数学模型：$R=U/I$，其中$R$为被测量的电阻值，$U$为测量的电压值，$I$为测量的电流值。不确定度来源识别：测量设备：电压表的精度等级为$0.5$级，量程为$0-10V$；电流表的精度等级为$0.5$级，量程为$0-1A$。测量环境：测量环境温度为$(25\pm2)^{\circ}C$，电阻的温度系数为$5\times10^{-4}/^{\circ}C$，温度变化会影响电阻的阻值。测量方法：采用伏安法测量电阻，由于电压表和电流表的内阻会引入测量误差。标准不确定度评定：A类评定：对电阻的阻值进行8次重复测量，得到的数据分别为$100.1\Omega,100.2\Omega,\cdots,100.8\Omega$，计算样本均值$\bar{R}=100.45\Omega$，样本标准差$s=0.25\Omega$，则A类标准不确定度$u_A=0.25\Omega$。B类评定：电压表的最大允许误差为$\DeltaU=0.5%\times10V=0.05V$，服从均匀分布，标准不确定度$u_U=0.05/\sqrt{3}\approx0.0289V$。电流表的最大允许误差为$\DeltaI=0.5%\times1A=0.005A$，服从均匀分布，标准不确定度$u_I=0.005/\sqrt{3}\approx0.00289A$。温度变化引入的不确定度：温度变化量$\DeltaT=2^{\circ}C$，电阻阻值的变化量$\DeltaR=R\times\alpha\times\DeltaT=100\times5\times10^{-4}\times2=0.1\Omega$，假设温度变化服从均匀分布，标准不确定度$u_T=0.1/\sqrt{3}\approx0.0577\Omega$。测量方法引入的不确定度：由于电压表内阻$R_V$和电流表内阻$R_A$的影响，测量值$R_{测}=U/I$与真实值$R_{真}$之间存在偏差。假设$R_V\ggR_{真}$，$R_A\llR_{真}$，则测量误差$\DeltaR/R_{真}\approxR_A/R_{真}$，已知电流表内阻$R_A=0.1\Omega$，则相对误差为$0.1/100=0.1%$，服从均匀分布，标准不确定度$u_{rel}=0.1%/\sqrt{3}\approx0.0577%$，标准不确定度$u_M=100\times0.0577%=0.0577\Omega$。合成标准不确定度计算：首先计算各输入量的相对标准不确定度：$u_{rel}(U)=u_U/U$，假设测量电压$U=5V$，则$u_{rel}(U)=0.0289/5=0.00578$。$u_{rel}(I)=u_I/I$，假设测量电流$I=0.05A$，则$u_{rel}(I)=0.00289/0.05=0.0578$。$u_{rel}(T)=u_T/R=0.0577/100=0.000577$。$u_{rel}(M)=0.000577$。$u_{rel}(A)=u_A/R=0.25/100=0.0025$。相对合成标准不确定度$u_{rel}(R)=\sqrt{u_{rel}(U)^2+u_{rel}(I)^2+u_{rel}(T)^2+u_{rel}(M)^2+u_{rel}(A)^2}=\sqrt{0.00578^2+0.0578^2+0.000577^2+0.000577^2+0.0025^2}\approx0.0582$。合成标准不确定度$u_c(R)=u_{rel}(R)\timesR=0.0582\times100=5.82\Omega$。6.扩展不确定度评定：假设被测量服从正态分布，置信水平为$95%$，包含因子$k=1.96$，则扩展不确定度$U=1.96\times5.82\approx11.41\Omega$，相对扩展不确定度$U_{rel}=U/R=11.41/100=11.41%$。7.测量不确定度报告：“被测量为电阻的阻值，测量结果为$100.45\Omega$，合成标准不确定度$u_c=5.82\Omega$，扩展不确定度$U=11.41\Omega$（$k=1.96$，置信水平$95%$）。测量不确定度主要来源于电压表和电流表的示值误差、温度变化、测量方法的不完善以及重复测量的随机性。”四、测量不确定度评定中的常见问题及解决方法（一）不确定度来源识别不充分在实际评定过程中，常常会出现不确定度来源识别不充分的情况，导致评定结果不准确。解决方法是对测量过程进行全面、细致的分析，从测量设备、测量环境、测量人员、测量方法、被测量本身等多个方面进行考虑，列出所有可能的不确定度来源。同时，可以参考相关的标准、规范和文献，借鉴他人的评定经验。（二）数学模型建立不合理数学模型建立不合理会直接影响测量不确定度的评定结果。解决方法是深入理解被测量的物理意义和测量过程，确保数学模型能够准确反映被测量与各输入量之间的关系。在建立数学模型后，要进行合理性验证，可以通过实验比对、理论分析等方法来验证数学模型的正确性。（三）标准不确定度评定方法选择不当A类评定和B类评定方法的选择应根据具体情况而定。如果可以通过重复测量得到足够多的数据，应优先采用A类评定方法；如果无法进行重复测量，或者测量成本过高，则采用B类评定方法。在进行B类评定时，要合理选择概率分布和包含因子，确保评定结果的准确性。（四）合成标准不确定度计算错误合成