The equivalence of systems plays an important role in 2D systems, which is closely related to the equivalence of bivariate polynomial matrices. In this paper, the Smith forms of several classes of bivariate polynomial matrices are investigated. Some new results and criteria on the equivalence of these matrices are obtained. These criteria are easily checked by computing the reduced Gr\"{o}bner basis of the ideal generated by the minors of lower order of the given matrix.