自适应差分脉码调制

1、 第第4章章 模拟信号的数字传输模拟信号的数字传输4.2自适应差分脉码调制自适应差分脉码调制（选修）（选修）4.2.1 概述概述n几十年来，人们一直致力于压缩数字化语音几十年来，人们一直致力于压缩数字化语音占用频带的工作，也就是在相同质量指标下，占用频带的工作，也就是在相同质量指标下，努力降低数字化语音比特率，以提高数字通努力降低数字化语音比特率，以提高数字通信系统的频带利用率信系统的频带利用率n语音压缩方法有语音压缩方法有dpcm sbc atc celp mplpc vq adpcmnadpcm在更低的码率下在更低的码率下(16-8kbps)，质量，质量明显下降明显下降 ( )( )( )

2、( )( )( )( )( )( )( )( )ereqrqd ks ks ks ks kd ke ks ks kd kd k22( )( )qsnre dke e k4.2.2 dpcm基本原理基本原理4.2.3 adpcm基本原理基本原理n极点预测器极点预测器n预测误差滤波器预测误差滤波器n重建滤波器重建滤波器1( )()neiris ka s ki1( )1niiid za z 11( ) 1( )1niiih zd zaz4.2.4 语音压缩编码语音压缩编码n原理框图原理框图 a/d 信号 变换 压缩编码 输出 4.2.4.1 语音编码技术概述语音编码技术概述n语音编码：语音编码：移动

3、通信数字化的基础移动通信数字化的基础 第第1/2代蜂窝系统根本区别代蜂窝系统根本区别n语音编码的意义：语音编码的意义：u提高通话质量 (数字化+信道编码纠错)u提高频谱利用率 (低码率编码)u提高系统容量 (低码率，语音激活技术)4.2.4.2 语音编码分类语音编码分类语音编码器波形编码器参量编码器频域时域非差分子带编码自适应变换域编码差分pcmdpcmm连续可变斜率madpcmapc线性预测编码声码器信道声码器共振峰声码器倒频谱声码器语音激励声码器多脉冲激励lpc码本激励lpc矢量和激励lpc4.2.4.3 语音编码分类语音编码分类n波形编码：波形编码：将时域模拟话音的波形信号进过采样、量化

4、和编将时域模拟话音的波形信号进过采样、量化和编码形成数字语音信号码形成数字语音信号-现代通信原理现代通信原理u编码速率较高，16k64ku包括：pcm，adpcm，m，cvsdm，apc等u占用较高带宽，适合有线n参量编码：参量编码：基于人类语音的产生机理建立数学模型，根据输基于人类语音的产生机理建立数学模型，根据输入语音得出模型参数并传输，在收端恢复。入语音得出模型参数并传输，在收端恢复。u编码速率较低，1.24.8 kbpsu包括各种线形预测编码(lpc)方法和余弦声码器u语音质量中等，不满足商用要求n混合编码：混合编码：波形编码波形编码+参量编码参量编码 (lpas)u包括gsm的rpe

5、-lpc编码和vselp编码4.2.4.4 语音编码的标准语音编码的标准ng.711 pcm (64k bps)ng.721 adpcm (32k bps)ng.722 7khz带宽带宽64k bps速率内的音频编码速率内的音频编码ng.723 6.3k/5.6k 双速率多媒体语音编码双速率多媒体语音编码ng.728 16k bps 语音编码语音编码 ld-celpng.729 8k bps多媒体语音编码多媒体语音编码4.2.4.5 语音编码的发展语音编码的发展n极低速率语音编码，极低速率语音编码，600bpsn高保真语音编码高保真语音编码n自适应多速率语音编码自适应多速率语音编码n新的编码分

6、析技术：新的编码分析技术：u非线形预测u多精度时频分析技术(子波分析技术)u高阶统计分析技术4.2.4.6 线形预测编码的原理线形预测编码的原理n原理：模型化人类语音信号产生的机制，提取原理：模型化人类语音信号产生的机制，提取模型参数，并且只传输模型的参数。模型参数，并且只传输模型的参数。n语音信号的产生模型：语音信号的产生模型：u语音的产生，声带和声道u不同语音产生的原因：声音激励源和声道不同u声音分类：浊音和清音n发声过程发声过程u口腔和鼻腔形成时变滤波器周期基音周期浊音声带振动语音波形三角形周期脉冲波，含丰富谐波清音声带不振动类似白噪声4.2.4.7 清音和浊音清音和浊音 清浊音频谱清浊

7、音频谱 清浊音信号清浊音信号4.2.4.8 语音产生模型语音产生模型n语音模型的建立：语音模型的建立：1. 产生激励，产生激励，2.响应响应n模型参数：基音，共振峰频率及强度，清浊音模型参数：基音，共振峰频率及强度，清浊音判决判决u低码率编码v码率降低有限码率降低有限v合成的语音波形失去了自然度和音质合成的语音波形失去了自然度和音质周期脉冲发生器随机噪声发生器g时变数字滤波器)(nu清/浊音开关声道参数)(ns基音周期4.2.4.9 lpc语音编码语音编码n使误差均方最小，可求得一组预测系数使误差均方最小，可求得一组预测系数ak，n传送：预测系数，基音周期合增益，清浊音判决传送：预测系数，基音

8、周期合增益，清浊音判决111()()()1( )()( )( )()( )( )( )nkkknkknkks zgh zu za zs na s nkg u ns na s nke ns ns n时变滤波器：预测器：误差：4.2.5 图像编码图像编码n变换：变换：n编码编码.261dftdcthjpegwthtslant静止二维预测变换利用象素间的相关特性活动三维离散傅立叶变换变换离散余弦变换最广，函数变换利用变换域沃尔什变换参数分布特征哈达码变换倾斜变换帧内编码局部更新编码帧间编码 帧间预测运动补偿4.2.5.1 图像压缩标准图像压缩标准n可压缩的原因可压缩的原因n已有的标准已有的标准统计冗

9、余度低频分量多视觉冗余度人眼的感觉冗余度空间冗余度 图像结构时间冗余度 图像活动性12,4mpegiso mpegdvd dvbmpeginternet.261.263.264hituhhn为什么要变换和量化为什么要变换和量化变换变换量化量化熵编码熵编码具有优良频率特性的滤波器作变换工具，有利于得到平稳的能量分布，有利于进一步的量化处理原始信号的信息损失主要发生在量化阶段，好的量化能以尽可能少的量化误差换取最小的码率+=高质量编码器4.2.5.2 编码与变换、量化编码与变换、量化4.2.6 dft与与dct120( ) ( )( ),01njnknnx kdft x nx n ekn1201(

10、 )( )( ),01njnknkx nidft x kx k ennnn dft变换变换 结果有虚部和实部结果有虚部和实部n 对对dft变换系数量化后，吉布斯效应明显变换系数量化后，吉布斯效应明显n 把把x(n)对称延拓，可以消除其中的虚部，得对称延拓，可以消除其中的虚部，得 到到dct变换变换4.2.6.1 正交变换的一个例子正交变换的一个例子假设图象的大小只有一个假设图象的大小只有一个44的块，即一共的块，即一共16个点，个点，其中每个象素点都有一定的值，我们把这个值叫灰其中每个象素点都有一定的值，我们把这个值叫灰度值。对图象而言，这些值为度值。对图象而言，这些值为0255之间的整数。之

11、间的整数。假设一幅假设一幅44的图象，它的灰度值如下：的图象，它的灰度值如下： 1281658182345214200168486424362204下面我们对这个下面我们对这个44的图象做变换，变换矩阵我的图象做变换，变换矩阵我们选：们选：111111111111121111t一般我们希望它是正交矩阵，即变换矩阵和它的一般我们希望它是正交矩阵，即变换矩阵和它的逆矩阵的乘积，为单位矩阵，即：逆矩阵的乘积，为单位矩阵，即：1ttii1000010000100001i对上述图像作变换，可以验证：对上述图像作变换，可以验证： t1111111111111211111281658182345214200

12、168486424362204370243318290781652181541661391322662286125212146我们希望这个过程是可逆的，即可以验证：我们希望这个过程是可逆的，即可以验证： = 11111111111112111137024331829078165218154166139132266228612521214612816581823452142001684864243622044.2.6.2 二维二维44 dct变换变换aaaaabacbccbacabaaaaacabcbbcabacty= axax1113,cos ,cos22828abcty= cxcex2222

13、2222aac aac11111 d11d11d111dac cac c1111111 daac aac1dd11d11ac cac c若若 表示矩阵对应元素相乘，表示矩阵对应元素相乘，d=b/c，那么上式可，那么上式可以表示为以表示为 时，为最优变换矩阵时，为最优变换矩阵1/( 2 1 ) 2.4142.d 在这里，在这里，d取取52，非常接近最优值，非常接近最优值 此时，此时， ，得到整数变换矩阵的正变换得到整数变换矩阵的正变换为：为： 1/( 2 1 ) 2.4142.d 1522,22 2929abc4.2.6.3 二维二维44 dct整数整数变换变换tfffy= c xcex2222

14、222211111511aac/2aac/251151115ac/2 c /4 ac/2 c /411111115aac/2aac/215511511ac/2 c /4 ac/2 c /44.2.6.4 二维二维44 dct变换与量化变换与量化对于对于9比特的预测误差数据做变换，变换矩阵每一行比特的预测误差数据做变换，变换矩阵每一行绝对值之和的最大值为绝对值之和的最大值为14，因此对于变换后的矩阵，因此对于变换后的矩阵元素需要元素需要 位来表示。现在大多数处理装置位来表示。现在大多数处理装置（如（如cpu，dsp）为）为32位宽，用上述变换进行处理位宽，用上述变换进行处理就没有数据溢出的危险。

15、就没有数据溢出的危险。 29+log14 14=17一旦得到变换后的系数，我们就可以对变换系数进一旦得到变换后的系数，我们就可以对变换系数进行量化操作。对于某个量化步长行量化操作。对于某个量化步长qstep，对应一个量，对应一个量化参数化参数q，量化参数，量化参数q的取值范围为的取值范围为0到到44整数，量整数，量化就是把变换得到的结果都除以化就是把变换得到的结果都除以qstep： tfffy/qstep= c xcex/qstep1/41/2321/41/23211111511511511151/2321/581/2321/58/qstep111111151/41/2321/41/23215

16、5115111/2321/581/2321/58tfff= c xcmqbitsq/(2)y2fqbitsfemqstep可以得到：可以得到：其中，其中， ，floor表示向负无穷取整。量表示向负无穷取整。量化参数化参数q与量化步长与量化步长qstep的对应关系如表的对应关系如表1。从。从表表1中可以看出，量化参数中可以看出，量化参数q每增加每增加5时，量化步时，量化步长长qstep就随之翻倍增加。就随之翻倍增加。15( /5)qbitsfloor q可以得到可以得到q01234567891011qstep 0.625 0.6875 0.8125 0.875 11.251.375 1.625

17、1.75 2 2.5 2.75 q1419 24 29 34 39 qstep 48 163264128 表表1 量化参数量化参数q与量化步长与量化步长qstep的对应关系的对应关系 nmoret首先提出了小波的概念首先提出了小波的概念(moret(moret小波函数是小波函数是grossmanngrossmann和和morletmorlet在在19841984年提出的年提出的) )，在很多学者的共同努力下，小，在很多学者的共同努力下，小波理论得到了完善和发展。特别是波理论得到了完善和发展。特别是mallat，他提出的多分，他提出的多分辨率分析辨率分析(multiresolution anal

18、ysis, mra)理论，对小波构理论，对小波构造理论起到了十分重要的作用。造理论起到了十分重要的作用。 n1976年，年，crochiere等人首先把子带编码应用于语音编码。等人首先把子带编码应用于语音编码。在语音编码中，用作频谱在语音编码中，用作频谱qmf (quadrature mirror filter)滤波器能够很好地解决过度带频率混叠问题，而且滤波器能够很好地解决过度带频率混叠问题，而且很自然地，人们把很自然地，人们把qmf引入二维信号处理。引入二维信号处理。n子带就是小波分解后，得到不同的频带分量的时域或者空子带就是小波分解后，得到不同的频带分量的时域或者空域域表现。表现。n子带

19、编码的基本思想是，把信号分成多个子带，然后对各子带编码的基本思想是，把信号分成多个子带，然后对各个子带进行编码。个子带进行编码。4.2.7 小波变换与小波变换与qmf4.2.7.1 qmf滤波器滤波器n 混叠问题混叠问题(aliasing)n 从信号处理角度看：从信号处理角度看：qmf滤波器能够很好地解决过度带频率混叠问题滤波器能够很好地解决过度带频率混叠问题。11( )( ) ( )x zh z x z1/21/21111( )()()2y zx zxz)()(211zyzu1111111( )( )( )( ) ( )() ()( )2t zk z u zh z x zhz xz k z2

20、222221( )( )( )( ) ( )() ()( )2t zk z u zh z x zhz xz k z假设：11( )( )k zh z221( )( )()k zh zhz22121( )( )( )( )2s zhzhzx z得到：同理：211( )()uzy z4.2.7.2 小波分解小波分解n 从函数正交分解角度看：小波分解从函数正交分解角度看：小波分解能够很好地解决信号的重构问题能够很好地解决信号的重构问题。n meyer于于1985-1986年度在年度在boubaki研讨会上提出了二进小波概念，完善了研讨会上提出了二进小波概念，完善了qmf滤波器的理论，为滤波器的理论，

21、为qmf把频谱以二分分解提供了理论依据。把频谱以二分分解提供了理论依据。( )(2)nkthtn1( )( 1)(2)nnkthtnn 对于具有紧支的标准正交函数系，任何连续可积函数的取样值为： 0( ), ()ncx ttnn 那么，x(t)可以写成： 0( )()nnx tctnn 由于 (t-n)对不同的n正交，对x(t)的小波分解算法可以写成： 11( ), (2)ncx tt n 0(),(2)kmkmctkhtnm 02n mmmch 1021( 1)mnn mmmdchn 可以证明， (t-n)对不同的n正交，等价于： 22( )()1hh mallat提出了多分辨率分析提出了多

22、分辨率分析(multiresolution analysis, mra)理论，希望得到图象理论，希望得到图象3 3个方向的纹个方向的纹理特征，以便机器识别理特征，以便机器识别4.2.7.3 小波和小波和qmf的比较的比较n 完备性概念完备性概念n 具有相同的信号完全重构的条件具有相同的信号完全重构的条件n 具有相同的信号分解和重构的算法具有相同的信号分解和重构的算法n 小波分解采用的是小波分解采用的是l2空间中完备的正交基空间中完备的正交基n 小波理论，能构造出相同抽头小波理论，能构造出相同抽头(taps)下，不同性能的滤波下，不同性能的滤波器器n 小波分析有更广泛的应用，除了编码，还能用在信

23、号检测、小波分析有更广泛的应用，除了编码，还能用在信号检测、微分方程、曲面拟合等领域。微分方程、曲面拟合等领域。4.2.7.4 正交小波的构造正交小波的构造22222 (0)(0)(1)(2) .( )2 (1)(0) (1) .(1) ( )2 (2)(0) (2) .(2) ( ).2 ( )(0) ( )ghhhh nghhh nh nghhh nh ng lhh n 根据根据mra的思想，的思想，daubechies设计了正交性能良设计了正交性能良好的高低通滤波器组，好的高低通滤波器组，但是它们不对称即不满足线但是它们不对称即不满足线性相位，因此不能完全重构图象的边界性相位，因此不能完

24、全重构图象的边界4.2.7.5 双正交小波的构造双正交小波的构造 1 2(0)(0) (0) 2 (1) (1)9162 (1)2 (0) (1) 2 (1) (0) 2 (1) (2)02 (2)2 (0) (2) 2 (1) (1)1162 (3)2 (1) (2)ghhhhghhhhhhghhhhghh3/5双正交小波双正交小波 -1/8 -1/8 3/5小波分解与重构中的对齐问题小波分解与重构中的对齐问题 x2 x1 x0 x1 x2 xn-4 xn-3 xn-2 xn-1 xn-2 ck -1/8 1/4 3/4 1/4 -1/8-1/8 1/4 3/4 1/4 -1/8 dk (-

25、1)n 1/4 1/2 1/4 1/4 1/2 1/4 分解取偶数点，重构时补分解取偶数点，重构时补0在奇数点，分解时在奇数点，补在奇数点，分解时在奇数点，补0在偶数点。然后在偶数点。然后以边界点为中心对称延拓，分解和重构时延拓方式是一样的。以边界点为中心对称延拓，分解和重构时延拓方式是一样的。4.2.7.6 矢量小波矢量小波(multiwavelets)som多重小波具有很好的低通、带通和多重小波具有很好的低通、带通和高通特性，具有提高编码质量的潜力高通特性，具有提高编码质量的潜力4.2.7.7 零树量化零树量化(ezw spiht)4.2.7.8 wavelet与与dct的比较的比较(y)

26、wavelet与dct的比较(u v)dct变换编码变换编码, 17.95 dbwavelet变换编码变换编码, 20.98 db小波变换可以有效的克服块效应小波变换可以有效的克服块效应(a) lenna原图原图 (b) 0.005bpp, psnr=20.76db (c) 0.01bpp, psnr=22.72db (d) 0.02bpp, psnr=24.21db (e) 0.05bpp, psnr=27.30db (f) 0.1bpp, psnr=28.81db 小波编码在高码率和低码率情况下都取得很很好的实验结果小波编码在高码率和低码率情况下都取得很很好的实验结果编码的尺寸编码的尺寸可

27、伸缩性可伸缩性编码的质量编码的质量可伸缩性可伸缩性4.2.7.9 小波变换编码具有可伸缩性小波变换编码具有可伸缩性小波编码的数据结构的量化手段十分有小波编码的数据结构的量化手段十分有利于可伸缩性编码和解码利于可伸缩性编码和解码 4.2.8 视频编码与视频编码与h.264随着数字技术的发展，欣赏数字影视成为大众娱随着数字技术的发展，欣赏数字影视成为大众娱乐生活的重要内容。数字影视以数字信号形式，乐生活的重要内容。数字影视以数字信号形式，或者说以或者说以bit的形式在各种介质（或者叫存储器）的形式在各种介质（或者叫存储器）存储和传输。我们希望在容量有限的存储器内存存储和传输。我们希望在容量有限的存

28、储器内存放更多的电影数量，或者希望在网络点播中看到放更多的电影数量，或者希望在网络点播中看到更流畅的电影作品。然而存储器的容量或者网络更流畅的电影作品。然而存储器的容量或者网络能提供的带宽总是不能满足人们不断增长的需要，能提供的带宽总是不能满足人们不断增长的需要，这就推动了视频压缩编码技术的发展。这就推动了视频压缩编码技术的发展。对于对于cif格式的视频流，它的画面大小是格式的视频流，它的画面大小是352288个象个象素，每个象素点用素，每个象素点用8个个bit来表示，对于黑白的画面来说，来表示，对于黑白的画面来说，一幅画面需要一幅画面需要3522888=811,008bit，如果是彩色画，如果是彩色画面，则需要面，则需要1,216,512。为了欣赏到连续的画面，根据人。为了欣赏到连续的画面，根据人类的视觉特性，每秒差不多要播放类的视觉特性，每秒差不多要播放25个这样的画面（或个这样的画面（或者叫帧），才不会有快镜头的感觉（比如象卓别林的老者叫帧），才不会有快镜头的感觉（比如象卓别林的老电影）。那么，对于一部电影）。那么，对于一部60分钟长的电影，它需要的存分钟长的电影，它需要的存储量为储量为1,216,512256060=109,486,080,000，约为，约为13,685,760,000个字节，即个字节，即13个个g的存储空间。市场上的的存储空间。市场上的