A majorization-minimization based solution to penalized nonnegative matrix factorization with orthogonal regularization