dc_dpcm.m

来自「用离散余弦变换实现JPEG编码」· M 代码 · 共 24 行

%dc系数的dpcm
% Differential Pulse Coded Modulation
function r = dc_dpcm(x)
%最简单的预测
a = [1];
% x: 输入列向量
% a: 预测滤波器（系数）列向量
% r: 取整后的差分编码输出
% before: 过去值序列
% now: 当前预测值
%获取信号长度和滤波器长度
m = length(x);  
p = length(a);
%与滤波器长度对应的前P个值直接输出
r = round( x(1:p) );
%预测序列初始值
before = r;
%预测滤波循环
for t=p+1:m,
   now(t)=a'*before(t-1:-1:t-p);
   r(t)=round(x(t)-now(t));
   before(t)=now(t)+r(t);
end