function?[msepsnr]?=?psnr_mse（imageimage_prime）

????%?convert?to?doubles
????image=double（image）;
????[MN]=size（image）;
????image_prime=double（image_prime）;

????%?avoid?divide?by?zero?nastiness
????if?（（sum（sum（（image-image_prime）.^2）））?=