一种基于历史数据的卷烟叶组配方优选方法、设备及存储介质组成比例

本发明专利技术公开了一种基于历史数据的卷烟叶组配方优选方法、设备及存储介质，在评价卷烟叶组配方稳定性和相似性中，基于误差统计及二维高斯分布的联合密度函数，通过对数据进行降维处理并引入误差椭圆约束方程构建计算待选配方在历史配方中偏离度指数，并根据偏离度指数从待选配方中进行优选，用于解决现有技术无法基于叶组配方历史数据指标的量化待用配方的稳定性和相似性的问题。本发明专利技术可应用于目前公司卷烟叶组配方设计中待用叶组配方与历史配方的相似性及符合度的评价工作，在提高卷烟叶组配方设计效率的同时，可推广到潜在烟叶资源挖掘、烟叶保障资源评估以及卷烟产品叶组配方设计的实务工作中。方设计的实务工作中。方设计的实务工作中。

【技术实现步骤摘要】
一种基于历史数据的卷烟叶组配方优选方法、设备及存储介质


[0001]本专利技术涉及一种基于历史数据的卷烟叶组配方优选方法、设备及存储介质，属于卷烟产品设计


技术介绍

[0002]卷烟叶组配方设计是形成卷烟产品风格特点的前提，而配方稳定性则是保障产品质量的基本要求。烟草作为一种农产品，品质波动变化较大，受到年份、气候、地域、醇化时间等因素的影响。叶组配方设计是不同年份、产地、香型、部位烟叶原料之间的优化组合。在实际的叶组配方设计中，配方设计同时兼顾成本、库存、产量等多种因素的考量，配方中烟叶的更替及配方结构的调整均要最大限度地保持配方的稳定性和相似性。
[0003]当前，卷烟叶组配方设计与维护仍存在评价指标量化难、设计效率低等问题，极大限制了烟叶资源的利用效率及潜在资源的挖掘。传统的评价方式主要参考客观定量指标，如均值、变异系数、极差等常规统计参数，同时参照少数指标的波动或单项指标阈值反映配方的变化，进而预判待用配方的可行性，缺乏一种客观、综合型的指标。与此同时，在卷烟设计过程中，大量的历史配方表征数据则是该配方的风格品质特征的集中体现。
[0004]因此，如何有效结合历史配方表征数据，提供一种衡量待用叶组配方在历史配方中的偏离度，并预测出待用叶组配方的方法是一个亟需解决的难题。

技术实现思路

[0005]目的：为了克服现有技术中存在的无法基于叶组配方历史表征数据量化待用配方稳定性和相似性的问题。本专利技术提供一种基于历史数据的卷烟叶组配方优选方法、设备及存储介质。
[0006]技术方案：为解决上述技术问题，本专利技术采用的技术方案为：
[0007]一种基于历史数据的卷烟叶组配方优选方法，包括如下步骤：
[0008]根据已有卷烟配方的历史表征数据，获取历史配方的第l对综合变量数据集。
[0009]根据历史配方的第l对综合变量数据集，计算历史配方的第l对综合变量数据集的误差椭圆长、短半轴和倾斜角。
[0010]根据历史配方的第l对综合变量数据集的误差椭圆长、短半轴，计算历史配方的第l对综合变量数据集误差椭圆面积。
[0011]根据待选配方的第l对综合变量，历史配方的第l对综合变量数据集的误差椭圆长、短半轴和倾斜角，计算待选配方的误差椭圆的长半轴及短半轴。
[0012]根据待选配方的误差椭圆的长半轴及短半轴，计算待选配方的第l对综合变量误差椭圆面积。
[0013]根据待选配方的第l对综合变量误差椭圆面积和历史配方的第l对综合变量数据集误差椭圆面积，计算待选配方的第l对综合变量偏离度值。
[0014]根据待选配方的第l对综合变量偏离度值，计算待选配方的偏离度指数。
[0015]计算所有待选配方的偏离度指数，依据偏离度指数数值进行排序，选取所有待选配方中偏离度指数较小的待选配方作为最终的卷烟配方优选方案。
[0016]进一步的，还包括：对已有卷烟配方的历史表征数据进行预处理，所述预处理包括：空值处理、异常值删除；空值即缺失值，采用平均值替换缺失值；异常值处理，即对采集的数据，采用Tukey
’
s test方法对异常值进行直接删除。
[0017]进一步的，还包括：将历史配方的第l对综合变量数据集的误差椭圆和待选配方的第l对综合变量误差椭圆通过R软件进行可视化处理。
[0018]进一步的，根据已有卷烟配方的历史表征数据，获取历史配方的第l对综合变量数据集，包括如下步骤：
[0019]根据已有卷烟配方的历史表征数据X，计算其协方差矩阵M。
[0020]其中,
[0021]其中，x
ij
为第i个历史配方的第j个表征数据指标，i＝1,2,3...n；j＝1,2,3...k，n为参数输入的历史配方表征数据个数，k为参数输入的评价指标个数。
[0022]计算协方差矩阵M的特征值F＝[λ1，λ2，λ3，...，λ
k
]和第j个特征值对应的特征向量U
j
＝[l
j1
，l
j2
，l
j3
，...，l
jk
]T
。
[0023]根据第j个特征值对应的特征向量U
j
，计算第i个历史配方综合变量。
[0024]计算所有历史配方的历史配方综合变量，得到历史配方的数据集N。
[0025]其中，
[0026]其中，N
il
为第i个历史配方所选取的第l对综合变量i＝1,2,3...n，l＝1,2,3...k/2。
[0027]根据历史配方的数据集N，选取计算第l对综合变量数据集[N
1l
，N
2l
，N3
l，
...，N
nl
]T
。
[0028][0029]其中，Dim
i2l
‑1，Dim
i2l
分别表示第l对综合变量中的第i个历史配方的第2l
‑
1，2l综合变量。
[0030]进一步的，所述根据历史配方的第l对综合变量数据集，计算历史配方的第l对综合变量数据集的误差椭圆长、短半轴和倾斜角，包括：
[0031]根据第l对综合变量数据集[N
1l
，N
2l
，N
3l
，...，N
n
l]T
计算二维协方差矩阵A
l
。
[0032]其中，
[0033]其中，σ
xx
σ
yy
分别表示历史配方的第l对综合变量数据集的第1、2列综合变量的方差。σ
xy
，σ
yx
分别表示历史配方的第l对综合变量数据集第1、2列综合变量的协方差。
[0034]根据二维协方差矩阵A
l
计算两个特征值为β1，β2及椭圆倾斜角θ。
[0035]其中，
[0036][0037][0038]根据特征值为β1，β2，计算历史配方的第l对综合变量数据集的误差椭圆长、短半轴a
l
、b
l
。
[0039]其中，a
l
＝sβ1，b
l
＝sβ2。
[0040]其中，s为常数。
[0041]进一步的，所述根据历史配方的第l对综合变量数据集的误差椭圆长、短半轴，计算历史配方的第l对综合变量数据集误差椭圆面积，包括：
[0042]S
Nl
＝πa
l
b
l
[0043]其中，S
Nl
表示历史配方的第l对综合变量数据集误差椭圆面积，π为圆周率。
[0044]进一步的，所述根据待选配方的第l对综合变量，历史配方的第l对综合变量数据集的误差椭圆长、短半轴和倾斜角，计算待选配方的误差椭圆的长半轴及短半轴，包括：...

【技术保护点】

【技术特征摘要】
1.一种基于历史数据的卷烟叶组配方优选方法，其特征在于：包括如下步骤：根据已有卷烟配方的历史表征数据，获取历史配方的第1对综合变量数据集；根据历史配方的第1对综合变量数据集，计算历史配方的第1对综合变量数据集的误差椭圆长、短半轴和倾斜角；根据历史配方的第1对综合变量数据集的误差椭圆长、短半轴，计算历史配方的第1对综合变量数据集误差椭圆面积；根据待选配方的第1对综合变量，历史配方的第1对综合变量数据集的误差椭圆长、短半轴和倾斜角，计算待选配方的误差椭圆的长半轴及短半轴；根据待选配方的误差椭圆的长半轴及短半轴，计算待选配方的第1对综合变量误差椭圆面积；根据待选配方的第1对综合变量误差椭圆面积和历史配方的第1对综合变量数据集误差椭圆面积，计算待选配方的第1对综合变量偏离度值；根据待选配方的第1对综合变量偏离度值，计算待选配方的偏离度指数；计算所有待选配方的偏离度指数，依据偏离度指数数值进行排序，选取所有待选配方中偏离度指数较小的待选配方作为最终的卷烟配方优选方案。2.根据权利要求1所述的一种基于历史数据的卷烟叶组配方优选方法，其特征在于：还包括：对已有卷烟配方的历史表征数据进行预处理，所述预处理包括：空值处理、异常值删除；空值即缺失值，采用平均值替换缺失值；异常值处理，即对采集的数据，采用Tukey
’
s test方法对异常值进行直接删除。3.根据权利要求1所述的一种基于历史数据的卷烟叶组配方优选方法，其特征在于：还包括：将历史配方的第1对综合变量数据集的误差椭圆和待选配方的第1对综合变量误差椭圆通过R软件进行可视化处理。4.根据权利要求1
‑
3任一项所述的一种基于历史数据的卷烟叶组配方优选方法，其特征在于：根据已有卷烟配方的历史表征数据，获取历史配方的第1对综合变量数据集，包括如下步骤：根据已有卷烟配方的历史表征数据X，计算其协方差矩阵M；其中，其中，x
ij
为第i个历史配方的第j个表征数据指标，i＝1,2，3...n；j＝1，2，3...k，n为参数输入的历史配方表征数据个数，k为参数输入的评价指标个数；计算协方差矩阵M的特征值F＝[λ1，λ2，λ3...，λ
k
]和第j个特征值对应的特征向量U
j
＝[l
j1
，l
j2
，l
j3
，...l
jk
]
T
；根据第j个特征值对应的特征向量U
j
，计算第i个历史配方综合变量；计算所有历史配方的历史配方综合变量，得到历史配方的数据集N；其中，
其中，N
i1
为第i个历史配方所选取的第1对综合变量i＝1，2，3...n，l＝1，2，3...k/2；根据历史配方的数据集N，选取计算第1对综合变量数据集[N
1l
，N
2l
，N
3l
，...，N
nl
]
T
；其中，Dim
i2l
‑1，Dim
i2l
分别表示第1对综合变量中的第i个历史配方的第21
‑
1，21综合变量。5.根据权利要求1
‑
3任一项所述的一种基于历史数据的卷烟叶组配方优选方法，其特征在于：所述根据历史配方的第1对综合变量数据集，计算历史配方的第1对综合变量数据集的误差椭圆长、短半轴和倾斜角，包括：根据第1对综合变量数据集[N
1l
，N
2l
，N
3l
，...，N
nl
]
T
计算二维协方差矩阵A
l
；其中，其中，σ
xx σ
yy
分别表示历史配方的第1对综合变量数据集的第1、2列综合变量的方差；σ
xy
，σ
yx
分别表示历史配方的第1对综合变量数据集第1、2列综合变量的协方差；根据二维协方差矩阵A
l
计算两个特征值为β1，β
２
及椭圆倾斜角θ；其中，其中，其中，根据特征值为β1，β2，计算历史配方的第1对综合变量数据集的误差椭圆长、短半轴a
l
、b
l
；其中，a
l
＝sβ1，b
l
＝sβ2；其中，s为常数。6.根据权利要求1
‑
3所述的一种基于历史数据的卷烟叶组配方优选方法，其特征在于：所述根据历史配方的第1对综合变量数据集的误差椭圆长、短半轴，计算历史配方的第1对综合变量数据集误差椭圆面积，包括：S
Nl
＝πa
l
b
l
其中，S

【专利技术属性】
技术研发人员：赵振杰，何文苗，毕一鸣，王辉，吴继忠，廖付，李海峰，宣润泉，李石头，李永生，李炳华，郝贤伟，田雨农，
申请(专利权)人：浙江中烟工业有限责任公司，
类型：发明
国别省市：