The proper way to do this would be to use the
kalman_transition_matrix.m
See [State space form of DSGE)