fwt_5scales.m

vector to matrix function

%This function performs the fast wavelet transform in 5 scales.

function [W1 W2 W3 W4 W5] = fwt_5scales (S)

[Lo_d,Hi_d]=wfilters ('db1','d');  %decomposition filters for the fwt algorithm ('gaus'?)

%a=1:1:128;

%S=zeros(128,19);
%S(:,1)=1:1:128;
%S(:,19)=1:1:128;

%perform fast wavelet transform

[N,M]=size(S); %dimensions of the features matrix

for m=1:1:M
[cA1,cD1]=dwt (S(:,m),Lo_d,Hi_d); %first scale
W1(:,m)=cA1';

[cA2,cD2]=dwt (cA1,Lo_d,Hi_d); %second scale
W2(:,m)=cA2';

[cA3,cD3]=dwt (cA2,Lo_d,Hi_d); %third scale
W3(:,m)=cA3';

[cA4,cD4]=dwt (cA3,Lo_d,Hi_d); %forth scale
W4(:,m)=cA4';

[cA5,cD5]=dwt (cA4,Lo_d,Hi_d); %fifth scale
W5(:,m)=cA5';
end