标准DPCM编码（Matlab实现）

　　以下代码中，f代表源码，f1代表预测码，e代表预测码与源码差值，en代表变换后的差分值（通过含e的公式计算得到，因为有特定对应值，也可以通过查表实现）。f2是通过变换后的差分值得到的预测码。

　　最终输出的结果是en，也就是变换后的差分值。

function en = my_dpcm(f)
[~,n] = size(f);
f1 = zeros(1,n);
f2 = zeros(1,n);
e = zeros(1,n);
en = zeros(1,n);

f1(1) = f(1);
f1(2) = f(1);
e(1) = 0;
en(1) = 0;
f2(1) = f(1);

for i = 2:n
if(i~=2)
f1(i) = ( f2(i-1)+f2(i-2) )/2;
end
e(i) = f(i) - f1(i);
en(i) = 16*trunc((255 + e(i))/16)-256+8;
f2(i) = en(i) + f1(i);
end

end

function a = trunc(b)
if(b>0)
a = floor(b);
elseif (b<0)
a = ceil(b);
else
a = 0;
end
end