function?[uv]=lk（img1img2wsize）;
[fxfyft]=deriv（img1img2）;
u=zeros（size（img1））;
v=zeros（size（img2））;
half_w=floor（wsize/2）;
for?i=half_w+1:size（fx1）-half_w
????for?j=half_w+1:size（fx2）-half_w
????????curx=fx（i-half_w:i+half_wj-half_w:j+half_w）;
????????cury=fy（i-half_w:i+half_wj-half_w:j+half_w）;
????????curt=ft（i-half_w:i+half_wj-half_w:j+half_w）;
????????curx=curx‘;
????????cury=cury‘;
????????curt=curt‘;
????????curx=curx（:）;
????????cury=cury（:）;
????????curt=-curt（:）;
????????a=[curxcury];
????????UV=pinv（a‘*a）*a‘*curt;
????????u（ij）=UV（1）;
????????v（ij）=UV（2）;
????end;
end;
u（isnan（u））=0;
v（isnan（v））=0;

function?[fxfyft]=deriv（img1img2）;
fx=conv2（img10.25*[-11;-11]）+conv2（img20.25*[-11;-11]）;
fy=conv2（img10.25*[-1-1;11]）+conv2（img20.25*[-1-1;11]）;
ft=conv2（img10.25*ones（2））+conv2（img2-0.25*ones（2））;

fx=fx（1:size（fx1）-11:size（fx2）-1）;
fy=fy（1:size（fy1）-11:size（fy2）-1）;
ft=ft（1:size（ft1）-11:size（ft2）-1）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1922??2018-04-10?17:56??optial_flow.m

?????文件???????1516??2018-04-13?10:19??optical_flow_detect.m

?????文件????????995??2012-05-25?16:22??lk.m

-----------?---------??----------?-----??----

?????????????????4433????????????????????3