关于matlab影像处理imhist直方图均衡化问题

我们得到一张影像，利用imhist函数我们可以获得一个直方图，横轴是灰度，纵轴是数量，当前问题怎么将一张直方图均衡化。

我的想法是通过函数计算出原最大最小值与0-255之间的映射关系，然后在利用双重循环遍历所有的像素， 将每个像素的灰度值都映射出来，最后的得到的图像就是均衡化后的直方图

I=imread('pout.tif');
for i=1:size(I,1)
    for j=1:size(I,2)
        a=I(i,j);
        I(i,j)=a*1.7-127.5;
    end
end
image(I)

但是我使用该段代码运行时，花了四十分钟只是完成了一般的均衡化，并且出现了黄色rgb色块，可恶