mat2lpc.m

%mat2lpc源程序
function y=mat2lpc(x,f)
%mat2lpc使用一维无损预测编码压缩一个矩阵
%y=mat2lpc(x,f)使用一个无损预测编码来编码矩阵y
%x的线性预测是建立在f的系数中，如果f省略，则f=1(先前像素编码)
%接着计算预测误差和编码矩阵y的输出。

%也看lpc2mat

error(nargchk(1,2,nargin));%检查输入的参数
if nargin<2%如果忽略，设置默认滤波器
    f=1;
end

x=double(x);%确认为double型计算
[m,n]=size(x);%获取输入矩阵的维数
p=zeros(m,n);%开始线性预测为0
xs=x;zc=zeros(m,1);%准备输入转换和填充

for j=1:length(f)%每个滤波器系数的建立
    xs=[zc xs(:,1:end-1)];%平移和零填充
    p=p+f(j)*xs;%形成预测和
end