模拟信号数字化PCM编码设计宿迁学院孩纸专用

1、仅供个人参考课程设计报告For personal use only in study and research; not for commercial use课程设计名称：通信系统原理系 部：三 系学生姓名：班 级：学 号：成 绩：指导教师：开课时间：学年_学期一.设计题目：模拟信号的数字化处理二.主要内容：a.熟悉模拟信号数字化的处理步骤：抽样、量化、编码；b.模拟信号的抽样过程，理解抽样频率的变化对抽样信号的影响；c.用MATLA喊其它EDAX具软件对PCMI码进彳T使用 A律和科律的压缩和扩张进行 软件仿真；d.PCM的 8 位编码 GGGGaaGGe.仿真实现增量调制的过程和并理解噪声

2、产生的原理。三.具体要求1 .课程设计的内容独立自主完成，课程设计报告内容完整、格式规范、排版整 洁美观；2 .设计选用的语言不限，推荐使用 MATLAB、C均可，编写出的程序，必须有 详细的注释说明；四.进度安排18周的时间安排：星期一 查阅资料，确定选题和软件，思考总体设计方案星期二熟悉软件的编程环境星期三总体设计方案的确定与设计星期四各部分的具体实现（程序调试并程序注释）星期五 整理完成设计报告的电子版，并答辩上交的材料：包含三样：（课程设计的电子稿、打印稿、程序）五.成绩评定考核方法：现场验收（占50%）,课程设计报告（占50%）。考核内容：学习态度（出勤情况，平时表现等）、方案合理性

3、、程序编制质 量、演示效果、设计报告质量。成绩评定：优，良，中，及格，不及格。特别说明：（1）如发现抄袭，按照不及格处理。（2）材料不齐的，考核等级降一级。电子稿件以文件夹的形式 上交的文件夹的命名为：（学号的最后两位+姓名）模拟信号数字化PCM编码设计1、 基本原理：本设计主要会用到的知识通信原理中的脉冲编码调制（ PCM）。本次课程设 计主要是做量化、编码、信道传输、译码四部分PCM系统的原理：脉码调制一将模拟调制信号的采样值变换为脉冲码组。PCM编码包括如三个过程。抽样：将模拟信号转换为时间离散的样本脉冲序列。量化：将离散时间连续幅度的抽样信号转换成为离散时间离散幅度的数字信 号。编码：

4、用一定位数的脉冲码组表示量化采样值。译码：把编码后的码流在通过信道传输后得到的码流译成离散时间连续幅度 的信号。其系统框图如下：图1-1 PCM系统原理图（一） 量化信噪比设压缩曲线则：噪声功率,1f (x) ln xB2."-2 v二 q = u 0f'(X) Px(X)dX6 022B2A v 2二60 X Px(X)dX不得用于商业用途仅供个人参考信号功率量化信噪比V 2S=20X px(x)dx3L2B2V2上式表明：压缩特性为对数特性时，量化器的信噪比始终保持常数，与输入信号的幅度无关，实现了最佳非均匀量化。国际上通用的两种对数压缩特性： A律和N律。本次设计采用的

5、是A律13 折线。（二）A律13折线的原理归一化信号（x/V）,过载电压为土 1, A律对数Ax1 ln A f (x)=,1 ln A其中A为压缩系数10 < X < A1x _ 1A,国际标准A=87.6。1 ln Ax图1-2为13折线特性表1-1列出了 13折线时的x值与计算x值的比较。表1-1y0182838485868781x01128160.6130.6115.417.7913.931 1.981按折线分段时的x011281641321161814121段落12345678斜率16168421124表1中第二行的x值是根据A=876时计算得到的，第三行的X值是13折线

6、 分段时的值。可见，13折线各段落的分界点与A=87.6曲线十分逼近，同时x按 2的幕次分割有利于数字化图中只有正幅度部分，共7折线，负幅度部分也有7折线。但正负部分第一 段折线斜率相等，为同一折线，共13折线。采用13折线近似后，加入正弦信号，通过测试，其量化信噪比与输入幅度的关图1-3 正弦输入SNR曲线（13折线近似）如图所示：采用折线近似方式后，加入正弦信号，信噪比的曲线会出现起 伏现象，这是因为在每段折线起始部分，量化间隔成倍增加，导致量化噪声增 加很快，而信号功率的增加却没有那么快，因而 SNR反而略有下降。但随信 号功率的增加，噪声功率基本保持不变，因此SNR有开始增加。这样共有

7、 6个起伏，7个峰值。（三）编码所谓编码就是把量化后的信号变换成代码， 其相反的过程称为译码。当然， 这里的编码和译码与差错控制编码和译码是完全不同的，前者是属于信源编码的范畴。在现有的编码方法中，若按编码的速度来分，大致可分为两大类：低 速编码和高速编码。通信中一般都采用第二类。编码器的种类大体上可以归结 为三类：逐次比较型、折叠级联型、混合型。在逐次比较型编码方式中，无论 采用几位码，一般均按极性码、段落码、段内码的顺序排列。下面结合13折线的量化来加以说明。在13折线法中，无论输入信号是正是负，均按 8段折 线（8个段落）进行编码。若用8位折叠二进制码来表示输入信号的抽样量化 值，其中用

8、第一位表示量化值的极性，其余七位（第二位至第八位）则表示抽 样量化值的绝对大小。具体的做法是：用第二至第四位表示段落码，它的8种 可能状态来分别代表8个段落的起点电平。其它四位表示段内码，它的 16种 可能状态来分别代表每一段落的16个均匀划分的量化级。这样处理的结果，8 个段落被划分成27= 128个量化级。段落码和8个段落之间的关系如表1-2所 示；段内码与16个量化级之间的关系见表1-3。表1-2段落码表1-3段内码段落序号段落码量化级段内码81111511111411107110131101121100610111101110101051009100181000401170111601

9、103010501014010020013001120010100010001000002、 仿真程序、程序编制、流程图、仿真结果:(1) 整个程序设计流程如图图1-4程序设计大致流程1)输入要进行编码的正弦信号，因为正弦信号的幅度是变化的所以定义了一个变量用于存放幅度，通过每次 的循环来改变幅度。由于幅度的的值最大为1 ,因此要对产生的正弦信号进行标准化。t = 0:0.01:10；m1 = -60:1:0; % 输入的正弦信号幅度不同m = 10A(m1/20); figure(1) for i = 1:length(m) for j = 1:2x1 = m(i)*sin(2*pi*t+2

10、*pi*rand(1); v = 1;x = x1/v;%normalize sxx = floor(x*4096);2)将采样后得到sxx进行A律PCM编码， y = pcm_encode(sxx);pcm_encode函数是按照附表A律PCM编码进行编程的。在话音信号的数字化 同际承准G. 711中，采用了折叠码编码。表 71是G. 711建议的语音信号的 4律编码规则。输入的语音信号经过抽样、量化后，每个抽样值编码成8个比特的二进制码组。量化时，A律中的每个区间义被均勺量化成 16个量化电平，其 编码规则为： 其中： b0为极性码,b0=0时对应输入为负。b0= 1时对应输入为正。 b1

11、b2b3为段落码，分别对应x的8个区间； b4b5b6b7为段内码，对应x区间中的l 6个量化电平值3)将经过PCM编码的信号送入信道题目要求的信道没有误码，所以将编码后的序列送入信道没有什么改变只是传输，没有产生误码。4)再将经过信道的信号y进行PCM解码y2 = pcm_decode(y,v);pcm_decode函数实际上是pcm_encode的逆过程。5)画出经过A律PCMg码、译码后得到的信号波形与未编码的信号波形。subplot(211)plot(t,x1);title('未编码信号波形');subplot(212)piot(t,y2)title('编译码后

12、得到的信号波形');6)画出不同幅度AH青况下，PCMt码后得到的信号量化信噪比曲线。对于一定Ac的情况下，量化信噪比的计算方法如下：根据公式SNR=： (6)CF覆行nq(j) = sum(x1-y2).*(x1-y2)/length(x1);sq(j) = mean(y2.A2);snr(j) = (sq(j)/nq(j);snrq(i)=10*log10(mean(snr);对于要绘制不同Ac情况下的量化信噪比曲线，只是进行循环，每次循环的Ac从而得到不同的量化信噪比。循环结束后，绘出量化信噪比曲线。plot(m1,snrq);axis(-60 0 0 60);title(

13、9;xulie 的量化信噪比') grid;下图为simulink搭建的PCM编码器框图：图1-5所示测试模型和仿真结果如图所示，其中以Saturation作为限幅器，将输入信号幅 度值限制在PCM®码的定义范围内，Relay模块的门限设置为0,其输出既可作 为PCM码出去的最高位一一极性码。阳值取绝对值后，以 Look-Up Table模块进行13折线压缩，并用增益模块将样值范围放大到0127,然后用间距为1的Quantizer进行四舍五入的取整，最后将整数编码为7位二进制序列，作为PCM编码的低7位。3、1、未编码的信号与经过PCM编码，PCM解码后的信号波形对比如下图结

14、论及其分析由图1-6可知，从图中看出图中两个信号波形几乎一样，说明未编码的 信号与把信号经过PCM编码、译码后得到的信号是一样的。与理论上经过 编译码后得到的波形与原信号的应该是相同的一致。图1-6是截取两个信号并将其放大得到的，从图1-6中可以看出两个波形大致一样，但是经编码译 码后得到的信号没有未编码信号的波形光滑。同时波形上也有一些细小的区别。2、PCM译码后的量化信噪比st =128将图1-8的量化信噪比与原理中的图3对比，发现两个图形大致一样，图 8 中的曲线也有6个起伏，7个波峰，说明程序得到的量化信噪比曲线是正确的。 因为题目中是对正弦信号进行 A律PCM编码、译码后得到的量化信

15、噪比曲线。 而在原理部分图3正是加入正弦信号，通过 A律13折线近似后得到的SNR曲 线。出现这样曲线是因为 在每段折线起始部分，量化间隔成倍增加，量化噪声功 率增加很快，而信号功率并没有很快增加，造成量化信噪比有所下降。但虽图 1-6与图1-3的波形基本吻合，却存在一定的不同。曲线的起伏和波峰没有理论 上的曲线的起伏明显。4、 心得体会一周的课程设计结束了，在这次的课程设计中不仅检验了我所学习的知识，也培养了我如何去把握一件事情，如何去做一件事情，又如何完成一件事情。学 会了合作，学会了运筹帷幄，学会了宽容，学会了理解，也学会了做人与处世。课程设计是我们专业课程知识综合应用的实践训练， 是我

16、们迈向社会，从事职业不得用于商业用途仅供个人参考工作前一个必不少的过程 ”千里之行始于足下”， 通过这次课程设计，我深深体会到这句千古名言的真正含义我今天认真的进行课程设计，学会脚踏实地迈开这一步，就是为明天能稳健地在社会大潮中奔跑打下坚实的基础在这次设计过程中，体现出自己单独设计的能力以及综合运用知识的能力，体会了学以致用、突出自己劳动成果的喜悦心情，从中发现自己平时学习的不足和薄弱环节，从而加以弥补。在此感谢我们的吴琼老师. ， 老师严谨细致、一丝不苟的作风一直是我工作、学习中的榜样；老师循循善诱的教导和不拘一格的思路给予我无尽的启迪；这次设计的每个实验细节和每个数据，都离不开老师您的细心

17、指导。而您开朗的个性和宽容的态度，帮助我能够很顺利的完成了这次课程设计。同时感谢对我帮助过的同学们，谢谢你们对我的帮助和支持，让我感受到同学的友谊。由于本人的设计能力有限，在设计过程中难免出现错误，恳请老师们多多指教 , 我十分乐意接受你们的批评与指正，本人将万分感谢。通过这次设计，掌握了 PCM®码的工彳原理及PCMR统的工作过程，学会了使用 matlab (通信系统的动态仿真软件)，并学会通过应用软件仿真来实现各种通信系统的设计，对以后的学习和工作都起到了一定的作用，加强了动手能力和学业技能。5、 参考资料1. 通信原理樊昌信曹丽娜国防工业出版社2. 现代通信系统分析与仿真 ma

18、tlab 通信工具箱李建新刘乃安西安电子科技大学出版社3. 数字信号处理教程 matlab 释义与实现陈怀琛电子工业出版社4. 现代通信系统一一使用 matlab约翰-G-普罗克斯西安交通大学出版社5. MATLAB 通信工程仿真张德丰机械工业出版社6. MATLAB/SIMULINK 通信系统建模与仿真实例分析邵玉斌 清华大学出版社附源程序代码：t = 0:0.01:10;m1 = -60:1:0; %输入的正弦信号幅度不同m = 10A(m1/20);figure(1)for i = 1:length(m)for j = 1:2x1 = m(i)*sin(2*pi*t+2*pi*rand(

19、1);v = 1;x = x1/v;%normalizesxx = floor(x*4096);y = pcm_encode(sxx);y2 = pcm_decode(y,v);nq(j) = sum(x1-y2).*(x1-y2)/length(x1);sq(j) = mean(y2.A2);snr(j) = (sq(j)/nq(j);drawnowsubplot(211)plot(t,x1);title('未编码信号波形');subplot(212)plot(t,y2)title('编译码后得到的信号波形');endsnrq(i)=10*log10(mean

20、(snr);endfigure(2)drawnowsubplot(211)plot(t,x1);axis(0 0.7 -1 1);title('未编码信号波形');subplot(212)plot(t,y2)axis(0 0.7 -1 1);title('编译码后得到的信号波形');figure(3)plot(m1,snrq);axis(-60 0 0 60);title('xulie 的量化信噪比')grid;PCM 编码函数程序：functionout=pcm_encode(x) %x encode to pcm code n=length(

21、x);%-4096<x<4096for i = 1:nif x(i)>0out(i,1)=1;elseout(i,1)=0;endif abs(x(i)>=0 & abs(x(i)<32out(i,2)=0;out(i,3)=0;out(i,4)=0;step=2;st=0;elseif 32<=abs(x(i)&abs(x(i)<64out(i,2)=0;out(i,3)=0;out(i,4)=1;step=2;st=32;elseif 64<=abs(x(i)&abs(x(i)<128out(i,2)=0;out

22、(i,3)=1;out(i,4)=0;step=4;st=64;elseif 128<=abs(x(i)&abs(x(i)<256out(i,2)=0;out(i,3)=1;out(i,4)=1;step=8;st=128elseif 256<=abs(x(i)&abs(x(i)<512out(i,2)=1;out(i,3)=0;out(i,4)=0;step=16;st=256;elseif 512<=abs(x(i)&abs(x(i)<1024out(i,2)=1;out(i,3)=0;out(i,4)=1;step=32;st=

23、512;elseif 1024<=abs(x(i)&abs(x(i)<2048out(i,2)=1;out(i,3)=1;out(i,4)=0;step=64;st=1024;elseif 2048<=abs(x(i)&abs(x(i)<4096out(i,2)=1;out(i,3)=1;out(i,4)=1;step=128;st=2048; elseout(i,2)=1;out(i,3)=1;out(i,4)=1;step=128;st=2048; endif(abs(x(i)>=4096)out(i,2:8)=1 1 1 1 1 1 1;elsetmp=floor(abs(x(i)-st)/step);t=dec2bin(tmp,4)-48;%函数 dec2bin输出的是 ASCII 字符串，48 对应 0 out(i,5:8)=t(1:4);end end out=reshape(out',1,8*n);PCM 解码函数程序：functionout= pcm_decode(in,v)%decode the input pcm code%in : input the pcm code 8 bits sample%v:quantized level n=length(in);in=reshape(in',