[signal,fs] = audioread('chakitham1.wav');
frame(i,:) = signal(temp + 1 : temp + f_size);
frame(n_f+1,:) = zeros(size(frame(1,:)));
save('framedchakitham1');
load('framedchakitham1','n_f','frame');
fr_wind(j,:) = frame(j,:)'.*hamming(length(frame(j,:)));
mel_low=1125*log(1+(f_low/700));
mel_high=1125*log(1+(f_high/700));
Mel = linspace(mel_low,mel_high,filt_num+2);
Freq = 700*((exp(Mel/1125))-1);
f(i) = floor((nfft+1)*Freq(i)/fs);
h(m-1,k) = formula(k,f,m)
function H=formula(k,f,m)
elseif k>=f(m-1)&&k<=f(m)
H = (k-f(m-1))/(f(m)-f(m-1));
elseif k<=f(m+1)&&k>=f(m)
H = (f(m+1)-k)/(f(m+1)-f(m));
load('periodogram','psdy');
load('mel_filter_bank','h','filt_num');
num(i,:) = psdy(j,:).*h(i,:);
energy_log(j,i) = log(sum(num(i,:)));
dct_energy(j,:) = dct(energy_log(j,:));
if (nargin < 1 || nargin > 2)
transpose = (rows (x) == 1);
if transpose, x = x (:); end
x = [ x ; zeros(n-nr,nc) ];
w = [ sqrt(1/4/n); sqrt(1/2/n)*exp((-1i*pi/2/n)*[1:n-1]') ] * ones (1, nc);
if ( realx && rem (n, 2) == 0 )
y = fft ([ x(1:2:n,:) ; x(n:-2:1,:) ]);
y = fft ([ x ; flipud(x) ]);
if (realx) y = real (y); end
if transpose, y = y.'; end