wt04fig06.m

来自「斯坦福大学的小波的工具箱 比matlab自带的工具箱更完善」· M 代码 · 共 38 行

% CAPTION
fprintf('\n');
disp('Figure 4.6')
disp('Left: graph of the Mexican hat wavelet, equal to minus the second')
disp('derivative of a Gaussian.')
disp('Right: graph of the Fourier transform of the Mexican hat wavelet.')

close all;
N = 1024;
t = -5 + 10*(1:N)/N;
om = -2*pi + 4*pi*(1:N)/N;
pi4 = pi^.25;
sqrt3 = sqrt(3);

psi = -2./(pi4.*sqrt3).*exp((-t.^2)./2).*(t.^2-1);
hatpsi = sqrt(8).*pi4./sqrt3.*(om.^2).*exp(-(om.^2)./2);

figure(1);clf
subplot(221)
	plot(t,psi);

subplot(222)
	plot(om,hatpsi);
	axis([-5 5 -.2 1.7])
	
% Written by Maureen Clerc and Jerome Kalifa, 1997
% [e-mail], [e-mail]
    
    
%   
% Part of WaveLab Version 802
% Built Sunday, October 3, 1999 8:52:27 AM
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%