The present invention discloses a digital halftoning method. The method comprises steps of: (a1) dividing an original image into non-overlapping blocks; (a2) obtaining a Least-Mean-Square trained (LMS-trained) filter by comparing at least a training image and a halftone result corresponding to the training image (a3) optimizing a class matrix with the LMS-trained filter, which involves the diffused area and the diffused weightings; and (a4) processing the non-overlapping blocks by performing a dot diffusion procedure with the optimized class matrix and the diffused weightings to generate a halftone image corresponding to the original image. A detailed class matrix optimizing method as in the above-mentioned step (a3) is also disclosed.