%?-----------------------------------------------------------------------－－
%???exa130202.m???例13.2.2
%??求沖激函數(shù)、階躍函數(shù)及三角函數(shù)的小波變換、模極大線并求它們的李氏指數(shù)
%??注：在該程序中，用到了子程序?RWT.m，
%??????該程序請(qǐng)讀者在如下的網(wǎng)站上下載:
%??????????????????????http://www-stat.stanford.edu/~wavelab/
%??????因?yàn)樵摼W(wǎng)站上的小波分析軟件并沒(méi)有列入MATLAB的工具箱，涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，
%??????因此不能將其列入國(guó)內(nèi)公開(kāi)出版物上，但讀者可以自由下載。
%------------------------------------------------------------------------
close?all;

nvoice?=?16;
%???nvoice表示每2倍s之間被分為多少點(diǎn)

wavelet?=?‘Sombrero‘;?
%?選擇小波：可選擇的小波包括：?‘Gauss‘‘DerGauss‘‘Morlet‘a(chǎn)nd‘Sombrero‘（mexh）
????N?=?64*8;
????n?=?N;
????oct?=?5;
????scale?=?8;
???	noctave?=?floor（log2（n））-oct;
	nscale??=?nvoice?.*?noctave;
????ytix???=?linspace（2+（oct-floor（log2（scale）））log2（N）+2-floor（log2（scale））nscale）;
%??ytix?從imageRWT函數(shù)得到?理論上應(yīng)等于log2（s）
%??ytix?=?linspace（log2（scale）-noctavelog2（scale）nscale）;
????xtix???=?linspace（0NN）;
for?count=1:3
???if?count==1
????????sig?=?zeros（1N）;sig（N/2）?=sig（N/2）?+1;%?沖激函數(shù)
????????ti=‘沖激函數(shù)‘;
????elseif?count==2
????????sig?=?zeros（1N）;sig（N/4:3*N/4）?=?sig（N/4:3*N/4）+1;?%?矩形窗
????????ti=‘矩形窗‘;
????else
????????sig（N/4:N/2）?=?linspace（01N/4+1）;
????????sig（N/2:3*N/4）?=?linspace（10N/4+1）;?%?三角窗
????????ti=‘三角窗‘;
????end
????sig?=?sig+1e-8*randn（1N）;
????r