稀疏分解matlab程序

load bat -ascii; signal=bat; iterative_number=30; %the matching pursuit processing iterative number [a,N]=size(signal); %the length of signal and the length of atoms;N signal_reconstruct=zeros(1,N); signal_r=signal; a_base=2; %scale: the scale is determined by a and j; j_min=0; j_max=log2(N); u_base=1/2; %the transmission or displacement u should determined as following v_base=pi; %the frequency v k_min=0; w_base=pi/6; %the phase w i_min=0; i_max=12; signal_reconstruct=(1/N)*sum(signal); signal_r=signal-signal_reconstruct; [atoms]=dic_a(signal_r,N,a_base,j_min,j_max,u_base,v_base,k_min,w_base,i_min,i_max); for n=1:iterative_number n [proj,scale,translation,freq,phase]=select_best(signal_r,N,a_base,j_min,j_max,u_base,v_base,k_min,w_base,i_min,i_max,atoms); t=0:N-1; t=(t-translation)/scale; g1=(1/sqrt(scale))*exp(-pi*t.*t).*cos(freq*t+phase); g=g1/sqrt(sum(g1.*g1)); signal_reconstruct =signal_reconstruct+proj*g; signal_r=signal_r-proj*g; subplot(221); plot(signal); %original signal subplot(222); plot(g); %the atom subplot(223); plot(signal_r); %the residual subplot(224); plot(signal_reconstruct); %reconstruct signal end