%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后向差分
???