谐波小波程序

% hwt.m: harmonic wavelet transform % whwt:windowed harmonic wavelet transform % Gabor_wt: Gabor wavlet transform

% test computational time of estimating phase with a simple model

clear
close all
clc
fs=100;% sampling frequency
N=10*fs;% data length
snr=1;
% % % error with different realization of noise
for p=1:1000 % realizations
    for i=1:N
        phase_true（i）=2*pi*5*（i-1）/fs+pi/4;
        x（i）=2*cos（phase_true（i））;% true signal
    end
    rx=awgn（x-mean（x）1）+mean（x）;

    tic;
    x_whwt=whwt（rx37fsN）;%extract by HWT
    phase_whwt=angle（x_whwt）;
    time_whwt（p）=toc;
    
    tic;
    x_Gwt=Gabor_wt（rx37fsN）;
    phase_Gwt=angle（x_Gwt）;
    time_Gwt（p）=toc;
    
    tic;
    x_filt=eegfilt（rxfs37）;% by FIR          
    phase_eegfilt=angle（hilbert（x_filt））;
    time_filt（p）=toc;
    
end
time_whwt=time_whwt（find（time_whwt<1））;
time_Gwt=time_Gwt（find（time_Gwt<1））;
time_filt=time_filt（find（time_filt<1））;

mean（time_whwt）std（time_whwt）
mean（time_Gwt）std（time_Gwt）
mean（time_filt）std（time_filt）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2009-11-04 15:15  HWT_PS\
     文件         961  2009-11-01 15:22  HWT_PS\compute_time.m
     文件        1583  2009-11-01 14:44  HWT_PS\Gabor_wt.m
     文件         813  2009-11-04 15:14  HWT_PS\help.m
     文件        1326  2009-11-01 14:44  HWT_PS\hwt.m
     文件        2358  2009-11-01 15:16  HWT_PS\test_PS.m
     文件        1512  2009-11-01 15:20  HWT_PS\test_simple.m
     文件        1601  2009-11-01 14:44  HWT_PS\whwt.m
     文件        1328  2009-11-04 02:18  license.txt