test_TV.rar_TV_去高斯噪声_高斯噪声

function [ obj, X1 ] = BTV( X0,Y,stepLength,lambda,P,a) %UNTITLED Summary of this function goes here % Detailed explanation goes here X0 = double(X0); Y = double(Y); BTV_term = 0; dBTV = zeros(size(X0,1),size(X0,2),'double'); for i=-P:P for j=0:P if(i+j >= 0) Xshift = circshift(X0,[j,i]); BTV_term = BTV_term + power(a,abs(i)+abs(j))*sum(sum(abs(X0-Xshift))); Xsign = sign(X0 - Xshift); dBTV_temp = Xsign; dBTV_temp = dBTV_temp - circshift(Xsign,[-j,-i]); dBTV = dBTV + power(a,abs(i)+abs(j))*dBTV_temp; end end end obj = sum(sum((Y-X0).*(Y-X0)))/2; obj = obj + lambda * BTV_term; X1 = X0 - stepLength*((X0 - Y) + lambda*dBTV); %X1 = X0 - stepLength*lambda*dBTV; end