Constrained energy minimization (CEM) has been extended to several generalized versions, iterative CEM (ICEM), Nystrom method-based kernel CEM (NKCEM), and iterative training sampling-based NKCEM (ITS-NKCEM) for hyperspectral image classification (HSIC). Since CEM is a subpixel target detector that specifies a target signature to detect its target abundance fractions present in data samples, this article takes the advantage of CEM’s ability in subpixel detection to consider NKCEM as a nonlinear mixed pixel classifier for hyperspectral mixed pixel classification (HMPC). Recently, a new concept of band sampling (BSam) was proposed by utilizing random signal sampling derived from compressive sensing (CS) to show its performance better than band selection (BSel) for HSIC. Thus, incorporating BSam into NKCEM and ITS-NKCEM yields two new versions for HPMC, called band sampling NKCEM (BSam-NKCEM) and band sampling ITS-NKCEM (BSam-ITS-NKCEM). Interestingly, despite that BSam-ITS-NKCEM uses sampled bands as well as training samples, extensive experiments demonstrate that it can perform better than all other CEM and KCEM versions using full bands and ground truth. Specifically, it also shows that BSam-ITS-NKCEM can do better than spectral–spatial HSIC techniques using BSel.