How to extract delta band (2-4 hz), theta band (4-8 hz), alpha band (8-12 hz), beta band (12-30 hz) and gama band (30-43 hz) from a single channel of EEG by wavelet decomposition?

hello friends.................

i m new in EEG analysis, i have written the code for wavelet decomposition of a single channel of eeg which is eeg1-

[C,L]=wavedec(s,6,'db4'); %Db4 wavelet multiscale decomposition of the signal (6 layers)

D6=detcoef(C,L,6);%Scale 64

D5=detcoef(C,L,5);%Scale 32

D4=detcoef(C,L,4);%Scale 16

D3=detcoef(C,L,3);%Scale 8

D2=detcoef(C,L,2);%Scale 4

D1=detcoef(C,L,1);

[C,L]=wavedec(s,6,'db4');%Scale 64

C6=appcoef(C,L,'db4',6);

[C,L]=wavedec(s,5,'db4');%Scale 32

C5=appcoef(C,L,'db4',5);

[C,L]=wavedec(s,4,'db4');%Scale 16

C4=appcoef(C,L,'db4',4);

[C,L]=wavedec(s,3,'db4');%Scale 8

C3=appcoef(C,L,'db4',3);

[C,L]=wavedec(s,2,'db4');%Scale 4

C2=appcoef(C,L,'db4',2);

[C,L]=wavedec(s,1,'db4');%Scale 2

C1=appcoef(C,L,'db4',1);

SRC6=wrcoef('a',C,L,'db4',6); % with the decomposition coeff reconstruction

SRC5=wrcoef('a',C,L,'db4',5);

SRC4=wrcoef('a',C,L,'db4',4);

SRC3=wrcoef('a',C,L,'db4',3);

SRC2=wrcoef('a',C,L,'db4',2);

SRC1=wrcoef('a',C,L,'db4',1);

SRD6=wrcoef('d',C,L,'db4',6);

SRD5=wrcoef('d',C,L,'db4',5);

SRD4=wrcoef('d',C,L,'db4',4);

SRD3=wrcoef('d',C,L,'db4',3);

SRD2=wrcoef('d',C,L,'db4',2);

SRD1=wrcoef('d',C,L,'db4',1);

am i correct or not?

[Wayne King]

What is the sampling frequency for your data?