Dekompozycja falkowa - matlab

[kowa19956 0]

Dzień dobry,

Proszę o przerobienie poniższego programu aby liczył dekompozycje z 3 poziomów ( poprzez użycie pętli)

q = 4;
%DWT 

[ca, cd] = dwt(y,'db1');% Współczynnik aproksymacji
ca = q*round(ca/q); %: zaokraglij do integera
cd = q*round(cd/q); 
s = symbols(ca); pa = s(:,2)/sum(s(:,2)); %: prawdopodobienstwa
s = symbols(cd); pd = s(:,2)/sum(s(:,2)); %: prawdopodobienstwa
ent_dwt1 = (entropy(pa) + entropy(pd))/2

xrec = idwt(ca, cd,'db1');
n = 1:length(ca);

plot(n, ca, n, cd); legend('Aproksymacja', 'Detal');title('Falka Haara (Falka Daubechies n=1)')

e_dwt1 = mean((xrec-y).^2)/A*100  %blad : wynik w procentach
figure(2)
plot(t,y,'b', t,xrec,'r')
legend('Sygnał oryginalny','Rekonstrukcja sygnału')