WO2023236842A1 - 一种卷烟燃烧外观的评价方法 - Google Patents

一种卷烟燃烧外观的评价方法 Download PDF

Info

Publication number
WO2023236842A1
WO2023236842A1 PCT/CN2023/097776 CN2023097776W WO2023236842A1 WO 2023236842 A1 WO2023236842 A1 WO 2023236842A1 CN 2023097776 W CN2023097776 W CN 2023097776W WO 2023236842 A1 WO2023236842 A1 WO 2023236842A1
Authority
WO
WIPO (PCT)
Prior art keywords
cigarette
factor
evaluation
matrix
factors
Prior art date
Application number
PCT/CN2023/097776
Other languages
English (en)
French (fr)
Inventor
郑晗
王浩
谢姣
王旭
余振华
詹建波
余婷婷
余江
张莹
王涛
Original Assignee
云南中烟工业有限责任公司
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 云南中烟工业有限责任公司 filed Critical 云南中烟工业有限责任公司
Publication of WO2023236842A1 publication Critical patent/WO2023236842A1/zh

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84—Systems specially adapted for particular applications
    • G01N21/88—Investigating the presence of flaws or contamination
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N31/00—Investigating or analysing non-biological materials by the use of the chemical methods specified in the subgroup; Apparatus specially adapted for such methods
    • G01N31/12—Investigating or analysing non-biological materials by the use of the chemical methods specified in the subgroup; Apparatus specially adapted for such methods using combustion

Definitions

  • the invention belongs to the field of research on the burning appearance of cigarettes, and specifically relates to a method for evaluating the burning appearance of cigarettes.
  • Cigarette research and development needs to pass the review procedures of clinical trial observation, pharmacology and toxicity testing and ingredient analysis by the health department. Cigarette manufacturing needs to pass quality appraisal and make quality judgments in accordance with the national standards GB5606 ⁇ 5610-85.
  • Quality inspection mainly focuses on evaluating the intrinsic quality, external quality and tar content of cigarettes.
  • the burning appearance of cigarettes is an important appearance that consumers intuitively feel during the burning process of cigarette products. Its quality directly affects consumers' judgment of the quality of cigarette products.
  • cigarettes frequently shed ash and fly ash, and in severe cases, the combustion cones fall off. This not only pollutes the environment, but also causes consumers and forced smokers to feel uncomfortable and disgusted with the environment.
  • the present invention aims to comprehensively evaluate the burning appearance of cigarettes during smoking, form a complete set of analysis and evaluation methods, provide effective cigarette quality information for the cigarette manufacturing industry, and provide information for the improvement of cigarette manufacturing processes. Good at providing effective feedback.
  • the present invention provides a cigarette combustion appearance analysis and evaluation method, establishes a set of fuzzy comprehensive evaluation models for combustion appearance applicable to all types of cigarettes, and uses a full vision camera to evaluate various external attributes of cigarettes. Recorded and input into the model, a comprehensive score is finally obtained to reflect the burning appearance of the tested cigarette.
  • a first aspect of the present invention provides a method for evaluating the burning appearance of cigarettes, which includes the following steps:
  • the external attributes of the cigarette include: crack rate, ash column deviation, ash column grayscale, ash column length, ash column Area, ash shrinkage rate, fly ash value, ash retention rate, carbon line width, carbon line uniformity and burning speed, a total of 11 attributes;
  • the robotic arm that simulates the human cigarette smoking process and smoking environment can refer to the application filed by the applicant on April 20, 2020.
  • the application number is 202010329623.9; the patent name is a machine that simulates the entire process of human cigarette smoking.
  • the method of using the robotic arm to collect the external attributes of cigarettes can refer to the application filed by the applicant on April 26, 2020.
  • the application number is 202010340747.7; the patent name is a robot-based method.
  • the technical solutions of the above two patents on the measurement method of the mechanical arm and cigarette pack ash performance are incorporated into this application.
  • the crack rate is the ratio of the area of ash pieces falling or cracking on the ash column to the overall surface area of the ash column;
  • the ash column deviation refers to the maximum distance between the ash column of the cigarette burning and the axis of the cigarette filter rod. Deviation angle;
  • ash column grayscale reflects the relative color difference of the remaining part of the cigarette except for the crack after burning under the same reference system. The larger the value, the whiter the ash column; the ash column length is the length of the ash column formed after the cigarette is burned.
  • the ash column area is the width of the ash column formed after the cigarette is burned; the ash shrinkage rate is the ratio of the ash column area after the cigarette is burned to the unburned area; the fly ash value is the area of ash scattered outside the fixed ash point; persistence
  • the ash rate is the ratio of the longest value of the ash column that can be maintained during the cigarette burning process to the length of the cigarette from burning to the butt;
  • the charcoal line width is the width of the charcoal line at the bottom of the cigarette's burning cone; the charcoal line uniformity refers to the bottom of the cigarette's burning cone
  • the maximum deviation angle between the carbon line and the cross section of the cigarette; the burning speed is the ratio of burning time and burning length.
  • Standardization of cigarette data Standardize the data to obtain standard data X through the following formula:
  • min(X og ) represents the minimum value of X og ;
  • max( X og ) represents the maximum value of X og ;
  • Definite attribute set u (u 1 , u 2 ,..., u n ), where u 1 , u 2 ,..., u n respectively represent different external attributes of cigarettes;
  • fuzzy sets of single-factor evaluations of all attributes together form a fuzzy comprehensive evaluation matrix R with a length of n ⁇ m;
  • each common factor is an independent normal random variable with a mean value of 0 and a variance of 1, and its covariance matrix is the unit matrix I p , that is, F ⁇ N(0,I p ); it is determined that the relationships between each special factor and the special factors are independent of each other and the common factors, that is,
  • ⁇ 1 2 , ⁇ 2 2 ,..., ⁇ n 2 represent the variance of special factors
  • the variance of the special factor is estimated by the diagonal element of R c -CC T , that is
  • the residual matrix can be represented by R c -CC T -Cov(U); when CC T +Cov(U) is close to the correlation coefficient matrix R c , that is, when the residual is less than 5%, the model is considered to have a better fitting degree. Continue after recovery;
  • the elements in the fuzzy vector B are the same as the evaluation set v
  • the elements in have a one-to-one correspondence, and the evaluation corresponding to the largest element is the fuzzy comprehensive final evaluation grade of the cigarette.
  • step S1 based on a robotic arm that simulates the entire process of human cigarette smoking, the real-time smoking status of the cigarette is captured through a camera system, and the external attributes of the cigarette are collected.
  • r ij is obtained by expert scoring.
  • the present invention has the following beneficial effects:
  • the present invention establishes a fuzzy comprehensive model of cigarette burning appearance, which can determine the data factor weight vectors of different attributes from different samples, thereby obtaining a graded evaluation of the cigarette burning appearance.
  • the evaluation method of the present invention is consistent with the actual appearance performance results of cigarette combustion.
  • the evaluation results can objectively reflect the combustion appearance performance of cigarettes. It is important to objectively analyze the performance of cigarette pack ash and fly ash and to compare the external attributes of different cigarettes. significance.
  • Figure 1 is an overall flow chart of the method of the present invention.
  • the invention provides a method for analyzing and evaluating the burning appearance of cigarettes, which includes the following steps:
  • the collected external attributes of cigarettes include: crack rate, ash column deviation, ash column grayscale, ash column length, ash column area, ash shrinkage rate, fly ash value, ash retention rate, carbon line width, and carbon line uniformity. and burning speed, a total of 11 attributes.
  • the meanings of these attributes are as follows: the crack rate is the ratio of the falling or cracked area of ash pieces on the ash column to the overall surface area of the ash column; the ash column deviation refers to the maximum deviation angle between the ash column of cigarette burning and the axis of the cigarette filter rod;
  • the gray scale of the ash column reflects the relative color difference of the remaining part of the cigarette except for the crack after burning under the same reference system.
  • min(X og ) represents the minimum value of X og ;
  • max(X og ) represents the maximum value of X og ;
  • n 11 in the present invention.
  • the variance of the special factor is estimated by the diagonal element of R c -CC T , that is

Abstract

本发明公开一种卷烟燃烧外观评价方法，其包括如下步骤：S1、卷烟外在属性数据采集；S2、卷烟数据标准化；S3、构建卷烟外在属性模糊综合评价矩阵；S4、通过因子分析确定每个属性的权向量以完成模糊综合评价模型计算，并根据计算得到的模糊向量B确定最终评价等级。针对卷烟燃烧外观分析与评定问题，本发明建立一套卷烟燃烧外观模糊综合模型，能够确定不同属性和来自不同样本的数据的因素权向量，结合对卷烟样本各外观属性的主观评判以获取对卷烟燃烧外观的等级评价。

Description

一种卷烟燃烧外观的评价方法 技术领域
本发明属于卷烟燃烧外观研究领域，具体涉及一种卷烟燃烧外观的评价方法。
背景技术
随着市场经济的高速发展，卷烟行业一直保持着良好的发展前景，逐步实现烟叶生产智能化、防控技术环保化、质量追溯链条化。随着人民群众消费水平的逐步上升，卷烟的生产质量标准不断提高，卷烟制造也随之持续进行着技术更新与迭代，探索减小吸入伤害、改善产品质量的生产工艺。
卷烟研创需要通过卫生部门的临床试验观察、药理和毒性检验和成分分析等审查程序，卷烟制造则需要通过质量鉴定，按照国家颁布的GB5606～5610-85标准进行质量判定。质检主要围绕卷烟的内在质量、外在质量、焦油含量进行评价。卷烟燃烧外观是卷烟产品在燃烧过程中消费者直观感受到的重要外观形态，其优劣直接影响着消费者对卷烟产品质量的判断。同时，卷烟在抽吸过程中，烟支的频繁落灰、飞灰，严重时发生燃烧锥脱落的现象，不仅污染环境，也是引起消费者及被迫吸烟者对所处环境不适、反感的重要影响因素之一，还会影响消费者对于卷烟品牌的忠诚度，甚至带来一定安全隐患。另一方面，随着卷烟技术的发展，消费水平和质量的不断提升，消费者对卷烟质量要求也不断提高，对卷烟产品的关注和选择由口味、价格、包装等方面逐渐向更多因素拓展，卷烟燃烧后的外观状态由于其直观性，越来越受到消费者的重视。从消费者的角度而言，卷烟燃烧和抽吸时带来的体验则是最重要的评价方向，如燃烧时间、烟灰可停留在烟支上的长度等。
近年来，烟草行业内针对卷烟燃烧外观形态开展了相关研究，主要包括卷烟包灰与飞灰性能研究和检测方法研究两个方面，但是没有一套完备的卷烟在抽吸过程中燃烧外观的评价方法，而卷烟燃烧外观的客观评价是对卷烟外观的客观分析以及不同卷烟之间燃烧外观形态进行比较具有重要意义。
为了解决以上问题，提出本发明。
发明内容
本发明旨在对卷烟在抽吸过程中的燃烧外观进行综合评测，形成一套完备的分析评定方法，为卷烟制造行业提供有效的卷烟质量信息，为卷烟制造工艺的改 善提供有效反馈。
为了填补上述评价项目的空白，本发明提供一种卷烟燃烧外观分析与评定方法，建立一套适用于所有类型卷烟的燃烧外观模糊综合评价模型，通过全视觉相机对卷烟的各项外在属性进行记录并输入模型中，最终得到一个综合评分从而反映所测卷烟的燃烧外观。
本发明第一方面提供一种卷烟燃烧外观评价方法，其包括如下步骤：
S1、卷烟外在属性数据采集：在卷烟实时抽吸状态时，采集卷烟外在属性，所述卷烟外在属性包括：裂口率、灰柱偏离度、灰柱灰度、灰柱长度、灰柱面积、缩灰率、飞灰值、持灰率、炭线宽度、炭线整齐度及燃烧速度，共11个属性；
采集卷烟外在属性，可以采用现有技术中常用的设备和方法，当然更优选地，基于模拟人体卷烟抽吸全过程的机械手臂，通过摄像系统捕捉卷烟实时抽吸状态，采集卷烟外在属性。其中模拟人体卷烟抽吸过程和抽吸环境的机械手臂，可以参考本申请人于2020年04月20日已经申请的，申请号为202010329623.9；专利名称为一种模拟人体卷烟抽吸全过程的机械手臂及其模拟方法，的专利申请，同时利用该机械手臂采集卷烟外在属性的方法可以参考本申请人于2020年04月26日已经申请的，申请号为202010340747.7；专利名称为一种基于机械手仿真人体卷烟抽吸全过程及全视觉测定表征卷烟燃烧包灰性能的方法，的专利申请。同时将上述两篇专利对于机械手臂和卷烟包灰性能的测定方法的技术方案表述并入本申请中。
上述卷烟外在属性的含义如下：裂口率是灰柱上灰片掉落或裂开面积与灰柱整体表面积的比值；灰柱偏离度指烟支燃烧灰柱与烟支滤棒轴心的最大偏离角度；灰柱灰度反映了卷烟燃烧后除去裂口以外剩余部分在同一参考系下的相对颜色差异，数值越大，灰柱越白；灰柱长度是烟支燃烧后形成的灰柱的长度；灰柱面积是烟支燃烧后形成的灰柱的宽度；缩灰率是卷烟燃烧后灰柱面积和未燃烧前面积的比值；飞灰值是散落在定点弹灰点以外的烟灰面积；持灰率是卷烟燃烧过程中能保持灰柱的最长值与烟支燃烧至烟蒂长度的比值；炭线宽度是烟支燃烧锥底部炭线的宽度；炭线整齐度是指烟支燃烧锥底部炭线与烟支横截面的最大偏离角度；燃烧速度是燃烧时间和燃烧长度的比值。
S2、卷烟数据标准化：通过下式对数据进行标准化处理得到标准数据X：
式中min(Xog)，表示Xog的最小值；max(Xog)，表示Xog的最大值；Xog＝(x1og,x2og,...,x11og)为步骤S1采集的卷烟数据，其中x1og,x2og,...,x11og表示烟支11个外在属性的测量计算值；n＝11；
S3、构建卷烟外在属性模糊综合评价矩阵：
定属性集u＝(u1,u2,...,un)，其中u1,u2,...,un分别表示不同的卷烟外在属性；
定评价集v＝(v1,v2,...,vm)，其中v1,v2,...,vm表示不同的评价结果；
定rij为属性集中第i个属性对评价集V中第j种评价的隶属度，则第i个属性进行单因素评价的结果表示为模糊集合Ri：Ri＝(ri1,ri2,...,rim)；
则全部属性单因素评价的模糊集合共同构成了长度为n×m的模糊综合评价矩阵R；
S4、通过因子分析确定每个属性的权向量以完成模糊综合评价模型计算，并根据计算得到的模糊向量B确定最终评价等级：
给每个属性建立如下形式的因子分析模型：xi＝ci1F1+ci2F2+...+cipFp+qiUi(i＝1,2,...,n)，
其中p代表公共因子数量(即变量数量)，F＝F1,...,Fp为每个变量都有的公共因子，Ui为仅与变量xi有关的特殊因子；cij,qi(i＝1,2,...,n,j＝1,2,...,p)为因子载荷，其中cij(i＝1,2,...,n,j＝1,2,...,p)构成了因子载荷矩阵C；
则上述因子分析模型表示为如下的矩阵形式：X＝CF+QU，
其中：Q＝diag(q1,q2,...,qn)，U＝(U1,U2,...,Un)；
认定各公共因子都是均值为0、方差为1的独立正态随机变量，其协方差矩阵为单位阵Ip，即F～N(0,Ip)；认定各特殊因子之间以及特殊因子与公共因子之间相互独立，即
其中σ1 2,σ2 2,...,σn 2表示特殊因子的方差；
则p个公共因子对第i个变量方差的贡献为：hi 2＝ai1 2+ai2 2+...+aip 2；
第i个变量的方差有如下分解：Var xi＝hi 2+σi 2,i＝1,2,...,n；
定λ1≥λ2≥...≥λn为样本相关系数Rc的特征值，η＝η1,η2,...,ηn为相应的标准正交化特征向量，此时样本相关系数矩阵Rc的主成分因子载荷矩阵C为：
特殊因子的方差用Rc-CCT的对角元来估计，即
其中残差矩阵可用Rc-CCT-Cov(U)来表示；CCT+Cov(U)与相关系数矩阵Rc比较接近时，即残差小于5％时，认为模型的拟合程度较好后继续进行；
通过因子逆时针旋转使变换后的因子具有清晰的可解释性及唯一性，当针对两个因子进行平面正交旋转时，因子载荷矩阵C具有如下形式：C＝(cij),i＝1,2,...,n,j＝1,2；
取正交矩阵如下：
则D＝CS＝(dij),i＝1,2,...,n,j＝1,2)为旋转因子载荷矩阵，模型变为如下形式：X＝D(STF)+QU；
此时公共因子由F变为QTF，当前旋转双因子的目的是将变量分为由不同因子说明的两个部分，故需要将(d11 2,d21 2,...,dn1 2)和(d12 2,d22 2,...,dn2 2)的方差最大化，二者的相对方差可以表示为：以消除dij符号的影响，消除各个属性变量对公共因子依赖程度不同的影响；
接下来将总方差V＝V1+V2最大化，令经计算，φ满足： 其中：
对于p个公共因子，将这些因子两两进行旋转；在一遍完整的旋转过程中，从p个因子中每次选两个进行逆时针旋转，共有p(p-1)/2次旋转，完成一遍旋转后，继续进行下一遍循环，在此重复过程中C各列的相对方差和总方差V会递增，直至某一次循环的总方差比前和一次循环的总方差相比差距不大于5％则停止旋转，此时得到最终旋转因子载荷矩阵D，和标准正交化特征向量相乘即得到模糊综合评价的因素权向量：A＝DηT，A实际上反映了u中属性的重要程度；
然后通过模糊变换将反映u的模糊向量A转换成反映v的模糊向量B，即B＝AR＝(b1,b2,...,bm)，模糊向量B中的元素与评价集v中的元素一一对应，最大的元素对应的评价即为卷烟的模糊综合最终评价等级。
优选地，步骤S1中，基于模拟人体卷烟抽吸全过程的机械手臂，通过摄像系统捕捉卷烟实时抽吸状态，采集卷烟外在属性。
优选地，步骤S3中，m＝4，按照优劣排序，v1是特等级，v2是优等级，v3是中等级，v4是次等级。
优选地，步骤S3中，rij由专家打分获得。
相对于现有技术，本发明具有以下有益效果：
1、针对卷烟燃烧外观分析与评定问题，本发明建立一套卷烟燃烧外观模糊综合模型，能够确定来自不同样本的不同属性的数据因素权向量，从而获取对卷烟燃烧外观的等级评价。
2、本发明评价方法与卷烟燃烧实际外观性能结果一致，评价结果可以客观反映出卷烟燃烧外观性能，对卷烟包灰及飞灰性能的客观分析以及不同卷烟之间外在属性能进行比较具有重要意义。
附图说明
图1为本发明方法的总体流程图。
具体实施方式
下面结合具体实施例对本发明进行说明，但本发明的实施方式不限于此。实施例中未注明具体条件的实验方法,通常按照常规条件以及手册中所述的条件,或按照制造厂商所建议的条件所用的通用设备、材料、试剂等,如无特殊说明,均可从商业途径得到。以下实施例和对比例中所需要的原料均为市售。
实施例1
本发明提供一种卷烟燃烧外观分析与评定方法，其包括如下步骤：
S1、卷烟外在属性数据采集。在机器人仿真抽吸状态和环境模拟条件下，通过多组摄像头在不同角度捕捉卷烟实时状态，包括裂口率、灰柱偏离度、灰柱灰度、灰柱长度、灰柱面积、缩灰率、飞灰值、持灰率、炭线宽度、炭线整齐度和燃烧速度。
使用两组多轴工业机械手夹持卷烟，以模拟吸烟者吸烟时的手臂及手腕动作，其中机械手的移动速度和角度在一定范围内可调。通过机械手指施力模拟弹烟灰动作，其中力度、间隔、次数在一定范围内可调。使用抽吸装置模拟吸烟过程，其中抽吸模式、抽吸时间、抽吸量、抽吸间隔在一定范围内可调。采用分辨率1600*1200、帧率50fps以上的高速精密智能相机组完成数据采集。
采集到的卷烟外在属性有：裂口率、灰柱偏离度、灰柱灰度、灰柱长度、灰柱面积、缩灰率、飞灰值、持灰率、炭线宽度、炭线整齐度及燃烧速度，共11个属性。这些属性的含义如下：裂口率是灰柱上灰片掉落或裂开面积与灰柱整体表面积的比值；灰柱偏离度指烟支燃烧灰柱与烟支滤棒轴心的最大偏离角度；灰柱灰度反映了卷烟燃烧后除去裂口以外剩余部分在同一参考系下的相对颜色差异，数值越大，灰柱越白；灰柱长度是烟支燃烧后形成的灰柱的长度；灰柱面积是烟支燃烧后形成的灰柱的宽度；缩灰率是卷烟燃烧后灰柱面积和未燃烧前面积的比值；飞灰值是散落在定点弹灰点以外的烟灰面积；持灰率是卷烟燃烧过程中能保持灰柱的最长值与烟支燃烧至烟蒂长度的比值；炭线宽度是烟支燃烧锥底部炭线的宽度；炭线整齐度是指烟支燃烧锥底部炭线与烟支横截面的最大偏离角度；燃烧速度是燃烧时间和燃烧长度的比值。
实施例采集了常规圆周卷烟样品390支，对每支样品采集如上所示的11个 外观指标，各样本各指标的采集值如表1所示，其中：x1代表裂口率，x2代表灰柱偏离度，x3代表灰柱灰度，x4代表炭线宽度，x5代表灰柱长度，x6代表灰柱面积，x7代表缩灰率，x8代表炭线整齐度，x9代表燃烧速度，x10代表飞灰值，x11代表持灰率。
表1实施例样本燃烧外观采集值









S2、卷烟数据标准化。设Xog＝(x1og,x2og,...,x11og)为步骤S1采集的卷烟数据，其中x1og,x2og,...,x11og表示烟支11个外在属性的测量计算值。通过下式对数据进行标准化处理得到具有零均值和单位方差的标准数据X：
式中min(Xog)，表示Xog的最小值；max(Xog)，表示Xog的最大值；本发明中n＝11。
以实施例的四支卷烟样本1、100、167、300为例，标准化后各指标数据如表2所示：
表2标准化后的样本数据示例
S3、构建卷烟外在属性模糊综合评价矩阵。设属性集u＝(u1,u2,...,un)，其中u1,u2,...,un分别表示不同的卷烟外在属性。设评价集v＝(v1,v2,...,vm)，其中v1,v2,...,vm表示不同的评价结果，本发明中m＝4。按照优劣排序，v1是最佳等级，v4是最差等级。实施例设定评价集为v＝(很好,好,较好,一般)。
设rij为属性集中第i个属性对评价集V中第j种评价的隶属度，由专家打分获得，那么对第i个属性进行单因素评价的结果表示为模糊集合Ri：
Ri＝(ri1,ri2,...,rim)
则全部属性单因素评价的模糊集合共同构成了长度为n×m的模糊综合评价矩阵R。
以实施例的四支卷烟样本1、100、167、300为例，每支样本的模糊综合评价矩阵如表3-6所示：
表3模糊综合评价矩阵(1号样本)
表4模糊综合评价矩阵(100号样本)
表5模糊综合评价矩阵(167号样本)

表6模糊综合评价矩阵(300号样本)
S4、通过因子分析确定每个属性的权向量以完成模糊综合评价模型计算。给每个属性建立如下形式的因子分析模型：
xi＝ci1F1+ci2F2+...+cipFp+qiUi(i＝1,2,...,n)
其中p代表公共因子数量(即变量数量)，F＝F1,...,Fp为每个变量都有的公共因子，Ui为仅与变量xi有关的特殊因子。cij,qi(i＝1,2,...,n,j＝1,2,...,p)为因子载荷，其中cij(i＝1,2,...,n,j＝1,2,...,p)构成了载荷矩阵C。此时因子分析模型可以表示为如下的矩阵形式：
X＝CF+QU
其中：Q＝diag(q1,q2,...,qn)，U＝(U1,U2,...,Un)
假设各公共因子都是均值为0、方差为1的独立正态随机变量，其协方差矩阵为单位阵Ip，即F～N(0,Ip)；假设各特殊因子之间以及特殊因子与公共因子之间相互独立，即
其中σ1 2,σ2 2,...,σn 2表示特殊因子的方差。则p个公共因子对第i个变量方差的贡献为：
hi 2＝ai1 2+ai2 2+...+aip 2
第i个变量的方差有如下分解：
Var xi＝hi 2+σi 2,i＝1,2,...,n
设λ1≥λ2≥...≥λn为样本相关系数Rc的特征值，η＝η1,η2,...,ηn为相应的标准正交化特征向量，此时样本相关系数矩阵Rc的主成分因子分析载荷矩阵C为：
至此，对于实施例，经由390支样本计算得到相关系数矩阵如表7所示，主成分因子分析载荷矩阵如表8所示：
表7相关系数矩阵
表8载荷矩阵

特殊因子的方差用Rc-CCT的对角元来估计，即
其中残差矩阵可用Rc-CCT-Cov(U)来表示。
主成分解的解释不是严格唯一的，因此继续通过因子逆时针旋转使变换后的因子具有清晰的可解释性及唯一性。当针对两个因子进行平面正交旋转时，因子载荷矩阵C具有如下形式：
C＝(cij),i＝1,2,...,n,j＝1,2
取正交矩阵如下：
则D＝CS＝(dij),i＝1,2,...,n,j＝1,2)为旋转因子载荷矩阵，模型变为如下形式：
X＝D(STF)+QU
此时公共因子由F变为QTF。当前需要将(d11 2,d21 2,...,dn1 2)和(d12 2,d22 2,...,dn2 2)的方差最大化。二者的相对方差可以表示为：
这样的形式可以消除dij符号的影响，消除各个属性变量对公共因子依赖程 度不同的影响。接下来将总方差V＝V1+V2最大化，令经计算，φ应该满足：
其中：
对于p个公共因子，将这些因子两两进行旋转。在一遍完整的旋转过程中，从p个因子中每次选两个进行逆时针旋转，共有p(p-1)/2次旋转。完成一遍旋转后，继续进行下一遍循环，在此重复过程中C各列的相对方差和总方差V会递增，直至某一次循环的总方差比前和一次循环的总方差相比差距不大于5％则停止旋转。此时得到最终旋转因子载荷矩阵D，和标准正交化特征向量相乘即得到模糊综合评价的因素权向量：
A＝DηT
对于实施例，经由390支样本计算得到旋转后因子载荷矩阵如表9所示：
表9旋转因子载荷矩阵
相应地，模糊综合评价因素权向量为：
A＝(0.12,0.08,0.10,0.14,0.06,0.10,0.07,0.13,0.09,0.04,0.07)
A实际上反映了u中属性的重要程度，现在通过模糊变换将反映u的模糊向量A转换成反映v的模糊向量B，即
B＝(0.337,0.296,0.233,0.134)
模糊向量B中的元素与评价集v中的元素一一对应，最大的元素对应的评价即为卷烟的模糊综合最终评价等级。
以实施例的四支卷烟样本1、100、167、300为例，样本1的模糊向量为B＝(0.153,0.238,0.356,0.253)，其卷烟燃烧外观评级为“较好”；样本100的模糊向量为B＝(0.199,0.348,0.327,0.126)，其卷烟燃烧外观评级为“好”；样本167的模糊向量为B＝(0.046,0.165,0.345,0.444)，其卷烟燃烧外观评级为“一般”；样本300的模糊向量为B＝(0.337,0.296,0.233,0.134)，其卷烟燃烧外观评级为“很好”。

Claims (4)

  1. 一种卷烟燃烧外观评价方法，其特征在于，其包括如下步骤：
    S1、卷烟外在属性数据采集：在卷烟实时抽吸状态时，采集卷烟外在属性，所述卷烟外在属性包括：裂口率、灰柱偏离度、灰柱灰度、灰柱长度、灰柱面积、缩灰率、飞灰值、持灰率、炭线宽度、炭线整齐度、燃烧速度，共11个属性；
    S2、卷烟数据标准化：通过下式对数据进行标准化处理得到标准数据X：
    式中min(Xog)，表示Xog的最小值；max(Xog)，表示Xog的最大值；Xog＝(x1og,x2og,...,x11og)为步骤S1采集的卷烟数据，其中x1og,x2og,...,x11og表示烟支11个外在属性的测量计算值；n＝11；
    S3、构建卷烟外在属性模糊综合评价矩阵：
    定属性集u＝(u1,u2,...,un)，其中u1,u2,...,un分别表示不同的卷烟外在属性；
    定评价集v＝(v1,v2,...,vm)，其中v1,v2,...,vm表示不同的评价结果；
    定rij为属性集中第i个属性对评价集V中第j种评价的隶属度，则第i个属性进行单因素评价的结果表示为模糊集合Ri：Ri＝(ri1,ri2,...,rim)；
    则全部属性单因素评价的模糊集合共同构成了长度为n×m的模糊综合评价矩阵R；
    S4、通过因子分析确定每个属性的权向量以完成模糊综合评价模型计算，并根据计算得到的模糊向量B确定最终评价等级：
    给每个属性建立如下形式的因子分析模型：xi＝ci1F1+ci2F2+...+cipFp+qiUi(i＝1,2,...,n)，
    其中p代表公共因子数量(即变量数量)，F＝F1,...,Fp为每个变量都有的公共因子，Ui为仅与变量xi有关的特殊因子；cij,qi(i＝1,2,...,n,j＝1,2,...,p)为因子载荷，其中cij(i＝1,2,...,n,j＝1,2,...,p)构成了因子载荷矩阵C；
    则上述因子分析模型表示为如下的矩阵形式：X＝CF+QU，
    其中：Q＝diag(q1,q2,...,qn)，U＝(U1,U2,...,Un)；
    认定各公共因子都是均值为0、方差为1的独立正态随机变量，其协方差矩阵为单位阵Ip，即F～N(0,Ip)；认定各特殊因子之间以及特殊因子与公共因子之间相互独立，即
    其中σ1 2,σ2 2,...,σn 2表示特殊因子的方差；
    则p个公共因子对第i个变量方差的贡献为：hi 2＝ai1 2+ai2 2+...+aip 2；
    第i个变量的方差有如下分解：Varxi＝hi 2+σi 2,i＝1,2,...,n；
    定λ1≥λ2≥...≥λn为样本相关系数Rc的特征值，η＝η1,η2,...,ηn为相应的标准正交化特征向量，此时样本相关系数矩阵Rc的主成分因子载荷矩阵C为：
    特殊因子的方差用Rc-CCT的对角元来估计，即
    其中残差矩阵可用Rc-CCT-Cov(U)来表示；CCT+Cov(U)与相关系数矩阵Rc残差小于5％时，认为模型的拟合程度较好后继续进行；
    通过因子逆时针旋转使变换后的因子具有清晰的可解释性及唯一性，当针对两个因子进行平面正交旋转时，因子载荷矩阵C具有如下形式：C＝(cij),i＝1,2,...,n,j＝1,2；
    取正交矩阵如下：
    则D＝CS＝(dij),i＝1,2,...,n,j＝1,2)为旋转因子载荷矩阵，模型变为如下形式：X＝D(STF)+QU；
    此时公共因子由F变为QTF，当前旋转双因子的目的是将变量分为由不同因子说明的两个部分，故需要将(d11 2,d21 2,...,dn1 2)和(d12 2,d22 2,...,dn2 2)的方差最大化，二者的相对方差可以表示为：以消除dij符号的影响，消除各个属性变量对公共因子依赖程度不同的影响；
    接下来将总方差V＝V1+V2最大化，令经计算，φ满足：其中：
    对于p个公共因子，将这些因子两两进行旋转；在一遍完整的旋转过程中，从p个因子中每次选两个进行逆时针旋转，共有p(p-1)/2次旋转，完成一遍旋转后，继续进行下一遍循环，在此重复过程中C各列的相对方差和总方差V会递增，直至某一次循环的总方差比前和一次循环的总方差相比差距不大于5％则停止旋转，此时得到最终旋转因子载荷矩阵D，和标准正交化特征向量相乘即得到模糊综合评价的因素权向量：A＝DηT，A实际上反映了u中属性的重要程度；
    然后通过模糊变换将反映u的模糊向量A转换成反映v的模糊向量B，即B＝AR＝(b1,b2,...,bm)，模糊向量B中的元素与评价集v中的元素一一对应，最大的元素对应的评价即为卷烟的模糊综合最终评价等级。
  2. 根据权利要求1所述的卷烟燃烧外观评价方法，其特征在于，步骤S1中，基于模拟人体卷烟抽吸全过程的机械手臂，通过摄像系统捕捉卷烟实时抽吸状态，采集卷烟外在属性。
  3. 根据权利要求1所述的卷烟燃烧外观评价方法，其特征在于，步骤S3中，m＝4，按照优劣排序，v1是特等级，v2是优等级，v3是中等级，v4是次等级。
  4. 根据权利要求1所述的卷烟燃烧外观评价方法，其特征在于，步骤S4中，rij由专家打分获得。
PCT/CN2023/097776 2022-06-10 2023-06-01 一种卷烟燃烧外观的评价方法 WO2023236842A1 (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202210655706.6A CN115078380A (zh) 2022-06-10 2022-06-10 一种卷烟燃烧外观的评价方法
CN202210655706.6 2022-06-10

Publications (1)

Publication Number Publication Date
WO2023236842A1 true WO2023236842A1 (zh) 2023-12-14

Family

ID=83250795

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2023/097776 WO2023236842A1 (zh) 2022-06-10 2023-06-01 一种卷烟燃烧外观的评价方法

Country Status (2)

Country Link
CN (1) CN115078380A (zh)
WO (1) WO2023236842A1 (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115078380A (zh) * 2022-06-10 2022-09-20 云南中烟工业有限责任公司 一种卷烟燃烧外观的评价方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103577681A (zh) * 2013-06-26 2014-02-12 长沙理工大学 基于因子分析锅炉效率影响指标的定量评价方法
KR101731465B1 (ko) * 2015-12-17 2017-04-28 태영산업 주식회사 담뱃재 낙하 측정장치
CN111426795A (zh) * 2020-04-26 2020-07-17 云南中烟工业有限责任公司 一种基于机械手仿真人体卷烟抽吸全过程及全视觉测定表征卷烟燃烧包灰性能的方法
CN112862356A (zh) * 2021-03-08 2021-05-28 东南大学 基于数理方法构建传统村落活态化保护评价体系的方法
CN114062589A (zh) * 2021-11-23 2022-02-18 云南中烟工业有限责任公司 一种卷烟灰烬特征的评价方法
CN115078380A (zh) * 2022-06-10 2022-09-20 云南中烟工业有限责任公司 一种卷烟燃烧外观的评价方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103577681A (zh) * 2013-06-26 2014-02-12 长沙理工大学 基于因子分析锅炉效率影响指标的定量评价方法
KR101731465B1 (ko) * 2015-12-17 2017-04-28 태영산업 주식회사 담뱃재 낙하 측정장치
CN111426795A (zh) * 2020-04-26 2020-07-17 云南中烟工业有限责任公司 一种基于机械手仿真人体卷烟抽吸全过程及全视觉测定表征卷烟燃烧包灰性能的方法
CN112862356A (zh) * 2021-03-08 2021-05-28 东南大学 基于数理方法构建传统村落活态化保护评价体系的方法
CN114062589A (zh) * 2021-11-23 2022-02-18 云南中烟工业有限责任公司 一种卷烟灰烬特征的评价方法
CN115078380A (zh) * 2022-06-10 2022-09-20 云南中烟工业有限责任公司 一种卷烟燃烧外观的评价方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
GU DELIANG, LIU XIN: "Research on Fuzzy Comprehensive Evaluation of Cigarette Sensory Quality", CHINA BUSINESS [ECONOMIC THEORY RESEARCH], no. 2, 28 February 2009 (2009-02-28), pages 33 - 34, XP009551227, ISSN: 1009-5292, DOI: 10.16834/j.cnki.issn1009-5292.2009.04.015 *
JIANG JIANGANG, QI SIHUA: "Application of Synthetic the Fuzzy Evaluating in Cigarettes Sampling Process", JILIANG-YU-CESHI-JISHU = METROLOGY AND MEASUREMENT TECHNIQUE, JILIANG YU CESHI JISHU ZASHISHE, CN, no. 2, 31 December 2002 (2002-12-31), CN , pages 32 - 33, XP009550952, ISSN: 1004-6941 *
SUN LI-CHENG, ZHOU DE-QUN, MEI QIANG: "The Evaluation of the Comprehensive Performance of Oil Industry Listed Companies", JOURNAL OF TECHNICAL ECONOMY AND MANAGEMENT, no. 7, 26 July 2012 (2012-07-26), pages 3 - 6, XP009550953, ISSN: 1004-292X *
ZHU, FENGMEI: "Research on Evaluation of Medical and Health System Reform Based on Fuzzy Comprehensive Evaluation Method", CHINESE JOURNAL OF HEALTH STATISTICS, vol. 33, no. 2, 25 April 2016 (2016-04-25), pages 267 - 270, XP009550954, ISSN: 1002-3674 *

Also Published As

Publication number Publication date
CN115078380A (zh) 2022-09-20

Similar Documents

Publication Publication Date Title
WO2023236842A1 (zh) 一种卷烟燃烧外观的评价方法
Zhang et al. Sound quality prediction of vehicle interior noise and mathematical modeling using a back propagation neural network (BPNN) based on particle swarm optimization (PSO)
CN110133049B (zh) 基于电子鼻和机器视觉的茶叶等级快速无损检测方法
CN109063366B (zh) 一种基于时间和空间加权的建筑性能数据在线预处理方法
CN105678332A (zh) 火焰图像cnn识别建模的转炉炼钢终点判断方法及系统
WO2024007781A1 (zh) 一种卷烟燃烧质量指标重要性评价综合分析方法
CN109784348A (zh) 一种红外电力设备识别和在线诊断方法及其系统
CN104713835A (zh) 一种烟叶颜色在线数值化识别方法
CN110727223A (zh) 基于井下采煤工作面的环式轨道智能巡检机器人及其应用
CN101685092B (zh) 用烟叶pH检测值判定烤烟醇化质量的方法
CN113252839B (zh) 基于卷烟纸参数的卷烟静态包灰综合评价指标测定方法
CN110309609A (zh) 一种基于粗糙集和小波神经网络的建筑室内空气品质评价方法
CN111474137A (zh) 一种基于卷积神经网络的柑橘糖度无损检测方法
CN116610998A (zh) 一种基于多模态数据融合的开关柜故障诊断方法和系统
CN111694827A (zh) 一种电力设备状态监测数据缺失值分类插补方法和系统
CN110082106A (zh) 一种基于Yu范数的深度度量学习的轴承故障诊断方法
Gao et al. Physical fitness clustering analysis based on self-organizing feature maps network
CN113222268B (zh) 一种基于多模式推理的烟草烘烤质量预测模型建立方法
CN114596296A (zh) 一种高灵敏度的热轧钢卷端面缺陷识别系统及方法
TW202238527A (zh) 基於頻譜影像之空氣汙染檢測方法
Cai et al. Composition analysis and identification of glass products based on Pearson correlation analysis
CN105699302A (zh) 一种利用光谱法测定苹果主要营养成分的方法
CN112485372A (zh) 一种评价烟气中杂气的方法
CN112488869B (zh) 一种适用于农业生产多因素管理决策系统
Cui et al. Using PCA and ANN to identify significant factors and modeling customer satisfaction for the complex service processes

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 23819007

Country of ref document: EP

Kind code of ref document: A1