PDE的一个除噪程序

clc;clear all;close all; I = imread('finger.bmp');I = double(I); figure,imshow(uint8(I)),title('Orignal image'); y = I; sigma = 1;% need some modification according image rol = 3;% need some modification according image dt = 0.1; % timestep, more smaller, more better result gamma = 0.000001; % lamda in the direction of gradient Cm = eps; % very small value Numiter = 100;% iteration num filter = fspecial('gaussian',3,1); y = imfilter(I,filter,'replicate');