Документ взят из кэша поисковой машины. Адрес оригинального документа : http://lnfm1.sai.msu.ru/grav/english/lecture/filtering/Lab4.m Дата изменения: Wed Feb 4 20:10:49 2009 Дата индексирования: Mon Oct 1 23:38:48 2012 Кодировка: Поисковые слова: m 44

%wavelet-transformation of the signal
% program is written 27.01.2009 by L.V. Zotov
clear;

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;
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);

%making a sum
signal=garm+ar;

plot(signal);


cd C:\work\course\filtr\eng\Lab4;

savefile = 'signal.mat';
save(savefile,'signal');

%------------------------------------------------
a_max=128;

c = cwt(signal,[1:a_max],'morl','plot');

clf;

%3d-plot
[A,B] = meshgrid([1:100],[1:a_max]);
for(i=1:1:a_max)
%morl
scale_ampl=sqrt(2*pi);
%mexh
% scale_ampl=sqrt(2*pi)*i;
for(j=1:1:100)
Z(i,j)=c(i,j*10)/scale_ampl;
end;
end;

%morl
scale_period=2*pi/5;
%mexh
%scale_period=pi*sqrt(2);




surface(B*scale_period,A,abs(Z))