%?img=imread（‘./lena.bmp‘）;
%?img=double（img）;
%?y=fspecial（‘gaussian‘size（img）11）;
%?w_f=（?fft2（img）.*fft2（y）?）?./?（?fft2（img）.*fft2（img）?）;
%?pred=ifft2（fft2（circshift（img[100100]））.*w_f）;
%?imshow（pred[]）????%超簡易版公式驗證

clear?all
clc
%將圖片放到一個文件夾下
names=dir（‘.\img‘）;?????
%獲取圖片名稱列表
names（1:2）=[];
img=double（rgb2gray（imread（[‘.\img\‘names（1）.name]）））;??????%讀取第一張圖
loc_len=[302008040];????????????????????????????
%目標框左上角的列坐標，行坐標，目標框列數(shù)，目標框的行數(shù)
patch=img（loc_len（2）:loc_len（2）+loc_len（4）loc_len（1）:loc_len（1）+loc_len（3））;??????????
%獲得目標
label=fspecial（‘gaussian‘size（patch）3）;????
%生成標簽
w_fft2_conj=fft2（patch）.*（fft2（label））./（（fft2（patch））.*fft2（patch））;????????????????????????
%訓練得到權(quán)重w的傅里葉變換（沒必要變到時域，因為一會兒直接用w的頻域矩陣）
for?i=1:size（names1）
????img=double（rgb2gray（imread（[‘.\img\‘names（i）.name]）））;??
????%讀圖，轉(zhuǎn)灰度，轉(zhuǎn)double
????tmp_patch=img（loc_len（2）:loc_len（2）+loc_len（4）loc_len（1）:loc_len（1）+loc_len（3））;??
????%獲得搜索范圍
????pred=ifft2（fft2（tmp_patch）.*（w_fft2_conj））;????????
????%使用w預(yù)測
????[rowcol]=find（pred==max（max（pred）））;??????????
????%找到相應(yīng)最大的位置
????loc_len（1）=loc_len（1）?+?col-loc_len（3）/2;??????
????%更新列坐標，原有列坐標=原有列坐標+響應(yīng)矩陣最高點列坐標-響應(yīng)矩陣中心點列坐標
????loc_len（2）=loc_len（2）?+?row-loc_len（4）/2;??????
????%更新行坐標，原有行坐標=原有行坐標+響應(yīng)矩陣最高點行坐標-響應(yīng)矩陣中心點行坐標
????imshow（img[]）
????rectangle（‘Position‘loc_len‘EdgeColor‘‘r‘）??
????%繪制矩形框
????pause（0.03）
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-05-13?22:07??img\
?????文件???????10547??2015-12-12?01:35??img\frame_101.jpg
?????文件???????10211??2015-12-12?01:35??img\frame_102.jpg
?????文件???????10204??2015-12-12?01:35??img\frame_103.jpg
?????文件???????10170??2015-12-12?01:35??img\frame_104.jpg
?????文件???????10079??2015-12-12?01:35??img\frame_105.jpg
?????文件????????9975??2015-12-12?01:35??img\frame_106.jpg
?????文件????????9990??2015-12-12?01:35??img\frame_107.jpg
?????文件????????9994??2015-12-12?01:35??img\frame_108.jpg
?????文件????????9974??2015-12-12?01:35??img\frame_109.jpg
?????文件???????10013??2015-12-12?01:35??img\frame_110.jpg
?????文件????????9967??2015-12-12?01:35??img\frame_111.jpg
?????文件???????10049??2015-12-12?01:35??img\frame_112.jpg
?????文件???????10022??2015-12-12?01:35??img\frame_113.jpg
?????文件????????9932??2015-12-12?01:35??img\frame_114.jpg
?????文件????????9891??2015-12-12?01:35??img\frame_115.jpg
?????文件????????9955??2015-12-12?01:35??img\frame_116.jpg
?????文件???????10102??2015-12-12?01:35??img\frame_117.jpg
?????文件???????10101??2015-12-12?01:35??img\frame_118.jpg
?????文件???????10212??2015-12-12?01:35??img\frame_119.jpg
?????文件???????10168??2015-12-12?01:35??img\frame_120.jpg
?????文件???????10251??2015-12-12?01:35??img\frame_121.jpg
?????文件???????10157??2015-12-12?01:35??img\frame_122.jpg
?????文件???????10207??2015-12-12?01:35??img\frame_123.jpg
?????文件???????10253??2015-12-12?01:35??img\frame_124.jpg
?????文件???????10376??2015-12-12?01:35??img\frame_125.jpg
?????文件???????10400??2015-12-12?01:35??img\frame_126.jpg
?????文件???????10307??2015-12-12?01:35??img\frame_127.jpg
?????文件???????10340??2015-12-12?01:35??img\frame_128.jpg
?????文件???????10369??2015-12-12?01:35??img\frame_129.jpg
?????文件???????10300??2015-12-12?01:35??img\frame_130.jpg
............此處省略55個文件信息