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 PDFInfo
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/148—Wavelet transforms
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M7/00—Conversion 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/30—Compression; 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)
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)
| 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)
| 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)
| 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 |
-
1999
- 1999-06-02 JP JP11154458A patent/JP2000349646A/ja active Pending
-
2000
- 2000-06-01 CA CA002339127A patent/CA2339127C/en not_active Expired - Fee Related
- 2000-06-01 EP EP00931637A patent/EP1146649A4/en not_active Withdrawn
- 2000-06-01 US US09/744,912 patent/US6772185B1/en not_active Expired - Fee Related
- 2000-06-01 WO PCT/JP2000/003545 patent/WO2000076074A1/ja not_active Ceased
Patent Citations (3)
| 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)
| 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 |