通信原理--13折线PCM量化编码解码

程序用matlab实现，全都运行过，完美实现PCM功能班级_______通信原理大作业--抽样量化编码译码的Matlab实现学院电子工程学院学号********学生姓名***授课老师***一、前言通信系统的信源有两大类：模拟信号和数字信号。例如：话筒输出的语音信号属于模拟信号；而文字、计算机数据属于数字信号。数字信号相比于模拟信号有抗干扰能力强、无噪声积累的优点。因此，若输入是模拟信号，则在数字通信系统的信源编码部分需对输入模拟信号进行数字化。数字化需要三个步骤：抽样、量化和编码。抽样是指用每隔一定时间的信号样值序列来代替原来在时间上连续的信号，也就是在时间上将模拟信号离散化。量化是用有限个幅度值近似原来连续变化的幅度值，把模拟信号的连续幅度变为有限数量的有一定间隔的离散值。编码则是按照一定的规律，把量化后的值用二进制数字表示，然后转换成二值或多值的数字信号流。三个步骤如下图所示：二、模拟信号的抽样（1）低通模拟信号的抽样原理抽样定理：设一个连续模拟信号m(t)中的最高频率Hf，则以间隔时间为T≤Hf2/1的周期性冲激脉冲对它抽样时，m(t)将被这些抽样值所完全确定。下面对这个定理进行证明。设有一个最高频率小于Hf的信号m(t)。将这个抽样信号抽样信号量化信号t011011011100100100100编码信号信号和周期性单位冲激脉冲(t)T相乘，乘积就是抽样信号，它是一系列间隔为T秒的强度不等的冲激脉冲。这些冲激脉冲的强度等于相应时刻上信号的抽样值。现用)()(kTmtmS表示此抽样信号序列。故有)()()(ttmtmTS令M(f)、)(f和)(fMS分别表示m(t)、(t)T和)(tms的频谱。计算可得：)(1)()(1)(snssnffMTnfffMTfM上式表明，抽样信号的频谱)(fMS是无数间隔频率为Snf的原信号频谱M(f)相叠加而成。m(t)、(t)T和)(tms的时域和频域波形如下图所示：信号m(t)的最高频率小于Hf，若频率间隔Hsff2，则)(fMS中包含的每个原信号频谱M(f)之间互不重叠。这样就能够从)(fMS中用一个低通滤波器分离出m(t)的频谱M(f)，也就是能从抽样信号中恢复原信号。（2）模拟信号抽样的Matlab实现及采样定理的验证（a）编程思想令输入的模拟信号m（t）=)140cos(2)400sin(35.2tt，由已学知识可知：m（t）的最高频率Hf=200Hz，由抽样定理知：当采样频率Hzfs400（采样周期小于1/400s）时，抽样信号可以完全确定原信号。由傅里叶变换知识得：模拟信号m（t）=)140cos(2)400sin(35.2tt的频谱)(fM=)70()200(5.1)70()200(5.1)(5.2fffff。分别令采样频率HzfHzfss400400＜和，绘制出两种情况下各自抽样信号的频谱，与原模拟信号的频谱进行比较。若只有当采样频率大于400Hz时，抽样信号的频谱才与原模拟信号相同，则可以验证采样定理的正确性。（b）Matlab输出波形图00.10.20.30.40.50.60.70.80.91-50510模拟信号时域波形-500-400-300-200-10001002003004005000100020003000模拟信号频域波形020406080100120-50510采样序列时域（不符合采样定理）-60-40-2002040600100200300采样序列频域（不符合采样定理）0100200300400500600700800-50510采样序列时域（符合采样定理）-400-300-200-10001002003004000500100015002000采样序列频域（符合采样定理）（c）结论由上图可知：当采样频率小于2Hf时，抽样信号的频谱发生了混叠。而采样频率大于2Hf时，抽样信号的频谱与原模拟信号频谱相同。（d）Matlab源代码%sampling.m对模拟信号采样产生序列信号，并验证采样定理t=0:0.001:1;%模拟信号时域和频域波形x=2.5+3.*sin(200*2*pi*t)+2.*cos(70*2*pi*t);figure(1)subplot(2,1,1)plot(t,x);title('模拟信号时域波形')fs=linspace(-1000/2,1000/2,length(t));xf=fftshift(fft(x));subplot(2,1,2)plot(fs,abs(xf))%可知，信号最大频率200Hztitle('模拟信号频域波形')T=0.009;%不满足采样定理，采样周期大于1/400n=0:1:1/T;xn=2.5+3.*sin(200*2*pi*n*T)+2.*cos(70*2*pi*n*T);figure(2)subplot(2,1,1)stem(n,xn)title('采样序列时域（不符合采样定理）')fns=linspace(-0.5/T,0.5/T,length(n));xnf=fftshift(fft(xn));subplot(2,1,2)plot(fns,abs(xnf))title('采样序列频域（不符合采样定理）')T=0.0013;%满足采样定理，采样周期小于1/400n=0:1:1/T;xn=2.5+3.*sin(200*2*pi*n*T)+2.*cos(70*2*pi*n*T);figure(3)subplot(2,1,1)stem(n,xn)title('采样序列时域（符合采样定理）')fns=linspace(-0.5/T,0.5/T,length(n));xnf=fftshift(fft(xn));subplot(2,1,2)plot(fns,abs(xnf))title('采样序列频域（符合采样定理）')三、抽样信号的非均匀量化编码（13折线法）（1）脉冲编码调制简介模拟信号抽样后变成时间离散的信号，经过量化后，此抽样信号才能成为数字信号。分析可知：最简单的均匀量化器对于小输入信号很不利。为了改善小信号时的信号量噪比，在实际应用中常采用非均匀量化。非均匀量化时，量化间隔随信号抽样值的不同而变化。信号抽样值小时，量化间隔v也小；信号抽样值大时，量化间隔v也变大。实际应用中，用13折线法近似A压缩律，来进行非均匀量化。图中横坐标x在0至1区间中分为不均匀的8段。1/2至1间的线段称为第8段；1/4至1/2间的线段称为第7段；1/8至1/4间的线段称为第6段；依此类推。图中纵坐标y则均匀地划分作8段。将与这8段相应的座标点(x,y)相连，就得到了一条折线。在语音通信中，通常采用8位的PCM编码就能够保证满意的通信质量。在13折线法中采用的折叠码有8位。第一位c1表示量化值的极性正负。后面的7位分为段落码和段内码两部分，用于表示量化值的绝对值。其中第2至4位(c2c3c4)是段落码，共计3位，可以表示8种斜率的段落；其他4位(c5--c8)为段内码，可以表示每一段落内的16种量化电平。段内码代表的16个量化电平是均匀划分的。所以，这7位码总共能表示72＝128种量化值。在下面的表中给出了段落码和段内码的编码规则。段落序号段落码c2c3c4段落范围（量化单位）81111024~20487110512~10246101256~5125100128~256401164~128301032~64200116~3210000~16在上述编码方法中，段内码是按量化间隔均匀编码的，但是因为各个段落的斜率不等，长度不等，故不同段落的量化间隔是不同的。其中第1和2段最短，斜率最大，其横坐标x的归一化动态范围只有1/128。再将其等分为16小段后，每一小段的动态范围只有(1/128)/(1/16)=1/2048。第8段最长，其横坐标x的动态范围为1/2。将其16等分后，每段长度为1/32。假若采用均匀量化而仍希望对于小电压保持有同样的动态范围1/2048，则需要用11位的码组才行。现在采用非均匀量化，只需要7位就够了。（2）13折线法编码的Matlab实现量化间隔段内码c5c6c7c815111114111014110112110011101110101091001810007011160110501014010030011200101000100000（a）编程思想上面“模拟信号的抽样”已经得到了时间离散、幅度连续的抽样信号。若将对抽样信号非均匀量化编码编出的8位码组用1c2c3c4c5c6c7c8c表示。（1）抽样信号是1×n的矩阵，例如采样周期T=0.0013;n=0:1:1/T;则有770个抽样值，13折现法编码后得到的则是770×8的矩阵，一行对应一个抽样值的8位PCM码组。每个抽样值的编码思想都是一样的，若要进行多个抽样值的编码，进行for循环即可。下面阐述单个抽样值的编码思路。（2）确定极性码1c：利用matlab的符号函数sign(x)：x0时，sign(x)=-1；x=0时，sign(x)=0；x0时，sign(x)=1。所以，若sign（x）大于等于0，1c=1；否则1c=0。例如：+1000，sign（1000）=1，所以1c=1。（3）确定段落码2c3c4c：PCM13折线编码的动态范围为-2048--2048，而上一步骤中已经求出了各抽样信号极性，于是只要对抽样信号的绝对值分析即可。故对抽样值依次进行取模、归一、乘以2048、取整的操作，可以将抽样值转化为0--2048之间的整数。根据段落码与段落范围的关系，使用if语句即可确定2c3c4c。例如：+1000，因为1000≥128，故2c=1；又1000≥512，故3c=1；又1000≤1024，故4c=0。对于其他取值情况，判断方法与此类似。（4）确定段内码5c6c7c8c：每一段落均被均匀地划分为16个量化间隔，不过不同段落的量化间隔是不同的。（3）中确定了段落编码，可以确定每段的起始值，再根据待编码值、所在段的起始值、所在段量化间隔的大小即可确定段内码。例如：+1000，2c3c4c=110，故1000处于第7段（二进制110转化为十进制得到6,6+1=7）sp=[0,16,32,64,128,256,512,1024];%每段起始值spmin=[1,1,2,4,8,16,32,64];%每段的最小量化间隔再得到段起始值sp（7）=512，段最小量化间隔spmin（7）=32。25.15325121000，向负无穷方向取整得到15，再将15十进制转化为二进制，得到1111，故5c6c7c8c=1111。令Matlab程序输入为S=1000，可得：S=1000抽样信号进行13折线编码后的码组为（每一行代表一个抽样值，共1个值）code=11101111（可知：程序运行结果与上述分析结果相同）（5）若码组矩阵初始化为全0阵，则只有当某个码为1时，才需要进行赋值。例如：-100，sign（-100）=-1，1c=0。而1c初始化值就为0，故可以不采取任何操作，保持初始值即可。（b）Matlab输出结果0100200300400500600700800-50510抽样得到的序列0100200300400500600700800-2000-1000010002000取模、归一、乘2048、取整后的序列抽样信号进行13折线编码后的码组为（每一行代表一个抽样值，共770个值）code=11110011111111101110100101001001111000001110110101000001010110101110111111111001111001001110001011111010111101110001111001000100111010101110000101011001110100111111101011110101……（c）源代码%pcmdepcmcomp.m对抽样信号进行量化编码译码，将译码输出值与原抽