基于MP的时频分析MATLAB代码

用MP对信号和图像进行时频分析MTLAB代码,详细地编写了时频分析过程

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%     子程序 dic_a 构建时频原子库       %%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%% 这部分的程序的主要思想是构造一个合适大小的矩阵，将范围内所有的时频原子填入矩
%%%% 阵中，即用矩阵的形式构造出原子库。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [atoms]=dic_a（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-3
%   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;             %% 计算原子个数
       end
      end
%   end
end
atoms=zeros（size_dicN）;                   %% 创建size_dic×N的矩阵存放该所有原
size_dic=0;                                %% 子库的所有原子并清零size_dic
for j=j_min:j_max-3
%   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;
            u=N/2;
            v=k*（1/s）*v_base;
            w=i*w_base;
            t=0:N-1;                                  %%把所有原子写入矩阵
            t=（t-u）/s;                                %%  atom中，构成原子库
        g1=（1/sqrt（s））*exp（-pi*t.*t）.*cos（v*t+w）;
        g=g1/sqrt（sum（g1.*g1））;
            atoms（size_dic:）=g;
       end
      end
%   end
end
 size_dic 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1456  2008-06-06 21:58  时频分析\dic_a.asv

     文件       1456  2008-06-06 22:03  时频分析\dic_a.m

     文件       3703  2008-06-06 21:58  时频分析\main.asv

     文件       3703  2008-06-06 22:02  时频分析\main.m

     文件       1712  2008-06-06 21:23  时频分析\select_best.m

     目录          0  2009-04-05 22:15  时频分析

----------- ---------  ---------- -----  ----

                12030                    6