WO2000076074A1 - Methode predictive pour series chronologiques au moyen de series numeriques d'ondelettes et dispositif y relatif - Google Patents

Methode predictive pour series chronologiques au moyen de series numeriques d'ondelettes et dispositif y relatif Download PDF

Info

Publication number
WO2000076074A1
WO2000076074A1 PCT/JP2000/003545 JP0003545W WO0076074A1 WO 2000076074 A1 WO2000076074 A1 WO 2000076074A1 JP 0003545 W JP0003545 W JP 0003545W WO 0076074 A1 WO0076074 A1 WO 0076074A1
Authority
WO
WIPO (PCT)
Prior art keywords
series
time
prediction
time series
wavelet
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.)
Ceased
Application number
PCT/JP2000/003545
Other languages
English (en)
French (fr)
Japanese (ja)
Inventor
Naoki Masuda
Kazuyuki Aihara
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.)
Japan Science and Technology Agency
Original Assignee
Japan Science and Technology Corp
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 Japan Science and Technology Corp filed Critical Japan Science and Technology Corp
Priority to US09/744,912 priority Critical patent/US6772185B1/en
Priority to EP00931637A priority patent/EP1146649A4/en
Priority to CA002339127A priority patent/CA2339127C/en
Publication of WO2000076074A1 publication Critical patent/WO2000076074A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/148Wavelet transforms
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M7/00Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
    • H03M7/30Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction

Definitions

  • the present invention relates to a time series prediction method using a wavelet series and an apparatus therefor.
  • time-series prediction has been performed only in the time domain or only in the frequency domain.
  • the prediction in the time domain is shown in Fig. 1, and many models such as the linear regression model belong to this.
  • the prediction in the frequency domain is as shown in Fig. 2, and the model for estimating the shape of the power spectrum belongs to this.
  • the prediction accuracy of a time series having a time structure only in a specific frequency band or a complex time series in which different time structures exist in different frequency bands is low. This is because the prediction method using only the time domain or only the frequency domain cannot handle information in the time domain and the frequency domain simultaneously.
  • the time series when the time series is subjected to the Exublet expansion, the time series is expressed as a sum of frequency components localized in a time domain.
  • the present invention performs a wavelet transform on a time series, decomposes the time series into a plurality of time series that are band-limited in a frequency domain, and reverses the time series of each frequency component from individual prediction to a single-letter inverse. By performing the conversion, an original predicted value of the time series can be obtained. It is an object of the present invention to provide a time series prediction method and a device using a series of Eulet series.
  • a time series prediction method using a wavelet series (a) the time series is subjected to a quadrature transform by a transform unit to form a plurality of time series band-limited in the frequency domain. (B) predict each of the decomposed frequency components with a prediction unit, and (c) reconstruct the predicted value of each of the frequency components with an inverse ⁇ ublet transform unit to obtain the original time series. This is to obtain the predicted value of
  • FIG. 1 is a diagram showing a conventional time domain prediction method.
  • FIG. 2 is a diagram illustrating a conventional prediction method in the frequency domain.
  • FIG. 3 is a diagram showing a method of performing prediction through a low-pass filter.
  • FIG. 4 is a configuration diagram of the present time-series prediction system using wavelets, showing an embodiment of the present invention.
  • FIG. 5 is a diagram showing sampling of a time frequency component by a wavelet transform (WT) showing the embodiment of the present invention.
  • WT wavelet transform
  • FIG. 6 is a diagram showing sampling of a time frequency component by steady wavelet transform (S WT) showing the embodiment of the present invention.
  • FIG. 7 is a diagram showing a prediction result with respect to a prediction time showing the embodiment of the present invention (part 1). 1).
  • FIG. 8 is a diagram (part 2) illustrating a prediction result with respect to a prediction time according to the embodiment of the present invention.
  • FIG. 4 is a configuration diagram of the present time-series prediction system using tablets, showing an embodiment of the present invention.
  • 1 is a wavelet transform unit
  • 2 i,..., 2 ⁇ - !, 2 n are prediction units of each frequency component
  • 3 is an inverse wavelet transform unit. Therefore, first, the time series X is decomposed for each frequency component by the wavelet transform unit 1, and each frequency component is predicted by the prediction units 2,,..., 2 relieve-,, 2 n , respectively. The predicted values of each frequency component are inversely wavelet-transformed by the inverse wavelet transform unit 3 to reconstruct the original time-series predicted values.
  • each frequency component represents one having a higher frequency in order from the top.
  • prediction is performed using prediction unit 2 !,..., 2n-l, 2shire
  • Each of the predicted frequency components is reconstructed by the inverse ⁇ ublet transform unit 3, and the predicted value of the original time series is obtained.
  • the wavelet transform / inverse tablet transform unit 1> 3 may be a normal single-bit (WT) or a stationary (stationsry) -single-bit transform (SWT).
  • time series is subjected to a ⁇ -Ublet transform and uniquely divided into the time series of each frequency component.
  • a method for solving the problem will be described.
  • steady-state wavelet transform necessary to observe time-series frequency components at the same time.
  • CR (k) ⁇ f (t), 0 R., (t)> (1) and is uniquely subjected to a single-blet transform (hereinafter abbreviated as WT).
  • WT a single-blet transform
  • f (t) obtained by projecting ⁇ z (11) ⁇ to the subspace spanned by ⁇ ⁇ 1 , 1 ⁇ (t) IkeZ ⁇ is transformed by the above equation (1). That is, f (t)- ⁇ di (k) ⁇ (t) + y, CR (k) 0R .. (t) (1) '
  • ⁇ (t) can be obtained by projection, and the stationary wavelet transform (hereinafter abbreviated as SWT) of the above equation (3) can be used. Note that J3 is displayed in place of bold d.
  • (k) ⁇ f (t), 2 j / 2 0 [2 j (t-2 " J k)]>, (R ⁇ j ⁇ J)
  • SWT sampling in the time frequency domain is as shown in FIG.
  • the original dynamical system is the sum of the deterministic chaotic dynamical system with high-frequency components and the deterministic force-ososdynamic system with low-frequency components. In some cases, each frequency component is dominated by only one dynamical system, so it can be predicted with high accuracy until a certain time. Therefore, the original time series can be predicted with relatively high accuracy.
  • high-frequency noise components decrease the prediction accuracy of low-frequency components, and the error increases exponentially.
  • the time series ( ⁇ ( ⁇ )) is obtained by adding the discrete-time observation value of the X coordinate of the Wrestler equation to y (n) and the time series ⁇ ( ⁇ ) that normalized the x n of the Ikeda map to mean 0.
  • y ( ⁇ ) represents low frequency components and ⁇ ( ⁇ ) represents high frequency noise.
  • is the variance of the time series
  • the low-frequency component y (n) is the discrete-time observation of the X component of the Lorentz equation
  • 1 (n) is the white noise
  • Figure 8 shows the results of estimating ⁇ up to 128 periods ahead in the same manner as in Example ⁇ ⁇ above. For comparison, a method of estimating after removing high-frequency components through a low-pass filter (see Fig. 3) is also listed.
  • a time series is decomposed into frequency components using a single tablet, and a predicted value of the original time series is obtained from the predicted value of each frequency component by inverse inverse transform.
  • SWT is advantageous from the point that the number of sampling points is large and the frequency components are observed at the same time.
  • redundant data is taken too much, data with close times are strongly associated with each other, and the prediction method deteriorates the prediction accuracy because the prediction method tries to forcibly construct a causal relationship between such data. It was shown that it is better to use WT to make predictions while avoiding WT is also advantageous from the viewpoint of short calculation time.
  • the prediction method of the present invention is particularly effective for predicting a time series in which the low frequency range and the high frequency range are governed by different rules. Furthermore, the forecast unit is set according to the nature of the time series. All kinds of time series can be predicted.
  • time series prediction system using the wave series of the present invention can be applied to prediction of a wide variety of time series regardless of whether the time series is deterministic or stochastic.
  • the original time series is obtained by decomposing the time series into frequency components using the wavelet, and obtaining the predicted value of the original time series from the forecasted value of each frequency component by inverse ⁇ -Urblet transform. Can be accurately predicted.
  • the present invention is based on low-frequency hydrodynamic time evolution. It can be applied to prediction fields such as added brain measurement data, low-frequency membrane potential fluctuations, and high-frequency channel noise plus neuronal membrane potential measurement data.
  • time series prediction method and apparatus using the wavelet series of the present invention are particularly effective for prediction of a time series in which the low frequency range and the high frequency range are governed by different rules. Further, the present invention can be applied to fields for predicting all kinds of time series by changing the prediction unit according to the nature of the time series.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)
  • Complex Calculations (AREA)
PCT/JP2000/003545 1999-06-02 2000-06-01 Methode predictive pour series chronologiques au moyen de series numeriques d'ondelettes et dispositif y relatif Ceased WO2000076074A1 (fr)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US09/744,912 US6772185B1 (en) 1999-06-02 2000-06-01 Time-series predicting method using wavelet number series and device thereof
EP00931637A EP1146649A4 (en) 1999-06-02 2000-06-01 METHOD FOR PREDICTING TIME SERIES WITH WAVELET NUMBER SERIES
CA002339127A CA2339127C (en) 1999-06-02 2000-06-01 Time-series prediction method and apparatus utilizing wavelet coefficient series

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP11154458A JP2000349646A (ja) 1999-06-02 1999-06-02 ウェーブレット系数列を用いた時系列予測方法及びその装置
JP11/154458 1999-06-02

Publications (1)

Publication Number Publication Date
WO2000076074A1 true WO2000076074A1 (fr) 2000-12-14

Family

ID=15584684

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2000/003545 Ceased WO2000076074A1 (fr) 1999-06-02 2000-06-01 Methode predictive pour series chronologiques au moyen de series numeriques d'ondelettes et dispositif y relatif

Country Status (5)

Country Link
US (1) US6772185B1 (https=)
EP (1) EP1146649A4 (https=)
JP (1) JP2000349646A (https=)
CA (1) CA2339127C (https=)
WO (1) WO2000076074A1 (https=)

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100004976A1 (en) * 2008-04-08 2010-01-07 Plan4Demand Solutions, Inc. Demand curve analysis method for analyzing demand patterns
US20100070426A1 (en) 2008-09-15 2010-03-18 Palantir Technologies, Inc. Object modeling for exploring large data sets
US8429194B2 (en) 2008-09-15 2013-04-23 Palantir Technologies, Inc. Document-based workflows
JP5018809B2 (ja) * 2009-03-04 2012-09-05 沖電気工業株式会社 時系列データ予測装置
US8732574B2 (en) 2011-08-25 2014-05-20 Palantir Technologies, Inc. System and method for parameterizing documents for automatic workflow generation
US9348677B2 (en) 2012-10-22 2016-05-24 Palantir Technologies Inc. System and method for batch evaluation programs
US8930897B2 (en) 2013-03-15 2015-01-06 Palantir Technologies Inc. Data integration tool
US8903717B2 (en) 2013-03-15 2014-12-02 Palantir Technologies Inc. Method and system for generating a parser and parsing complex data
US8855999B1 (en) 2013-03-15 2014-10-07 Palantir Technologies Inc. Method and system for generating a parser and parsing complex data
US8868486B2 (en) 2013-03-15 2014-10-21 Palantir Technologies Inc. Time-sensitive cube
US8909656B2 (en) 2013-03-15 2014-12-09 Palantir Technologies Inc. Filter chains with associated multipath views for exploring large data sets
US8938686B1 (en) 2013-10-03 2015-01-20 Palantir Technologies Inc. Systems and methods for analyzing performance of an entity
US9105000B1 (en) 2013-12-10 2015-08-11 Palantir Technologies Inc. Aggregating data from a plurality of data sources
US8924429B1 (en) 2014-03-18 2014-12-30 Palantir Technologies Inc. Determining and extracting changed data from a data source
CN109992608B (zh) * 2019-03-26 2022-03-04 浙江大学 一种基于频域的多模型融合预测方法和系统
CN110414442B (zh) * 2019-07-31 2022-04-01 广东省智能机器人研究院 一种压力时序数据分段特征值预测方法
CN110839016B (zh) * 2019-10-18 2022-07-15 平安科技(深圳)有限公司 异常流量监测方法、装置、设备及存储介质
EP3896625A1 (en) * 2020-04-17 2021-10-20 Tata Consultancy Services Limited An adaptive filter based learning model for time series sensor signal classification on edge devices
US11681914B2 (en) 2020-05-08 2023-06-20 International Business Machines Corporation Determining multivariate time series data dependencies
CN114944831A (zh) * 2022-05-12 2022-08-26 中国科学技术大学先进技术研究院 多周期时间序列数据分解方法、装置、设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09307903A (ja) * 1996-05-16 1997-11-28 Ricoh Co Ltd 多階調画像符号化方法
US5883978A (en) * 1995-12-22 1999-03-16 Matsushita Electric Industrial Co., Ltd Image compressing method and image compressing apparatus
JPH1188183A (ja) * 1997-09-11 1999-03-30 Seiko Epson Corp ウェーブレット変換装置およびその方法並びにウェーブレット逆変換装置およびその方法並びに画像符号化装置およびその方法並びに画像復号化装置およびその方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2997355B2 (ja) * 1991-12-04 2000-01-11 シャープ株式会社 復号装置
GB2362772B (en) * 1997-06-10 2002-01-16 Sony Uk Ltd Concealing errors

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5883978A (en) * 1995-12-22 1999-03-16 Matsushita Electric Industrial Co., Ltd Image compressing method and image compressing apparatus
JPH09307903A (ja) * 1996-05-16 1997-11-28 Ricoh Co Ltd 多階調画像符号化方法
JPH1188183A (ja) * 1997-09-11 1999-03-30 Seiko Epson Corp ウェーブレット変換装置およびその方法並びにウェーブレット逆変換装置およびその方法並びに画像符号化装置およびその方法並びに画像復号化装置およびその方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP1146649A4 *

Also Published As

Publication number Publication date
CA2339127C (en) 2005-03-29
CA2339127A1 (en) 2000-12-14
EP1146649A1 (en) 2001-10-17
US6772185B1 (en) 2004-08-03
JP2000349646A (ja) 2000-12-15
EP1146649A4 (en) 2010-06-02

Similar Documents

Publication Publication Date Title
WO2000076074A1 (fr) Methode predictive pour series chronologiques au moyen de series numeriques d'ondelettes et dispositif y relatif
So et al. Comparison of various periodograms for sinusoid detection and frequency estimation
CN102891770A (zh) 一种利用短相关模型预测长相关序列的方法
CN105320995A (zh) 一种基于频域分解的风力发电短期负荷预测方法及装置
Conroy et al. The limits of extended Kalman filtering for pulse train deinterleaving
Flornes et al. A direct interpolation method for irregular sampling
CN111753681A (zh) 基于Legendre小波多分辨率分析的分数阶系统辨识方法
Amburgey et al. Time-series models for event counts
EP3080640A1 (en) Phase retrieval algorithm for generation of constant time envelope with prescribed fourier transform magnitude signal
Rupniewski Period and signal reconstruction from the curve of trains of samples
Thomson et al. Estimation of second-order properties from jittered time series
CN107729868A (zh) 一种基于小波分析的信号处理方法
Boche et al. On the calculation of the Hilbert transform from interpolated data
Cevri et al. Performance analysis of Gibbs sampling for Bayesian extracting sinusoids
Saha et al. Empirical mode decomposition and stationary wavelet transformation in internet traffic prediction
US7489665B1 (en) Enhanced time resolution using multiple receptions
CN112034253A (zh) 一种moa在线监测方法
Jayawardena Stochastic vs. deterministic approaches of modelling hydrological time series
CN101603525B (zh) 一种压缩机工况在线监控方法及装置
Cutler Computing pointwise fractal dimension by conditioning in multivariate distributions and time series
Xia et al. Wavelet-based signal processing of large SHM data
CN120602310B (zh) 一种用于电能表通信故障的诊断方法及系统
Wan et al. Inverse Differential Equation Modeling of ENSO Prediction Based on Memory Kernel Functions
Rui WPANFIS: combine fuzzy neural network with multiresolution for network traffic prediction
Masry et al. Field estimation from randomly located binary noisy sensors

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): CA US

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE

ENP Entry into the national phase

Ref document number: 2339127

Country of ref document: CA

Ref country code: CA

Ref document number: 2339127

Kind code of ref document: A

Format of ref document f/p: F

WWE Wipo information: entry into national phase

Ref document number: 2000931637

Country of ref document: EP

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 09744912

Country of ref document: US

WWP Wipo information: published in national office

Ref document number: 2000931637

Country of ref document: EP