CN109961315A - 一种基于非线性组合模型的卷烟月度销量预测方法 - Google Patents
一种基于非线性组合模型的卷烟月度销量预测方法 Download PDFInfo
- Publication number
- CN109961315A CN109961315A CN201910085257.4A CN201910085257A CN109961315A CN 109961315 A CN109961315 A CN 109961315A CN 201910085257 A CN201910085257 A CN 201910085257A CN 109961315 A CN109961315 A CN 109961315A
- Authority
- CN
- China
- Prior art keywords
- model
- cigarette
- data
- monthly
- neural network
- 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
- 235000019504 cigarettes Nutrition 0.000 title claims abstract description 83
- 238000000034 method Methods 0.000 title claims abstract description 56
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 51
- 238000003062 neural network model Methods 0.000 claims abstract description 37
- 238000012549 training Methods 0.000 claims abstract description 30
- YHXISWVBGDMDLQ-UHFFFAOYSA-N moclobemide Chemical compound C1=CC(Cl)=CC=C1C(=O)NCCN1CCOCC1 YHXISWVBGDMDLQ-UHFFFAOYSA-N 0.000 claims abstract description 23
- 238000013528 artificial neural network Methods 0.000 claims description 19
- 210000002569 neuron Anatomy 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 9
- 238000010276 construction Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 7
- 230000001537 neural effect Effects 0.000 claims description 6
- 238000000714 time series forecasting Methods 0.000 claims description 6
- 210000005036 nerve Anatomy 0.000 claims description 5
- 230000001932 seasonal effect Effects 0.000 claims description 3
- 230000006978 adaptation Effects 0.000 claims 1
- 230000006872 improvement Effects 0.000 abstract description 10
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 abstract description 6
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 abstract description 6
- 239000000284 extract Substances 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 abstract description 3
- 239000010410 layer Substances 0.000 description 37
- 230000006870 function Effects 0.000 description 11
- 230000000694 effects Effects 0.000 description 9
- 238000012360 testing method Methods 0.000 description 5
- 238000011478 gradient descent method Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 230000004913 activation Effects 0.000 description 2
- 230000006399 behavior Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- FNMKZDDKPDBYJM-UHFFFAOYSA-N 3-(1,3-benzodioxol-5-yl)-7-(3-methylbut-2-enoxy)chromen-4-one Chemical compound C1=C2OCOC2=CC(C2=COC=3C(C2=O)=CC=C(C=3)OCC=C(C)C)=C1 FNMKZDDKPDBYJM-UHFFFAOYSA-N 0.000 description 1
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 210000004218 nerve net Anatomy 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
- 238000012892 rational function Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/044—Recurrent networks, e.g. Hopfield networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q30/00—Commerce
- G06Q30/02—Marketing; Price estimation or determination; Fundraising
- G06Q30/0201—Market modelling; Market analysis; Collecting market data
- G06Q30/0202—Market predictions or forecasting for commercial activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/04—Manufacturing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- Finance (AREA)
- Development Economics (AREA)
- General Health & Medical Sciences (AREA)
- Accounting & Taxation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computational Linguistics (AREA)
- Biomedical Technology (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Economics (AREA)
- Biophysics (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Computing Systems (AREA)
- Evolutionary Computation (AREA)
- General Business, Economics & Management (AREA)
- Manufacturing & Machinery (AREA)
- Tourism & Hospitality (AREA)
- Primary Health Care (AREA)
- Human Resources & Organizations (AREA)
- Game Theory and Decision Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于非线性组合模型的卷烟月度销量预测方法,包括构建一阶模型以及构建组合模型,通过构建的组合模型进行准确预测。本发明通过利用ARIMA模型、基于梯度下降算法的BP神经网络模型、基于Levenberg‑Marquardt算法改进的BP神经网络模型三个单项模型对各月卷烟销量数据进行拟合并预测;然后采用集成学习结合策略中的学习法,提出建立ARIMA模型、基于梯度下降算法的BP神经网络模型、基于Levenberg‑Marquardt算法改进的BP神经网络模型的非线性组合模型,对历年各月卷烟销量数据进行仿真训练,在大量数据中找到销量数据之间的关系,以期预测未来几个月的卷烟销量。该模型充分利用上述3种单个模型的优点,有效提取数据各方面信息,其预测误差最为平稳,适应性最强,更加贴合卷烟销售实际情况。
Description
技术领域
本发明涉及卷烟生产预测技术领域,尤其涉及一种基于非线性组合模型的卷烟月度销量预测方法。
背景技术
目前,传统卷烟预测模式立足于市场实际,把握市场的能力不断加强,但由于受人为因素影响较大,预测准确率难以有效提升。目前卷烟预测研究分为四类:第一类为以年度数据为基础,以宏观经济数据为自变量、卷烟消费数据为因变量,建立回归方程,分析各指标的相关性,进而对卷烟销量进行预测,多见于卷烟需求预测方法初探阶段,预测误差较大;第二类为线性时间序列模型,通过建立单一的线性时间序列模型来预测卷烟销量,由于时间序列模型仅利用数据自身的历史时间序列信息,获得过去行为的有关结论,进而对未来时间序列行为进行推断,突出了时间因素在预测中的作用,暂不考虑外界具体因素的影响,具有获取成本低、建模时间快等优点,但预测精度有待进一步提升;第三类为单个神经网络预测模型,通过对数据集的建模训练,构建网络,进而预测未来月份的数据,该方法根据选择算法不同,其预测精度也会不同;第四类为多个单一预测模型进行组合,给予每个预测模型一定的权重,从而完成组合模型的构建,有效提升了模型的预测精度,但各模型的权重没有固定的模式,进而影响预测结果。
发明内容
本发明的目的是提供一种基于非线性组合模型的卷烟月度销量预测方法,能够有效弥补上述模型弊端,准确预测各月卷烟销量。
本发明采用的技术方案为:
一种基于非线性组合模型的卷烟月度销量预测方法,包括如下步骤:
A、根据卷烟N个月度销量数据受节日因素较大、季节性较强的规律特征,选择多种适应性较强的时间序列预测方法;
B、通过输入卷烟N个月度销售数据端分别通过多种适应性较强的时间序列预测方法对卷烟N个月度销售数据进行建模和训练,分别得到对应的模型,记为一阶模型;
C:选取一阶模型中其中一个模型作为二阶模型,以一阶模型的输出数据作为二阶模型的输入数据,然后对一阶模型输入N个月度销量数据,一阶模型得到的对应的多个拟合值同时作为输入向量,以N月度相对应的真实销量作为输出向量,对二阶模型进行训练,从而得到一阶模型和二阶模型的组合模型;
D:对组合模型输入新的月度销售数据,即可输出得到最新的月度销售预测值。所述的一阶模型分别包括有ARIMA模型、基于梯度下降算法的BP神经网络模型和基于Levenberg-Marquardt的BP神经网络模型,二阶模型为基于Levenberg-Marquardt的BP神经网络模型。
所述的步骤B分别包括有如下步骤:
B1,选择ARIMA模型对卷烟月度销售数据进行差分处理、模型识别与定阶、仿真等过程,进而预测未来月份的卷烟销量;
B2,选择基于梯度下降算法的BP神经网络模型,对卷烟月度销售数据进行处理,确定输入层、输出层,创建基于梯度下降算法的网络,编制神经元测试程序,确定神经元数量;设定网络参数,经过训练,完成模型构建;选择预测输入层数据,从而预测未来月份卷烟销量;
B3,选择基于Levenberg-Marquardt的BP神经网络模型,对卷烟月度销售数据进行处理,确定输入层、输出层,创建基于LM算法的神经网络;编制神经元测试程序,确定神经元数量;设定网络参数,经过训练,完成模型构建;选择预测输入层数据,从而预测未来月份卷烟销量。
本发明通过利用ARIMA模型、基于梯度下降算法的BP神经网络模型、基于Levenberg-Marquardt算法改进的BP神经网络模型三个单项模型对各月卷烟销量数据进行拟合并预测;然后采用集成学习结合策略中的学习法,提出建立ARIMA模型、基于梯度下降算法的BP神经网络模型、基于Levenberg-Marquardt算法改进的BP神经网络模型的非线性组合模型,对历年各月卷烟销量数据进行仿真训练,在大量数据中找到销量数据之间的关系,以期预测未来几个月的卷烟销量。该模型充分利用上述3种单个模型的优点,有效提取数据各方面信息,其预测误差最为平稳,适应性最强,更加贴合卷烟销售实际情况。
附图说明
图1本发明的建模流程图。
图2ARIMA模型对卷烟月度销售数据的拟合及预测图。
图3基于梯度下降算法的BP神经网络模型对卷烟月度销售数据的拟合及预测图。
图4基于LM算法的BP神经网络模型对卷烟月度销售数据的拟合及预测图。
图5非线性组合预测模型对卷烟月度销售数据的拟合及预测图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1、2和3所示,本发明包括如下步骤:
A、根据卷烟N个月度销量数据受节日因素较大、季节性较强的规律特征,选择多种适应性较强的时间序列预测方法;
B、通过输入卷烟N个月度销售数据端分别通过多种适应性较强的时间序列预测方法对卷烟N个月度销售数据进行建模和训练,分别得到对应的模型,记为一阶模型;所述的一阶模型分别包括有ARIMA模型、基于梯度下降算法的BP神经网络模型和基于Levenberg-Marquardt的BP神经网络模型,二阶模型为基于Levenberg-Marquardt的BP神经网络模型。所述的步骤B分别包括有如下步骤:
B1,选择ARIMA模型对卷烟月度销售数据进行差分处理、模型识别与定阶、仿真等过程,进而预测未来月份的卷烟销量;
B2,选择基于梯度下降算法的BP神经网络模型,对卷烟月度销售数据进行处理,确定输入层、输出层,创建基于梯度下降算法的网络,编制神经元测试程序,确定神经元数量;设定网络参数,经过训练,完成模型构建;选择预测输入层数据,从而预测未来月份卷烟销量;
B3,选择基于Levenberg-Marquardt的BP神经网络模型,对卷烟月度销售数据进行处理,确定输入层、输出层,创建基于LM算法的神经网络;编制神经元测试程序,确定神经元数量;设定网络参数,经过训练,完成模型构建;选择预测输入层数据,从而预测未来月份卷烟销量。
C:选取一阶模型中其中一个模型作为二阶模型,以一阶模型的输出数据作为二阶模型的输入数据,然后对一阶模型输入N个月度销量数据,一阶模型得到的对应的多个拟合值同时作为输入向量,以N月度相对应的真实销量作为输出向量,对二阶模型进行训练,从而得到一阶模型和二阶模型的组合模型;
具体的,采用集成学习结合策略中的学习法,以ARIMA模型拟合值、基于梯度下降算法的BP神经网络拟合值、基于LM算法的BP神经网络拟合值作为输入向量,即选取相同月份的3种模型拟合值作为输入向量,当月真实销量作为输出向量,选择基于LM算法的BP神经网络构造非线性组合模型,编制神经元测试程序,确定神经元数量;设定网络参数,进行训练后进行仿真拟合;输入上述3种模型对预测月份卷烟销量的预测值,得到组合预测模型对销量的预测结果。通过测试,其拟合及预测结果如图5所示。
D:对组合模型输入新的月度销售数据,即可输出得到最新的月度销售预测值。
具体的以实际数据进行举例说明:ARIMA模型是由Box和Jenkins提出的一种时间预测方法,该方法将原始时间序列数据看作一个随机序列,并用一定的数学模型来近似描述,一旦被识别后,就可以利用该模型根据历史数据来预测未来。本发明以2006~2017年共144个月的全国卷烟销量数据(万箱)作为基础数据,如图2所示,应用MATLAB对其进行ARIMA拟合,并预测2018年1~4月的销量。
从销量趋势图中可以看出,该时间序列非平稳,对数列先进行1阶差分处理,将其平稳化,其公式为:
Δxt=xt+1-xt 式(1)
一阶差分的时间序列的均值和方差均基本平稳,接近一个白噪声序列。通过平稳性检验可知,即将建立的ARIMA(p,d,q)模型,其中d=1,其基本公式为:
xt=θ+α1xt-1+α2xt-2+...+αpxt-p+ut+β1ut-1+β2ut-2+...+βqut-q 式(2)
其中,x为时间序列,θ是常数项,α为自回归系数,β为移动平均系数,u为白噪声,p为自回归项,q为移动平均项。
为确认模型中p和q的值,以AIC准则为评价该时间序列数据的准则,利用MATLAB进行测试。本次基础数据为144个,故将p和q的范围定为1~29之间,通过测试确定AIC=7.06,p=25,q=7,即选择的模型为ARIMA(25,1,7)。
将p和q的值代入模型,计算得到2006~2017年各月卷烟销量的拟合值(如图2所示),模型拟合度为76.07%,平均绝对误差为14.82,平均绝对误差率(平均绝对误差/平均真实销量,下同)为3.81%,拟合精度整体良好。
将预测长度设定为4个月,利用该模型对2018年1~4月全国卷烟销量进行预测,如表1所示,可以看出预测误差均较大,除1月份外,其余3个月的绝对误差值均大于当月的日均销量。
表1 ARIMA模型对2018年1~4月全国卷烟销量预测表(万箱,%,天)
注:2018年,1月共22个工作日、2月共17个工作日、3月共22个工作日、4月共20个工作日。
BP神经网络是Rumelhart和McClellend于1986年提出一种多层网络的“逆推”学习算法,是一种单向传播的多层前向网络,具有3层或3层以上,包括输入层、隐含层和输出层,其中隐含层可以有多个层次,上下层之间的神经元全部连接,但每层神经元之间相互独立没有连接。当学习样本进入神经网络后,神经元的激活值从输入层经过各隐含层传向输出层,若达不到期望,则进入误差反向传播过程,逐次调整网络各层的权值和偏置,最后回到输入层,再重复计算,通过不断调整各层的权值和阈值,从而达到网络误差最小,学习过程结束;再利用训练好的网络,进行预测。
BP神经网络通过训练误差来逐步调整各层间的输入权重和偏置,这个调整过程的算法依据为最基本的梯度下降法(Gradient Descent),由于训练误差是关于输入权重和偏置的二次函数,分别对权重和偏置求偏导数,即梯度向量,沿着梯度向量的方向,是训练误差增加最快的地方,而沿着梯度向量相反的方向,梯度减少最快,在这个方向上更容易找到训练误差函数(损失函数)的最小值。本节利用MATLAB软件应用梯度下降算法的BP神经网络对全国卷烟历史销量数据进行训练并预测,其处理步骤如下:
1.样本数据确定。
由于卷烟销量执行的是年度计划,为更好的反应卷烟销量的时间特点,本文将2006至2017年每12个月销量数据作为一组输入向量,形成12*129的输入矩阵;以之后4个月的销量数据作为网络的输出,形成4*129的输出矩阵,用以神经网络进行训练。
2.样本数据处理。
数据的预处理和后处理是训练神经网络的关键步骤,它直接影响到训练后神经网络的性能,常见的方法是将原始数据进行归一化处理,即通过一定的线性变换将输入数据和输出数据统一限制在[0,1]或[-1,1]区间内。本节选择最大最小法进行处理,其函数形式是:
其中,xmin为数据序列的最小值,xmax为数据序列的最大值。最大最小值法可以通过MATLAB的mapminmax函数对数据进行处理,让其落入在[-1,1]区间内。
3.创建网络及参数设定。
本模型设计主要涉及到网络的层数、每层中的神经元个数和激活函数、初始值以及学习速率等几个方面。
(1)网络层数确定。理论证明,具有偏差和至少一个S型隐层加上一个线性输出层的网络,能够逼近任何有理函数,增加层数可以进一步降低误差,提高精度,但同时也是网络复杂化。本节模型选择3层网络,即1个输入层、1个隐含层、1个输出层。
(2)创建网络。在MATLAB中创建神经网络要确定以下函数:①输入层和隐含层的传递函数采用正切S型函数,以tansig表示,即:
②输出层的传递函数采用线性函数,以purelin表示;③网络训练函数采用动量反转和动态自适应学习率的梯度下降算法训练函数tianingd;④确定隐含层神经元个数,理论上没有计算出多少隐含层节点最合适,在具体设计时,要通过对不同神经元数进行训练比较,以网络最小误差为依据确定神经元数。通过循环计算,确定此模型的神经元个数为9个。
(3)网络训练参数设定。①根据梯度下降法的特点,训练次数选取为80000;②误差性能目标值选取为0.01;③学习速率一般取值范围为0.01-0.8之间,学习速率太大可能导致系统的不稳定,但学习速率太小导致收敛太慢,需要较长的训练时间,为达到目标误差精度,此模型选取0.01。
利用基于梯度下降算法的BP神经网络预测模型对卷烟销量数据进行训练和仿真,经过79976次训练之后达到目标误差0.01的要求,模型拟合效果如图3所示,其平均绝对误差为17.87,平均绝对误差率为4.59%,其拟合精度不及ARIMA模型。
输入2017年12个月的卷烟销量数据,模型得到2018年1~4月销量的预测值,如表2所示,可以看出虽然该模型拟合效果不及ARIMA模型,但预测结果要明显好于ARIMA模型,除2月份由于春节因素导致预测效果较差外,其余月份预测结果十分逼近真实销量,均不足半天的当月卷烟销量。
表2基于梯度下降算法的BP神经网络模型对2018年1~4月全国卷烟销量预测表(万箱,%,天)
注:2018年,1月共22个工作日、2月共17个工作日、3月共22个工作日、4月共20个工作日。
从上述基于梯度下降算法的BP神经网络模型训练过程发现,虽然该模型能较好的拟合卷烟销量并比较准确的预测未来4个月的销量,但其学习算法的收敛速度较慢,训练次数达到了近8万次,另外其预测误差精度较低,仅能达到0.01。针对该模型存在的不足,出现了几种基于该模型的改进算法,如拟牛顿法、Levenberg-Marquardt算法等。通过实验发现,对中小规模的神经网络,使用Levenberg-Marquardt算法(简称LM算法)的收敛速度最快,且计算精度较高。作为BP神经网络的改进算法,LM算法实际上是梯度下降法和高斯牛顿法的结合,该算法期望在不计算Hessian矩阵的情况下获得高阶的训练速度,其公式表达为:
xk+1=xk-[JTJ+μI]-1JTe 式(5)
其中,JT为雅克比矩阵,e是网络误差向量。随着μ越来越小,该算法越来越接近高斯牛顿法;随着μ越来越大,该算法越来越接近梯度下降法。由于高斯牛顿法收敛速度更快更准确,因此LM算法的目的就是尽快接近高斯牛顿法,如果某次迭代成功,误差减小,则减小μ的值;如果迭代失败,则增加μ的值,从而使得误差性能函数随着迭代的进行而下降到极小值。
基于LM算法改进BP神经网络模型构建过程:
基于LM算法改进BP神经网络模型的输入向量和输出向量与基于梯度下降算法的BP神经网络模型相同,并进行归一化的数据处理,确定网络层数,选择传递函数,设定训练次数、训练精度、学习速率等参数。但与前者不同的是,基于LM算法改进的BP神经网络模型需要将模型的训练函数更改成LM算法,在MATLAB中以trianlm表示。鉴于LM算法快速、准确收敛的特点,将模型训练次数取为3000,误差性能目标值取为0.0005,学习速率仍为0.01,经过测试神经元数量取为30个。
利用LM算法改进BP神经网络模型对卷烟销量数据进行训练和仿真,经过32次即可达到目标误差0.0005的要求,模型拟合效果明显更好,如图4所示。其平均绝对误差为3.99,平均绝对误差率为1.03%,其拟合精度明显高于ARIMA模型和基于梯度下降算法的BP神经网络模型。
输入2017年12个月的卷烟销量数据,从而得到2018年1~4月销量的预测值,如表3所示,可以看出经过改进后的BP神经网络模型,预测误差更加均衡,其预测效果更加平稳,适应性更强。
表3基于LM算法的BP神经网络模型对2018年1~4月全国卷烟销量预测表(万箱,%,天)
注:2018年,1月共22个工作日、2月共17个工作日、3月共22个工作日、4月共20个工作日。
上述三种模型均在一定程度上解释了卷烟销量数据之间的联系,并对未来4个月的卷烟销量进行预测,但单项预测方法存在未能充分利用历史数据的缺点,从而导致预测结果不够精准,为此,通过构建非线性组合预测模型,吸收各模型的优点,提高预测精度。
本发明通过基于LM算法改进的BP神经网络构造非线性组合模型,非线性组合模型构建过程:①以ARIMA模型拟合值、基于梯度下降算法的BP神经网络拟合值、基于LM算法的BP神经网络拟合值作为输入向量,即选取相同月份的3种模型拟合值作为输入向量,当月真实销量作为输出向量,进行训练后进行仿真拟合;②输入上述3种模型对2018年1~4月全国卷烟销量的预测值,得到组合预测模型对销量的预测结果。构建构成中对销量数据的训练参数设定为:训练次数取为3000,误差性能目标值取为0.0005,学习速率为0.01,经过测试神经元数量取为30个。通过训练,组合预测模型训练的销量数据经过4次达到目标误差精度0.0005的要求,如图5所示,模型拟合效果明显更好,其平均绝对误差为3.3,平均绝对误差率为0.85%,其拟合精度均高于上述3种模型。
输入上述3种模型对2018年1~4月销量的预测值,利用训练好的模型得到组合预测模型对2018年1~4月销量的预测值,如表4所示,可以看出组合预测模型的预测绝对误差较为均衡,均不足1天的当月日均销量。
表4非线性组合预测模型对2018年1~4月全国卷烟销量预测表(万箱,%,天)
注:2018年,1月共22个工作日、2月共17个工作日、3月共22个工作日、4月共20个工作日。
模型预测效果对比如表5所示表5四种模型对2018年1~4月全国卷烟销量误差值对比(万箱,天)
注:2018年,1月共22个工作日、2月共17个工作日、3月共22个工作日、4月共20个工作日。
对比4种模型的预测值,发现各预测模型呈现不同的特点:ARIMA模型预测近期销量相对较准,但较远月份的预测值误差较大;基于梯度下降算法的BP神经网络模型可以较为准确的预测大部分月份的卷烟销量,但对于波动性较大的月份预测误差较大;基于LM算法的BP神经网络模型作为改进后的神经网络模型,预测结果精度更高,预测误差更为平稳;非线性组合预测模型综合利用上述3种模型的优点,有效提取数据各方面信息,预测绝对误差均小于当月的日均量,其预测结果最为平稳,适应性最强,更加贴合卷烟销售实际情况。
以上利用ARIMA模型、基于梯度下降算法的BP神经网络模型、基于Levenberg-Marquardt算法改进的BP神经网络模型构建的非线性组合模型对2006~2017年144个月的全国卷烟销量进行仿真训练,并对2018年1~4月的卷烟销量进行预测。通过对比,表明非线性组合模型有更好的预测效果,对卷烟调拨计划制定,科学安排生产计划等具有一定的参考意义。此外,该模型的构建思路还可用于卷烟各价类的销售预测。
Claims (3)
1.一种基于非线性组合模型的卷烟月度销量预测方法,其特征在于:包括如下步骤:
A、根据卷烟N个月度销量数据受节日因素较大、季节性较强的规律特征,选择多种适应性较强的时间序列预测方法;
B、通过输入卷烟N个月度销售数据端分别通过多种适应性较强的时间序列预测方法对卷烟N个月度销售数据进行建模和训练,分别得到对应的模型,记为一阶模型;
C:选取一阶模型中其中一个模型作为二阶模型,以一阶模型的输出数据作为二阶模型的输入数据,然后对一阶模型输入N个月度销量数据,一阶模型得到的对应的多个拟合值同时作为输入向量,以N月度相对应的真实销量作为输出向量,对二阶模型进行训练,从而得到一阶模型和二阶模型的组合模型;
D:对组合模型输入新的月度销售数据,即可输出得到最新的月度销售预测值。
2.根据权利要求1所述的基于非线性组合模型的卷烟月度销量预测方法,其特征在于:所述的一阶模型分别包括有ARIMA模型、基于梯度下降算法的BP神经网络模型和基于Levenberg-Marquardt的BP神经网络模型,二阶模型为基于Levenberg-Marquardt的BP神经网络模型。
3.根据权利要求2所述的基于非线性组合模型的卷烟月度销量预测方法,其特征在于:所述的步骤B分别包括有如下步骤:
B1,选择ARIMA模型对卷烟月度销售数据进行差分处理、模型识别与定阶、仿真等过程,进而预测未来月份的卷烟销量;
B2,选择基于梯度下降算法的BP神经网络模型,对卷烟月度销售数据进行处理,确定输入层、输出层,创建基于梯度下降算法的网络,编制神经元测试程序,确定神经元数量;设定网络参数,经过训练,完成模型构建;选择预测输入层数据,从而预测未来月份卷烟销量;
B3,选择基于Levenberg-Marquardt的BP神经网络模型,对卷烟月度销售数据进行处理,确定输入层、输出层,创建基于LM算法的神经网络;编制神经元测试程序,确定神经元数量;设定网络参数,经过训练,完成模型构建;选择预测输入层数据,从而预测未来月份卷烟销量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910085257.4A CN109961315A (zh) | 2019-01-29 | 2019-01-29 | 一种基于非线性组合模型的卷烟月度销量预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910085257.4A CN109961315A (zh) | 2019-01-29 | 2019-01-29 | 一种基于非线性组合模型的卷烟月度销量预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109961315A true CN109961315A (zh) | 2019-07-02 |
Family
ID=67023519
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910085257.4A Pending CN109961315A (zh) | 2019-01-29 | 2019-01-29 | 一种基于非线性组合模型的卷烟月度销量预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109961315A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111047101A (zh) * | 2019-12-18 | 2020-04-21 | 陕西中烟工业有限责任公司 | 基于烟支轴向密度预测燃烧锥落头的方法 |
CN112001740A (zh) * | 2020-06-19 | 2020-11-27 | 南京理工大学 | 一种基于自适应神经网络的组合预测方法 |
CN112633576A (zh) * | 2020-12-22 | 2021-04-09 | 华中科技大学 | 应用于卷烟工厂生产排程的两阶段排程优化方法与系统 |
WO2023168916A1 (zh) * | 2022-03-08 | 2023-09-14 | 太原理工大学 | 一种基于不锈钢极薄带退火工艺神经网络模型优化方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103903071A (zh) * | 2014-04-17 | 2014-07-02 | 上海电机学院 | 一种风电功率预测组合方法和系统 |
CN105956699A (zh) * | 2016-04-29 | 2016-09-21 | 连云港天马网络发展有限公司 | 基于电商销售数据的商品分类提取及在该分类提取下的销量预测方法 |
CN105976199A (zh) * | 2016-04-26 | 2016-09-28 | 重庆大学 | 一种基于混合模型的医药销量的预测方法及系统 |
CN106779859A (zh) * | 2016-12-27 | 2017-05-31 | 南京安讯科技有限责任公司 | 一种移动终端产品实时销量预测方法 |
CN106845686A (zh) * | 2016-12-26 | 2017-06-13 | 南通同心烟智道文化传播有限公司 | 基于因素评分的烟草销量预测方法 |
CN107705157A (zh) * | 2017-10-19 | 2018-02-16 | 大连理工大学 | 基于统一动态集成模型和元启发式算法的汽车备件销量预测方法及系统 |
CN108416636A (zh) * | 2018-03-30 | 2018-08-17 | 大连理工大学 | 一种预测药品销量的arima和支持向量回归融合方法 |
CN108764460A (zh) * | 2018-05-16 | 2018-11-06 | 华中科技大学 | 一种基于时间卷积和lstm的时间序列预测方法 |
CN109242265A (zh) * | 2018-08-15 | 2019-01-18 | 杭州电子科技大学 | 基于误差平方和最小的城市需水量组合预测方法 |
-
2019
- 2019-01-29 CN CN201910085257.4A patent/CN109961315A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103903071A (zh) * | 2014-04-17 | 2014-07-02 | 上海电机学院 | 一种风电功率预测组合方法和系统 |
CN105976199A (zh) * | 2016-04-26 | 2016-09-28 | 重庆大学 | 一种基于混合模型的医药销量的预测方法及系统 |
CN105956699A (zh) * | 2016-04-29 | 2016-09-21 | 连云港天马网络发展有限公司 | 基于电商销售数据的商品分类提取及在该分类提取下的销量预测方法 |
CN106845686A (zh) * | 2016-12-26 | 2017-06-13 | 南通同心烟智道文化传播有限公司 | 基于因素评分的烟草销量预测方法 |
CN106779859A (zh) * | 2016-12-27 | 2017-05-31 | 南京安讯科技有限责任公司 | 一种移动终端产品实时销量预测方法 |
CN107705157A (zh) * | 2017-10-19 | 2018-02-16 | 大连理工大学 | 基于统一动态集成模型和元启发式算法的汽车备件销量预测方法及系统 |
CN108416636A (zh) * | 2018-03-30 | 2018-08-17 | 大连理工大学 | 一种预测药品销量的arima和支持向量回归融合方法 |
CN108764460A (zh) * | 2018-05-16 | 2018-11-06 | 华中科技大学 | 一种基于时间卷积和lstm的时间序列预测方法 |
CN109242265A (zh) * | 2018-08-15 | 2019-01-18 | 杭州电子科技大学 | 基于误差平方和最小的城市需水量组合预测方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111047101A (zh) * | 2019-12-18 | 2020-04-21 | 陕西中烟工业有限责任公司 | 基于烟支轴向密度预测燃烧锥落头的方法 |
CN112001740A (zh) * | 2020-06-19 | 2020-11-27 | 南京理工大学 | 一种基于自适应神经网络的组合预测方法 |
CN112633576A (zh) * | 2020-12-22 | 2021-04-09 | 华中科技大学 | 应用于卷烟工厂生产排程的两阶段排程优化方法与系统 |
CN112633576B (zh) * | 2020-12-22 | 2022-04-29 | 华中科技大学 | 应用于卷烟工厂生产排程的两阶段排程优化方法与系统 |
WO2023168916A1 (zh) * | 2022-03-08 | 2023-09-14 | 太原理工大学 | 一种基于不锈钢极薄带退火工艺神经网络模型优化方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109961315A (zh) | 一种基于非线性组合模型的卷烟月度销量预测方法 | |
CN103745273B (zh) | 一种半导体制造过程的多性能预测方法 | |
El Hamzaoui et al. | Application of artificial neural networks to predict the selling price in the real estate valuation process | |
CN112990500B (zh) | 基于改进加权灰色关联分析的台区线损分析方法及系统 | |
CN110598929B (zh) | 一种风电功率非参数概率区间超短期预测方法 | |
CN107729684B (zh) | 基于故障机理综合损伤累积规则的多阶段任务系统可靠性分层建模方法 | |
Yan et al. | Investment estimation of prefabricated concrete buildings based on XGBoost machine learning algorithm | |
CN107886160A (zh) | 一种bp神经网络区间需水预测方法 | |
CN101004593A (zh) | 基于AR(p)模型的过程工业生产过程动态成本控制方法 | |
CN105550747A (zh) | 一种新型卷积神经网络的样本训练方法 | |
Mithiya et al. | Time series analysis and forecasting of oilseeds production in India: using autoregressive integrated moving average and group method of data handling–neural network | |
CN103903070B (zh) | 应用系统资源需求测算系统 | |
Akhbari | Project time and cost forecasting using Monte Carlo simulation and artificial neural networks | |
Ming-Tao et al. | The research on stock price forecast model based on data mining of BP neural networks | |
Mohammadi et al. | Identifying and prioritizing criteria for selecting sustainable façade materials of high-rise buildings | |
CN113537540B (zh) | 一种基于自动特征工程的供热用气量预测模型 | |
Zvyagin | System analysis in optimization and decision making | |
CN102955978B (zh) | 一种基于周期性灰色系统的时装生产控制方法 | |
Ren et al. | Research of LSTM-RNN Model and Its Application Evaluation on Agricultural Products Circulation | |
Abrishami et al. | Survey of a rule based expert system for gas price forecasting | |
CN109325861A (zh) | 采用经验回放机制的目标股票选择方法、装置及存储介质 | |
Abdullah | Forecasting Iraqi oil production using artificial neural networks | |
CN103116272A (zh) | 一种氢氰酸生产过程的在线自适应建模方法 | |
Chen | Optimization Control of Construction Engineering Management Projects Based on Deep Learning Algorithm | |
Dimitrova | Application of the Analytic Network Process (ANP) in a framework of ERP systems implementation success |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190702 |