%下面的程序使用一個子程序來選擇最佳的原子；
%項目，比例尺，平移，頻率，相位
function?[projscaletranslationfreqphasesize_dic]=select_best（signal_rNa_basej_minj_maxu_basep_minv_basek_minw_basei_mini_max）;
proj_trans=0;
proj=0;
%%是一個參數，用于顯示所用過全字典的大小?
size_dic=0;

for?j=j_min:j_max
????for?p=p_min:N*2^（-j+1）
????????for?k=k_min:2^（j+1）
????????????for?i=i_min:i_max
????????????????
??????????????size_dic=size_dic+1;
??????????????s=a_base^j;??%尺度因子??r?=?（suw?）?=?（a^j??pa^jΔuka^-jΔwiΔ?）
??????????????u=p*s*u_base;%平移因子
??????????????v=k*（1/s）*v_base;%頻率
??????????????w=i*w_base;
??????????????t=0:N-1;
???????????????%?t=1366.28:0.1:1446.18;
??????????????t=（t-u）/s;
??????????????
??????????????g=（1/sqrt（s））*exp（-pi*t.*t）.*cos（v*t+w）;?%?獲得原子
??????????????g=g/sqrt（sum（g.*g））;?%?歸一化??
??????????????proj_trans=sum（signal_r.*g）;?%?計算內積
??????????????
????????????????if?abs（proj_trans）>abs（proj）%?判斷是否為最佳原子
?????????????????proj=proj_trans;
?????????????????scale=s;
?????????????????translation=u;
?????????????????freq=v;
?????????????????phase=w;
????????????????end
?????????????end
?????????end
????end
end
?%disp（‘the?size?of?the?dictionary=‘）disp（size_dic）;??%?size_dic:?the?size?of?the?dictionary
end

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

?????文件??????16897??2020-03-29?17:39??gabor_mp2\R192C056.picotd

?????文件???????1324??2020-03-29?23:31??gabor_mp2\select_best.m

?????文件???????3005??2020-03-30?09:52??gabor_mp2\Untitled.m

?????目錄??????????0??2020-04-05?20:42??gabor_mp2

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

????????????????21226????????????????????4