%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%?????子程序?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