# How to implement the Daubechies wavelet?

http://en.wikipedia.org/wiki/Daubechies_wavelet#Transform.2C_D4

I find it is difficult to understand the pseudo-code on this Wiki page.

N = length(S);
s1 = S(1:2:N-1) + sqrt(3)*S(2:2:N);
d1 = S(2:2:N) - sqrt(3)/4*s1 - (sqrt(3)-2)/4*[s1(N/2); s1(1:N/2-1)];
s2 = s1 - [d1(2:N/2); d1(1)];<<<<<<<<<<<<<<<<<<<<<<<
s = (sqrt(3)-1)/sqrt(2) * s2;
d = (sqrt(3)+1)/sqrt(2) * d1;


I do not understand what does the square-bracket means.

-
I think the best way to understand the code is to implement it in Matlab. Can you do that? –  Nana Jun 20 '11 at 6:41