Page 185 - Algorithm Collections for Digital Signal Processing Applications using MATLAB

174 Chapter 4

assigncenno.m

function [res]=assigncenno(x)
x=double(x);
load C1
load C2
load I
x=reshape(x,1,64);
p=lpc(x,16);;
temp=p(2:1:17);
a=[ ];
for n=1:1:size(I,1)
a= [a sum(temp.*I(n,:)) ];
end
d1=sum((C1-a).^2);
d2=sum((C2-a).^2);
if(d1<d2)
res=0;
else
res=1;
end
___________________________________________________________

8.3 Program Illustration