最大残差法的蒙特卡洛模拟研讨.pdfVIP

最大残差法的蒙特卡洛模拟研究 张海滨1，王中宇1，刘智敏2 (1．北京航空航天大学仪器科学与光电工程学院，北京100083：2．中国计量科学研究院。北京100013) 摘疆：本文研究了评定测量不确定度的最大残差法，它可以简单、迅速、可靠地计算出评定不确定度所需的标准 差。文中研究了最大残差的分布及其分布函数与数字特征的计算，并采用蒙特卡洛模拟法求出测量次数小于50时，各 分，甫函数、均值、标准差和自由度的具体数值，最后用实例验证了理论分析的正确性。 关键词：测量不确定度最大残差法蒙特卡洛模拟 l引言 在科学技术、生产、商业、贸易、医疗、环保等各个领域，进行着大量的测量工作，用以认识客 观事物。测量结果的质量如何?传统误差理论用“精度”一词来评价。精度是精密度和准确度的统称。 随着误差理论的发展，在计量测试领域，世界各国已广泛地采用测量不确定度来定量表示测量结果的质 量。 测量不确度的定义为：表征合理地赋予被测量之值的分散性，与测量结果相关联的参数。它是定量 说明测量结果的质量的一个参数。不确定度越小，测量结果的质量越高，其使用价值也越高；反之则 越低。不确定度在质量管理和质量保证体系中也极为重要，IS09001《质量体系一设计／开发、生产、 Guide25 安装和服务的质量保证模式》中规定，使用时应保证所用设备的不确定度为已知；ISO／IEC 《校准实验室和测试实验室能力的通用要求》中指出，实验室每个证书或报告，必须包含有关评定校 准或测试结果不确定度的说明。 按照BIPM建议INC—ln屯1，不确定度以标准差表示。标准差的评定方法口1常用贝塞尔法、最大残差 法、彼得斯法、和最大误差法等。其中最大残差法是由刘智敏于1979年提出的H1，并对行≤10，以及 门=15，20，25，30进行了研究，得出了以取以上值时的残差法系数表。该方法只需算出残差： u，：=，。一，，再找出最大残差maxlvl，以求标准差仃，计算过程简化了许多，也更易于掌握。另外， 通过下文分析可得，由最大残差法算标准差仃具有较高的精度。 本文使用MatLab软件进行蒙特卡洛模拟，得出以从2至50时的最大残差的分布函数、均值、标 准差及自由度。 2最大残差的分布函数与数字特征计算 2．1最大残差的分布 若测量值‘独立服从均值为∥标准差为仃的正态分布，即‘～Ⅳ(∥，仃)i=1，2…门 则残差为： q=‘一厂。‘一去喜‘ ·——337．．．． 因正态分布变量的多个线性函数仍服从正态分布，故随机向量(％，％，…，‰)服从正态分布，但 ∑q=o，故(q，呸，…，％)所服从的正态分布为退化的，我们研究了(q，呸，…，％一1)的分布，注意 t=l 到q的均值E(q)=o 故(q，呸，…，％一1)分布密度函数为： 八％b…稚扣南n-IeXp·羽磊n-I心h)zl朋协1 (2) (2万)了√iMI 式中M为(q，％，…，％一，)协分方差矩阵，l咏I为IMl中元素以的代数余子式。 由(1)式知，q的标准差q与q、q间相关系数岛为： 岛2—1／旷1 (3) 吼2 (4) 故(玎一1)×(刀一1)阶协方差矩阵为 刀一1． 1 ， 1 2 ——口‘ 一一盯‘ 一一仃‘ 刀 刀 以 1 ， 玎一1， 1 2 一一仃‘ ——盯‘ 一一盯‘ M=乃 玎