Skip to content
Snippets Groups Projects
Commit b63e4e7e authored by Frédéric Karamé's avatar Frédéric Karamé
Browse files

Add the possibility of proposal approximated with Monte Carlo.

(cherry picked from commit ee6eaa84)
parent f6c9a0f7
Branches
No related tags found
No related merge requests found
...@@ -71,7 +71,11 @@ if isempty(init_flag2) ...@@ -71,7 +71,11 @@ if isempty(init_flag2)
init_flag2 = 1; init_flag2 = 1;
end end
if ParticleOptions.proposal_approximation.cubature || ParticleOptions.proposal_approximation.montecarlo if ParticleOptions.proposal_approximation.montecarlo
nodes = randn(ParticleOptions.number_of_particles,number_of_state_variables+number_of_structural_innovations) ;
weights = 1/ParticleOptions.number_of_particles ;
weights_c = weights ;
elseif ParticleOptions.proposal_approximation.cubature
[nodes,weights] = spherical_radial_sigma_points(number_of_state_variables+number_of_structural_innovations) ; [nodes,weights] = spherical_radial_sigma_points(number_of_state_variables+number_of_structural_innovations) ;
weights_c = weights ; weights_c = weights ;
elseif ParticleOptions.proposal_approximation.unscented elseif ParticleOptions.proposal_approximation.unscented
...@@ -118,5 +122,5 @@ else ...@@ -118,5 +122,5 @@ else
StateVectorMean = PredictedStateMean + KalmanFilterGain*PredictionError; StateVectorMean = PredictedStateMean + KalmanFilterGain*PredictionError;
StateVectorVariance = PredictedStateVariance - KalmanFilterGain*PredictedObservedVariance*KalmanFilterGain'; StateVectorVariance = PredictedStateVariance - KalmanFilterGain*PredictedObservedVariance*KalmanFilterGain';
StateVectorVariance = .5*(StateVectorVariance+StateVectorVariance'); StateVectorVariance = .5*(StateVectorVariance+StateVectorVariance');
StateVectorVarianceSquareRoot = chol(StateVectorVariance)'; %reduced_rank_cholesky(StateVectorVariance)'; StateVectorVarianceSquareRoot = reduced_rank_cholesky(StateVectorVariance)';
end end
\ No newline at end of file
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment