holdOffTime=holdOFF,
fix_init_Q=1,
r=1,
evalMetrics='F',
ignoreUp2=ignoreUp2)
res_new['Alg'] = 'My Implimentation of FRAUST Version 3.4'
pltSummary2(res_new, data, (e_high, e_low))
if v3_3:
'''My previous version'''
res_v3_3 = FRAHST_V3_3(data,
alpha=alpha,
e_low=e_low,
e_high=e_high,
holdOffTime=holdOFF,
fix_init_Q=1,
r=1,
evalMetrics='F',
ignoreUp2=ignoreUp2)
res_v3_3[
'Alg'] = 'My Implimentation of FRAUST V3.3 Best simple version'
pltSummary2(res_v3_3, data, (e_high, e_low))
if v3_1:
'''My older version'''
res_v3_1 = FRAHST_V3_1(data,
alpha=alpha,
e_low=e_low,
e_high=e_high,
'alg': algInfo
},
'data': res_fr_n
})
# Plot
pltSummary2(
res_fr_n, streams, (el, eh))
if run_frahst_previous == 1:
# My version of Frahst
res_fr_p = FRAHST_V3_3(
streams,
alpha=a,
e_low=el,
e_high=eh,
holdOffTime=h,
fix_init_Q=1,
r=1,
evalMetrics='F',
r_upper_bound=upper_bound,
ignoreUp2=ignoreUp2)
res_fr_p[
'Alg'] = 'MyFrahst-OLD: alpha = ' + str(
a
) + ' ,E_Thresh = (' + str(
el) + ',' + str(eh) + ')'
# Store
algInfo = 'FRAHST-OLD: E_Thresh = (' + str(el) + ',' + str(eh) + ')\n' + \
'alpha = ' + str(a) + '\nHoldOff = ' + str(h)
results.append({
# -*- coding: utf-8 -*-
"""
Created on Sat Mar 05 23:19:11 2011
THIS IS THE VERSION TO RUN ON MAC, NOT WINDOWS
@author: musselle
"""
from numpy import eye, zeros, dot, sqrt, log10, trace, arccos, nan, arange, ones
import numpy as np
import scipy as sp
import numpy.random as npr
from numpy.linalg import qr, eig, norm, solve
from matplotlib.pyplot import plot, figure, title, step, ylim
from artSigs import genCosSignals_no_rand , genCosSignals
import scipy.io as sio
from utils import analysis, QRsolveA, pltSummary2
from PedrosFrahst import frahst_pedro_original
from Frahst_v3_1 import FRAHST_V3_1
from Frahst_v3_3 import FRAHST_V3_3
from Frahst_v3_4 import FRAHST_V3_4
from Frahst_v4_0 import FRAHST_V4_0
from load_syn_ping_data import load_n_store
from QR_eig_solve import QRsolve_eigV
from create_Abilene_links_data import create_abilene_links_data
from MAfunctions import MA_over_window
def FRAHST_V6_2(data, r=1, alpha=0.96, EW_mean_alpha = 0.1, EWMA_filter_alpha = 0.3,
residual_thresh = 0.1, epsilon = 0.05,
e_low = 0.96, e_high = 0.98, fix_init_Q = 0, holdOffTime=0,
evalMetrics = 'F', static_r = 0, r_upper_bound = None, L = 5, ignoreUp2 = 0,
data_norm_window = 50):
"""
# -*- coding: utf-8 -*-
"""
Created on Sat Mar 05 23:19:11 2011
THIS IS THE VERSION TO RUN ON MAC, NOT WINDOWS
@author: musselle
"""
from numpy import eye, zeros, dot, sqrt, log10, trace, arccos, nan, arange, ones
import numpy as np
from numpy.random import rand
from numpy.linalg import qr, eig, norm, solve
from matplotlib.pyplot import plot, figure, title, step
from artSigs import genCosSignals_no_rand , genCosSignals
import scipy.io as sio
from utils import analysis, QRsolveA, pltSummary2
from PedrosFrahst import frahst_pedro_original
from Frahst_v3_1 import FRAHST_V3_1
from Frahst_v3_3 import FRAHST_V3_3
from Frahst_v3_4 import FRAHST_V3_4
from load_syn_ping_data import load_n_store
from QR_eig_solve import QRsolve_eigV
def FRAHST_V4_0(data, r=1, alpha=0.96, e_low=0.96, e_high=0.98, fix_init_Q = 0, holdOffTime=0,
evalMetrics = 'F', static_r = 0, r_upper_bound = None, L = 5, ignoreUp2 = 0):
"""
Fast Rank Adaptive Householder Subspace Tracking Algorithm (FRAHST)
Version 4.0 - Now also approximates eigenvalues for the approximated tracked basis for the eignevectors
- Approach uses an orthogonal iteration arround X.T
- Note, only a good approximation if alpha ~< 1. Used as its the fastest method
as X.T b --> b must be solved anyway.
- New entries in res
# Use Active Normal
data, ant, time = prime_data('SYN', 'AN', useCol = range(2,8), use2time = 7000)
# Use Passive Normal
#data, ant, time = prime_data('SYN', 'PN', useCol = range(2,8), use2time = 7000)
"""PARAMETERS"""
e_high = 0.98
e_low = 0.85
alpha = 0.9
holdOFF = 0
'''My latest version'''
res_new = FRAHST_V3_3(data, alpha=alpha, e_low=e_low, e_high=e_high,
holdOffTime=holdOFF, fix_init_Q = 1, r = 1, evalMetrics = 'F')
res_new['Alg'] = 'My Latest Implimentation of FRAUST '
pltSummary2(res_new, data, (e_high, e_low))
'''My older version'''
res_old = FRAHST_V3_1(data, alpha=alpha, e_low=e_low, e_high=e_high,
holdOffTime=holdOFF, fix_init_Q = 1, r = 1, evalMetrics = 'F')
res_old['Alg'] = 'My Older Implimentation of FRAUST '
pltSummary2(res_old, data, (e_high, e_low))
'''Pedros Version'''
res_ped = frahst_pedro_original(data, r=1, alpha=alpha, e_low=e_low, e_high=e_high,
holdOffTime=holdOFF, evalMetrics='F')