标准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 endfunction a = trunc(b)if(b>0)    a = floor(b);elseif (b<0)    a = ceil(b);else    a = 0;endend