我有一个kalman滤波器的实现,我正在做矩阵运算。在某个点上,我应该减去两个1x1矩阵。我有个错误,我不知道它是从哪里来的。
在python中进行矩阵运算的最佳方法是什么?在import numpy as np
import pylab as pl
import scipy as Sci
import scipy.linalg as linalg
class GetPos(object):
def __init__(self):
self.Posp = 0
self.Velp = 80
self.z = np.matrix(0)
def __repr__(self):
return "from GetPos.__repr__ z=%s" % (self.z)
def __call__(self):
self.dt = 0.1
self.w = 0 + 10*np.random.random()
self.v = 0 + 10*np.random.random()
self.z = self.Posp + self.Velp*self.dt + self.v
self.Posp = self.z - self.v
self.Velp = 80 + self.w
print 'from GetPos.__call__ z = %s' % self.z
return self.z
class DvKalman(object):
def __init__(self):
self.dt