一种工频电场校准装置测量不确定度的评定方法及系统与流程

1.本发明涉及电学参数测试技术领域，更具体地，涉及一种工频电场校准装置测量不确定度的评定方法及系统。


背景技术：

2.工频电场参数作为输变电工程重要的环境影响因子，在输变电工程建设运行前后的工程环境影响评价中不可或缺。其测量结果直接应用于导线选型，端点选址，路径规划等工程设计领域，同时也是缓解公众焦虑和环保纠纷评判的重要法律依据，工频电场测量结果的准确与否，会直接对工程造价控制的难易程度，环境生态的友好水平，人民群众切身权益的保障质量产生影响。同时在电力行业当中，工频电场测量作为研究分析电场的一种十分重要的技术手段，也在电气电缆故障定位、电气绝缘性能优化、气象灾害预警，电场生物效应等方面有着诸多运用。
3.而对工频电场参数的准确校准是上述技术得以实际运用的前提条件，工频电场参数的测量结果质量，很大程度由溯源链上计量学等级更高的校准系统水平高低直接决定。合理评定的工频电场校准装置不确定度是判断工频电场参数的测量结果质量的重要参考依据，对科学、有效、正确的利用测量结果具有重要的科研意义与工程价值。
4.自1993年iso第四技术咨询工作组(tag4)中的测量不确定度表示工作组(wg3)以7个国际组织的名义发布《测量不确定度第3部分：测量不确定度表示指南》(uncertaintyofmeasurement
‑
part3:guidetotheexpressionofuncertaintyinmeasurement，1993)以来，经过多年的研究与发展，各个国家和一部分专业学科的国际组织在指南的通用规则基础上，先后发布了相应的测量不确定度评定指南文件。在测量不确定度评定领域，国际上已形成了较为统一规范的评定方法共识。
5.我国由中国计量科学研究院，中国法制计量管理计量技术委员会等机构组织，依据《测量不确定度第3部分：测量不确定度表示指南》(uncertaintyofmeasurement
‑
part3:guidetotheexpressionofuncertaintyinmeasurement)，同样先后制定发布了一系列测量不确定度指南和技术规范文件，包括《测量不确定度规范》，jjf1059
‑
2012《测量不确定度评定与表示》等。这些指南和技术规范文件构成了我国测量不确定度评定的框架，也使得我国在计量、测试领域对测量结果不确定度的表述与国际接轨。
6.现有技术中，国内外在测量结果的不确定度评定上主要有gum法和蒙特卡洛法(mcm)，其中gum法基于中心极限定理，通过应用不确定度传播律传播输入量的最佳估计值和其标准不确定度，最终获得输出量的最佳估计值和其测量不确定度；mcm法利用对概率分布进行随机抽样，以分布传播的方式，最终确定输出量的概率密度函数(pdf)。目前国内外公开报道中，通常使用gum法对工频电场校准装置测量不确定度的进行评定，未见使用蒙特卡洛法评定的应用方法与实例。
7.因此，需要一种技术，以基于蒙特卡洛法对工频电场校准装置测量不确定度进行评定。


技术实现要素：

8.本发明技术方案提供一种工频电场校准装置测量不确定度的评定方法及系统，以解决如何基于蒙特卡洛法工频电场校准装置测量不确定度的评定问题。
9.为了解决上述问题，本发明提供了一种工频电场校准装置测量不确定度的评定方法，所述方法包括：
10.基于预先建立的测量模型，确定工频电场校准装置产生工频电场名义值的测量函数；
11.为所述测量函数的每一个输入量设定概率密度函数，并根据所述概率密度函数，确定每一个输入量的概率分布；
12.根据每一个输入量的概率分布，基于蒙特卡洛抽样生成最终输出量的离散分布样本，基于输出量的离散分布样本，确定蒙特卡洛抽样次数m；
13.基于所述蒙特卡洛抽样次数m，进行蒙特卡洛抽样，获取所述测量函数m个输出量；
14.对所述测量函数m个输出量进行非递减排序，获取排序后的测量函数m个输出量，并基于排序后的测量函数m个输出量，获取排序后的测量函数m个输出量的分布函数的离散表达；
15.基于所述分布函数的离散表达，计算预设概率下的m个输出量的包含区间。
16.优选地，所述测量函数为：
[0017][0018]
d＝d
t
[1+2α(t
‑
t0)]，
[0019]
其中，e为工频电场校准装置产生的工频电场名义值；n为电压互感器变比；v为电压表示值；e0为背景电场；d为实验室温度下的平行极板间距；d
t
为溯源温度下的平行极板间距；rep为测量重复性；t为实验室温度；t0为溯源时的温度；α为平行极板材质的线膨胀系数。
[0020]
优选地，所述确定蒙特卡洛抽样次数m，包括：
[0021]
确定m1的值大于等于1/(1
‑
p)的104倍：
[0022]
m1≥1/(1
‑
p)
×
104，
[0023]
其中，m1为在保证不确定度提供期望的包含区间时的最小抽样次数，p为期望的包含概率；
[0024]
确定m2的值大于等于为不确定度保留足够的有效数字而所需的自由度：
[0025][0026]
其中，m2为在保证不确定度保留足够有效数字时的最小抽样次数，v为所需的自由度，u(y)为标准不确定度，σ[u(y)]为标准不确定度u(y)的标准差；
[0027]
确定m取m1和m2中的较大值：
[0028]
m≥max(m1,m2)。
[0029]
优选地，所述基于所述蒙特卡洛抽样次数m，进行蒙特卡洛抽样，包括：
[0030]
将所述测量模型记为y＝f(x1,
…
x
n
)，从n个输入量的x
i
的概率密度函数gx
i
(ξ
i
)中
抽取m个向量x
r
＝(x
1,r
,
…
x
n,r
)，通过计算获得m个输出量计算获得m个输出量y
r
＝f(x
r
)，r＝1，
…
，m，进一步可得简化后测量模型y的估计值和标准不确定u(y)，输入矩阵x
r
，输出矩阵y
r
分别表示：
[0031][0032][0033]
其中，(x1,
…
x
n
)为测量模型中的n个输入量，ξ
i
为描述输入量x
i
可能值的变量，(x
1,r
,
…
x
n,r
)为输入量的抽取样本，(y1,
…
,y
m
)。
[0034]
优选地，所述预设概率为50％、90％、95％、99％或99.73％。
[0035]
根据本发明的另一方面，本发明提供一种工频电场校准装置测量不确定度的评定系统，所述系统包括：
[0036]
初始单元，用于基于预先建立的测量模型，确定工频电场校准装置产生工频电场名义值的测量函数；
[0037]
第一确定单元，用于为所述测量函数的每一个输入量设定概率密度函数，并根据所述概率密度函数，确定每一个输入量的概率分布；
[0038]
第二确定单元，用于根据每一个输入量的概率分布，基于蒙特卡洛抽样生成最终输出量的离散分布样本，基于输出量的离散分布样本，确定蒙特卡洛抽样次数m；
[0039]
第一获取单元，用于基于所述蒙特卡洛抽样次数m，进行蒙特卡洛抽样，获取所述测量函数m个输出量；
[0040]
第二获取单元，用于对所述测量函数m个输出量进行非递减排序，获取排序后的测量函数m个输出量，并基于排序后的测量函数m个输出量，获取排序后的测量函数m个输出量的分布函数的离散表达；
[0041]
计算单元，用于基于所述分布函数的离散表达，计算预设概率下的m个输出量的包含区间。
[0042]
优选地，所述测量函数为：
[0043][0044]
d＝d
t
[1+2α(t
‑
t0)]，
[0045]
其中，e为工频电场校准装置产生的工频电场名义值；n为电压互感器变比；v为电压表示值；e0为背景电场；d为实验室温度下的平行极板间距；d
t
为溯源温度下的平行极板间距；rep为测量重复性；t为实验室温度；t0为溯源时的温度；α为平行极板材质的线膨胀系数。
[0046]
优选地，所述确定蒙特卡洛抽样次数m，包括：
[0047]
确定m1的值大于等于1/(1
‑
p)的104倍：
[0048]
m1≥1/(1
‑
p)
×
104，
[0049]
其中，m1为在保证不确定度提供期望的包含区间时的最小抽样次数，p为期望的包含概率；
[0050]
确定m2的值大于等于为不确定度保留足够的有效数字而所需的自由度：
[0051][0052]
其中，m2为在保证不确定度保留足够有效数字时的最小抽样次数，v为所需的自由度，u(y)为标准不确定度，σ[u(y)]为标准不确定度u(y)的标准差；
[0053]
确定m取m1和m2中的较大值：
[0054]
m≥max(m1,m2)。
[0055]
优选地，所述
[0056]
基于所述蒙特卡洛抽样次数m，进行蒙特卡洛抽样，包括：
[0057]
将所述测量模型记为y＝f(x1,
…
x
n
)，从n个输入量的x
i
的概率密度函数gx
i
(ξ
i
)中抽取m个向量x
r
＝(x
1,r
,
…
x
n,r
)，通过计算获得m个输出量计算获得m个输出量y
r
＝f(x
r
)，r＝1，
…
，m，进一步可得简化后测量模型y的估计值和标准不确定u(y)，输入矩阵x
r
，输出矩阵y
r
分别表示：
[0058][0059][0060]
其中，(x1,
…
x
n
)为测量模型中的n个输入量，ξ
i
为描述输入量x
i
可能值的变量，(x
1,r
,
…
x
n,r
)为输入量的抽取样本，(y1,
…
,y
m
)。
[0061]
优选地，所述预设概率为50％、90％、95％、99％或99.73％。
[0062]
本发明技术方案提供一种工频电场校准装置测量不确定度的评定方法及系统，根
据工频电场校准的实施过程建立测量模型，基于可用信息为输入量设定概率密度函数，再对输入量的概率密度进行离散抽样，并通过不确定度评定的测量模型传播输入量的分布，获得输出量分布函数的离散表达；基于分布函数的离散表达，计算预设概率下的m个输出量的包含区间。本发明技术方案提出了一种基于蒙特卡洛法工频电场校准装置测量不确定度的评定方法，在精度，可靠性和计算效率上相较于现有方法更具优势，能够产生一定的理论意义和社会经济效益。
附图说明
[0063]
通过参考下面的附图，可以更为完整地理解本发明的示例性实施方式：
[0064]
图1为根据本发明优选实施方式的一种工频电场校准装置测量不确定度的评定方法100的流程图；
[0065]
图2为根据本发明优选实施方式的工频电场校准装置结构示意图；
[0066]
图3为根据本发明优选实施方式的评定流程图；
[0067]
图4为根据本发明优选实施方式的输入量的概率分布传播示意图；
[0068]
图5为根据本发明优选实施方式的各输入量的概率分布直方图；
[0069]
图6为根据本发明优选实施方式的输出量的概率分布直方图；
[0070]
图7为根据本发明优选实施方式的95％包含概率下的输出量的包含区间示意图；
[0071]
图8为根据本发明优选实施方式的一种工频电场校准装置测量不确定度的评定系统800的结构示意图。
具体实施方式
[0072]
现在参考附图介绍本发明的示例性实施方式，然而，本发明可以用许多不同的形式来实施，并且不局限于此处描述的实施例，提供这些实施例是为了详尽地且完全地公开本发明，并且向所属技术领域的技术人员充分传达本发明的范围。对于表示在附图中的示例性实施方式中的术语并不是对本发明的限定。在附图中，相同的单元/元件使用相同的附图标记。
[0073]
除非另有说明，此处使用的术语(包括科技术语)对所属技术领域的技术人员具有通常的理解含义。另外，可以理解的是，以通常使用的词典限定的术语，应当被理解为与其相关领域的语境具有一致的含义，而不应该被理解为理想化的或过于正式的意义。
[0074]
图1为根据本发明优选实施方式的一种工频电场校准装置测量不确定度的评定方法100的流程图。本发明涉及工频电场参数的计量水平评估，具体涉及一种基于蒙特卡洛法的工频电场校准装置测量不确定度的评定方法。本发明提出了一种基于蒙特卡洛法工频电场校准装置测量不确定度的评定方法，方法根据工频电场校准的实施过程建立测量模型，基于可用信息为输入量设定概率密度函数，再对输入量的概率密度进行离散抽样，并通过不确定度评定的测量模型传播输入量的分布，获得输出量分布函数的离散表达及其特征值，在精度，可靠性和计算效率上相较于现有方法更具优势。
[0075]
如图1所示，本发明提供的工频电场校准装置测量不确定度的评定方法100，包括：
[0076]
步骤101：基于预先建立的测量模型，确定工频电场校准装置产生工频电场名义值的测量函数。
[0077]
优选地，测量函数为：
[0078][0079]
d＝d
t
[1+2α(t
‑
t0)]，
[0080]
其中，e为工频电场校准装置产生的工频电场名义值；n为电压互感器变比；v为电压表示值；e0为背景电场；d为实验室温度下的平行极板间距；d
t
为溯源温度下的平行极板间距；rep为测量重复性；t为实验室温度；t0为溯源时的温度；α为平行极板材质的线膨胀系数。
[0081]
本发明首先建立测量模型。在工频电场校准的实施过程中，将工频电场探头放置在平行极板结构的几何中心，对上极板施加一定电压，可以在平行极板间产生相应的工频电场场强水平。同时电压信号经由电压互感器传输至电压表，通过记录电压数值，从而可以计算获得校准装置产生的工频电场名义值。其中，工频电场校准的系统如图2所示。
[0082]
可知工频电场场强水平取决于输入电压，电压互感器变比，和平行极板间距
[0083][0084]
其中，e为校准系统产生的工频电场名义值，单位为v/m；v为电压表示值，单位为v；n为电压互感器变比，无量纲；d为平行极板间距，单位为m。
[0085]
由于实验环境背景电场的存在，式(1)变为
[0086][0087]
其中e0为背景电场，单位为v/m。
[0088]
引入校准过程的测量重复性，式(2)变为
[0089][0090]
其中rep为测量重复性，其估计值rep＝1.0。
[0091]
测量函数中的部分输入量本身又是其它量的函数，应表示成由基本量所组成。进一步的,平行极板间距d是与实验室环境温度相关的函数：
[0092]
d＝d
t
[1+2α(t
‑
t0)]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
[0093]
其中，t为实验室温度，单位为℃；t0为溯源时的温度，单位为℃；α为平行极板材质的线膨胀系数，单位为℃
‑1。
[0094]
最终工频电场校准装置的测量函数变为：
[0095][0096]
步骤102：为测量函数的每一个输入量设定概率密度函数，根据概率密度函数，确定每一个输入量的概率分布。
[0097]
本发明设定输入量的概率密度函数。为了上述测量函数的每一个输入量设定一个适当的概率密度函数，概率密度函数基于可用信息推导得到。输入量包括电压互感器变比、
电压表示值、平行极板间距、平行极板材质的线膨胀系数、实验室温度、溯源时的温度、背景电场、测量重复性。
[0098]
步骤103：根据每一个输入量的概率分布，基于蒙特卡洛抽样生成最终输出量的离散分布样本，基于输出量的离散分布样本，确定蒙特卡洛抽样次数m。
[0099]
优选地，确定蒙特卡洛抽样次数m，包括：
[0100]
确定m1的值大于等于1/(1
‑
p)的104倍：
[0101]
m1≥1/(1
‑
p)
×
104，
[0102]
其中，m1为在保证不确定度提供期望的包含区间时的最小抽样次数，p为期望的包含概率；
[0103]
确定m2的值大于等于为不确定度保留足够的有效数字而所需的自由度：
[0104][0105]
其中，m2为在保证不确定度保留足够有效数字时的最小抽样次数，v为所需的自由度，u(y)为标准不确定度，σ[u(y)]为标准不确定度u(y)的标准差；
[0106]
确定m取m1和m2中的较大值：
[0107]
m≥max(m1,m2)。
[0108]
本发明确定蒙特卡洛抽样次数，依据第二步骤为输入量设定的概率分布，使用蒙特卡洛抽样可以生成最终输出量的离散分布样本，该分布样本与实际样本分布的接近程度取决于蒙特卡洛抽样次数m，也即样本容量m。抽样次数越多，样本容量越大，统计规律越接近输出量的真实情况，但所需的计算时间也越长。抽样次数越少，样本容量越小，统计规律越远离输出量的真实情况。因此m的取值对于工频电场的测量不确定度评定十分关键。
[0109]
确定的具体步骤如下：
[0110]
1)m1的值应大于等于1/(1
‑
p)的104倍，即：
[0111]
m1≥1/(1
‑
p)
×
104ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0112]
2)m2的值应大于等于为不确定度保留足够的有效数字而所需的自由度，即：
[0113][0114]
3)m至少取m1和m2中的较大值，即：
[0115]
m≥max(m1,m2)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
[0116]
步骤104：基于蒙特卡洛抽样次数m，进行蒙特卡洛抽样，获取测量函数m个输出量。
[0117]
优选地，基于所述蒙特卡洛抽样次数m，进行蒙特卡洛抽样，包括：
[0118]
将所述测量模型记为y＝f(x1,
…
x
n
)，从n个输入量的x
i
的概率密度函数gx
i
(ξ
i
)中抽取m个向量x
r
＝(x
1,r
,
…
x
n,r
)，通过计算获得m个输出量计算获得m个输出量y
r
＝f(x
r
)，r＝1，
…
，m，进一步可得简化后测量模型y的估计值和标准不确定u(y)，输入矩阵x
r
，输出矩阵y
r
分别表示：
[0119][0120][0121]
其中，(x1,
…
x
n
)为测量模型中的n个输入量，ξ
i
为描述输入量x
i
可能值的变量，(x
1,r
,
…
x
n,r
)为输入量的抽取样本，(y1,
…
,y
m
)。
[0122]
第四步骤：进行蒙特卡洛抽样
[0123]
将简化的测量模型记为y＝f(x1,
…
x
n
),从n个输出量的x
i
的概率密度函数gx
i
(ξ
i
)中抽取m个向量x
r
＝(x
1,r
,
…
x
n,r
)，通过计算获得m个输出量y
r
＝f(x
r
)，r＝1，
…
，m，进一步可得测量模型y的估计值和标准不确定u(y)。输入矩阵x
r
，输出矩阵y
r
分别表示为式(9)和式(10)：
[0124][0125][0126]
步骤105：对所述测量函数m个输出量进行非递减排序，获取排序后的测量函数m个输出量，并基于排序后的测量函数m个输出量，获取排序后的测量函数m个输出量的分布函数的离散表达。
[0127]
本发明对输出量做非递减排序，将第四步骤中的y
r
＝f(x
r
)，r＝1，
…
，m按照非递减顺序进行排序，得到排序后的输出量y
(r)
，r＝1，
…
，m，进一步可得输出量的分布函数g
y
(η)的离散表达g。
[0128]
步骤106：基于分布函数的离散表达，计算预设概率下的m个输出量的包含区间。
[0129]
优选地，预设概率为50％、90％、95％、99％或99.73％。
[0130]
本发明确定包含区间，应用离散表达g计算约定包含概率下的输出量y的包含区间如图3所示。
[0131]
本发明有效解决了工频电场校准装置测量结果因影响因素较多，导致测量模型规律性差、非线性显著，进而在使用gum法评定测量不确定度时，结果不可靠的问题。
[0132]
本发明有效解决了现有方法在面对工频电场校准装置非线性模型时，求导计算困难的问题，同时不存在模型函数的近似转换，实现了更为精准的测量不确定评定。
[0133]
本发明计算处理过程明晰，避免了现有方法中的主观假设条件的引入，便于使用计算机语言执行，为工频电场校准装置实现高精度高效率高可靠的测量不确定度评定工作提供了一种新方法和一定的有益参考。
[0134]
在本发明实施例中，本发明提供的一种基于蒙特卡洛法工频电场校准装置测量不确定度的评定方法，工频电场校准在t＝23.0℃下进行，校准工频电场点位为1000v/m，包括以下步骤：
[0135]
第一步骤：建立测量模型
[0136]
工频电场校准装置的测量函数由式(5)给出：
[0137][0138]
第二步骤：设定输入量的概率密度函数
[0139]
基于可用信息推导为上述测量函数的每一个输入量设定一个适当的概率密度函数，输入量有电压互感器变比n、电压表示值v、平行极板间距d
t
、平行极板材质的线膨胀系数α、实验室温度t、溯源时的温度t0、背景电场e0、测量重复性rep，表1将列出这些信息。
[0140]
表1工频电场校准不确定度设定输入量的概率密度函数
[0141]
[0142]
第三步骤：确定蒙特卡洛抽样次数
[0143]
这里取p＝0.95，即以95％的概率包含实际样本分布，由下式计算可得：
[0144]
m1≥1/(1
‑
p)
×
104＝2.0
×
105，
[0145]
选择为包含区间长度保留至2位有效数字，由下式计算可得：
[0146][0147]
蒙特卡洛抽样次数m至少取m1和m2中的较大值，为了便于描述以及为抽样结果留有一定的冗余，设定m＝1.0
×
106。
[0148]
第四步骤：进行蒙特卡洛抽样
[0149]
依据电压互感器变比n、电压表示值v、平行极板间距d
t
、平行极板材质的线膨胀系数α、实验室温度t、溯源时的温度t0、背景电场e0、测量重复性rep的概率密度函数，分别进行m＝1.0
×
106次离散抽样后，各输入量的概率分布直方图如图4所示，同时获得m＝1.0
×
106个向量的输入矩阵和输出矩阵。
[0150][0151][0152]
进一步计算得的测量模型y的估计值和标准不确定u(y)，
[0153][0154]
u(y)＝22.82v/m。
[0155]
第五步骤：对输出量做非递减排序
[0156]
将第四步骤中的y
r
按照非递减顺序进行排序，得到排序后的输出量y
(r)
概率分布，进一步可得输出量的分布函数g
y
(η)的离散表达g，如图6所示。
[0157]
第六步骤：确定包含区间
[0158]
应用离散表达g计算在95％包含概率下的输出量y的包含区间如图7所示。
[0159]
图8为根据本发明优选实施方式的一种工频电场校准装置测量不确定度的评定系统结构图。如图8所示，本发明提供一种工频电场校准装置测量不确定度的评定系统，系统包括：
[0160]
初始单元801，用于基于预先建立的测量模型，确定工频电场校准装置产生工频电场名义值的测量函数。
[0161]
优选地，测量函数为：
[0162][0163]
d＝d
t
[1+2α(t
‑
t0)]，
[0164]
其中，e为工频电场校准装置产生的工频电场名义值；n为电压互感器变比；v为电压表示值；e0为背景电场；d为实验室温度下的平行极板间距；d
t
为溯源温度下的平行极板间距；rep为测量重复性；t为实验室温度；t0为溯源时的温度；α为平行极板材质的线膨胀系数。
[0165]
第一确定单元802，用于为测量函数的每一个输入量设定概率密度函数，根据概率密度函数，确定每一个输入量的概率分布。
[0166]
第二确定单元803，用于根据每一个输入量的概率分布，基于蒙特卡洛抽样生成最终输出量的离散分布样本，基于输出量的离散分布样本，确定蒙特卡洛抽样次数m。
[0167]
优选地，确定蒙特卡洛抽样次数m，包括：
[0168]
确定m1的值大于等于1/(1
‑
p)的104倍：
[0169]
m1≥1/(1
‑
p)
×
104，
[0170]
其中，m1为在保证不确定度提供期望的包含区间时的最小抽样次数，p为期望的包含概率；
[0171]
确定m2的值大于等于为不确定度保留足够的有效数字而所需的自由度：
[0172][0173]
其中，m2为在保证不确定度保留足够有效数字时的最小抽样次数，v为所需的自由度，u(y)为标准不确定度，σ[u(y)]为标准不确定度u(y)的标准差；
[0174]
确定m取m1和m2中的较大值：
[0175]
m≥max(m1,m2)。
[0176]
第一获取单元804，用于基于蒙特卡洛抽样次数m，进行蒙特卡洛抽样，获取测量函数m个输出量。
[0177]
优选地，
[0178]
基于所述蒙特卡洛抽样次数m，进行蒙特卡洛抽样，包括：
[0179]
将所述测量模型记为y＝f(x1,
…
x
n
)，从n个输入量的x
i
的概率密度函数gx
i
(ξ
i
)中抽取m个向量x
r
＝(x
1,r
,
…
x
n,r
)，通过计算获得m个输出量y
r
＝f(x
r
)，r＝1，
…
，m，进一步可得简化后测量模型y的估计值和标准不确定u(y)，输入矩阵x
r
，输出矩阵y
r
分别表示：
[0180]
[0181][0182]
其中，(x1,
…
x
n
)为测量模型中的n个输入量，ξ
i
为描述输入量x
i
可能值的变量，(x
1,r
,
…
x
n,r
)为输入量的抽取样本，(y1,
…
,y
m
)。
[0183]
第二获取单元805，用于对所述测量函数m个输出量进行非递减排序，获取排序后的测量函数m个输出量，并基于排序后的测量函数m个输出量，获取排序后的测量函数m个输出量的分布函数的离散表达。
[0184]
计算单元806，用于基于分布函数的离散表达，计算预设概率下的m个输出量的包含区间。
[0185]
优选地，预设概率为50％、90％、95％、99％或99.73％。
[0186]
本发明优选实施方式的一种工频电场校准装置测量不确定度的评定系统800与本发明优选实施方式的一种工频电场校准装置测量不确定度的评定方法100相对应，在此不再进行赘述。
[0187]
已经通过参考少量实施方式描述了本发明。然而，本领域技术人员所公知的，正如附带的专利权利要求所限定的，除了本发明以上公开的其他的实施例等同地落在本发明的范围内。
[0188]
通常地，在权利要求中使用的所有术语都根据他们在技术领域的通常含义被解释，除非在其中被另外明确地定义。所有的参考“一个/所述/该[装置、组件等]”都被开放地解释为所述装置、组件等中的至少一个实例，除非另外明确地说明。这里公开的任何方法的步骤都没必要以公开的准确的顺序运行，除非明确地说明。