In this paper we study behavioral systems whose trajectories are given as solutions of
quaternionic difference equations. As happens in the commutative case, it turns out that
quaternionic polynomial matrices play an important role in this context. Therefore we focus
our attention on such matrices and derive new results concerning their Smith form. Based
on these results, we obtain characterizations of system theoretic properties of quaternionic
behaviors