多输出量测量模型的不确定度评定方法

'2017年6月宇航计测技术Jun．，2017第37卷第3期JoumalofAstronauticMetr010胛andMeasurementV01．37．No．3文章编号：1000—7202(2017)03—0001一04DOI：10．12060／j．issn．1000一7202．2017．03．0l多输出量测量模型的不确定度评定方法谌贝1龚鹏伟1谢文1姜河1马红梅1杨春涛1’2(1．北京无线电计量测试研究所；2．计量与校准技术重点实验室，北京100039)摘要本论文基于GuM不确定度评定方法，提出了多输出量测量模型的测量不确定度评定方法，并对如何获得多输出量测量模型的包含范围进行了举例说明。关键词CuM多输出量不确定度包含范围中图分类号：TB9文献标识码：AAMethodofMeasurementUncertaintyEValuationforMultiVariateTⅥeasurementModelsCHENBeilGONGPeng·weilXIEWenlJIANGHelMAHong—meilYANGChun-ta01"2(1．BeijingInstituteofRadioMetmlogyandMeasurement；2．NationalKeyLaboratoryofMetrologyandCalibrationTechn0109y，Beijing，100039，China)AbstractAmethodbasedonGUMuncertaintyframeworkisproposedformeasurementuncenaintyevalu—ationofmultiv撕atemeasurementmodels．Theway0fspecifyingthecoverageregionofmultivariatemeasurementmodelsisa180illustrated．KeywOrdsGUMMultivariateUncertaintyCoverageRegion1引言测量不确定度表示指南中推荐的评定方法(简称GuM)‘13是最为常用的测量不确定度评定方法，经过几十年的发展和完善，已经为国际上广泛使用，涉及到标准建立、量值比对、技术文件编制、计量服务等诸多领域。GuM方法所针对的测量模型是单输出量测量模型，即测量模型包含若干实数输入量墨，墨，．．．，x。和一个实数输出量y，表示为：l，=厂皤)(1)式中：卜实数的输入量矢量，x=瞄，，墨，．．．，xⅣ)71。然而在日常评定中，其测量模型可能是包含多个输出量的，GUM方法并不适用于该类测量模型。例如测量反射系数的过程中，校准后的复反射系数r和未修正的复反射系数形之间的关系表示为r=躺(2)式中：n，6，c——复校准系数{2"3|。复数量r可视为双输出量(rR，rI)，使用GuM方法进行不确定度评定时，需要首先根据复数的实部和虚部将该测量模型分解为两个实数模型。收稿日期：2017一03一01，修回日期：2017—06—14作者简介：谌贝(1985．12一)，男，高级工程师，主要研究方向：无线电计量技术。 ·2宇航计测技术若该测量过程最终是为了获得校准系数，那么式(2)转化为隐函数：(c形+1)0一@％+6)=o(3)此时，三个校准系数口，6，c就是测量模型的输出量。当考虑到校准系数为复数时，该模型的输出量个数为6。此时利用GUM方法对每个输出量进行不确定度评定是无法实现的。虽然GUM的相关标准文件中并没有直接考虑多输出量测量模型，但是从GuM的评定原理和不确定度传播律中，并未发现该方法的应用会受到输出量数量的限制。利用矩阵的表示方法，能够将多输出量和多输入量同时纳入至单次计算过程中，使GUM方法能够适用于多输出量测量模型。本文中将对GUM如何应用于多输出量测量模型的不确定度评定进行讨论。2不确定度评定步骤无论针对何种测量模型，使用何种不确定度评定方法，其评定步骤是基本一致的，包括建模、传播和总结三部分。建模部分是不确定度评定中最基础的环节，其中的主要步骤包括定义被测的输出量y(被测量矢量)，确定与】，相关的输入量x，定义与x和l，相关的测量模型，以及根据x中各个量的已知信息来给出各自的概率密度函数(PDF)+。常用的概率密度函数包括高斯(正态)分布、矩形(均匀)分布等，但是当x中的两个量之间不独立时需要给出联合PDF。传播部分是指根据测量模型，由x中每个量的PDF获得l，的(联合)PDF。根据y的PDF，将l，的期望值作为Y的估计值y，将y的协定差矩阵作为y的协方差矩阵嘶，并根据给定的包含概率p获得包含范围。2．1建模1)定义被测的输出量y(被测量矢量)；2)确定与】，相关的输入量x；3)定义与x和l，相关的测量模型；4)根据x中各个量的已知信息来给出各自的概率密度函数(PDF)，例如高斯(正态)分布、矩形(均匀)分布等；或者当x中的两个量之间不独立时给出联合PDF。2．2传播根据测量模型，由x中每个量的PDF获得y的(联合)PDF。2．3总结利用】，的PDF，可以获得：1)y的期望值作为】，的估计值y；2)】，的协方差矩阵作为y的协方差矩阵以；3)根据给定的包含概率p获得y的包含范围。其中，包含范围是指基于可获得的信息确定的包含被测量矢量的一组值的范围，被测量矢量的值以一定概率落在该范围中。该定义与“包含区间”的定义类似，不同的是，包含区间在数学上只是一维的概念，而包含范围不再局限在一维空间中。3多变量的概率分布在JJF1059．2—2012《用蒙特卡洛法评定测量不确定度》H?中，给出了在不确定度评定的建模过程中一些常用的输入量PDF类型，其中涉及的多变量分布只有多变量正态分布。除了多变量正态分布，本文中还会介绍另一种常用的多变量分布——多变量￡分布。，3．1多变量正态分布(多元正态分布)如果仅知Ⅳ维量x=(瓦x扩．．，xⅣ)71的最佳估计值石=(z，，戈扩．．，戈Ⅳ)7’及严格正定不确定度矩阵，那么x的分布为多变量正态分布，表示为g互(国2ii；赫exp(、一吉(孝一肛)7y。1(f一肛)j式中：肚——期望；卜x的正定协方差矩阵。3．2多变量f分布假设有n组示值戈，，z，，．．．，戈。，每组示值都为Ⅳ×1维，其中n>Ⅳ，它们都由多变量高斯分布Ⅳ(肛，渤独立获得，但期望¨和Ⅳ×Ⅳ维协方差∑未知。假设Ⅳ×1维的输人量x等于斗，那么基于联合先验分布和贝叶斯原理，x的边缘(联合)分布为多变量￡分布f，，伍，s／n)，自由度为Ⅳ=n一Ⅳ，其中：牙=丢@。以，．．≯。)，s=孝[@．一司@，一动1+⋯+@。一习伍。一动7]贝0，X的PDF为：gx(9=ii；乏：嚣第×[det(s／功]‘1佗l+吉(f一功7(要)一1(亭一羽一州2其中，r亿)是变量z的伽马函数为：r(名)=F厂1e一‘d￡，z>o 第3期多输出量测量模型的不确定度评定方法此时，X的期望和方差分别为：E(柳=牙，y(柳=老鲁4基于GUM的多输出量测量模型的不确定度评定测量模型根据其表达式形式可分为显函数和隐函数两种，不同形式模型的不确定度评定过程会有一定差异，但原理基本相同。4．1显函数模型如果测量模型中具有多个输出量K，一，．．．，y卅，且测量模型可表示为l，可(柳(4)其中y=(y。，y2，．．．，ym)7，则该模型为显函数形式，，代表了测量函数。任意具体的函数，(X)只与x的某个子集相关，而每个输入量分量咒出现在至少一个函数中。当X的估计值为菇时，则y的估计值为y=／(功。y的协方差矩阵表示为l“◇。，y，)⋯M◇。，)，。)Iu，=l；。．；I“(y。，y。)⋯M◇。，，，。)J其中，协方差cov(yj，yi)=u2◇，)，并有U=c，u。c：(5)式中：e——x=戈时的m×Ⅳ维灵敏系数矩阵。灵敏系数矩阵表示为Cx=嘶aXl弧aX，(6)以式(2)所表示的测量过程对本方法进行举例说明。当使用上文的符号进行相关表示时，有Ⅳ-8，m=2，以及x=@R，o，，6R，6，，cR，c，，％，肜∥】，=叮R，r∥其中的下标R和1分别表示实部和虚部。那么反射系数的估计值y=矿。，ry可根据式(2)计算获得。当X：戈，以为z的8×8维协方差矩阵时，y的2×2维的协方差矩阵以可由式(5)计算，其中2×8维的灵敏系数e由下式给出arRa11RarRar月arRar—a口Rao，a6Ra6，acR，ac，a％a盯l『ar，ar，ar，la口Ra口，a6Ra6，acRac，a％a啊l代人具体数据后，就能够获得所有输出量的不确定度。4．2隐函数模型某些情况下，测量模型只能以隐函数形式表示^(y，X)=0(7)矗=(危，，^：，⋯，^。)7代表了若干测量表达式。当x的估计值为戈，y的估计值为y时，有^◇，功=0(8)通常式(8)的解很难直接获得，需要利用例如牛顿法‘51或牛顿法的衍生方法，通过迭代获得最终解。y的m×m维协方差矩阵以可由式(9)获得c，u，c：=c。虬c：(9)c’是m×m维灵敏系数矩阵，其中包括偏导数a^f／ar，z=1，2，⋯，m，．i=1，2，⋯，m；e是m×Ⅳ维灵敏系数矩阵，其中包括偏导数a危，／axi，Z=1，2，⋯，m，江1，2，⋯，Ⅳ。所有导数都在x=戈和】，=y时获得。为了形式上与式(4)相比较，式(9)可进行变形，得到u，=c以C7(10)式中：C=c：1c。——m×Ⅳ维的灵敏系数矩阵。将式(3)的例子用上文的符号进行表示。因为要获得三个复校准系数的解必须具有至少三次测量过程，所以有』V=12，m=6，以及x=孵印％．『，r1∽r1，，，畈一畈，，，r2一r2，，"耽一职_，，F3∽r3．，)』y=@R，o，，6R，6，，cR，c∥令危：H∞x)=O和^：，∞柳=O分别代表式(3)中涉及实部或虚部计算的隐函数形式方程，其中有，=1，2，3。将由此得到的6个方程联立，形和D的估计值代人后计算得到校准系数的估计值y=@R，a，，6R，6，，6R，6，)』。在x=戈，l，=y时，根据式(9)计算得到y的6×6维协方差矩阵以，其中G是包含跳?／aK的6×6维灵敏系数矩阵，G是包含娩，／axi的6×12维灵敏系数矩阵，f_1，2，⋯，6，尼=1，2，⋯，6，扣1，2，⋯，12，玑是戈的12×12维协方差矩阵。在本示例中，虽然有12个输入量，但是由测量鞔一呱；吮一％ 宇航计测技术2017年过程可以分析得出，6个方程的任意一个只涉及4个输入量形∽形，，，1．。，1．，，这会使原计算变得更为简单。5矢量输出量的包含范围包含范围是不确定度的重要评价指标。对多输出量测量模型而言，输冉量y=(y，，y2，．，ym)7。的估计值为y，其协方差矩阵为以。当要求的包含概率为p时，则需要在m维空间中确定一个包含范罔R，，使l，以概率p落在该范围中。通常，当y的PDF确定后，无论是某个具体包含范阿的包含概率，还是具有某个包含概率的包含范周，都是已经确定的。当测量模型中包含多输出量时，对其直观地表示将很闲难。而双输出量测量模型是最简单的多输出量测量模型，足以说明多输出量的包含范同与单输出量的包含区间之间的区别。本文将利用双输出量测量模型对包含范围进行阐述。假设如下情况：在直角坐标系中，】，=(y．，y")71分别对应坐标轴K(横坐标)和y2(纵坐标)。y的已知信息包括估计值y．和y2，对应的标准不确定度u(y．)和u(y：)，以及由相同仪器获得估计值所产生的协方差u(y，，y2)。在GuM不确定度框架中，由于没有其它信息，表征l，的联合PDFg，，(叼。，叩：)为双变量高斯PDFⅣ(y，U、)，其中：y=眵l以2I菇2，警爱翻当y，=yz=0，M2(y，)=2．0，u2@：)=1．0，“◇-，yz)=M眵z，y。)=0．0时，在该分布下进行1000次随机抽样，如图1(a)所示。网中矩形表示的是95％包含概率的包含范围，面积27．9。矩形包含范围是通过分别计算每个输出最的边缘PDF获得的gK("71)=J一。gH，y，(77I，卵2)d叩2gH("72)=J-。gkr，("7l，叼2)d771当y-=y：=0，“2(yt)=u2(y：)=2．0，“(y-，yz)=u(弛，)，．)=1．9时，1000次随机抽样如图1(b)所示，95％包含概率的矩形的包含范围，面积35．6。需要指出，矩形包含范同只是包含范围众多表示形式巾的一种。另一种常用的包含范围是椭网形式。6J：(叩一y)～、“∽一y)=七：式中：矗，——已知的常数；p——取该PDF时椭网范围中的概率。该形式能够更好地表现输出量之间的关系。例如对于图1(b)中的随机分布来说，两个输出量是相关的，当矩形的边长与坐标轴平行时，其表示的包含范同中大部分区域中并无数据，无法很好地反映数据分布，而相同包含概率下，椭圆形式的包含范围会小得多，所以更加合理。(a)双输¨5量不相关时(b)双输量相关时图1双输出量的测量模型随机抽样后的包含范围(包含概率95％)若输出量多于两个，那么该测虽模型很难通过直观的图形方式进行分析，但寻找包含范围的原理及计算方法与双输出量模型中是相同的，只不过将其中的矩形、椭圆等二维图形换为超矩形、超椭球体等多维体。(下转第24页) 宇航计测技术一50．OO崮100．0jj酣坚一150．O图12优化后的环路相位噪声测试结果这是电路优化后的PLL相位噪声测试电路仿真结果。从优化后的仿真结果可以看到，系统频率在900MHz基本稳定，即达到了锁定状态，锁定时间为252．5“s，相位噪声相对减小，由集中变得分散，且偏离了lokHz，减小了对系统锁定状态的影响，使电路达到了设计目标的要求。5结束语通过插人对电路的优化目标，设置额外的仿真器，对电路元器件的下层结构运行情况进行监视与调整，再次进行仿真，得到的结果与改进前相比有了很明显的变化，电路震荡特性减小，很快达到了锁定状态，锁定时间缩短，且相位噪声减小并分散，偏离目标合成频率，减少杂散频率，提高了电路的稳定性和频率合成效果。参考文献[1]ChenakinA．FrequencySynthesizers：ConcepttoProd-uct[M]．Boston：ArtechHouse，2011．[2]李伟，凌伟，白轶荣，等．一种频率合成模块的设计和实现[J]．微波学报，2014，6：177～181．[3]李智鹏，郭伟，刘永智，鲍景富．自偏置锁相环的相位噪声分析[J]．微波学报，2015，lO：59～63．[4]徐兴福_ADs2011射频电路设计与仿真实例[E]．北京：电子工业出版社，2叭4，(05)．[5]魏福立，张华鹏．具有快速锁定时间的新型低相噪频率合成技术[J]．半导体技术，2006，5：321～324．[6]牛吉凌，康弘俊．一种基于ADs仿真的射频放大器设计方法[J]．通信技术，2015，4：489—494．(上接第4页)6结束语根据GUM不确定度评定方法，提出了针对显函数和隐函数类型的多输出量测量模型的测量不确定度评定方法，并利用双输出量测量模型，举例说明了如何获得给定概率的矩形或椭网包含范围。利用本论文的结果，能够将GUM方法的应用拓展到多输出量测量模型领域，为测量不确定度的评定提供新的思路。参考文献[1]国家质量监督检验检疫总局．JJFL059．卜2012，测量不确定度评定与表示[s]．北京：中国质检出版社，2012．12JEngenGF．MicrowavecircuittheoryandfoundationsofmicrowavemetrologylMj．London：PeterPeregrinus，1992．[3]KernsDM，BeattyRW．BasictheoryofwaveguideJunc【ionsandintroductorymicrowavenetworkanalysis[MJ．London：PergamonPress，1967．[4]国家质量监督检验检疫总局．JJF1059．2—2叭2，用蒙特卡洛法评定测量不确定度[S]．北京：中国质检出版社，2012．15JGillPE，Murrayw，Wr培htMH．PracticalOptimiza-tionlMJ．London：AcademicPress，1981．[6]MardiaKV，KentJT，BibbyJM．MultivariateAnaly—sis『M]．London：AcademicPress，1979．加∞∞加∞¨B¨"'