mp_least_square_gpr_2.m

基于morlet小波的最小二乘匹配追踪算法

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%   mp-least square algorithm
%import GPR signal
receiver_1=load('e:\gpr_test\case 1\receiver_1.txt');
t=(receiver_1(2:(receiver_1(1,1)+1),1))';    %sample time
s=(receiver_1(2:(receiver_1(1,2)+1),2))'; %signal
global n atom_num 
[m,n]=size(t);
atom_num=0;
figure_num=1;
figure (figure_num)
plot(t,s,'k-');
%hold on
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%   ST time-frequency distribution
%[st,t1,f] = st(real(s),0,3*1.0e6,Ts,1);
%st_abs=abs(st);
%figure_num=figure_num+1;
%figure (figure_num)
%contourf(t,f,st_abs);
%grid on
%xlabel('t(s)');
%ylabel('f(Hz)');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%   WVD distribution
%[tfr,t1,f1] = tfrwv(signal_real);
%tfr=fftshift(tfr,1);
%f=fs*linspace(0,1,nb);
%figure_num=figure_num+1;
%figure (figure_num)
%contourf(t,f,tfr);
%grid on
%xlabel('t(s)');
%ylabel('f(Hz)');
%----------------------------------------function 1
t_scan_num=3; 
f_scan_num=3;
threshold=0.0001;
amp=1;  
f_width=0.09; 
delt=0.5*1.0e-19;   
k=0.02;     
ebs=0;    
%----Run Time test
tic
[t_array,f_array,a_array,phi_array,select_array,residual] = LS(t,s,t_scan_num,f_scan_num,threshold,amp,f_width,delt,k,ebs);
toc
run_time=toc;
%compute error & draw the restruction
[atom_num,n]=size(select_array);
restruction=zeros(1,n);
for ii=1:atom_num
    restruction=restruction+select_array(ii,1:n);
end
error=s-restruction;
figure (figure_num)
plot(t,real(restruction),'b-');
hold on
plot(t,error,'r-');
figure_num=figure_num+1;
figure (figure_num)
error_dB=20*log10((real(error)./real(s)));
plot(t,error_dB,'b-');