function?Sxy=jnd（A）
%最小可覺差Jnd計算，及人眼刺激圖計算。
%A=imread（‘D:\我的文檔\MATLAB\新建文件夾\評價方法\CCD.bmp‘）;
%計算背景灰度圖I
[m?n]=size（A）;
%A=rgb2gray（A）;
A=double（A）;
for?i=1:m
????for?j=1:n
???????flag=2????????if?flag
????????????for?a=1:5
????????????????for?b=1:5
????????????????????C（ab）=A（i-3+aj-3+b）;
????????????????end
????????????end
????????????C=double（C）;
????????????B=[11111;12221;12021;12221;11111];?
????????????II=（1/32）*C.*B;
????????????I（ij）=sum（sum（II））;
????????else
????????????I（ij）=A（ij）;
????????end
????end
end
I=double（I）;
%計算某像素點處，其背景灰度值所對應的jnd值
for?i=1:m
????for?j=1:n
????????if?0<=I（ij）&I（ij）<=60
????????????E（ij）=log2（1+（A（ij）-I（ij））/I（ij））;
????????????E0（ij）=log2（1+（-0.009*I（ij）+0.575））;
????????????Sxy（