标准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