A general approach is developed for discriminating strong and hereditary symmetric operators. The recursion operator of the Blaszak-Marciniak (BM) equation hierarchy is proved to be strong and hereditary symmetric. As an example of discrete soliton equations related to 3×3 matrix spectral problems, the τ-symmetries and Lie algebra structure of the BM equation are built firstly.