The evolution of both quantum and classical ensembles may be described via the probability density P on configuration space, its canonical conjugate S, and an ensemble Hamiltonian
. For quantum ensembles this evolution is, of course, equivalent to the Schrödinger equation for the wavefunction, which is linear. However, quite simple constraints on the canonical fields P and S correspond to nonlinear constraints on the wavefunction. Such constraints act to prevent certain superpositions of wavefunctions from being realized, leading to superselection-type rules. Examples leading to superselection for energy, spin direction and 'classicality' are given. The canonical formulation of the equations of motion, in terms of a probability density and its conjugate, provides a universal language for describing classical and quantum ensembles on both continuous and discrete configuration spaces, and is briefly reviewed in an appendix.