def method_ERA(self, _measSignal, _simSignal):
'''
_measSignal as output
_simSignal as input
'''
'''TODO: match sampletime from meas with sim '''
# if (_measSignal.get_sampleTime()!= _simSignal.get_sampleTime()):
timeSignal= _measSignal.get_sampleTime()
if _measSignal!= None and _simSignal!= None:
outSignal= _measSignal.get_signalMag()
inSignal= _simSignal.get_signalReal()
else:
if _measSignal!= None:
outSignal= _measSignal.get_signalMag()
inSignal= _measSignal.get_signalMag()
if _simSignal!= None:
outSignal= _simSignal.get_signalReal()
inSignal= _simSignal.get_signalReal()
num_states = 2
# a,b,c = mr.compute_ERA_model([timeSignal,outSignal,inSignal], num_states)
a,b,c = mr.compute_ERA_model(numpy.array(outSignal[0:1000]), num_states)
print 'Measurements: '
print 'A= ', a
print 'B= ', b
print 'C= ', c
a,b,c = mr.compute_ERA_model(numpy.array(inSignal), num_states)
print 'Simulation: '
print 'A= ', a
print 'B= ', b
print 'C= ', c
def calculate_eigenvalues(self, oneSignal):
''' eigenvalues and eigen vectors '''
print len(oneSignal)
#first era method
self.__A, self.__B, self.__C = mr.compute_ERA_model(np.array(oneSignal), 3)
# second, eigenvalues and eigenvectors
self.__elambda, self.__vlambda = linalg.eig(self.__A)
'''
Created on Jul 30, 2015
@author: ekj05
'''
import numpy as N
import modred as MR
'''----num_vecs------------------------------column number'''
num_vecs = 4
'''test case---------------------generating arbitrary data'''
vecs = N.random.random((10, num_vecs))
print vecs
'''what is eigenvalues?'''
a,b,c =MR.compute_ERA_model(vecs,4)
print 'printing a matrix'
print a
print 'printing b matrix'
print b
print 'printing c matrix'
print c
# """ opening the h5 file """
# #File1=h5py.File('Simu.h5','r')
# ''' using solution from load_pandaSource '''
# # h5Data= h5py.File(_sourceH5,'r')
# senyal = self.ioh5.get_senyal('block0')
# ''' format of the signal (sampletime, real/magnintude, imag/polar) '''
# """ getting the data set from the h5 file """
# #d1=File1[u'subgroup']
# # d1= h5Data[u'df']
