%f1?and?f2?are??two?consecutive?frames?of?a?video?sequence
clc;
clear;
close?all



f1=imread（‘CT0010.jpg‘）;
f40=imread（‘CT0011.jpg‘）;

figure;
imshow（f1）
figure
imshow（f40）


fr_f1=rgb2gray（f1）;
fr_f40=rgb2gray（f40）;?????
Xn=double（fr_f1）;
Xnp1=double（fr_f40）;
?
%get?image?size?and?adjust?for?border??獲取圖像對邊界進行調整
[hw]=size（fr_f1）;??
hm5=h-5;?wm5=w-5;
z=zeros（hw）;?v1=z;?v2=z;

%initialize????????初始化
dsx2=v1;?dsx1=v1;?dst=v1;
alpha2=625;
imax=20;

%Calculate?gradients??計算梯度
dst（5:hm55:wm5）?=?（?Xnp1（6:hm5+16:wm5+1）-Xn（6:hm5+16:wm5+1）?+?Xnp1（6:hm5+15:wm5）-Xn（6:hm5+15:wm5）+?Xnp1（5:hm56:wm5+1）-Xn（5:hm56:wm5+1）?+Xnp1（5:hm55:wm5）-Xn（5:hm55:wm5））/4;
dsx2（5:hm55:wm5）?=?（?Xnp1（6:hm5+16:wm5+1）-Xnp1（5:hm56:wm5+1）?+?Xnp1（6:hm5+15:wm5）-Xnp1（5:hm55:wm5）+?Xn

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2016-03-09?14:50??matlab\
?????文件????????2363??2016-03-09?15:55??matlab\aa.asv
?????文件????????2363??2016-03-09?14:40??matlab\aa.m
?????文件??????190100??2007-11-06?22:57??matlab\car0.png
?????文件??????191767??2007-11-06?22:57??matlab\car1.png
?????文件???????15776??2012-05-17?21:25??matlab\CT0010.jpg
?????文件???????15843??2012-05-17?21:25??matlab\CT0011.jpg
?????文件???????15913??2012-05-17?21:25??matlab\CT0012.jpg
?????文件???????15919??2012-05-17?21:25??matlab\CT0013.jpg
?????文件???????15584??2012-05-17?21:25??matlab\CT0014.jpg
?????文件???????24064??2016-03-04?10:33??matlab\Thumbs.db
?????文件?????????138??2013-03-09?16:46??matlab\www.hslogic.com.txt