Tikhonov-.zip_tikhonov_tikhonov 图像

clear; clc; close all f=imread('lenaNoise.bmp'); f=double(f); u=f; %输入初值可设为f或0 lamda=2; %参数手动输入 [m,n]=size(f); imshow(uint8(f)); for i=1:10 u0j=circshift(u,[0 -1]); u0j(:,n) = u0j(:,n-1); ui0=circshift(u,[1 0]); ui0(1,:) = ui0(2,:); u1j=circshift(u,[0 1]); u1j(:,1) = u1j(:,2); ui1=circshift(u,[-1 0]); ui1(m,:) = ui1(m-1,:); u=(f+lamda*(u0j+ui0+u1j+ui1))/(1+4*lamda); %半隐式迭代 end figure imshow(uint8(u));