A new method of calculating nonequilibrium density matrices with the aid of the quantum integrals of motion is proposed. The method is shown to be very effective in the case of systems described by means of quadratic Hamiltonians. The possibility of constructing phenomenological nonstationary Hamiltonians for a wide class of dissipative systems is discussed. The exact formulas for nonequilibrium density matrices of arbitrary quadratic systems are obtained. The quantum problem of the motion of a charged particle in uniform electric and magnetic fields in the presence of a frictional force proportional to the velocity is solved exactly by means of introducing the new phenomenological Hamiltonian.