function?[ImMAEPSNRRJ]=diffusion_PM（IoInNtauK）


[NxNy]=size（In）;?
a=5;
a1=5;
b=5;
b1=5;

J=In;

I_original_loc=Io（a:（Nx-a1）b:（Ny-b1））;


t=0;
for?i=1:N;???
??i
tic
??
%???%%%div的模式1
??%?coeffice
????cn?=?J（[1?1:end-1]:）?-?J;
????cs?=J（[2:end?end]:）?-?J;
????cw?=?J（:?[1?1:end-1]）?-J;
????ce?=J（:?[2:end?end]）?-?J;
?%
????div?=?1./（1+（abs（cn）/K）.^2）?.*?cn...
??????????+?1./（1+（abs（cs）/K）.^2）?.*?cs...
??????????+?1./（1+（abs（cw）/K）.^2）?.*??cw...
??????????+?1./（1+（abs（ce）/K）.^2）?.*?ce;

??????
??????
????J=?J+?tau?*?div;
%?????pause（0.01）
%?????hold?on
%?????figure（111）imshow（J[]）
%?????hold?off
toc;

t=t+toc;
????%%取圖像的一部分計算SNR
J_loc=J（a:（Nx-a1）b:（Ny-b1））;
???
??????PSNRR（i）=psnr（J_locI_original_loc）;
???NowPSNR=PSNRR（i）
???NowSNR=snr（J_locI_original_loc）
???NowMAE=mae（J_locI_original_loc）
???ImMAE（i）=NowMAE;
end;?%?for?i

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2016-07-04?21:40??PM-code\
?????文件?????????921??2016-05-09?08:55??PM-code\diffusion_PM.m
?????目錄???????????0??2016-06-18?18:54??PM-code\Image\
?????文件??????196662??2011-02-25?10:26??PM-code\Image\Aerial_Image.bmp
?????文件??????270054??2008-03-04?15:31??PM-code\Image\Article3.bmp
?????文件??????263222??2009-12-06?09:08??PM-code\Image\barb.bmp
?????文件???????66614??2009-12-06?09:42??PM-code\Image\barb256.bmp
?????文件??????786486??2010-11-03?21:08??PM-code\Image\Barbara.bmp
?????文件??????750054??2011-04-06?15:53??PM-code\Image\Barbara500.bmp
?????文件??????263222??2009-12-06?09:11??PM-code\Image\boat.bmp
?????文件??????196662??2009-12-06?09:12??PM-code\Image\cameraman.bmp
?????文件???????44982??2011-03-09?10:51??PM-code\Image\f1.bmp
?????文件??????786486??2009-12-06?09:11??PM-code\Image\fingerprint.bmp
?????文件??????121458??2011-03-07?10:11??PM-code\Image\fingerprint1.bmp
?????文件??????196662??2009-12-06?09:43??PM-code\Image\fingerprint256.bmp
?????文件?????????144??2015-11-12?13:18??PM-code\Image\Gnoise2.m
?????文件???????17462??2009-08-24?22:06??PM-code\Image\GZC01.bmp
?????文件?????????245??2015-11-12?12:47??PM-code\Image\jiazao.m
?????文件?????????364??2015-11-13?22:34??PM-code\Image\jiazao1.m
?????文件?????????280??2015-11-12?14:11??PM-code\Image\LaplaceDenoise.m
?????文件??????263222??2009-12-06?09:11??PM-code\Image\lena.bmp
?????文件???????66614??2009-12-06?09:42??PM-code\Image\lena256.bmp
?????文件???????91078??2011-03-30?20:11??PM-code\Image\lena300.bmp
?????文件???????13602??2008-03-08?13:57??PM-code\Image\lenapart.bmp
?????文件??????270054??2010-03-14?00:02??PM-code\Image\mdn4.bmp
?????文件??????270054??2010-02-03?11:53??PM-code\Image\Multi_N2.bmp
?????文件?????????599??2015-11-12?14:36??PM-code\Image\noise.m
?????文件???????36906??2011-10-25?08:06??PM-code\Image\PDE_seismic1.bmp
?????文件???????60774??2011-10-25?08:06??PM-code\Image\PDE_seismic2.bmp
?????文件???????35246??2011-10-25?08:06??PM-code\Image\PDE_seismic3.bmp
?????文件??????196662??2003-11-06?00:25??PM-code\Image\pepper.bmp
............此處省略26個文件信息