数字通信原_3：信源编码.pptVIP

第三章 信源编码 (1)极点预测器（续前） xe(n) = ?i=1N aixr(n-i), {ai}为预测系数； 预测：利用过去值来估计（当前）未来值。 dq(n)≈d(n) = x(n) - xe(n) = x(n)- ?i=1N aixr(n-i) ≈ xr(n)- ?i=1N aixr(n-i) 等式两边取 Z 变换： dq(Z) ≈ （1- ?i=1N aiZ－i）Xr(Z) 若定义：H(Z) = Xr(Z)/dq(Z) H(Z) = 1/（1- ?i=1N aiZ－i） H(Z)只有极点－－－“极点预测器”。 第三章 信源编码 (2) 零点预测器 若取： xe(n) = ?i=1N bidq(n-i), {bi}为预测系数； xr(n) = dq(n) ＋ xe(n) = dq(n) ＋ ?i=1N bidq(n-i) 等式两边取 Z 变换： Xr(Z)＝ （1+ ?i=1N biZ－i）dq(Z) 若定义：H(Z) = Xr(Z)/dq(Z) H(Z) = （1 + ?i=1N biZ－i） H(Z)只有零点－－－“零点预测器”。 第三章 信源编码 (2) 零点预测器（续前） 基于零点预测器的DPCM编码解码系统 － 量化器 （1+?i=1N biZ－i） 编码器 x(n) xe(n) d(n) dq(n) I(n) 解码 ＋ I(n) dq(n) xe(n) xr(n) （1+?i=1N biZ－i） 预测器 第三章 信源编码 (3) 零极点预测器 若取： xe(n)= ?i=1N aixr(n-i)＋?j=1M bjdq(n-j), {ai,bj}为预测系数； 由：xr(n) = xe(n) + dq(n) - xe(n) = xr(n) - dq(n) 即：xr(n)-dq(n)= xe(n) = ?i=1N aixr(n-i)＋?j=1M bjdq(n-j) 又由：H(Z) = Xr(Z)/dq(Z) = [1+?j=1M bjZ-j]/[1-?j=1N aiZ-i] H(Z)包含零点和极点－－－“零极点预测器”。 第三章 信源编码 (3) 零极点预测器（续前） 编码器 － ＋ 量化器 零点 预测器 编码 x(n) xe(n) d(n) dq(n) xr(n) I(n) ＋ 极点 预测器 ?j=0N bjZ－j 零点 预测器 ＝ 极点 预测器 ?i=0N aiZ－i ＝ 第三章 信源编码 (3) 零极点预测器（续前） 解码器 ＋ 零点 预测器 解码 xe(n) dq(n) xr(n) I(n) ＋ 极点 预测器 第三章 信源编码 极点预测器系数{ai}的确定 (1) E{d2}最小(最佳预测)条件下极点预测器系数的求解 E[d2(n)] = E{[x(n)-xe(n)]2} = E{[x(n)-?i=1N aixr(n-i)]2} ≈ E{[x(n)-?i=1N aix(n-i)]2} 令 ? E[d2]/?am= -2E{[x(n)-?i=1N aix(n-i)]x(n-m)}=0 (*) m = 1,2,3,…,N 设x(n)广义平稳的随机序列，则相关函数R(n,n-i)满足 R(n,n-i) = E[x(n)x(n-i)] = R(i) (*)式变为： R(1)=a1R(0)+a2R(1)+…+aNR(N-1) R(2)=a1R(1)+a2R(0)+…+aNR(N-2) … …… R(N)=a1R(N-1)+a2R(N-2)+…+aNR(0) (*1) 第三章 信源编码 （1） E{d2}最小条件下极点预测器系数的求解(续前) (*1)式的 R(1) R(0) R(1) … R(N-1) a1 矩阵形式： R(2) = R(1) R(0) … R(N-2) a2 … …… … R(N) R(N-1)R(N-2)… R(0) aN 解得（假定[R(i)]为非奇异矩阵）：