Документ взят из кэша поисковой машины. Адрес оригинального документа : http://lnfm1.sai.msu.ru/grav/russian/lecture/filtr/labs/Lab5.m Дата изменения: Thu Mar 24 18:08:42 2011 Дата индексирования: Mon Oct 1 23:41:10 2012 Кодировка: Поисковые слова: п п п п п п п п п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п р п

clear;
% install Wavelab http://www-stat.stanford.edu/~wavelab/Wavelab_850/download.html
%and put the path

cd '/home/leonid/lectures/labs/Wavelab850'%path to wavelab
WavePath

N_signal=1024;
% generating two-sin signal
for (k=1:1:N_signal)
if k<(N_signal/2)
garm(k)=5*sin(2*pi/20*(k-1));
else
garm(k)=5*sin(2*pi/50*(k-1));
end;
chirp(k)=5*sin(2*pi/5000*(k-1)^2);

end;

plot(garm);

% ARMA process generating
ar(1)=0.5*randn(1);
ar(2)=-0.2*ar(1)+0.5*randn(1);

for (i=3:1:N_signal)
ar(i)=0.9*ar(i-1)-0.7*ar(i-2)+0.5*randn(1);
end;

plot(ar);

signal=garm+ar+chirp;

plot(signal);

plot(chirp);
cd 'Orthogonal';

QMF_Filter = MakeONFilter('Coiflet',2);
wcNaive = FWT_PO(signal, 5, QMF_Filter);

%PlotWaveCoeff(wcNaive,1,0.5)
DisplayMultiRes(wcNaive,1,0.5,QMF_Filter,'Mesh')
%--------------------------------------------------------


specgm = WindowFT(signal,100,20,'Hamming','wft');
%--------------------------------------------------------

rwt = RWT(signal,7,'Morlet');

c = CWT_Wavelab(signal,5,'Morlet');
colors='20';
ImageCWT(c,'Individual',colors,'lin')