%計算圖像的峰態
clear?all;
I=imread（‘image47.bmp‘）;
[mn]=size（I）;
mean=mean2（I）;
for?i=1:1:m
???for?j=1:1:n
????????sum（ij）=（I（ij）-mean）.^4;
???end
end
meanx=mean2（sum）;
stdx=std2（I）;%std2計算矩陣元素的標準差
ratio=meanx/stdx

%PSNR計算公式


A?=?imread（‘image47.bmp‘）;