基于PCM编码的研究与仿真—课程设计论文

1、 课程设计评分标准课程设计评分标准 评分项目得分 具有题目、摘要、目录、正文、参考文献（5 分） 报告书写 及格式 正文格式，图、表、参考文献引用等正确，排版美观（5 分） 报告中是否体现被仿真系统的原理以及原理框图（5 分） 基础原理 仿真目的，仿真方法，仿真结果的意义表述清楚（5 分） 做出信源，调制信号，解调信号波形（10 分） m 文件仿真 仿真参量丰富（如对频谱，信噪比，误码率等的分析） ，仿真波形直 观。 （10 分） 是否实现设计功能，各个模块的设计参数是否清晰（10 分） 框图直观，有对不同参数条件下的仿真对比及结论（10 分）simulink 仿真 仿真参量丰富（如对频谱，信

2、噪比，误码率等的分析） ，仿真波形直 观。 （10 分） 是否存在抄袭（10 分） 对所仿真系统原理的提问回答情况（10 分）答 辩 对仿真过程提问的回答情况（10 分） 总分 课程设计评分标准.- 1 - 基于 matlab/simulink 的 pcm 编码的研究与仿真.- 3 - 摘 要.- 3 - abstract .- 3 - 1 背景知识.- 4 - 1.1 pcm 原理及仿真.- 4 - 1.1.1 脉冲编码调制.- 4 - 1.1.2 pcm 编码原理.- 5 - 2 m 文件仿真.- 10 - 3 simulink 仿真.- 15 - 3.1 原始模拟信号电路图及仿真图.-

3、15 - 3.2 pcm 编码器电路设计.- 17 - 3.3 pcm 解码器电路设计.- 24 - 4 心得体会.- 28 - 参考文献.- 28 - 基于 matlab/simulink 的 pcm 编码的研究与仿真 摘 要 本设计结合 pcm 的抽样、量化、编码原理，利用 matlab 软件编程和绘图功能，完成 了对脉冲编码调制（pcm）系统的建模与仿真分析。即学习通过利用计算机建立通信系统 模型的基本方法和基本技能，学习会利用仿真的手段对通信系统的基本理论和基本算法 进行验证。学习现有流行的通信系统仿真软件的基本使用方法，利用 matlab 软件解决通 信中存在的问题。 关键词：脉冲编

4、码调制（pcm） 均匀与非均匀量化matlab 仿真 abstract combined with the sampling, quantization, coding theory of pcm, using matlab software programming and graphics, the completion of the pulse code modulation (pcm) modeling and simulation analysis system.the basic method is to establish the model of communication sy

5、stem of learning through the use of computer and basic skills, learning by means of simulation of communication system of basic theory and basic algorithm are verified.the basic method of using communication system simulation software to learn of the existing, to solve communication problems using m

6、atlab software. keywords: pulse code modulation (pcm) with uniform and non-uniform quantization in matlab simulation 1 背景知识 1.1 pcm 原理及仿真 脉冲编码调制就是把一个时间，取值连续的模拟信号变换成时间离散，取值离散的 数字信号后在信道中传输。脉冲编码调制就是对模拟信号先抽样，再对样值幅度量化， 编码的过程。 1.1.1 脉冲编码调制 脉冲编码调制(pulse code modulation，pcm)是概念上最简单、理论上最完善的编 码系统，是最早研制成功、使用最为

7、广泛的编码系统，但也是数据量最大的编码系统。 pcm 的实现主要包括三个步骤完成：抽样、量化、编码。分别完成时间上离散、幅度上离 散、及量化信号的二进制表示。根据 ccitt 的建议，为改善小信号量化性能，采用压扩 非均匀量化，有两种建议方式，分别为 a 律和 律方式，我国采用了 a 律方式，由于 a 律压缩实现复杂，常使用 13 折线法编码,下图为 pcm 系统的原理框图： 图中，输入的模拟信号 m(t)经抽样、量化、编码后变成了数字信号(pcm 信号)，经 信道传输到达接收端，由译码器恢复出抽样值序列，再由低通滤波器滤出模拟基带信号 m(t)。通常，将量化与编码的组合称为模/数变换器（a/

8、d 变换器)；而译码与低通滤波的 组合称为数/模变换器(d/a 变换器)。前者完成由模拟信号到数字信号的变换，后者则相 反，即完成数字信号到模拟信号的变换。 pcm 在通信系统中完成将语音信号数字化功能，它的实现主要包括三个步骤完成：抽 样、量化、编码。分别完成时间上离散、幅度上离散、及量化信号的二进制表示。根据 ccitt 的建议，为改善小信号量化性能，采用压扩非均匀量化，有两种建议方式，分别为 a 律和 律方式，我国采用了 a 律方式，由于 a 律压缩实现复杂，常使用 13 折线法编 码,采用非均匀量化 pcm 编码。 1.1.2 pcm 编码原理 1)抽样 所谓抽样，就是对模拟信号进行周

9、期性扫描，把时间上连续的信号变成时间上离散 的信号。该模拟信号经过抽样后还应当包含原信号中所有信息，也就是说能无失真的恢 复原模拟信号。它的抽样速率的下限是由抽样定理确定的。 在一个频带限制在（0，f h）内的时间连续信号 f（t） ，如果以 1/2 f h 的时间间隔 对它进行抽样，那么根据这些抽样值就能完全恢复原信号。或者说，如果一个连续信号 f（t）的频谱中最高频率不超过 f h，当抽样频率 f s2 f h 时，抽样后的信号就包含 原连续的全部信息。抽样定理在实际应用中应注意在抽样前后模拟信号进行滤波，把高 于二分之一抽样频率的频率滤掉。这是抽样中必不可少的步骤。 2) 量化 量化，就

10、是把经过抽样得到的瞬时值将其幅度离散，即用一组规定的电平，把瞬时 抽样值用最接近的电平值来表示。 从数学上来看，量化就是把一个连续幅度值的无限数集合映射成一个离散幅度值的 有限数集合。一个模拟信号经过抽样量化后，得到已量化的脉冲幅度调制信号，它仅为 有限个数值。 如下图所示，量化器输出l个量化值，k=1，2，3，l。常称为重建电平或量 k y k y 化电平。当量化器输入信号幅度x落在与之间时，量化器输出电平为。这个量化 k x 1k x k y 过程可以表达为： 1 ( ),1,2,3, kkk yq xq xxxykl 这里称为分层电平或判决阈值。通常称为量化间隔。 k x kkk xx

11、1 图1-1 模拟信号的量化分为均匀量化和非均匀量化。 均匀量化：用这种方法量化输入信号时，无论对大的输入信号还是小的输入信号一 律都采用相同的量化间隔。为了适应幅度大的输入信号，同时又要满足精度要求，就需 要增加样本的位数。但是，对话音信号来说，大信号出现的机会并不多，增加的样本位 数就没有充分利用。为了克服这个不足，就出现了非均匀量化的方法。 非均匀量化：非均匀量化是根据信号的不同区间来确定量化间隔的。对于信号取值 小的区间，其量化间隔也小；反之，量化间隔就大。它与均匀量化相比，有两个突出 v 的优点。首先，当输入量化器的信号具有非均匀分布的概率密度（实际中常常是这样） 时，非均匀量化器的

12、输出端可以得到较高的平均信号量化噪声功率比；其次，非均匀量 化时，量化噪声功率的均方根值基本上与信号抽样值成比例。因此量化噪声对大、小信 号的影响大致相同，即改善了小信号时的量化信噪比。 实际中，非均匀量化的实际方法通常是将抽样值通过压缩再进行均匀量化。通常使 用的压缩器中，大多采用对数式压缩。广泛采用的两种对数压缩律是压缩律和 a 压缩 律。美国采用压缩律，我国和欧洲各国均采用 a 压缩律，所谓 a 压缩律也就是压缩器 具有如下特性的压缩律： ， 。 a x a ax y 1 0 , ln1 1 1 , ln1 ln1 x aa ax y 由于 a 律压缩实现复杂，常使用 13 折线法编码,

13、 压扩特性图如下图所示： 模拟入 y x 量化器 量化值 图 1-2 a 律函数 13 折线压扩特性图 这样，它基本上保持了连续压扩特性曲线的优点，又便于用数字电路实现，本设计 中所用到的 pcm 编码正是采用这种压扩特性来进行编码的。 表 1-1 13 折线时的值与计算值的比较xx y 0 8 1 8 2 8 3 8 4 8 5 8 6 8 7 1 x0 128 1 6 . 60 1 6 . 30 1 4 . 15 1 79. 7 1 93 . 3 1 98 . 1 1 1 按折线 分段时的x 0 128 1 64 1 32 1 16 1 8 1 4 1 2 1 1 段落12345678 斜

14、率16168421 2 1 4 1 表 1 中第二行的 值是根据 时计算得到的，第三行的 值是 13 折线分段时的值。可 见，13 折线各段落的分界点与 曲线十分逼近，同时 按 2 的幂次分割有利于数字化。 3) 编码 所谓编码就是把量化后的信号变换成代码，其相反的过程称为译码。当然，这里的 编码和译码与差错控制编码和译码是完全不同的，前者是属于信源编码的范畴。 在现有的编码方法中，若按编码的速度来分，大致可分为两大类：低速编码和高速 编码。通信中一般都采用第二类。编码器的种类大体上可以归结为三类：逐次比较型、 折叠级联型、混合型。在逐次比较型编码方式中，无论采用几位码，一般均按极性码、 段落

15、码、段内码的顺序排列。下面结合 13 折线的量化来加以说明。 表 1-2 段落码 表 1-3 段内码 段落码量化级段内码 151111 8111 141110 131101 7110 1100 111011 6101 101010 91001 5100 81000 70111 4011 60110 50101 3010 40100 30011 2001 20010 10001 1000 00000 pcm 编译码器的实现可以借鉴单片 pcm 编码器集成芯片，如：tp3067a、cd22357 等。 单芯片工作时只需给出外围的时序电路即可实现，考虑到实现细节，仿真时将 pcm 编译 码器分为编码

16、器和译码器模块分别实现，在 13 折线法中，无论输入信号是正是负，均按 8 段折线（8 个段落）进行编码。若用 8 位折叠二进制码来表示输入信号的抽样量化值， 其中用第一位表示量化值的极性，其余七位（第二位至第八位）则表示抽样量化值的绝 对大小。 具体的做法是：用第二至第四位表示段落码，它的 8 种可能状态来分别代表 8 个段 落的起点电平。其它四位表示段内码，它的 16 种可能状态来分别代表每一段落的 16 个 均匀划分的量化级。这样处理的结果，8 个段落被划分成 27128 个量化级。段落码和 8 个段落之间的关系如表 2 所示；段内码与 16 个量化级之间的关系见表 3。 4）译码 pc

17、m 译码器是实现 pcm 编码的逆系统。其中各模块功能如下： d/a 转换器：用来实现与 a/d 转换相反的过程，实现数字量转化为模拟量，从而达到 译码最基本的要求，也就是最起码的步骤。 瞬时扩张器：实现与瞬时压缩器相反的功能，由于采用 a 律压缩，扩张也必须采用 a 律 瞬时扩张器。 低通滤波器：由于采样脉冲不可能是理想冲激函数会引入孔径失真，量化时也会带 来量化噪声，及信号再生时引入的定时抖动失真，需要对再生信号进行幅度及相位的补 偿，同时滤除高频分量，在这里使用与编码模块中相同的低通滤波器。 2 m 文件仿真 程序如下： clear; t = -0.1:0.001:0.1; %该参数用于

18、画原信号图形 f = sin(3*pi*90*t)+cos(3*pi*37*t);%原函数, 由 t 的取值可得 f 有 201 个值 subplot(3,1,1) %matlab 矩阵区域设置 plot(t, f); %画出采原函数序列图 title(原信号); xlabel(时间 t(s); t= 1/500; %抽样周期，500 是抽样频率，可以调整抽样频率 gs = -0.1:t:0.1; fg = sin(2*pi*60*gs)+cos(2*pi*25*gs); %对信号进行以 t 周期抽样 subplot(3,1,2) stem(gs, fg) %画图 title(采样信号); x

19、label(时间 t(s)； 2-1 clear all; close all; %建立原信号 t=0.002; %取时间间隔为 0.01 t=-0.1:t:0.1; %时域间隔 dt 为间隔从 0 到 10 画图 xt=sin(3*pi*90*t)+cos(3*pi*37*t);%xt 方程%采样：时间连续信号变为时间离散模拟信号 fs=800; %抽样 fs=2fc，每秒钟内的抽样点数目将等于或大于 2fc 个 sdt=1/fs; %频域采样间隔 0.002 t1=-0.1:sdt:0.1; %以 sdt 为间隔从-0.1 到 0.1 画图 st=sin(2*pi*60*t1)+cos(2

20、*pi*25*t1); % 离散的抽样函数 figure(1); subplot(3,1,1); plot(t1,st); title(原始信号); %画出原始的信号图,以好对比 grid on %画背景 subplot(3,1,2); stem(t1,st,.); %这里画出来的是抽样后的离散图 title(抽样信号); grid on %画背景%量化过程 n=length(st); %取 st 的长度为 n m=max(st); a=(st/m)*2048; %a1(极性码) a2a3a4（段落码）a5a6a7a8（段内电平码） code=zeros(i,8); %产生 i*8 的零矩阵

21、%极性码 a1 for i=1:n %if 循环语句 if a(i)=0 code(i,1)=1; %代表正值 else code(i,1)=0; %代表负值 end if abs(a(i)=0code(i,3)=0;code(i,4)=0;step=1;start=0; elseif 16=abs(a(i)code(i,3)=0;code(i,4)=1;step=1;start=16; elseif 32=abs(a(i)code(i,3)=1;code(i,4)=0;step=2;start=32; elseif 64=abs(a(i)code(i,3)=1;code(i,4)=1;ste

22、p=4;start=64; elseif 128=abs(a(i)code(i,3)=0;code(i,4)=0;step=8;start=128; elseif 256=abs(a(i)code(i,3)=0;code(i,4)=1;step=16;start=256; elseif 512=abs(a(i)code(i,3)=1;code(i,4)=0;step=32;start=512; elseif 1024=abs(a(i)code(i,3)=1;code(i,4)=1;step=64;start=1024; end b=floor(abs(a(i)-start)/step); %段

23、内码编码 floor 取整(四舍五入) t=dec2bin(b,4)-48; %dec2bin 定义将 b 变为 4 位 2 进制码，-48 改变格式 code(i,5:8)=t(1:4); %输出段内码 end code=reshape(code,1,8*n); %reshape 代表从新塑形 code subplot(3,1,3); stem(code,.);axis(1 64 0 1); %这里我们先取前面八个点编码输出,输出时候有 64 个点 title(编码信号); grid on 2-2 3 3 simulinksimulink 仿真仿真 3.1 原始模拟信号电路图及仿真图 3-1

24、 原始模拟信号电路图 两个正弦波的参数设置分别为： 3-2 正弦波参数设置 3-3 正弦波参数设置 所得波形为 3-4 3.2 pcm 编码器电路设计 3-4 13 折线近似的 pcm 编码器测试模型和仿真结果 测试模型和仿真结果如上图所示。其中以 saturation 作为限幅器，将输入信号幅度值限 制在 pcm 编码的定义范围内，以 a-law compressor 作压缩器，relay 模块的门限值设置为 0， 其输出即可作为 pcm 编码输出的最高位极性码。样值取值绝对值后，用增益模块将样值放 大到 0-127，然后用间隔为 1 的 quantizer 进行四舍五入取整，最后将整数编码

25、为 7 位二进制 序列，作为 pcm 编码的低 7 位。可以将上图中 constant 和 display（不含）之间的模块封装 一个 pcm 编码子系统备用。 其中各模块的具体参数设置如下： 3-5 constant 3-6 saturation 3-7 abs 3-8 a-low compressor 3-9 gain 3-10 quantizer 3-11 integer to bit converter 3-12 display 3-13 relay 将该系统进行封装 3-14 封装之后的 pcm 编码子系统 3-15 封装之后的 pcm 编码子系统图标 3.3 pcm 解码器电路设计 3-16 13 折线近似的 pcm 解码