CN111160651A - 一种基于stl-lstm的地铁客流预测方法 - Google Patents

一种基于stl-lstm的地铁客流预测方法 Download PDF

Info

Publication number
CN111160651A
CN111160651A CN201911403973.9A CN201911403973A CN111160651A CN 111160651 A CN111160651 A CN 111160651A CN 201911403973 A CN201911403973 A CN 201911403973A CN 111160651 A CN111160651 A CN 111160651A
Authority
CN
China
Prior art keywords
time
data
lstm
prediction
subway
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.)
Granted
Application number
CN201911403973.9A
Other languages
English (en)
Other versions
CN111160651B (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.)
Fuzhou University
Original Assignee
Fuzhou University
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 Fuzhou University filed Critical Fuzhou University
Priority to CN201911403973.9A priority Critical patent/CN111160651B/zh
Publication of CN111160651A publication Critical patent/CN111160651A/zh
Application granted granted Critical
Publication of CN111160651B publication Critical patent/CN111160651B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • 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/045—Combinations of 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/04—Architecture, e.g. interconnection topology
    • G06N3/049—Temporal neural networks, e.g. delay elements, oscillating neurons or pulsed inputs
    • 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
    • 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/40—Business processes related to the transportation industry
    • 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
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Software Systems (AREA)
  • Economics (AREA)
  • Evolutionary Computation (AREA)
  • Computational Linguistics (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Biophysics (AREA)
  • Mathematical Physics (AREA)
  • Biomedical Technology (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Strategic Management (AREA)
  • Data Mining & Analysis (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种基于STL‑LSTM的地铁客流预测方法，首先，计算统计周期内的地铁客运量数据；使用STL加法模型对地铁客运量数据进行分解，得到3个分解序列；接着，设定LSTM参数，建立LSTM预测模型。对3个分解序列使用相同参数的LSTM模型进行训练，并进行测试集预测。然后，把3个分解序列的测试集预测结果相加得到地铁数据测试集的预测结果，计算平均相对误差；修改分解周期及时间步长，重复上述过程，平均相对误差基本稳定时停止修改；最后，用测试集平均相对误差最小的模型对得到的3个序列的预测步长后的客运量进行预测，将结果相加，得到最终预测结果。本发明可以提升对地铁客流预测的精度。

Description

一种基于STL-LSTM的地铁客流预测方法
技术领域
本发明涉及地铁客流预测领域，特别是一种基于STL-LSTM的地铁客流预测方法。
背景技术
地铁客流是地铁运营管理的重要依据，对地铁未来客流进行准确预测对地铁运营部门调整行车计划有着重要的意义。考虑到地铁客流有明显的周期性、趋势性和随机性，在建立预测模型时应予以体现。
对时间序列进行分解能剔除周期因素的影响，反映时间序列真正的客观规律和趋势。近些年已经有很多研究人员对其展开深入研究。时间序列分解模型根据采取的策略不同，可以分为X11系列、X12系列、SABL、小波分析和EMD等。这些模型可以较好的对时间序列进行分解，对提升预测精度也有一定帮助，但上述模型并不能充分发挥地铁的客流特性。STL分解模型将时间序列分解为周期序列、趋势序列和随机序列，与地铁客流的周期性、趋势性和随机性相对应。
目前已有一些学者对STL分解模型进行研究，也取得了一定的成果，但并没有将模型应用于地铁客流预测领域。STL将客流分解为周期序列、趋势序列和随机序列，但是地铁客流数据的趋势性不是很明显，地铁客流数据的周期性也是以周为主，STL对地铁客流数据的有效性还有待验证。
发明内容
有鉴于此，本发明的目的是提供一种基于STL-LSTM的地铁客流预测方法，可以提升对地铁客流预测的精度。
本发明采用以下方案实现：一种基于STL-LSTM的地铁客流预测方法，包括以下步骤：
步骤S1：从AFC地铁刷卡数据中获取统计周期内的地铁客运量，得到原始时间序列；
步骤S2：确定分解周期，使用STL加法模型对S1的地铁客运量数据进行分解，得到3个分解序列；
步骤S3：设定神经元个数、学习率、误差值计算公式、训练次数、激活函数、抽样数量、输入维度、时间步长、预测步长以建立LSTM预测模型，设置三个分解序列的训练集和测试集的比例均为x:y，所述x，y皆为大于1且小于输入地铁客运量数据个数的整数；分别对S2得到的3个分解序列使用相同参数的LSTM模型输入训练集进行训练，并分别对它们的测试集进行预测；
步骤S4：把3个分解序列的测试集的预测结果相加得到步骤S1地铁数据测试集的预测结果，计算平均相对误差；
步骤S5：参数调整，修改步骤S2中的分解周期及步骤S3中的时间步长，重复步骤S2到步骤S4的过程，当平均相对误差前后相差0.01％时停止修改，用以寻找精度最高的参数组合；
步骤S6：用测试集平均相对误差最小的模型对步骤S2得到的3个序列的预测步长后的客运量进行预测，将结果相加得到最终预测结果。
进一步地，所述步骤S2具体包括以下步骤：
步骤S21：确定分解周期：
分解周期T为大于1的正整数，T∈(2,3,...,N)，N＜l，l为时间序列长度；
步骤S22：使用STL加法模型分解出趋势项Tt：
yt＝Seasonalt+Tt+Randomt (1)
其中，t为时刻，yt为t时刻的地铁客运量，Seasonalt、Tt、Randomt分别是t时刻的周期成分、趋势成分、随机成分；
Figure BDA0002347937810000031
Figure BDA0002347937810000032
其中，Tt为趋势项，f为时间频率，l为时间序列长度，Tt是长度为l的趋势序列，当时刻t超出定义域时，值为NA，f为奇数时用公式(2)计算，f为偶数时用公式(3)计算；
步骤S23：计算季节项：
将原始序列减去趋势项：
St＝yt-Tt (4)
将各个周期内相同频率下的值平均化，得到季节项figure：
Figure BDA0002347937810000033
将figure中心化，得到中心化的长度为f的季节项figure；
将长度为f的季节项扩展成长度为l的季节项Seasonal：
Seasonalt＝figuret％％f,其中％％为取余运算；(6)步骤S24：计算随机项Randomt：
Randomt＝yt-Tt-Seasonalt (7)
其中，Xt表示t时刻的地铁客运量，Randomt表示t时刻的随机成分，Tt表示t时刻的趋势成分，Seasonalt表示t时刻的季节成分。
通过步骤S22至步骤S24将地铁客运量y，分解成趋势序列T，周期序列Seasonal和随机序列Random。
进一步地，所述步骤S3具体包括以下步骤：
步骤S31：分别对3个序列的数据进行归一化，将数据变成0到1之间的小数；
Figure BDA0002347937810000041
步骤S32：设定LSTM隐藏神经元个数，取值范围为大于1且小于训练集数据个数的整数；设定学习率，取值为0到1；设置抽样数量，所述抽样数量大小必须大于或等于1且小于或等于训练数据集中的样本数；训练次数为1和无穷大之间的整数值；设置误差计算公式和激活函数；设置输入维度即特征个数，时间步长即输入的长度和预测步长，输入时间步长个数据得到预测步长个的数据；
步骤S33：设置训练集与测试集的比例为x:y，每个序列使用相同参数的LSTM进行时间序列曲线拟合，用训练集拟合的曲线预测测试集的地铁客运量，即输入测试集时间步长个数据，得到预测步长个地铁客运量，用滑动窗口的形式即每次往后移动一个数据预测测试集的所有数据；
步骤S34：修改时间步长，重复步骤S32和步骤S33。
进一步地，所述步骤S4涉及的测试集的平均相对误差计算公式如下：
Figure BDA0002347937810000051
其中，表示Qi'是第i个预测客运量，Qi表示第i个实际客运量，n为测试集数据个数。
进一步地，步骤S32中所述误差值计算包括选取均方误差MSE来计算：
Figure BDA0002347937810000052
或选取平均绝对误差MAE来计算：
Figure BDA0002347937810000053
进一步地，步骤S32中所述激活函数计算包括选择Logistic函数：
Figure BDA0002347937810000054
或Tanh函数：
Figure BDA0002347937810000061
与现有技术相比，本发明具有以下有益效果：
本发明一方面是用STL分解模型将原始客流数据分解成周期序列、趋势序列和随机序列，提升了预测精度；另一方面是将STL分解与LSTM预测模型相结合，将地铁客运量数据分解后的数据，分别作为LSTM的输入，将3个序列的结果相加得到最后预测结果，拓宽了STL分解的应用领域。
附图说明
图1为本发明实施例的流程图。
图2为本发明实施例的分解周期为2-4的分解结果图。
图3为本发明实施例的分解周期与测试集平均误差变化情况图。
图4为本发明实施例的测试集误差对比图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
如图1所示，本实施例一种基于STL-LSTM的地铁客流预测方法，包括以下步骤：
步骤S1：通过数据库或编写代码从AFC地铁刷卡数据中获取统计周期内的地铁客运量(统计周期可以是天、小时等)，得到原始时间序列；
步骤S2：确定分解周期，使用STL加法模型对S1的地铁客运量数据进行分解，得到3个分解序列；
步骤S3：设定神经元个数、学习率、误差值计算公式、训练次数、激活函数、抽样数量、输入维度、时间步长、预测步长以建立LSTM预测模型，设置三个分解序列的训练集和测试集的比例均x:y，所述x，y皆为大于1且小于输入地铁客运量数据个数的整数，通常x取7，y取3；分别对S2得到的3个分解序列使用相同参数的LSTM模型输入训练集进行训练，并分别对它们的测试集进行预测；
步骤S4：把3个分解序列的测试集的预测结果相加得到S1地铁数据测试集的预测结果，计算平均相对误差；
步骤S5：参数调整，修改步骤S2中的分解周期及步骤S3中的时间步长，重复步骤S2到步骤S4的过程，当平均相对误差前后相差0.01％时停止修改，用以寻找精度最高的参数组合；
步骤S6：用测试集平均相对误差最小的模型对步骤S2得到的3个序列的预测步长后的客运量进行预测，将结果相加得到最终预测结果。
在本实施例中，所述步骤S2具体包括以下步骤：
步骤S21：确定分解周期：
分解周期T为大于1的正整数，T∈(2,3,...,N)，N＜l，l为时间序列长度；
步骤S22：使用STL加法模型分解出趋势项Tt：
yt＝Seasonalt+Tt+Randomt (1)
其中，t为时刻，yt为t时刻的地铁客运量，Seasonalt、Tt、Randomt分别是t时刻的周期成分、趋势成分、随机成分；
Figure BDA0002347937810000081
Figure BDA0002347937810000082
其中，Tt为趋势项，f为时间频率，l为时间序列长度，Tt是长度为l的趋势序列，当时刻t超出定义域时，值为NA，f为奇数时用公式(2)计算，f为偶数时用公式(3)计算；
步骤S23：计算季节项：
将原始序列减去趋势项：
St＝yt-Tt (4)
将各个周期内相同频率下的值平均化，得到季节项figure：
Figure BDA0002347937810000083
将figure中心化，得到中心化的长度为f的季节项figure；
将长度为f的季节项扩展成长度为l的季节项Seasonal：
Seasonalt＝figuret％％f,其中％％为取余运算；(6)步骤S24：计算随机项Randomt：
Randomt＝yt-Tt-Seasonalt
(7)
其中，Xt表示t时刻的地铁客运量，Randomt表示t时刻的随机成分，Tt表示t时刻的趋势成分，Seasonalt表示t时刻的季节成分。
通过步骤S22至步骤S24将地铁客运量y，分解成趋势序列T，周期序列Seasonal和随机序列Random。
本实施例中，所述步骤S3具体包括以下步骤：
步骤S31：分别对3个序列的数据进行归一化，将数据变成0到1之间的小数；
Figure BDA0002347937810000091
步骤S32：设定LSTM隐藏神经元个数units，取值范围为大于1且小于训练集数据个数的整数；即用多少个神经元进行函数拟合，一个神经元可对应一个特征，一般序列输入长度为多少就设置多少个；
设定学习率learningrate，取值为0到1；进行一轮训练后得到了一个拟合结果，该拟合结果与真实值有一定的差距，这个差距通过误差值计算公式来计算，如果要减小这个差距，需要进行参数调整，学习率是参数调整的幅度，根据loss的变化情况来调整learningrate的变化情况，loss减小时，可以适当调大learningrate，loss增大时，减小learningrate；
设定误差值loss计算，可以选均方误差MSE来计算：
Figure BDA0002347937810000101
或平均绝对误差MAE来计算：
Figure BDA0002347937810000102
等；可用来计算误差来评价此次训练结果的好坏；
设置抽样数量batch_size，输入样本过多时，不能每次都用所有数据进行训练，可以用抽样的方法选择一定数量的数据进行训练，batch_size大小必须大于或等于1且小于或等于训练数据集中的样本数；
设置训练次数epochs，epochs为1和无穷大之间的整数值，可根据误差变化来调整；
设置激活函数，可以选择Logistic函数：
Figure BDA0002347937810000103
Tanh函数：
Figure BDA0002347937810000104
等；不同的激活函数会有不同的输出，激活函数可以把输入转成一定形式的输出，方便模型的训练；
设置输入维度input_dim(特征个数)，时间步长timesteps(输入的长度)和预测步长predictsteps，输入时间步长个数据(每个数据有输入维度个数据)得到预测步长个的数据(输出为预测步长个数据，每个数据为1维)；
步骤S33：设置训练集与测试集的比例为x:y(一般为7:3)，每个序列使用相同参数(包括时间步长、预测步长、学习率、神经元个数等)的LSTM进行时间序列曲线拟合，用训练集拟合的曲线预测测试集的地铁客运量，即输入测试集时间步长个数据，得到该序列预测步长个地铁客运量，用滑动窗口的形式即每次往后移动一个数据预测测试集的所有数据；
步骤S34：修改时间步长，重复步骤S32和步骤S33。
在本实施例中，设置LSTM的神经元个数为100，初始学习率0.01，误差计算公式为均方误差计算公式，训练次数50次，激活函数为tanh，抽样数量50，预测步长为1，输入维度为1，训练集和测试集的比例为7:3。
本实施例中，所述步骤S4涉及的测试集的平均相对误差计算公式如下：
Figure BDA0002347937810000111
其中，表示Qi'是第i个预测客运量，Qi表示第i个实际客运量，n为测试集数据个数。
图2是用STL加法模型分解法对周期分别为2-4天的分解结果示例。
图3是分解周期为2-22，时间步长1-30(由于当时间步长大于30以后，预测的误差在10％以上，所以不再进行研究)时在测试集上的平均相对误差，可以发现周期为2天或者3天时，经过STL加法模型分解后再用LSTM模型预测时，其在测试集上的误差相对较低，分解周期大于22天时，变化已经稳定，所以不再对分解周期大于10的序列进行研究。
图4是不同分解周期2-10在测试集上的误差，可以发现，STL-LSTM预测结果在分解周期在2-10时，预测精度高于直接对原序列进行LSTM预测，且精度高于最近2年提出类似模型EMD-LSTM。以上所述仅为本发明的较佳实施例，凡依本发明申请专利范围所做的均等变化与修饰，皆应属本发明的涵盖范围。

Claims (6)

1.一种基于STL-LSTM的地铁客流预测方法，其特征在于：包括以下步骤：
步骤S1：从AFC地铁刷卡数据中获取统计周期内的地铁客运量，得到原始时间序列；所述统计周期包括天或小时；
步骤S2：确定分解周期，使用STL加法模型对步骤S1的地铁客运量数据进行分解，得到3个分解序列；
步骤S3：设定神经元个数、学习率、误差值计算公式、训练次数、激活函数、抽样数量、输入维度、时间步长、预测步长以建立LSTM预测模型，设置三个分解序列的训练集和测试集的比例均为x:y，所述x，y皆为大于1且小于输入地铁客运量数据个数的整数；分别对步骤S2得到的3个分解序列使用相同参数的LSTM模型输入训练集进行训练，并分别对它们的测试集进行预测；
步骤S4：把3个分解序列的测试集的预测结果相加得到步骤S1地铁数据测试集的预测结果，计算平均相对误差；
步骤S5：参数调整，修改步骤S2中的分解周期及步骤S3中的时间步长，重复步骤S2到步骤S4的过程，当平均相对误差前后相差0.01％时停止修改，用以寻找精度最高的参数组合；
步骤S6：用测试集平均相对误差最小的模型对步骤S2得到的3个序列的预测步长后的客运量进行预测，将结果相加得到最终预测结果。
2.根据权利要求1所述的一种基于STL-LSTM的地铁客流预测方法，其特征在于：所述步骤S2具体包括以下步骤：
步骤S21：确定分解周期：
分解周期T为大于1的正整数，T∈(2,3,...,N)，N＜l，l为时间序列长度；
步骤S22：使用STL加法模型分解出趋势项Tt：
yt＝Seasonalt+Tt+Randomt (1)
其中，t为时刻，yt为t时刻的地铁客运量，Seasonalt、Tt、Randomt
分别是t时刻的周期成分、趋势成分、随机成分；
Figure FDA0002347937800000021
Figure FDA0002347937800000022
其中，Tt为趋势项，f为时间频率，l为时间序列长度，Tt是长度为l的趋势序列，当时刻t超出定义域时，值为NA，f为奇数时用公式(2)计算，f为偶数时用公式(3)计算；
步骤S23：计算季节项：
将原始序列减去趋势项：
St＝yt-Tt (4)
将各个周期内相同频率下的值平均化，得到季节项figure：
Figure FDA0002347937800000023
将figure中心化，得到中心化的长度为f的季节项figure；
将长度为f的季节项扩展成长度为l的季节项Seasonal：
Seasonalt＝figuret％％f,其中％％为取余运算； (6)
步骤S24：计算随机项Randomt：
Randomt＝yt-Tt-Seasonalt (7)
其中，Xt表示t时刻的地铁客运量，Randomt表示t时刻的随机成分，Tt表示t时刻的趋势成分，Seasonalt表示t时刻的季节成分。
通过步骤S22至步骤S24将地铁客运量y，分解成趋势序列T，周期序列Seasonal和随机序列Random。
3.根据权利要求1所述的一种基于STL-LSTM的地铁客流预测方法，其特征在于：所述步骤S3具体包括以下步骤：
步骤S31：分别对3个序列的数据进行归一化，将数据变成0到1之间的小数；
Figure FDA0002347937800000031
步骤S32：设定LSTM隐藏神经元个数，取值范围为大于1且小于训练集数据个数的整数；设定学习率，取值为0到1；设置抽样数量，所述抽样数量大小必须大于或等于1且小于或等于训练数据集中的样本数；训练次数为1和无穷大之间的整数值；设置误差计算公式和激活函数；设置输入维度即特征个数，时间步长即输入的长度和预测步长，输入时间步长个数据得到预测步长个的数据；
步骤S33：设置训练集与测试集比例为x:y，每个序列使用相同参数的LSTM进行时间序列曲线拟合，用训练集拟合的曲线预测测试集的地铁客运量，即输入测试集时间步长个数据，得到预测步长个地铁客运量，用滑动窗口的形式即每次往后移动一个数据预测测试集的所有数据；
步骤S34：修改时间步长，重复步骤S32和步骤S33。
4.根据权利要求1所述的一种基于STL-LSTM的地铁客流预测方法，其特征在于：步骤S4中所述测试集的平均相对误差MRE计算公式如下：
Figure FDA0002347937800000041
其中，表示Qi'是第i个预测客运量，Qi表示第i个实际客运量，n为测试集数据个数。
5.根据权利要求3所述的一种基于STL-LSTM的地铁客流预测方法，其特征在于：步骤S32中所述误差值计算包括选取均方误差MSE来计算：
Figure FDA0002347937800000042
或选取平均绝对误差MAE来计算：
Figure FDA0002347937800000043
6.根据权利要求3所述的一种基于STL-LSTM的地铁客流预测方法，其特征在于：步骤S32中所述激活函数计算包括选择Logistic函数：
Figure FDA0002347937800000051
或Tanh函数：
Figure FDA0002347937800000052
CN201911403973.9A 2019-12-31 2019-12-31 一种基于stl-lstm的地铁客流预测方法 Expired - Fee Related CN111160651B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911403973.9A CN111160651B (zh) 2019-12-31 2019-12-31 一种基于stl-lstm的地铁客流预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911403973.9A CN111160651B (zh) 2019-12-31 2019-12-31 一种基于stl-lstm的地铁客流预测方法

Publications (2)

Publication Number Publication Date
CN111160651A true CN111160651A (zh) 2020-05-15
CN111160651B CN111160651B (zh) 2022-07-08

Family

ID=70559656

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911403973.9A Expired - Fee Related CN111160651B (zh) 2019-12-31 2019-12-31 一种基于stl-lstm的地铁客流预测方法

Country Status (1)

Country Link
CN (1) CN111160651B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112215422A (zh) * 2020-10-13 2021-01-12 北京工业大学 基于季节性分解的长短时记忆网络水质动态预警方法
CN112288158A (zh) * 2020-10-28 2021-01-29 税友软件集团股份有限公司 一种业务数据预测方法及相关装置
CN113139643A (zh) * 2021-03-09 2021-07-20 卓望数码技术（深圳）有限公司 网卡流量模型构建方法、流量预测方法、设备及存储介质
CN113705910A (zh) * 2021-08-31 2021-11-26 深信服科技股份有限公司 一种数据样本的扩充方法、装置、设备和介质
CN113792931A (zh) * 2021-09-18 2021-12-14 北京京东振世信息技术有限公司 数据预测方法、装置、物流货量预测方法、介质和设备
CN114037138A (zh) * 2021-11-04 2022-02-11 南京地铁建设有限责任公司 基于双层分解和深度学习的地铁短时进站客流预测系统及实施方法
CN115130786A (zh) * 2022-08-11 2022-09-30 华侨大学 一种景区游客饱和量的判断方法和装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110160927A1 (en) * 2009-12-30 2011-06-30 Wilson Kevin W Method for Prediction for Nonlinear Seasonal Time Series
CN106250999A (zh) * 2015-06-03 2016-12-21 阿里巴巴集团控股有限公司 预测流失率的方法、装置和系统
CN107067283A (zh) * 2017-04-21 2017-08-18 重庆邮电大学 基于历史商家记录及用户行为的电商消费客流量预测方法
CN107274001A (zh) * 2017-04-28 2017-10-20 国网浙江省电力公司 一种基于数据挖掘的用电客户流失电量预测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110160927A1 (en) * 2009-12-30 2011-06-30 Wilson Kevin W Method for Prediction for Nonlinear Seasonal Time Series
CN106250999A (zh) * 2015-06-03 2016-12-21 阿里巴巴集团控股有限公司 预测流失率的方法、装置和系统
CN107067283A (zh) * 2017-04-21 2017-08-18 重庆邮电大学 基于历史商家记录及用户行为的电商消费客流量预测方法
CN107274001A (zh) * 2017-04-28 2017-10-20 国网浙江省电力公司 一种基于数据挖掘的用电客户流失电量预测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李应兵 等: "基于VMD-BP神经网络模型的铁路车站月度客流发送量预测研究", 《交通运输研究》 *
武创: "基于BP神经网络的重庆市轨道交通客流预测分析", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112215422A (zh) * 2020-10-13 2021-01-12 北京工业大学 基于季节性分解的长短时记忆网络水质动态预警方法
CN112288158A (zh) * 2020-10-28 2021-01-29 税友软件集团股份有限公司 一种业务数据预测方法及相关装置
CN113139643A (zh) * 2021-03-09 2021-07-20 卓望数码技术（深圳）有限公司 网卡流量模型构建方法、流量预测方法、设备及存储介质
CN113705910A (zh) * 2021-08-31 2021-11-26 深信服科技股份有限公司 一种数据样本的扩充方法、装置、设备和介质
CN113792931A (zh) * 2021-09-18 2021-12-14 北京京东振世信息技术有限公司 数据预测方法、装置、物流货量预测方法、介质和设备
CN114037138A (zh) * 2021-11-04 2022-02-11 南京地铁建设有限责任公司 基于双层分解和深度学习的地铁短时进站客流预测系统及实施方法
CN115130786A (zh) * 2022-08-11 2022-09-30 华侨大学 一种景区游客饱和量的判断方法和装置
CN115130786B (zh) * 2022-08-11 2024-08-06 华侨大学 一种景区游客饱和量的判断方法和装置

Also Published As

Publication number Publication date
CN111160651B (zh) 2022-07-08

Similar Documents

Publication Publication Date Title
CN111160651B (zh) 一种基于stl-lstm的地铁客流预测方法
CN115329682B (zh) 一种基于多周期特征的lstm-svr地铁车站温度预测方法
CN111563706A (zh) 一种基于lstm网络的多变量物流货运量预测方法
CN107102969A (zh) 一种时间序列数据的预测方法和系统
CN111967688A (zh) 一种基于卡尔曼滤波器与卷积神经网络的电力负荷预测方法
CN105825040A (zh) 短期电力负荷预测方法
CN101739819A (zh) 预测交通流的方法和装置
CN110633859B (zh) 一种两阶段分解集成的水文序列预测方法
CN110909928B (zh) 能源负荷短期预测方法、装置、计算机设备及存储介质
CN112434891A (zh) 基于wcnn-alstm的太阳辐照度时间序列的预测方法
CN115128427B (zh) Mos器件寿命预测方法、装置、电子设备、介质及程序产品
CN112736902B (zh) 一种基于stl分解的时间序列短期电力负荷预测方法
CN114912666A (zh) 一种基于ceemdan算法和注意力机制的短时客流量预测方法
CN112434890A (zh) 基于CEEMDAN-BiLSTM的隧道沉降时间序列的预测方法
CN114694379B (zh) 一种基于自适应动态图卷积的交通流预测方法及系统
CN116911419A (zh) 一种基于趋势相关性特征学习的长时序预测方法
CN114091768A (zh) 基于stl和带有注意力机制的lstm的旅游需求预测方法
CN107704944B (zh) 一种基于信息论学习的股市波动区间预测模型的构建方法
CN105354644A (zh) 一种基于集成经验模态分解和1-范数支持向量机分位数回归的金融时间序列预测方法
CN117035155A (zh) 一种水质预测方法
CN116739168A (zh) 一种基于灰色理论与编解码器的径流量预测方法
CN111815458A (zh) 基于细粒度量化标注和集成方法的动态投资组合配置方法
CN116957639A (zh) 基于arima算法自适应参数选择的方法及系统
CN117094431A (zh) 一种多尺度熵门控的DWTformer气象数据时序预测方法及设备
CN108134687B (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220708

CF01 Termination of patent right due to non-payment of annual fee