gabor字典

gabor字典

%下面的程序使用一个子程序来选择最佳的原子；
%项目，比例尺，平移，频率，相位
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