I=imread('D:\学习\22.jpg');
imshow(I)
title('原始图象')
B=imhist(I,256) ;
for i=1:256
if i<64
N(i)=8+(-3/32)*i;
Fd(i)=0;
else if i>192
N(i)=(1/64)*i-1;
Fd(i)=1;
else if 64<i<192
N(i)=2;
Fd(i)=i/128-0.5;
end
end
end
Fr(i)=N(i)/8;
Q(i)=Fr(i)*Fd(i);
end
P=Q';
for i=1:256
A(i)=B(i)*P(i);
end
imshow(A)
title('处理后的图象')