%?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）

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????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