From dc2ff00938136debec6f24d4bfac4c9d44dd24e1 Mon Sep 17 00:00:00 2001 From: mu Date: Tue, 4 Feb 2014 13:17:00 +0000 Subject: [PATCH] xt --- GPy/models/state_space_xt_sep.py | 12 +++++++----- 1 file changed, 7 insertions(+), 5 deletions(-) diff --git a/GPy/models/state_space_xt_sep.py b/GPy/models/state_space_xt_sep.py index 9938c4e5..f263bd51 100644 --- a/GPy/models/state_space_xt_sep.py +++ b/GPy/models/state_space_xt_sep.py @@ -20,8 +20,8 @@ # The spatial coordinate of data point do not change over time # Kernel structure: separatable kernel # -# Spatial kernel : Matern32 -# Temporal kernel : state space of of rbf +# Spatial kernel : rbf +# Temporal kernel : state space of of Matern32 import numpy as np from scipy import linalg @@ -56,8 +56,10 @@ class StateSpace_1(Model): # Default kernel if kernel is None: - self.kern = kern.Matern32(1,lengthscale=0.1) - self.spacekern = kern.rbf(1,lengthscale=0.1) + self.kern = kern.Matern32(1,lengthscale=0.5) + #self.spacekern = kern.rbf(1,lengthscale=0.1) + self.spacekern = kern.exponential(1,lengthscale=0.4) + #self.spacekern = kern.Matern32(1,lengthscale=0.5) else: self.kern = kernel @@ -167,7 +169,7 @@ class StateSpace_1(Model): count = count+1 (M, P) = self.kalman_filter(F1,L1,Qc1,H1,self.sigma2,Pinf1,NX.T,NY) - stop + #stop # Run the Rauch-Tung-Striebel smoother #if not filter: #(M, P) = self.rts_smoother(F,L,Qc,X.T,M,P)