% Run the Kalman filter x_est = zeros(size(x_true)); P_est = zeros(size(t)); for i = 1:length(t) % Prediction step x_pred = A * x_est(:,i-1); P_pred = A * P_est(:,i-1) * A' + Q; % Update step K = P_pred * H' / (H * P_pred * H' + R); x_est(:,i) = x_pred + K * (y(i) - H * x_pred); P_est(:,i) = (eye(2) - K * H) * P_pred; end