The training algorithm based on variational approximation for separable 2D-HMM