IIUC, complex numbers are a number system that supports rotations -- one representation is as an angle and a magnitude. As such they work well at describing systems that have rotational components. This makes them useful for working with waves like in QM (light, etc.) and Fourier transformations/analysis (sine waves) which is why they are used in QM.
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.
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.