lpc2mat.m

%lpc2mat源程序
function x=lpc2mat(y,f)
%lpc2mat解压一维无损预测编码矩阵
% x=lpc2mat(y,f)在f中的线性预测系数下和一维无损预测编码下，解码输入矩阵y
%f忽略，则f=1(先前像素编码)

%也看mat2lpc

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

f=f(end:-1:1);%逆滤波器系数
[m,n]=size(y);%得到输出矩阵的维数
order=length(f);%得到线性预测的顺序
f=repmat(f,m,1);%赋值滤波器向量
x=zeros(m,n+order);%填充一列解码

%每次解码输出x的一列
%计算建立在先前顺序的元素和加它去预测误差
%结构是附上输出矩阵的

for j=1:n;
    jj=j+order;
    x(:,jj)=y(:,j)+round(sum(f(:,order:-1:1)).*...
        x(:,(jj-1):-1:(jj-order)),));
end

x=x(:,order+1:end);%删除左侧填充