compute Kalman gain
D = H@B@H.T + R
K = B @ H @ np.linalg.inv(D)