% this file creates a matrix with each sound vector in a row. it then mdwt each row, store it in another matrix, finally return a vector with the average of the coeff's.

% run adj_three.m before running this prog

temp1(1,:) = O3_1';
temp1(2,:) = O3_2';
temp1(3,:) = O3_3';
temp1(4,:) = M3_1';
temp1(5,:) = M3_2';
temp1(6,:) = M3_3';
temp1(7,:) = N3_1';
temp1(8,:) = N3_2';
temp1(9,:) = N3_3';
temp1(10,:) = J3_1';
temp1(11,:) = J3_2';
temp1(12,:) = J3_3';
temp1(13,:) = B3_1';
temp1(14,:) = B3_2';
temp1(15,:) = B3_3';
temp1(16,:) = E3_1';
temp1(17,:) = E3_2';
temp1(18,:) = E3_3';
temp1(19,:) = S3_1';
temp1(20,:) = S3_2';
temp1(21,:) = S3_3';
% insert more here

h=daubcqf(le,'min');

% Delete the sound input files to save space
clear O3* M3* N3* J3* B3* E3* S3*;

for(i = 1: size(temp1,1))
	[temp1(i,:), L] = mdwt(temp1(i,:), h, L); 
	%temp1(i,1:3) = zeros(1,1:3);
        %temp1(i,:) = temp1(i,:)/max(temp1(i,:));
end

THREE = mean(temp1);

% if you want to see the template

%figure
%plot(THREE)
%title('THREE')