?? wt07fig08.m
字號:
% CAPTION
fprintf('\n');
disp('Figure 7.8')
disp('The left graph shows a Meyer wavelet.')
disp('The modulus of its Fourier transform is on the right.')
close all;
n = 8192;
M = 2*pi;
omega = ((1:n) - n/2)* M / n;
hath = Meyerhath(omega);
hatg = exp(-i*pi) .* Meyerhath(omega + pi);
M = 8*pi*20;
omega = ((1:n) - n/2)* M / n;
hatpsi = zeros(1,n);
supp1 = find((abs(omega) > 2*pi/3) & (abs(omega) <= 4*pi/3));
hatpsi(supp1) = exp(-i.*(omega(supp1)./2)).*Meyerhath(omega(supp1)./2 + pi)./sqrt(2);
supp2 = find((abs(omega) > 4*pi/3) & (abs(omega) <= 8*pi/3));
hatpsi(supp2) = exp(-i.*omega(supp2)./2).*Meyerhath(omega(supp2)./4)./sqrt(2);
figure(1);clf
subplot(222)
plot(omega,abs(hatpsi))
axis([-4*pi 4*pi 0 1.1]);
psi = zeros(1,n);
[def,psi] = NormedIFFT(hatpsi,M);
supp = find(abs(psi) > max(abs(psi))./10000);
subplot(221)
plot(def(supp),real(psi(supp)));
axis([-8 8 -1 1.3]);
% Written by Maureen Clerc and Jerome Kalifa, 1997
% clerc@cmapx.polytechnique.fr, kalifa@cmapx.polytechnique.fr
%
% 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 wavelab@stat.stanford.edu
%
?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -