CN112163022A - 一种面向时序数据的预测方法 - Google Patents
一种面向时序数据的预测方法 Download PDFInfo
- Publication number
- CN112163022A CN112163022A CN202011095091.3A CN202011095091A CN112163022A CN 112163022 A CN112163022 A CN 112163022A CN 202011095091 A CN202011095091 A CN 202011095091A CN 112163022 A CN112163022 A CN 112163022A
- Authority
- CN
- China
- Prior art keywords
- value
- time point
- time
- matrix
- reaching
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/24—Querying
- G06F16/245—Query processing
- G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
- G06F16/2477—Temporal data queries
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Accounting & Taxation (AREA)
- Finance (AREA)
- Development Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Probability & Statistics with Applications (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Marketing (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Economics (AREA)
- Evolutionary Computation (AREA)
- General Business, Economics & Management (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Computational Linguistics (AREA)
- Databases & Information Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种面向时序数据的预测方法,属于数据分析技术领域。本发明的方法包括以下步骤:采集历史时序数据;根据历史时序数据设定若干个时间点,计算每个时间点相对前一个时间点的变化率;记录达到峰值或谷值的时间点,选取达到峰值或谷值之后的N个时间点数据进行计算得到累积变化率矩阵;再根据累积变化率矩阵计算得到新的矩阵;之后利用聚类算法对新的矩阵进行聚类得到聚类结果,并根据聚类结果进行预测。本发明克服了现有技术中,时序数据达到峰值或者谷值后的变化趋势预测不准确的不足,本发明可以精确判断时序数据在未来达到某一峰值或谷值后后续的变化趋势,进而提高了时序数据达到峰值或者谷值后后续变化趋势的预测精度。
Description
技术领域
本发明属于数据分析技术领域,更具体地说,涉及一种面向时序数据的预测方法。
背景技术
时序数据是基于时间的一系列数据在有时间的坐标中将这些数据点连成线,往过去看可以做成多维度报表,揭示其规律性、趋势性、异常性;往未来看可以做大数据分析,机器学习,实现预测和预警。时间序列数据是同一统一指标按时间顺序记录的数据列。在同一数据列中的各个数据必须是同口径的,要求具有可比性。时序数据可以是时期数,也可以是时点数。时间序列分析的目的是通过找出样本内时间序列的统计特性和发展规律性,构建时间序列模型,进行样本外预测。
目前使用较广泛的是运用统计原理对历史数据进行统计分析,找出其内在规律,建立映射关系模型用于预测。常用的统计方法有多元线性回归算法、人工神经网络算法和灰色理论算法等。使用差异较大的一组样本训练模型,得到的参数进行校验时往往产生较大的误差,即使重新训练,也得不到很好的效果。例如发明创造名称为:一种自适应时间序列数据预测方法(申请日:2020年5月25日;申请号:202010451280.3),该方案公开了一种自适应时间序列数据预测方法,包括以下步骤:收集并清理原始时间序列数据,增强数据维度;构建数据集,评定数据集特性;构建数据集的预测模型;对预测模型进行训练和参数优化;利用预测模型计算数据预测值。该方案解决了多时间序列特征条件下数据预测的问题,提供一种自适应时间序列数据预测方法,可以监测时间序列数据的特性,并根据特性选择对应的预测模型实现预测。但是,该方案的不足之处在于:当条件变化剧烈时,模型预测精度还有待提高,甚至有可能会失效。
综上所述,如何提高时序数据达到峰值或者谷值后后续变化趋势的预测精度,是现有技术亟需解决的问题。
1.要解决的问题
本发明克服了现有技术中,时序数据达到峰值或者谷值后的变化趋势预测不准确的不足,提出了一种面向时序数据的预测方法,可以分析时序数据达到峰值或谷值后的短期动量规律,从而可以精确判断时序数据在未来达到某一峰值或谷值后后续的变化趋势,进而提高了时序数据达到峰值或者谷值后后续变化趋势的预测精度。
2.技术方案
为了解决上述问题,本发明所采用的技术方案如下:
本发明的一种面向时序数据的预测方法,包括以下步骤:
采集历史时序数据;根据历史时序数据设定若干个时间点,计算每个时间点相对前一个时间点的变化率;然后记录达到峰值或谷值的时间点,选取达到峰值或谷值之后的N个时间点数据进行计算得到累积变化率矩阵;再根据累积变化率矩阵计算得到新的矩阵;之后利用聚类算法对新的矩阵进行聚类得到聚类结果,并根据聚类结果进行预测。
更进一步地,利用下列公式计算每个时间点相对前一个时间点的变化率:
Vt是当前时间点的值;Vt-1是当前时间点的前一个时间点的值;P表示变化率。
更进一步地,选取达到峰值或谷值之后的N个时间点数据的具体过程为:当某个时间点为峰值或谷值时,若在该时间点之前的5个时间点之内没有时间点值大于峰值或者小于谷值,则以该时间点为零点,并记录该时间点之后的N个连续时间点数据。
更进一步地,计算得到累积变化率矩阵的具体过程为:
设定若干个时间点有M个零点,根据M个零点得到(N+1)*M的矩阵X;
Y表示累积变化率矩阵;X[i,:]表示矩阵X的每一行数据,i为正整数;X[0,:]表示矩阵X第一行的时间点数据。
更进一步地,根据累积变化率矩阵计算得到新的矩阵的具体过程为:
先根据累积变化率矩阵计算每一行元素的正增长率,平均值,中值,最大值,最小值和T_Value;
根据每一行元素的正增长率,平均值,中值,最大值,最小值和T_Value组成新的矩阵。
更进一步地,利用下述公式计算正增长率,平均值,中值,最大值,最小值和T_Value:
Positive%表示在达到峰值或谷值后的第t个时间点的正增长率;P0表示第t个时间点Y值大于0的个数;P表示第t个时间点Y值的总共个数;
Median=Median(Yt)–1
Median表示在达到峰值或谷值后的第t个时间点的中值;
Max=Max(Yt)–1
Max表示在达到峰值或谷值后的第t个时间点的最大值;
Min=Min(Yt)–1
Min表示在达到峰值或谷值后的第t个时间点的最小值;
T_Value表示正增长率,平均值,中值,最大值,最小值的指标计算结果的可信度;Std表示每一行元素的方差。
更进一步地,利用k-means算法对新的矩阵进行聚类。
更进一步地,对新的矩阵进行聚类的具体过程为:
(a)根据新的矩阵得到k个聚类并对每个聚类确定一个初始聚类中心;
(b)将样本集中的样本xi按照最小距离原则分配到最邻近聚类uk;
dij=min(||xi-uk||),xi∈S,uk∈Z
其中,S为样本集,Z为整数集;
(c)使用每个聚类中的样本均值作为新的聚类中心;
重复步骤(b),(c)直到聚类中心不再变化。
3.有益效果
相比于现有技术,本发明的有益效果为:
本发明的一种面向时序数据的预测方法,通过记录峰值或谷值后的时间点数据,从而可以分析时序数据达到峰值或谷值后的短期动量规律,进一步通过聚类算法对新的矩阵进行聚类分析,进而可以精确判断时序数据在未来达到某一峰值或谷值后后续的变化趋势,进一步提高了时序数据达到峰值或者谷值后后续变化趋势的预测精度。
附图说明
图1为本发明的方法流程示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例;而且,各个实施例之间不是相对独立的,根据需要可以相互组合,从而达到更优的效果。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为进一步了解本发明的内容,结合附图和实施例对本发明作详细描述。
实施例1
结合图1所示,本发明的一种面向时序数据的预测方法,包括以下步骤:
(1)采集数据
采集历史时序数据;本实施例中历史时序数据为多年的商品销售数据。
(2)计算变化率
根据历史时序数据设定若干个时间点,值得说明的是,时间点可以是天数或者小时数,本实施例设定近6年的时间点,即时间点个数为365*6个。进一步地,计算每个时间点相对前一个时间点的变化率,具体地,利用下列公式计算每个时间点相对前一个时间点的变化率:
P表示变化率,Vt是当前时间点的值;Vt-1是当前时间点的前一个时间点的值。
(3)计算累积变化率矩阵
记录达到峰值或谷值的时间点,值得说明的是,峰值和谷值指的是单位时间段内的极值点,峰值为极大值点,谷值为极小值点,其中,单位时间段可任意设置。之后选取达到峰值或谷值之后的N个时间点数据进行计算得到累积变化率矩阵;具体地,选取达到峰值或谷值之后的N个时间点数据的具体过程如下:
当某个时间点为峰值或谷值时,若在该时间点之前的5个时间点之内没有时间点值大于峰值或者小于谷值,则以该时间点为零点,并记录该时间点之后的N个连续时间点数据,N为整数,且N不小于20。值得说明的是,通过记录峰值或谷值后的时间点数据,从而可以对时序数据达到峰值或谷值后的短期动量规律进行分析。
计算得到累积变化率矩阵的具体过程为:设定若干个时间点有M个零点,根据M个零点得到(N+1)*M的矩阵X;
Y表示累积变化率矩阵;X[i,:]表示矩阵X的每一行数据,i为正整数;X[0,:]表示矩阵X第一行的时间点数据。
(4)计算得到新的矩阵
根据累积变化率矩阵计算得到新的矩阵,具体过程如下:
先根据累积变化率矩阵计算每一行元素的正增长率,平均值,中值,最大值,最小值和T_Value;具体计算公式如下:
Positive%表示在达到峰值或谷值后的第t个时间点的正增长率;P0表示第t个时间点Y值大于0的个数;P表示第t个时间点Y值的总共个数;
Median=Median(Yt)–1
Median表示在达到峰值或谷值后的第t个时间点的中值;
Max=Max(Yt)–1
Max表示在达到峰值或谷值后的第t个时间点的最大值;
Min=Min(Yt)–1
Min表示在达到峰值或谷值后的第t个时间点的最小值;
T_Value表示正增长率,平均值,中值,最大值,最小值的指标计算结果的可信度;Std表示每一行元素的方差。
之后根据每一行元素的正增长率,平均值,中值,最大值,最小值和T_Value组成新的矩阵。
(5)对新的矩阵进行聚类
利用聚类算法对新的矩阵进行聚类得到聚类结果,本发明利用k-means算法对新的矩阵进行聚类;具体地,聚类过程如下;
(a)根据新的矩阵得到k个聚类并对每个聚类确定一个初始聚类中心;
(b)将样本集中的样本xi按照最小距离原则分配到最邻近聚类uk;
dij=min(||xi-uk||),xi∈S,uk∈Z
其中,S为样本集,Z为整数集;
(c)使用每个聚类中的样本均值作为新的聚类中心;
重复步骤(b),(c)直到聚类中心不再变化,共得到k个聚类,其中,k为整数,本实施例中k=4。值得说明的是,通过聚类可以将新的矩阵分为具有明显特征的类,从而便于判断时序数据在未来达到某一峰值或谷值后后续的变化趋势,进一步提高了时序数据达到峰值或者谷值后后续变化趋势的预测精度。
进一步地,计算聚类的误差平方和得到误差最小的聚类,具体公式如下:
其中,J表示误差平方和。值得说明的是,结合最小二乘法和拉格朗日原理,聚类中心为对应类别中个数据点的平均值,同时为了使算法收敛,在迭代的过程中,应使得最终的聚类中心尽可能的不变。
(6)预测结果
根据聚类结果进行预测,矩阵每行计算所得正增长率和平均值为后续每一天的正增长率和平均值。根据聚类分析结果,可得到四种特征,分别为(1)达到峰值时,正增长率超过平均值;(2)达到峰值时,正增长率低于平均值;(3)达到谷值后,正增长率超过平均值;(4)达到谷值后,正增长率低于平均值。本实施例以商家销售某个商品为例,当达到峰值或谷值后,观察其后的20天销售量。达到峰值时,当正增长率超过平均值时,判断商品每天销售量将继续增加,应扩大生产量,防止造成供不应求的情况;当正增长率低于平均值时,考虑销售量短期内将逐步下降,应适量减少生产,避免供过于求;达到谷值后,当正增长率超过平均值时,判断商品每天销售量将逐步增加,应适量扩大生产量;当正增长率低于平均值时,考虑销售量短期内将继续下降,应减少生产量或停止生产,并及时找出原因。
本发明的一种面向时序数据的预测方法,通过记录峰值或谷值后的时间点数据,从而可以分析时序数据达到峰值或谷值后的短期动量规律,进一步通过聚类算法对新的矩阵进行聚类分析,进而可以精确判断时序数据在未来达到某一峰值或谷值后后续的变化趋势,进一步提高了时序数据达到峰值或者谷值后后续变化趋势的预测精度。
在上文中结合具体的示例性实施例详细描述了本发明。但是,应当理解,可在不脱离由所附权利要求限定的本发明的范围的情况下进行各种修改和变型。详细的描述和附图应仅被认为是说明性的,而不是限制性的,如果存在任何这样的修改和变型,那么它们都将落入在此描述的本发明的范围内。此外,背景技术旨在为了说明本技术的研发现状和意义,并不旨在限制本发明或本申请和本发明的应用领域。
Claims (8)
1.一种面向时序数据的预测方法,其特征在于,包括以下步骤:
采集历史时序数据;
根据历史时序数据设定若干个时间点,计算每个时间点相对前一个时间点的变化率;
记录达到峰值或谷值的时间点,选取达到峰值或谷值之后的N个时间点数据进行计算得到累积变化率矩阵;再根据累积变化率矩阵计算得到新的矩阵;
利用聚类算法对新的矩阵进行聚类得到聚类结果,并根据聚类结果进行预测。
3.根据权利要求1所述的一种面向时序数据的预测方法,其特征在于,选取达到峰值或谷值之后的N个时间点数据的具体过程为:当某个时间点为峰值或谷值时,若在该时间点之前的5个时间点之内没有时间点值大于峰值或者小于谷值,则以该时间点为零点,并记录该时间点之后的N个连续时间点数据。
5.根据权利要求4所述的一种面向时序数据的预测方法,其特征在于,根据累积变化率矩阵计算得到新的矩阵的具体过程为:
先根据累积变化率矩阵计算每一行元素的正增长率,平均值,中值,最大值,最小值和T_Value;
根据每一行元素的正增长率,平均值,中值,最大值,最小值和T_Value组成新的矩阵。
6.根据权利要求5所述的一种面向时序数据的预测方法,其特征在于,利用下述公式计算正增长率,平均值,中值,最大值,最小值和T_Value:
Positive%表示在达到峰值或谷值后的第t个时间点的正增长率;P0表示第t个时间点Y值大于0的个数;P表示第t个时间点Y值的总共个数;
Median=Median(Yt)–1
Median表示在达到峰值或谷值后的第t个时间点的中值;
Max=Max(Yt)–1
Max表示在达到峰值或谷值后的第t个时间点的最大值;
Min=Min(Yt)–1
Min表示在达到峰值或谷值后的第t个时间点的最小值;
T_Value表示正增长率,平均值,中值,最大值,最小值的指标计算结果的可信度;Std表示每一行元素的方差。
7.根据权利要求1~6任一项所述的一种面向时序数据的预测方法,其特征在于,利用k-means算法对新的矩阵进行聚类。
8.根据权利要求7所述的一种面向时序数据的预测方法,其特征在于,对新的矩阵进行聚类的具体过程为:
(a)根据新的矩阵得到k个聚类并对每个聚类确定一个初始聚类中心;
(b)将样本集中的样本xi按照最小距离原则分配到最邻近聚类uk;
dij=min(||xi-uk||),xi∈S,uk∈Z
其中,S为样本集,Z为整数集;
(c)使用每个聚类中的样本均值作为新的聚类中心;
重复步骤(b),(c)直到聚类中心不再变化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011095091.3A CN112163022A (zh) | 2020-10-14 | 2020-10-14 | 一种面向时序数据的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011095091.3A CN112163022A (zh) | 2020-10-14 | 2020-10-14 | 一种面向时序数据的预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112163022A true CN112163022A (zh) | 2021-01-01 |
Family
ID=73866822
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011095091.3A Pending CN112163022A (zh) | 2020-10-14 | 2020-10-14 | 一种面向时序数据的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112163022A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113670970A (zh) * | 2021-07-22 | 2021-11-19 | 四川省建材工业科学研究院有限公司 | 石膏胶凝材料凝结硬化时间检测方法 |
-
2020
- 2020-10-14 CN CN202011095091.3A patent/CN112163022A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113670970A (zh) * | 2021-07-22 | 2021-11-19 | 四川省建材工业科学研究院有限公司 | 石膏胶凝材料凝结硬化时间检测方法 |
CN113670970B (zh) * | 2021-07-22 | 2023-12-15 | 四川省建材工业科学研究院有限公司 | 石膏胶凝材料凝结硬化时间检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2020063689A1 (zh) | 一种电力系统热负荷预测的方法和装置 | |
CN107992976B (zh) | 热点话题早期发展趋势预测系统及预测方法 | |
CN108762228A (zh) | 一种基于分布式pca的多工况故障监测方法 | |
CN111709454B (zh) | 一种基于最优copula模型的多风电场出力聚类评估方法 | |
CN111008726B (zh) | 一种电力负荷预测中类图片转换方法 | |
CN110297469B (zh) | 基于重采样的集成特征选择算法的生产线故障判断方法 | |
CN113515512A (zh) | 一种工业互联网平台数据的质量治理及提升方法 | |
CN116380445B (zh) | 基于振动波形的设备状态诊断方法及相关装置 | |
CN109670549B (zh) | 火电机组的数据筛选方法、装置以及计算机设备 | |
CN116128141A (zh) | 风暴潮预测方法、装置、存储介质及电子设备 | |
CN115096627A (zh) | 一种液压成形智能装备制造过程故障诊断与运维方法及系统 | |
CN112163022A (zh) | 一种面向时序数据的预测方法 | |
CN114997841A (zh) | 一种基于区块链的低碳行为数据管理系统 | |
CN108134687B (zh) | 一种基于马尔可夫链的灰色模型局域网峰值流量预测方法 | |
CN111832730A (zh) | 一种油液不确定状态的信度表征与状态识别的方法 | |
CN116894687A (zh) | 一种基于机器学习的用电量分析方法及系统、电子设备 | |
CN115794805A (zh) | 一种中低压配网量测数据补齐方法 | |
CN115983113A (zh) | 一种全连接神经网络模型下的丝印工艺优化方法 | |
CN115618708A (zh) | 基于增量式Informer算法的设备健康状态预测方法 | |
CN114662056A (zh) | 镀膜厚度控制方法、装置和存储介质 | |
CN115049136A (zh) | 一种变压器负荷预测方法 | |
JP7458183B2 (ja) | エネルギー需要変動パターンの分類方法及びシステム | |
CN110175705B (zh) | 一种负荷预测方法及包含该方法的存储器、系统 | |
CN113947904A (zh) | 基于s-g滤波和深信度网络的多尺度短时交通流预测方法 | |
CN113361750A (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 |