CN103063490A - 一种制备标准样品的方法 - Google Patents
一种制备标准样品的方法 Download PDFInfo
- Publication number
- CN103063490A CN103063490A CN2012105744339A CN201210574433A CN103063490A CN 103063490 A CN103063490 A CN 103063490A CN 2012105744339 A CN2012105744339 A CN 2012105744339A CN 201210574433 A CN201210574433 A CN 201210574433A CN 103063490 A CN103063490 A CN 103063490A
- Authority
- CN
- China
- Prior art keywords
- value
- data
- group
- standard
- unusual
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明提供了一种制备标准样品的方法。所述方法通过依次进行的均匀性检测、夏皮罗一威尔克正态数据检验、柯克仑检验、格拉布斯检验、狄克逊检验等步骤来确定准标准样品的待确定参数的标准值,从而根据标准值确定标准样品,并能够给出不确定度值。本发明的标准样品制备方法能够同时进行均匀性检验和定值分析,而且使用方便、灵活。
Description
技术领域
本发明涉及标准样品制备技术领域,更具体的讲,涉及一种制备标准样品(尤其是,冶金标准样品)的方法。
背景技术
通常,标准样品的研制必需遵照相关国家标准的要求,其中均匀性检验、定值分析涉及到非常复杂的结果计算,往往有数千个数据,较为复杂。
发明内容
本发明的目的在于解决上述现有技术存在的不足。
例如,本发明的目的之一在于提供一种能够同时进行均匀性检验和定值分析的标准样品制备方法。
本发明提供了一种制备标准样品的方法,所述标准样品具有待定值的参数,所述方法通过依次进行的以下步骤来确定所述参数的标准值:对m个准标准样品的待定值参数均进行n次测量,以形成m组数据,每组数据中均有n个数值;根据GB/T15000.3采用方差值Fa对上述m×n个数值进行统计检验,以确定所述m×n个数值是否满足均匀性要求;在所述m×n个数值满足均匀性要求的情况下,对所述m×n个数值进行夏皮罗一威尔克正态数据检验,当正态性检验统计量W大于临界值W0.05时,认为所述m×n个数值呈正态分布;在所述m×n个数值呈正态分布的情况下,对所述m组数据中的每组数据进行柯克仑检验,当该组数据的统计量C小于临界值C0.05,表明该组数据非异常;对所述m组数据中的经柯克仑检验非异常的数据组进行格拉布斯检验,当该组数据的统计量G小于临界值G0.05,表明该组数据非异常;对所述m组数据中的经格拉布斯检验非异常的数据组进行狄克逊检验,当该组数据的统计量D小于临界值D0.05,表明该组数据非异常;对所述m组数据中的经狄克逊检验非异常的数据组进行几何平均,得出待定值参数的标准值,将具有该标准值的样品作为标准样品;对所述m组数据中的经狄克逊检验非异常的数据组进行不确定度计算,给出不确定度值。
在本发明的一个示例性实施例中,所述狄克逊检验步骤的经格拉布斯检验非异常的数据组为8至10组。
在本发明的一个示例性实施例中,所述标准样品具有多个待定值的参数,并且对所述多个待定值的参数中的每个均确定其标准值。
在本发明的一个示例性实施例中,所述标准样品为冶金标准样品。
与现有技术相比,本发明的标准样品制备方法能够同时进行均匀性检验和定值分析,而且使用方便、灵活。
附图说明
通过下面结合附图进行的描述,本发明的上述和其他目的及特点将会变得更加清楚,其中:
图1示出了根据本发明的一个示例性实施例的制备标准样品的方法的流程示意图;
图2示出了根据本发明的一个示例性实施例的制备标准样品的方法对钛精矿中的S进行均匀性检验的示例;
图3示出了根据本发明的一个示例性实施例的制备标准样品的方法对不同单位的S元素进行定值分析检验的示例。
具体实施方式
在下文中,将参照附图来详细说明本发明的示例性实施例。
图1示出了根据本发明的一个示例性实施例的制备标准样品的方法的流程示意图。
如图1所示,在本发明的一个示例性实施例中,制备标准样品的方法通过从准标准样品的待确定参数的测量值中得出该待确定参数的标准值来得到标准样品。
具体来讲,所述方法通过依次进行的以下步骤来确定所述待确定参数的标准值:对m个准标准样品的待定值参数均进行n次测量,以形成m组数据,每组数据中均有n个数值;根据GB/T15000.3采用方差值Fa对上述m×n个数值进行统计检验,以确定所述m×n个数值是否满足均匀性要求,如果不满足均匀性要求,则表明该所述m×n个数值不能得出待确定参数的标准值,需哟啊重新采集数据;在所述m×n个数值满足均匀性要求的情况下,对所述m×n个数值进行夏皮罗-威尔克正态数据检验,当正态性检验统计量W大于临界值W0.05时,认为所述m×n个数值呈正态分布,如果所述m×n个数值呈非正态分布,则需要舍弃或替换数据;在所述m×n个数值呈正态分布的情况下,对所述m组数据中的每组数据进行柯克仑检验,当该组数据的统计量C小于临界值C0.05,表明该组数据非异常,并舍弃经柯克仑检验呈异常的数据组;对所述m组数据中的经柯克仑检验非异常的数据组进行格拉布斯检验,当该组数据的统计量G小于临界值G0.05,表明该组数据非异常,并舍弃经格拉布斯检验呈异常的数据组;对所述m组数据中的经格拉布斯检验非异常的数据组进行狄克逊检验,当该组数据的统计量D小于临界值D0.05,表明该组数据非异常,并舍弃经狄克逊检验呈异常的数据组;对所述m组数据中的经狄克逊检验非异常的数据组进行几何平均,得出待定值参数的标准值,将具有该标准值的样品作为标准样品;对所述m组数据中的经狄克逊检验非异常的数据组进行不确定度计算,给出不确定度值。
下面分别给出相关检测的计算公式。
I、标样的均匀性统计计算:
1、均匀性检验F值计算:
式中:
m-单元间的测定值组数;
n-单元内的测定组数;
Xij-各单元内测定值。
2、单元间均匀性检验
(1)计算公式为:
(2)测试方法的重复性欠佳时,用MSwithin代替单元内均匀性研究所用测量的重复性方差,此时sbb计算公式为:
v-自由度,在各测试结果等权重的情况下等于m(n-1);
MSamong-单元间均匀性研究中测量的重复性均方;
MSwithin-单元内均匀性研究中测量的重复性均方;
sbb-单元间不均匀性标准偏差;
ubb-因单元间不均匀性引起的标准不确定度;
n为观测数,n0为有效子(组)单元数,当没有缺失数据时等于n。
II、标样的定值结果统计计算:
1、测试数据的正态性检验
对全部组份的定值数据用夏皮罗-威尔克(Shapiro--Wilk)方法考查检验分析数据的正态性,当分析数据的正态性检验统计量W大于临界值(W>W0.05),表明该分析数据呈正态性分布,否则为非正态性分布。夏皮罗-威尔克检验W值用以下公式进行计算:
S=∑ak[x(n+1-k)-x(k)]
ak-为计算检验统计量W而用的系数;
xi-测定值;
-平均值;
2、异常值检验
(1)柯克仑(COCHRAN)检验
对单元间的定值数据用柯克仑检验(COCHRAN)方法考查检验分析单元数据是否异常,当分析数据的统计量C小于临界值(C<C0.05),表明该单元分析数据为正确值,否则为异常。柯克仑检验C值用以下公式进行计算:
C-柯克仑检验C值;
smax-单元间最大标准偏差;
si-单元标准偏差;
(2)格拉布斯(Grubbs)统计方法检验
实验室间分析数据用格拉布斯(Grubbs)统计方法检验,检查是否存在异常值,当分析数据的统计量G小于临界值(G<G0.05),表明该单元分析数据为正确值,否则为异常。
G-格拉布斯检验G值;
xmax、xmin-最大、最小测定值;
(3)狄克逊检验统计方法检验
实验室间分析数据用狄克逊统计方法检验,当检测数据为8~10组时,检查是否存在异常值,当分析数据的统计量D小于临界值(D<D0.05),表明该单元分析数据为正确值,否则为异常。
D-格拉布斯检验D值;
xmax、x1、x2-最大、最小、次小测定值;
III、不确定度计算:
特性值的标准不确定度计算公式为:
其中
特性值的扩展不确定度计算公式为:
U=kuCRM
不确定度分量ubb 2就是均匀性统计表中的Sbb 2,由于冶金标准物质的稳定性很好,ults(长期稳定性标准不确定度)和usts(短期稳定性标准不确定度)可忽略不计。为定值统计的单次测量标准偏差,n为数据组数。uchar为因测定引起的标准不确定度,k为包含因子,95%的置信区间时k=2。
此外,本发明的方法也可利用EXCEL电子表格进行计算。
图2示出了根据本发明的一个示例性实施例的制备标准样品的方法对钛精矿中的S进行均匀性检验的示例。
图3示出了根据本发明的一个示例性实施例的制备标准样品的方法对不同单位的样品中S元素进行定值分析检验的示例。
尽管上面已经结合附图和示例性实施例描述了本发明,但是本领域普通技术人员应该清楚,在不脱离权利要求的精神和范围的情况下,可以对上述实施例进行各种修改。
Claims (4)
1.一种制备标准样品的方法,所述标准样品具有待定值的参数,其特征在于,所述方法通过依次进行的以下步骤来确定所述参数的标准值:
对m个准标准样品的待定值参数均进行n次测量,以形成m组数据,每组数据中均有n个数值;
根据GB/T15000.3采用方差值Fa对上述m×n个数值进行统计检验,以确定所述m×n个数值是否满足均匀性要求;
在所述m×n个数值满足均匀性要求的情况下,对所述m×n个数值进行夏皮罗-威尔克正态数据检验,当正态性检验统计量W大于临界值W0.05时,认为所述m×n个数值呈正态分布;
在所述m×n个数值呈正态分布的情况下,对所述m组数据中的每组数据进行柯克仑检验,当该组数据的统计量C小于临界值C0.05,表明该组数据非异常;
对所述m组数据中的经柯克仑检验非异常的数据组进行格拉布斯检验,当该组数据的统计量G小于临界值G0.05,表明该组数据非异常;
对所述m组数据中的经格拉布斯检验非异常的数据组进行狄克逊检验,当该组数据的统计量D小于临界值D0.05,表明该组数据非异常;
对所述m组数据中的经狄克逊检验非异常的数据组进行几何平均,得出待定值参数的标准值,将具有该标准值的样品作为标准样品;
对所述m组数据中的经狄克逊检验非异常的数据组进行不确定度计算,给出不确定度值。
2.根据权利要求1所述的制备标准样品的方法,其特征在于,所述狄克逊检验步骤中的经格拉布斯检验非异常的数据组为8至10组。
3.根据权利要求1所述的制备标准样品的方法,其特征在于,所述标准样品具有多个待定值的参数,并且对所述多个待定值的参数中的每个均确定其标准值。
4.根据权利要求1所述的制备标准样品的方法,其特征在于,所述标准样品为冶金标准样品。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012105744339A CN103063490A (zh) | 2012-12-26 | 2012-12-26 | 一种制备标准样品的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012105744339A CN103063490A (zh) | 2012-12-26 | 2012-12-26 | 一种制备标准样品的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103063490A true CN103063490A (zh) | 2013-04-24 |
Family
ID=48106204
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012105744339A Pending CN103063490A (zh) | 2012-12-26 | 2012-12-26 | 一种制备标准样品的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103063490A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103674659A (zh) * | 2013-12-27 | 2014-03-26 | 株洲硬质合金集团有限公司 | 一种碳化钨粉游离碳标准样品制备方法 |
CN108037148A (zh) * | 2017-12-22 | 2018-05-15 | 江西瑞林稀贵金属科技有限公司 | X射线荧光光谱法测定铜冶炼炉渣多元素标准样品的研制 |
CN108335233A (zh) * | 2018-03-01 | 2018-07-27 | 深圳凯达通光电科技有限公司 | 一种基于大数据技术的智能电网量测数据处理系统 |
CN112945785A (zh) * | 2021-02-04 | 2021-06-11 | 华润怡宝饮料(中国)有限公司 | 一种利用铝箔片检验耐破度仪性能的方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102539208A (zh) * | 2011-12-28 | 2012-07-04 | 东北轻合金有限责任公司 | 2d70铝合金光谱与化学标准样品的制备方法 |
-
2012
- 2012-12-26 CN CN2012105744339A patent/CN103063490A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102539208A (zh) * | 2011-12-28 | 2012-07-04 | 东北轻合金有限责任公司 | 2d70铝合金光谱与化学标准样品的制备方法 |
Non-Patent Citations (3)
Title |
---|
全国标准样品技术委员会秘书处: "《中华人民共和国国家标准 GB/T15000.3-2008》", 7 March 2008 * |
田衎等: "松花江哈尔滨段水系沉积物中无机元素环境标准样品的研制", 《岩矿测试》 * |
苏仲鸣等: "碳素钢埋弧焊用焊剂标准样品研制", 《焊管》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103674659A (zh) * | 2013-12-27 | 2014-03-26 | 株洲硬质合金集团有限公司 | 一种碳化钨粉游离碳标准样品制备方法 |
CN103674659B (zh) * | 2013-12-27 | 2016-03-02 | 株洲硬质合金集团有限公司 | 一种碳化钨粉游离碳标准样品制备方法 |
CN108037148A (zh) * | 2017-12-22 | 2018-05-15 | 江西瑞林稀贵金属科技有限公司 | X射线荧光光谱法测定铜冶炼炉渣多元素标准样品的研制 |
CN108335233A (zh) * | 2018-03-01 | 2018-07-27 | 深圳凯达通光电科技有限公司 | 一种基于大数据技术的智能电网量测数据处理系统 |
CN112945785A (zh) * | 2021-02-04 | 2021-06-11 | 华润怡宝饮料(中国)有限公司 | 一种利用铝箔片检验耐破度仪性能的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103063490A (zh) | 一种制备标准样品的方法 | |
CN103246821B (zh) | 一种基于仿真的多应力小样本加速寿命试验方案设计优化方法 | |
CN101226094A (zh) | 一种用于六维力传感器标定装置的标定方法 | |
CN101556229B (zh) | 弹簧片刚度精密测试仪 | |
CN105572191B (zh) | 一种电化学气体传感器的压力补偿方法 | |
CN111829639B (zh) | 组合式汽车衡的标定方法 | |
CN102081143B (zh) | 电池电量的估测方法与系统 | |
CN105629052A (zh) | 芯片功耗实时检测方法 | |
CN103983179B (zh) | 电池厚度变化检测装置、电池安全检测及判断方法 | |
CN106053727A (zh) | 一种标准曲线的校正方法和系统 | |
CN102879061A (zh) | 一种基于拟合方程的水表误差校正方法 | |
CN110441727A (zh) | 一种对电能表现校仪进行状态评价的方法及装置 | |
CN106557652A (zh) | 判断样品检测数据可疑值的方法 | |
CN208206384U (zh) | 一种横向摩擦力校准装置 | |
CN105550383B (zh) | 一种非定常气动力测量试验系统的设计方法 | |
CN203811126U (zh) | 一种应力环变形检测装置 | |
CN204405306U (zh) | 疲劳寿命计标定装置 | |
CN101592692B (zh) | 测量机台的评估方法 | |
CN105260617B (zh) | 基于可靠性的仪表测量风险评估方法 | |
CN104101418A (zh) | 一种电子分析天平微量加载与漂移判别方法 | |
CN105101031B (zh) | 麦克风的配置方法 | |
CN202928993U (zh) | 一种钢屋架损伤诊断器 | |
CN202041488U (zh) | 生物产品的离子测量装置 | |
CN103336059B (zh) | 一种判别支柱瓷绝缘子机械强度的方法 | |
CN203534740U (zh) | 用于汽车后桥测试系统的三维力传感器 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20130424 |