# -*- coding: utf-8 -*-
"""
Created on Wed May 18 11:23:36 2011
Test SPIRIT and Frahst on Contol Time Series: Cyclic 
@author: - Musselle
"""
from ControlCharts import Tseries
from CUSUM import cusum
import numpy as np
import matplotlib.pyplot as plt
from Frahst_v3_1 import FRAHST_V3_1
from SPIRIT import SPIRIT
from utils import analysis, QRsolveA, pltSummary
#===============================================================================
# Runscript
#===============================================================================
#============
# Initialise
#============
'Create time series string'
series_1 = Tseries(0)
series_2 = Tseries(0)
series_3 = Tseries(0)
series_4 = Tseries(0)
series_5 = Tseries(0)
# S 1
# CYCLIC - Normal - Cyclic
series_1.cyclicEt(100, base=6, noise=0.1, amp=2, period=5, phase=0, \
                    noise_type = 'gauss')
series_1.normalEt(100, base=6, noise=0.00000000000000001, noise_type='gauss')
# ============
# Initialise
# ============

"Create time series string"

series_1 = Tseries(0)
series_2 = Tseries(0)
series_3 = Tseries(0)
series_4 = Tseries(0)
series_5 = Tseries(0)

# S 1
# CYCLIC - Normal - Cyclic
series_1.cyclicEt(100, base=6, noise=0.1, amp=2, period=5, phase=0, noise_type="gauss")
series_1.normalEt(100, base=6, noise=0.00000000000000001, noise_type="gauss")
series_1.cyclicEt(100, base=6, noise=0.1, amp=2, period=5, phase=0, noise_type="gauss")
# S 2
# Cyclic
series_2.cyclicEt(300, base=6, noise=0.1, amp=2, period=5, phase=1.5, noise_type="gauss")
# S 3
# Cyclic
series_3.cyclicEt(300, base=6, noise=0.1, amp=2, period=5, phase=0.5, noise_type="gauss")

# S 4
# Cyclic
series_4.cyclicEt(300, base=1, noise=0.1, amp=2, period=5, phase=0.2, noise_type="gauss")

# S 5
# Cyclic
series_5.cyclicEt(300, base=5, noise=0.1, amp=2, period=5, phase=0.3, noise_type="gauss")
Exemple #3
0
import matplotlib.pyplot as plt

#===============================================================================
# Runscript
#===============================================================================

#============
# Initialise
#============

'Create time series string'

series_1 = Tseries(0)

# NORMAL - normalEt(self, size, base=0, noise, noise_type = 'gauss')
series_1.normalEt(50, 0, 1, 'gauss')
series_1.normalEt(50, 3, 1, 'gauss')

# CYCLIC - cyclicEt(self, size, base=0, noise=2, amp=10, period=25, \
#                 phase = 0, noise_type = 'gauss'):
series_1.cyclicEt(100, 3, 1, 5, 50)

# NORMAL -
series_1.normalEt(50, 3, 1, 'gauss')
series_1.normalEt(50, -3, 1, 'gauss')

# UP - upEt(self,size,base=0,noise=1,gradient=0.2, noise_type = 'gauss')
series_1.upEt(100, -3, 1, 0.2, 'gauss')

# Normal
series_1.normalEt(50, 17, 1, 'gauss')
# -*- coding: utf-8 -*-
"""
Created on Wed May 18 11:23:36 2011
Test SPIRIT and Frahst on Contol Time Series: Cyclic 
@author: - Musselle
"""
from ControlCharts import Tseries
from CUSUM import cusum 
import numpy as np
import matplotlib.pyplot as plt
from Frahst_v3_1 import FRAHST_V3_1
from SPIRIT import SPIRIT
from utils import analysis, QRsolveA, pltSummary
#===============================================================================
# Runscript 
#===============================================================================
#============
# Initialise
#============
'Create time series string'
series_1 = Tseries(0)
series_2 = Tseries(0)
series_3 = Tseries(0)
series_4 = Tseries(0)
series_5 = Tseries(0)
# S 1
# CYCLIC - Normal - Cyclic
series_1.cyclicEt(100, base=6, noise=0.1, amp=2, period=5, phase=0, \
                    noise_type = 'gauss')
series_1.normalEt(100, base=6, noise = 0.00000000000000001, noise_type ='gauss')
Exemple #5
0
import matplotlib.pyplot as plt

#===============================================================================
# Runscript 
#===============================================================================

#============
# Initialise
#============

'Create time series string'

series_1 = Tseries(0)

# NORMAL - normalEt(self, size, base=0, noise, noise_type = 'gauss')
series_1.normalEt(50, 0, 1,'gauss')
series_1.normalEt(50, 3, 1,'gauss')

# CYCLIC - cyclicEt(self, size, base=0, noise=2, amp=10, period=25, \
#                 noise_type = 'gauss'):
series_1.cyclicEt(100, 3, 1, 5, 50, 'gauss')

# NORMAL - 
series_1.normalEt(50, 3, 1,'gauss')
series_1.normalEt(50, -3, 1,'gauss')

# UP - upEt(self,size,base=0,noise=1,gradient=0.2, noise_type = 'gauss')
series_1.upEt(100, -3, 1, 0.2, 'gauss')

# Normal 
series_1.normalEt(50, 17, 1,'gauss')