% -------------------------------------------------------------------
% exa090603.m, for fig 9.6.3
% 產生 mexhat 小波的時域波形及頻譜 
% 注：在該程序中，用到了子程序 mexhat.m，
%      該程序請讀者在如下的網站上下載:
%                      http://crttsn.univ-nantes.fr/~auger/tftb.html
%      因為該網站上的時頻分析軟件并沒有列入MATLAB的工具箱，涉及到知識產權問題，
%      因此不能將其列入國內公開出版物上，但讀者可以自由下載。
%--------------------------------------------------------------------
clear;

[wavelet]=mexhat(0.1);
subplot(2,2,1)
plot(wavelet);grid on
%title(' mexhat wavelet: Psi')

a=[1];
   fs=1;
[h,w]=freqz(wavelet,a,256,'whole',fs);
hr=abs(h);
subplot(2,2,2)
plot(w(1:128),hr(1:128));grid on
%title('The FT of Psi')
[fm,B]=locfreq(wavelet')