From f37094b20a574706a1c75620e3364d08aaf8d7e9 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Fr=C3=A9d=C3=A9ric=20Karam=C3=A9?= Date: Wed, 7 Oct 2015 15:25:16 +0200 Subject: [PATCH] Minor modifications in the inputs of the function. --- src/gaussian_densities.m | 3 --- src/probability2.m | 3 +-- 2 files changed, 1 insertion(+), 5 deletions(-) diff --git a/src/gaussian_densities.m b/src/gaussian_densities.m index 4f2b21438..73a731403 100644 --- a/src/gaussian_densities.m +++ b/src/gaussian_densities.m @@ -43,9 +43,6 @@ prior = probability2(st_t_1,sqr_Pss_t_t_1,particles) ; % likelihood yt_t_1_i = measurement_equations(particles,ReducedForm,ThreadsOptions) ; eta_t_i = bsxfun(@minus,obs,yt_t_1_i)' ; -%yt_t_1 = sum(yt_t_1_i*weigths1,2) ; -%tmp = bsxfun(@minus,yt_t_1_i,yt_t_1) ; -%Pyy = bsxfun(@times,weigths2',tmp)*tmp' + H ; Pyy = H ; sqr_det = sqrt(det(Pyy)) ; foo = (eta_t_i/Pyy).*eta_t_i ; diff --git a/src/probability2.m b/src/probability2.m index 57ab65f8f..a7486d1b8 100644 --- a/src/probability2.m +++ b/src/probability2.m @@ -34,5 +34,4 @@ function [density] = probability2(mu,S,X) dim = size(X,1) ; normfact = bsxfun(@power,(2*pi),(dim/2)) ; foo = S\(bsxfun(@minus,X,mu)) ; -density = exp(-0.5*sum(foo.*foo)')./abs((normfact*prod(diag(S)))) + 1e-99 ; - +density = exp(-0.5*sum(foo.*foo)')./abs((normfact*prod(diag(S)))) + 1e-99 ; \ No newline at end of file