If you exclude non-real operations and states you are removing part of the system such that it becomes impossible to work with certain cases -- like handling non-real roots of ax^2 + bx + c polynomials.
It is possible to represent complex numbers as 2x2 matrices as those can encode 2D rotations. With the matrix formulation you are not dealing with imaginary numbers -- or you are, but they are not encoded with i = sqrt(-1) but as a 45deg rotation. IIRC, there is a formulation of Dirac's QED (Quantum ElectroDynamics) using matrices.