clear all close all f = imread('./res/test.jpg'); [m, n] = size(f); pf = zeros(255, 1);%概率函数 maxpixels = size(f, 1) * size(f, 2); for i = 1:maxpixels pf(f(i)) = pf(f(i)) + 1; end pf = pf / maxpixels; c = zeros(255, 1);%对应于pf的累计概率函数 for i = 1:255 for j = 1:i c(i) = c(i) +pf(j); end end g = zeros(size( f ));%目标图像 for i = 1:maxpixels g(i) = c(f(i)); end imshow(f, []); figure imshow(g, []);