PM模型代码Matlab

PM模型的matlab代码

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个文件信息