Update EnKF.m

Revert a change

Author: JoonasJL

plugins/Kalman/m/EnKF.m

@@ -1,33 +1,27 @@
-function [z_inverse] = EnKF(m, A, P, Q, L, R, timeSteps, number_of_frames, n_ensembles, q_given)
+function [z_inverse] = EnKF(m, A, P, Q, L, R, timeSteps, number_of_frames, n_ensembles)
 %ENKF Summary of this function goes here
 %x_ensemble = mvnrnd(zeros(size(m)), 100* ones(size(m,1)), n_ensembles)';
 x_ensemble = mvnrnd(zeros(size(m)), P, n_ensembles)';
 z_inverse = cell(0);
-if not(q_given)
-    q_values = Q;
-end
-h = zef_waitbar(0, 'EnKF Filtering');
+h = zef_waitbar(0,1, 'EnKF Filtering');
 for f_ind = 1:number_of_frames
-    zef_waitbar(f_ind/number_of_frames,h,...
+    zef_waitbar(f_ind,number_of_frames,h,...
         ['EnKF Filtering ' int2str(f_ind) ' of ' int2str(number_of_frames) '.']);
     f = timeSteps{f_ind};
-    if not(q_given)
-        Q = diag(q_values(:,f_ind));
-    end
     w = mvnrnd(zeros(size(m)), Q, n_ensembles)';
     % Forecasts
 
     x_f = A * x_ensemble + w;
     C = cov(x_f');
     correlationLocalization = true;
     if correlationLocalization
-    T = corrcoef(x_f');
-    % explain How to find 0.05
-    T(abs(T) < 0.05) = 0;
-    C = C .* T;
+        T = corrcoef(x_f');
+        % explain How to find 0.05
+        T(abs(T) < 0.05) = 0;
+        C = C .* T;
     end
     v = mvnrnd(zeros(size(R,1),1), R, n_ensembles);
-    
+
     % method to calculate resolution D
     method = '3';
     if(method == '1')