CN107025378B - 一种基于标偏分位百分比的均匀性评价方法 - Google Patents

一种基于标偏分位百分比的均匀性评价方法 Download PDF

Info

Publication number
CN107025378B
CN107025378B CN201710225417.1A CN201710225417A CN107025378B CN 107025378 B CN107025378 B CN 107025378B CN 201710225417 A CN201710225417 A CN 201710225417A CN 107025378 B CN107025378 B CN 107025378B
Authority
CN
China
Prior art keywords
standard deviation
tobacco
data
value
chemical
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.)
Active
Application number
CN201710225417.1A
Other languages
English (en)
Other versions
CN107025378A (zh
Inventor
张军
薛庆逾
石超
王筑临
李勇
杨泽会
尚关兰
秦家文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
SHANGHAI MICRO-VISION TECHNOLOGY Ltd
Original Assignee
SHANGHAI MICRO-VISION TECHNOLOGY Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by SHANGHAI MICRO-VISION TECHNOLOGY Ltd filed Critical SHANGHAI MICRO-VISION TECHNOLOGY Ltd
Priority to CN201710225417.1A priority Critical patent/CN107025378B/zh
Publication of CN107025378A publication Critical patent/CN107025378A/zh
Application granted granted Critical
Publication of CN107025378B publication Critical patent/CN107025378B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/359Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass

Abstract

本发明提供一种基于标偏分位百分比的均匀性评价方法,包括步骤:一:烟叶取样;二:获取烟叶评价指标的化学值x;三:求解x的均值xm;计算置信概率的上下限[Up,Lo];四:计算不在[Lo,Up]内的化学值,计算[Lo,Up]内化学值的标准偏差;五:对化学值x进行降序排列得到y;六:计算所有数据的标偏;七:设y的最后一个数据的序号为n,剔除第一个数据计算标偏s1,剔除第n个数据得到标偏s2,经过k次得到PSD向量;八:选取数值区间,以p为数据间隔,设置预设标偏;九:计算PSD小于给定区间数值的数据对应的最少剔除数据个数及此个数占总体样本个数的百分比。本方法的发明评价方法较现有的评价方法更为客观。

Description

一种基于标偏分位百分比的均匀性评价方法
技术领域
本发明涉及一种基于标偏分位百分比的均匀性评价方法,属于烟草质量控制领域。
背景技术
随着近年来,人们对烟叶质量的要求的逐渐提高,均质化生产越来越多的应用于打叶复烤领域,片烟的化学成分的调控多以尼古丁为主,片烟的均匀性指标一般采用化学值的标偏与化学值的均值的比值作为变异系数来衡量,目前绝大多数打叶复烤的均质化控制水平都一般控制在5%的变异系数,但是由于烟叶作为一种特殊的农产品,其化学值均值水平,批次规模,样本量,都很难一致,不同复烤厂之间的片烟的波动水平很难由变异系数来衡量,因此在实际比较的过程中一般会采用标偏的大小,但是在不同复烤厂,不同批次标偏比较的过程中,又会导致标偏0.07与0.08差异不大,而在同等均值水平下变异系数3.5%与4%又差异很大的情况出现;有些由于取样量非常少,且集中取样,导致标偏很低,这样评价又有失公允;在深入研究打叶复烤均匀性评价体系中,华环国际烟草有限公司上海烟草集团有限责任公司发表的专利申请号为201610280750.8,发明名称为一种测试在打叶复烤中基于均质化控制的效果的方法中对打叶复烤的均质化,基于过程与整体,全局与局部对均质化进行了分三部分评价,整体变异系数,局部变异系数以及协变量变异系数的评价体系,上述方法是探索均质化评价体系中比较详实的尝试,能够很快的分析出整体的稳定性水平以及过程中是否出现异常,调控所发挥作用的来源,但是在上述评价烟叶化学值均质化的过程中忽略了异常数据对整个均质化评价的影响,如有一批非常稳定的化学值数据,但是偶尔出现一些离群的化学值,这种情况会容易把标偏与变异系数都拉高,在实际的执行的过程中又很难区分哪些的确是化学成分的异常,其异常数据的评判标准尚不明朗。显然,对于上述的数据是非常稳定的数据只不过是出现了异常,而非不稳定。
发明内容
鉴于上述方法的缺点,本发明的目的在于提供一种基于标偏分位百分比的均匀性评价方法,在获取烟叶化学值后根据化学值的分布范围求出化学值的置信分布区间,同时运用满足给定区间标偏数据的最大占比即标偏分位百分比来衡量比较不同复烤厂、不同批次,不同类型烟叶,不同异常数据类型下的各种数据的均匀性。
本发明采用了如下技术方案:
一种基于标偏分位百分比的均匀性评价方法,其特征在于,包括如下步骤:
步骤一:烟叶取样;
步骤二:获取烟叶的化学值并选择评价指标,形成该评价指标的化学值x;
步骤三:求解x的均值xm,以及标偏xstd;并求解在显著性水平alpha下的临界值t;
则xm的1-alpha置信概率的上限Up与下线Lo的计算公式如下:
Lo=xm-xstd×t (1)
Up=xm+xstd×t (2)
步骤四:计算不在[Lo,Up]区间内的化学值,并计算在[Lo,Up]区间内的化学值的1-alpha置信概率下的标准偏差;形成每个显著性水平下的标偏CPSD;
步骤五:对化学值x进行降序排列得到排序化学值y;
步骤六:计算所有数据的标偏STD0;
步骤七:设y的第一个数据序号为1,最后一个数据的序号为n,
只剔除第一个数据,计算标偏得到s1,只剔除第n个数据得到标偏s2,则
Figure BDA0001265063800000031
删除上述计算PSD取值对应序号的y值,然后重新计算,经过k次得到PSD向量;
步骤八:选取[0,S]数值区间,以p为数据间隔,每个数据点代表预设标偏;
步骤九:计算PSD小于给定区间数值的数据对应的最少剔除数据个数及此个数占总体样本个数的百分比PPSD。
进一步,本发明的基于标偏分位百分比的均匀性评价方法,还可以具有这样的特征:步骤二中,使用近红外检测法获取烟叶的化学值。
进一步,本发明的基于标偏分位百分比的均匀性评价方法,还可以具有这样的特征:其中,步骤二中,处理烟叶前应先进行预处理,预处理的方法如下:去除烟叶上的细土和泥沙,抽去主脉及直径超过2mm的支脉,将烟叶剪成片或者切成丝放入烘箱中,在40℃的温度下烘4h左右,直至可用手指捻碎,将烘好的烟片或者烟丝取出,立即用粉碎机粉碎,进行实验室近红外定量检测。
进一步,本发明的基于标偏分位百分比的均匀性评价方法,还可以具有这样的特征:步骤二中,选择尼古丁作为评价指标。
进一步,本发明的基于标偏分位百分比的均匀性评价方法,还可以具有这样的特征:步骤七中,还包括判断s1是否小于等于s2的步骤。
进一步,本发明的基于标偏分位百分比的均匀性评价方法,还可以具有这样的特征:步骤七中,还包括判断K是否小于n的步骤。
发明的有益效果
1、从均匀性的评价的角度出发,本发明方法规避了传统运用标偏或者变异系数评价带来的数据基准不一致,标偏差异不明显,尤其是异常数据对均匀性评价带来的影响。往往在很多时候一批数据很稳定出现极个别的数据很容易把数据的整体标偏给拉高,从而得出错误的结论,本方法的发明评价方法较现有的评价方法更为客观。
2、在均质化推进的过程中,往往有很多的结论模糊,含糊不清为探讨均质化调控方法带来了错误的引导,使得多数实验推进不下去,本方法避免了打叶复烤均质化调控过程中的各种调控方式尝试以及简单运用标偏评价带来的误导,为构建”均质化调控体系”与”均质化评价体系“做了有益的尝试”。
附图说明
图1是CHF批次尼古丁的置信区间;
图2是不同批次尼古丁的PSD值;
图3是不同批次的原始数据的百分比;
图4是基于标偏分位百分比的均匀性评价方法计算流程图。
具体实施方式
以下结合附图来说明本发明的具体实施方式。
基于标偏分位百分比的均匀性评价方法计算流程图见图4。
步骤11:在原烟或者片烟处取样,获取不同批次或者不同复烤厂的烟叶样本;
步骤12:选取烟叶,去除烟叶上的细土和泥沙,抽去主脉及直径超过2mm的支脉,将烟叶剪成片或者切成丝放入烘箱中,在40℃的温度下烘4h左右,直至可用手指捻碎。将烘好的烟片或者烟丝取出,立即用粉碎机粉碎,进行实验室近红外定量模型的检测;
步骤13:获取烟叶的化学值并选择评价指标,形成该评价指标的化学值x;
步骤14:求解x的均值,xm,以及标偏xstd;并求解在显著性水平alpha的临界值t;
则xm的1-alpha置信概率的上限Up与下限Lo的计算公式如下:
Lo=xm-xstd×t (1)
Up=xm+xstd×t (2)
步骤15:计算不在[Lo,Up]区间内的化学值;
步骤16:计算在[Lo,Up]区间内的化学值的1-alpha置信概率下的标准偏差;形成每个显著性水平下的标偏CPSD,并形成每个置信概率下的CPSD表;
步骤17:对化学值x进行降序排列得到排序化学值y;
步骤18:计算所有数据的标偏STD0;
步骤19:设y的第一个数据序号为1,最后一个数据的序号为n,
只剔除第一个数据,计算标偏得到s1,只剔除第n个数据得到标偏s2,
步骤20,判断s1是否小于等于s2,则有下列公式
Figure BDA0001265063800000061
若判断结果为是,则进入步骤21,PSD(k)=s1;
若判断结果为否,则进入步骤22,PSD(k)=s2。
步骤23,删除上述计算PSD取值对应序号的y值,然后重新计算,经过k次得到PSD向量;步骤24中判断K是否小于n,若是,则再进行一次上述求PSD的过程,若否,则进入步骤25,生成各组数据的PSD。
步骤26:选取[0,S]数值区间,以p为数据间隔,每个数据点代表预设标偏;
步骤27:计算PSD小于给定区间数值的数据对应的最少剔除数据个数及此个数占总体样本个数的百分比PPSD。
实际验证:
分别从某复烤厂CGY,CHE,CHF,CHB四个批次的烤机出口各取156,201,97,196个烟叶样本;
根据步骤12和步骤13对上述烟叶样本进行预处理后检测化学成分,并选择尼古丁作为均匀性的评价指标;
表1:不同批次的化学值
批次代码 样本个数 均值 标偏 CV值
CGY 156 2.250 0.0787 3.54%
CHE 201 2.388 0.0868 3.63%
CHF 97 2.867 0.0807 2.81%
CHB 196 2.091 0.0830 3.97%
根据步骤14到步骤16求解得到每个批次的上限与下限以及t值,并形成每个置信概率下的标偏CPSD,如表2所示。其中在0.05的显著性水平下,四个批次的化学值的上限与下限分别为[2.122,2.382],[2.27,2.50],[2.761,2.973],[1.982,2.2],t值分别为1.96,t值是按照norminv(1-alpha/2,0,1)来计算的。
CHF批次尼古丁的置信区间如图1所示。此图是典型的数据均匀性较好,但是在60-70序号附近出现了明显的异常数据,而造成不稳定性的表征。
表2:不同批次不同置信概率的标偏
Figure BDA0001265063800000071
Figure BDA0001265063800000081
在99%的置信概率下,标偏的顺序为CHF<CGY<CHB<CHE,CHF批次调控的均匀性最好;
在95%的置信概率下,标偏的顺序为CGY<CHF<CHB<CHE,CGY,CHF批次调控的均匀性最好,CGY批次略优于CHF;
根据步骤17到步骤19分别设置四个批次的K为126,171,67,166,每个批次的初始的标偏分别为0.0787、0.0868、0.0807、0.0830。
根据步骤20到步骤24进行运算,得到不同批次的PSD。不同批次尼古丁的PSD值见图2。
根据步骤26选取预设的标偏区间为S_SD=[0.045:0.0025:0.09];则PSD对应小于区间内数值对应的原始数据删除的最小个数以及百分比见表3:
表3:不同批次符合给定标偏的PPSD计算结果
Figure BDA0001265063800000082
Figure BDA0001265063800000091
并求得原始数据符合给定区间数量的百分比,见图3,从1-PPSD的分布图上可以看出数据在0.0725之前YN2116CGY与YN8116CHF一直无差异,并且远好于YN5116CHE与YN9116CHB批次。
本发明方法规避了传统运用标偏或者变异系数评价带来的数据基准不一致,标偏差异不明显,尤其是异常数据对均匀性评价带来的影响。往往在很多时候一批数据很稳定出现极个别的数据很容易把数据的整体标偏给拉高,从而得出错误的结论,本方法的发明评价方法较现有的评价方法更为客观。

Claims (5)

1.一种基于标偏分位百分比的烟叶近红外检测均匀性评价方法,其特征在于,包括如下步骤:
步骤一:烟叶取样;
步骤二:获取烟叶的化学值并选择评价指标,形成该评价指标的化学值x;
步骤三:求解x的均值xm,以及标准偏差xstd;并求解在显著性水平alpha下的临界值t;
则xm的1-alpha置信概率的上限Up与下限Lo的计算公式如下:
Lo=xm-xstd×t (1)
Up=xm+xstd×t (2)
步骤四:计算不在[Lo,Up]区间内的化学值,并计算在[Lo,Up]区间内的化学值在1-alpha置信概率下的标准偏差;
步骤五:对化学值x进行降序排列得到排序化学值y;
步骤六:计算化学值y的标准偏差STD0;
步骤七:设y的第一个数据序号为1,最后一个数据的序号为n,
只剔除第一个数据,计算标准偏差得到s1,只剔除第n个数据得到标准偏差s2,则
Figure FDA0002752653750000011
删除上述计算PSD取值对应序号的化学值,然后重新计算,经过k次得到PSD向量;
步骤八:选取[0,S]数值区间,以p为数据间隔,每个数据点代表预设标准偏差;
步骤九:计算PSD向量中小于给定数值区间的数据点的标准偏差所对应的最少剔除数据个数及此个数占总体样本个数的百分比PPSD。
2.如权利要求1所述的基于标偏分位百分比的烟叶近红外检测均匀性评价方法,其特征在于:
步骤二中,使用近红外检测法获取烟叶的化学值。
3.如权利要求1所述的基于标偏分位百分比的烟叶近红外检测均匀性评价方法,其特征在于:
其中,步骤二中,处理烟叶前应先进行预处理,预处理的方法如下:去除烟叶上的细土和泥沙,抽去主脉及直径超过2mm的支脉,将烟叶剪成片或者切成丝放入烘箱中,在40℃的温度下烘4h,直至可用手指捻碎,将烘好的烟片或者烟丝取出,立即用粉碎机粉碎,进行实验室近红外定量检测。
4.如权利要求1所述的基于标偏分位百分比的烟叶近红外检测均匀性评价方法,其特征在于:
步骤二中,选择尼古丁作为评价指标。
5.如权利要求1所述的基于标偏分位百分比的烟叶近红外检测均匀性评价方法,其特征在于:
步骤七中,还包括,判断k是否小于n的步骤。
CN201710225417.1A 2017-04-07 2017-04-07 一种基于标偏分位百分比的均匀性评价方法 Active CN107025378B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710225417.1A CN107025378B (zh) 2017-04-07 2017-04-07 一种基于标偏分位百分比的均匀性评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710225417.1A CN107025378B (zh) 2017-04-07 2017-04-07 一种基于标偏分位百分比的均匀性评价方法

Publications (2)

Publication Number Publication Date
CN107025378A CN107025378A (zh) 2017-08-08
CN107025378B true CN107025378B (zh) 2020-12-22

Family

ID=59527298

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710225417.1A Active CN107025378B (zh) 2017-04-07 2017-04-07 一种基于标偏分位百分比的均匀性评价方法

Country Status (1)

Country Link
CN (1) CN107025378B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110264031B (zh) * 2019-05-05 2024-03-01 贵州中烟工业有限责任公司 一种标准出片率的表征方法
CN110411978B (zh) * 2019-07-26 2020-05-22 中国农业大学 一种物料均匀性检测设备及检测方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288602A (zh) * 2011-07-17 2011-12-21 红云红河烟草(集团)有限责任公司 一种快速检测卷烟加料均匀性的方法
CN102680632A (zh) * 2012-05-24 2012-09-19 安徽中烟工业有限责任公司 一种基于木糖醇检测白肋烟处理加料均匀性评价方法
CN102707013A (zh) * 2012-05-24 2012-10-03 安徽中烟工业有限责任公司 一种白肋烟处理加料均匀性评价方法
CN105124731A (zh) * 2015-09-14 2015-12-09 上海创和亿电子科技发展有限公司 烟叶均质化调控区间划分方法
CA2955573A1 (en) * 2014-09-30 2016-04-07 Philip Morris Products S.A. Homogenized tobacco material and method of production of homogenized tobacco material
CN105717190A (zh) * 2016-04-20 2016-06-29 浙江中烟工业有限责任公司 一种烟叶质控样品砷含量的定值方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288602A (zh) * 2011-07-17 2011-12-21 红云红河烟草(集团)有限责任公司 一种快速检测卷烟加料均匀性的方法
CN102680632A (zh) * 2012-05-24 2012-09-19 安徽中烟工业有限责任公司 一种基于木糖醇检测白肋烟处理加料均匀性评价方法
CN102707013A (zh) * 2012-05-24 2012-10-03 安徽中烟工业有限责任公司 一种白肋烟处理加料均匀性评价方法
CA2955573A1 (en) * 2014-09-30 2016-04-07 Philip Morris Products S.A. Homogenized tobacco material and method of production of homogenized tobacco material
CN105124731A (zh) * 2015-09-14 2015-12-09 上海创和亿电子科技发展有限公司 烟叶均质化调控区间划分方法
CN105717190A (zh) * 2016-04-20 2016-06-29 浙江中烟工业有限责任公司 一种烟叶质控样品砷含量的定值方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Evaluation of tobacco mixing uniformity based on chemical composition;Li Junli 等;《Proceedings of the 31st Chinese Control Conference》;20121224;第7552-7555页 *
烟叶配方打叶均匀性的研究;肖明礼 等;《现代食品科技》;20111231;第27卷(第6期);第684-686页 *

Also Published As

Publication number Publication date
CN107025378A (zh) 2017-08-08

Similar Documents

Publication Publication Date Title
Bader et al. Automated threshold selection for extreme value analysis via ordered goodness-of-fit tests with adjustment for false discovery rate
CN107038292B (zh) 一种基于自适应多变量非参数核密度估计的多风电场出力相关性建模方法
CN111597729A (zh) 加工设备工艺参数优化方法及系统
CN107025378B (zh) 一种基于标偏分位百分比的均匀性评价方法
CN111680820A (zh) 分布式光伏电站故障诊断方法和装置
US10918014B2 (en) Fertilization precision control method for water and fertilizer integrated equipment and control system thereof
CN110334406B (zh) 一种考虑风速特大值的极值风速重现期确定方法和装置
CN115453064B (zh) 一种细颗粒物空气污染成因分析方法及系统
Watt et al. Use of advanced modelling methods to estimate radiata pine productivity indices
CN103868957B (zh) 一种评价打叶复烤工序片烟感官品质稳定性的方法
CN108376148A (zh) 一种基于数据隐结构分析的基站聚类检测方法
CN112651173B (zh) 一种基于跨域光谱信息的农产品品质无损检测方法及可泛化系统
CN116469466B (zh) 一种高效预测菊花耐涝性的方法及其应用
CN105787296B (zh) 一种宏基因组和宏转录组样本相异度的比较方法
CN116502943A (zh) 一种熔模铸造产品质量追溯方法
CN110232951A (zh) 判断测序数据饱和的方法、计算机可读介质和应用
CN110853711B (zh) 一种预测烟草果糖含量的全基因组选择模型及其应用
CN114330478A (zh) 针对电网风速预报的风速分类订正方法
Qi et al. Genetic structure and diversity of parental cultivars involved in China mainland sugarcane breeding programs as inferred from DNA microsatellites
CN110853710B (zh) 一种预测烟草淀粉含量的全基因组选择模型及其应用
CN110618458B (zh) 一种基于ica-ra的钻孔应变数据的分频段级联校正方法
CN113946621A (zh) 一种基于关联规则的制丝车间数据波动关系的挖掘方法
JP7025216B2 (ja) トランスクリプトーム解析装置及び解析方法
CN111178627A (zh) 一种基于spca的神经网络混合优化预测方法
CN117332358B (zh) 一种玉米浸泡水处理方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant