I1=wavread?（‘a(chǎn)rctic2.wav‘）‘;?
I1=I1（6021:6276）;

bat=I1;??
iterative_number=2;
[aN]=size（bat）;
signal_reconstruct=zeros（1N）;
signal_r=bat;
a_base=2;
j_min=1;
j_max=log2（N）;
u_base=1/2;
p_min=0;
v_base=pi;
k_min=0;
w_base=pi/6;
i_min=0;
i_max=12;



for?n=1:iterative_number
?

[projscaletranslationfreqphase]=select_best（signal_rNa_basej_minj_maxu_basep_minv_basek_minw_basei_mini_max）;


t=0:N-1;

t=（t-translation）/scale;

?g=（1/sqrt（scale））*exp（-pi*t.*t）.*cos（freq*t+phase）;
????????????
?g=g/sqrt（sum（g.*g））;
?
?figure（3）
?plot（g）
?hold?on

signal_reconstruct=signal_reconstruct+proj*g;

signal_r=signal_r-proj*g;


figure（2）
subplot（311）;plot（I1）;title（‘信號I1‘）
subplot（312）;plot（g）;title（‘Gabor字典中的最佳原子‘）
subplot（313）;plot（signal_reconstruct）;title（‘I1MP重構(gòu)‘）
y=signal_reconstruct;
end

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

?????文件????????873??2014-03-17?19:24??mp\MP1.m

?????文件????????985??2014-03-15?15:52??mp\select_best.m

?????目錄??????????0??2014-03-17?19:24??mp

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

?????????????????1858????????????????????3