图像去噪 TV模型程序

%TV去噪模型程序————2018年4月14日（4月25日程序修改完成可以实现）
clear all       %清空工作空间
I=imread（‘D:\Program Files\CX\图像库\灰度图像\128_128\I13.bmp‘）;   %读入图像I
[mn]=size（I）;       %将图片大小赋值给mn

B=imnoise（I‘gaussian‘0.08）;       %高斯白噪声，方差0.08
figure（1）imshow（B）     %显示噪声图像
B=double（B）;    %双精度
I=double（I）;    %双精度
Q=norm（I-B）;    %原图像u 与噪声图像u0 之差的范数
U0=zeros（m+2n+2）;  %将原图像增大一圈
for i=2:m+1
    for j=2:n+1
       U0（ij）=B（i-1j-1）;      %将图像B的值赋值给U0
    end
end
U=U0;   
h=1;t=0.2;  %步长为1；
r=40;

Ux=zeros（m+2n+2）;  %差分初值都为0
Uy=zeros（m+2n+2）;
Ux1=zeros（m+2n+2）;
Uy1=zeros（m+2n+2）;
d=zeros（m+2n+2）;
d1=zeros（m+2n+2）;
d2=zeros（m+2n+2）;
q=zeros（m+2n+2）;
p=zeros（m+2n+2）;

while （r>0）         %求 U表达式中的一些式子
  for i=2:m+1       
        for j=2:n+1
         Ux（ij）=U（i+1j）-U（ij）;       %Ux前向差分
         Uy（ij）=U（ij+1）-U（ij）;       %Uy前向差分
         Ux1（ij）=U（ij）-U（i-1j）;      %Ux后向差分
         Uy1（ij）=U（ij）-U（ij-1）;      %Uy后向差分