Jd=diffusion(J,method,N,K) J :待扩散的原图像 method : lin : 线性扩散 (常数 c=1). pm1 : perona-malik, c=exp{-(|grad(J)|/K)^2} pm2 : perona-malik, c=1/{1+(|grad(J)|/K)^2} rmp : 复数 K :边缘门限参数 N :迭代次数 dt:时间步长 (0 < dt <= 0.25, 缺省值为0.2)