%先用循环求出Q(i)
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
%再进行灰度变换
I=imread('D:\学习\7.bmp');
for i=0:255
for j=0:255
M=I(i,j)*Q(I(i,j))
Index into matrix is negative or zero. See release notes on changes to
logical indices.