mp重构一维信号（可运行）

本人最近研究MP，这个经过修改可以运行的。一维信号是语音信号的

I1=wavread （‘arctic2.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重构‘）
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