This article deals with a regularized version of the split gradient method (SGM), leading to multiplicative algorithms. The proposed algorithm is available for the optimization of any divergence depending on two data fields under positivity constraint. The SGM-based algorithm is derived to solve the nonnegative matrix factorization (NMF) problem. An example with a Frobenius norm on both the data consistency and the penalty term is developped and applied to hyperspectral data unmixing.