Statistically averaged rate equations for intense non-neutral beam propagation through a periodic solenoidal focusing field based on the nonlinear Vlasov–Maxwell equations