function x = mat2lpc( x,f )
error(nargchk(1,2,nargin));
if nargin<2
f=1;
end
x=double(x);
[m,n]=size(x);
xs=x; zc=zeros(m,1);
p=zeros(m,n);
for j=1:length(f)
xs=[zc xs(:,1:end-1)];
p=p+f(j)*xs;
end
y=x-round(p);
function x = mat2lpc( x,f )
error(nargchk(1,2,nargin));
if nargin<2
f=1;
end
x=double(x);
[m,n]=size(x);
xs=x; zc=zeros(m,1);
p=zeros(m,n);
for j=1:length(f)
xs=[zc xs(:,1:end-1)];
p=p+f(j)*xs;
end
y=x-round(p);