2022年增量调制仿真设计终版 .pdf

沈阳理工大学通信系统课程设计报告1 增量调制仿真设计1.课程设计目的（1）加深理解增量调制编译码的基本原理。（2）培养独立开展科研的能力和编程能力。（3）了解系统的过载特性,动态范围以及最大量化信噪比等三大指标的测试方法。2.课程设计要求（1）掌握课程设计的相关知识、概念清晰。（2）程序设计合理、能够正确运行。3.相关知识3.1 增量调制简介增量调制简称 M或增量脉码调制方式 （DM ），它是继 PCM 后出现的又一种模拟信号数字化的方法。1946年由法国工程师 De Loraine 提出，目的在于简化模拟信号的数字化方法。 主要在军事通信和卫星通信中广泛使用，有时也作为高速大规模集成电路中的A/D 转换器使用。对模拟信号采样，并用每个样值与它的预测值的差值对周期脉冲序列进行调制，简称墹M或 DM 。已调脉冲序列以脉冲的有、无来表征差值的正负号，也就是差值只编成一位二进制码。增量调制的基本原理是于1946 年提出的，它是一种最简单的差值脉冲编码。早期的语言增量调制编码器是由分立元件组成的。随着模拟集成电路技术的发展，70 年代末出现了音节压扩增量调制集成单片，80 年代出现了瞬时压扩集成单片，单片内包括了开关电容滤波器与开关电容积分器，集成度不断提高，使增量调制的编码器的体积减小，功耗降低。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 15 页 - - - - - - - - - 沈阳理工大学通信系统课程设计报告2 3.2 基本概念在 PCM 系统中， 为了得到二进制数字序列， 要对量化后的数字信号进行编码，每个抽样量化值用一个码组（码字）表示其大小。码长一般为7 位或 8 位，码长越大，可表示的量化级数越多，但编、解码设备就越复杂。那么能否找到其它更为简单的方法完成信号的模/ 数转换呢？我们看一下图 1。图中在模拟信号 f ( t ) 的曲线附近， 有一条阶梯状的变化曲线 f ( t ), f ( t )与 f (t ) 的形状相似。显然，只要阶梯“台阶” 和时间间隔 t 足够小，则 f ( t )与 f (t )的相似程度就会提高。对f ( t )进行滤波处理，去掉高频波动，所得到的曲线将会很好地与原曲线重合，这意味着f(t) 可以携带 f ( t ) 的全部信息（这一点很重要）。因此，f ( t )可以看成是用一个给定的“台阶” 对 f (t )进行抽样与量化后的曲线。我们把“台阶”的高度称为增量，用“ 1”表示正增量，代表向上增加一个；用“0”表示负增量，代表向下减少一个 。则这种阶梯状曲线就可用一个“0”、“1”数字序列来表示（如图（1）所示），也就是说，对f (t ) 的编码只用一位二进制码即可。此时的二进制码序列不是代表某一时刻的抽样值，每一位码值反映的是曲线向上或向下的变化趋势。这种只用一位二进制编码将模拟信号变为数字序列的方法（过程）就称为增量调制（ Delt a Modulat ion ），缩写为 DM或M调制。增量调制最早由法国人De Loraine 于 1946年提出， 目的是简化模拟信号的数字化方法。其主要特点是：(1)在比特率较低的场合，量化信噪比高于PCM 。(2)抗误码性能好。能工作在误比特率为102103 的信道中，而 PCM 则要求信道的误比特率为104106。(3)设备简单、制造容易。它与 PCM 的本质区别是只用一位二进制码进行编码，但这一位码不表示信号抽样值的大小，而是表示抽样时刻信号曲线的变化趋向。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 15 页 - - - - - - - - - 沈阳理工大学通信系统课程设计报告3 图 1 增量调制波形示意示4 课程设计分析4.1 M 的调制原理如何在发送端形成f (t ) 信号并编制成相应的二元码序列呢？仔细分析一上图（ 1），比较在每个抽样时刻 t 处的 f ( t )和 f (t ) 的值可以发现 , 当 f ( i t )f (i t_ )时，上升一个 ，发“1”码；当 f ( i t )0 ，则 Po(0)=1 (1) t =t 时， e( t )=f ( t )- f ( t_ )0，则 Po(t )=1 (2) t =2t 时，e(2 t )=f (2t )- f (2 t_ )0，则 Po(3t )=1; (4) t =4t 时，e(4 t )=f (4t )- f (4 t_ )0，则 Po(5t )=1; (6) t =6t 时，e(6 t )=f (6t )- f (6 t_ )0，则 Po(6t )=1; (7) 以此类推，即可得到如图3 所示的波形。会发现图3 中的 f ( t ) 和图 1 的波形不一样。 其实，图 1 的阶梯波只是为了形象地说明增量调制原理，而实际积分器的输出波形如图3d 所示。f (t)抽样判 决积分器f (t)T(t)e(t)Po(t)RC(a) 增量调 制器(编 码器 )框图(b) RC积分器名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 15 页 - - - - - - - - - 沈阳理工大学通信系统课程设计报告5 (C) 图 3 增量调制过程示意图4.2 M的解调原理为了完成整个通信过程， 发送端调制出的信号必须在接收端通过解调恢复出原始模拟信号。 M信号的解调比较简单，用一个和本地解码器一样的积分器即可。在接收端和发送端的积分器一般都是一个RC积分器。解调过程就是图43中的积分过程。当积分器输入“1”码时，积分器输出产生一个正斜变的电压并上升一个量化台阶 ；而当输入“ 0”码时，积分器输出电压就下降一个量化台tt02t 3t 4t 5 t 6t 7t 8t 9t 10t11t12t13t14tT(t)(a) 抽样脉冲(b) 样值信号tt02t3t4t5 t 6t 7t 8t 9t 10t11t 12t13t 14tf (t)f (t)(d) 积 分器 输出信 号t110101111110000t02 t 3 t 4 t 5t 6t 7 t 8 t 9 t 10 t11t12t13t14tPo(t)t0f (t)名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 15 页 - - - - - - - - - 沈阳理工大学通信系统课程设计报告6 阶。为了保证解调质量，对解码器有两个要求：（1） 每次上升或下降的大小要一致，即正负斜率大小一样。（2） （2）解码器应具有“记忆”功能， 即输入为连续“ 1”或“0”码时，输出能连续上升或下降。对积分器的输出信号进行低通滤波，滤除波形中的高频成分， 即可得到与原始模拟信号十分近似的解调信号，如图4 所示图 4 增量调制译码 ( 解调 ) 示意示5 增量调制 MATLAB 的仿真5.1 程序代码 : Ts=1e-3; t=0:Ts:20*Ts; x=sin(2*pi*50*t)+0.5*sin(2*pi*150*t); delta=0.4; D(1+length(t)=0; for k=1:length(t) 积分器低通滤波Po(t)f (t)fo(t)(a) 增量解调器(译码器)框图t0f (t)(b) 各点波形fo(t)名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 15 页 - - - - - - - - - 沈阳理工大学通信系统课程设计报告7 e(k)=x(k)-D(k); e_q(k)=delta*(2*(e(k)=0)-1); D(k+1)=e_q(k)+D(k); codeout(k)=(e_q(k)0); end subplot(3,1,1);plot(t,x,-o);axis(0 20*Ts,-2 2);hold on; subplot(3,1,2);stairs(t,codeout);axis(0 20*Ts,-2 2); Dr(1+length(t)=0; for k=1:length(t) eq(k)=delta*(2*codeout(k)-1); xr(k)=eq(k)+Dr(k); Dr(k+1)=xr(k); end subplot(3,1,3);stairs(t,xr);hold on; subplot(3,1,3);plot(t,x); 5.2 增量调制 Simulink 仿真实现采用 Simulink基本模块实现和采用DPCM 编解码模块实现。仿真测试模型如图（6）所示。仿真步进设置为0.001s ，模型中所有需要设置采样时间的地方均设置采样时间为0.001s 。在增量调制部分， Relay 模块作为量化器适应， 其门限设置为 0，输出值分别设置为0.4 和-0.4 ；Relay 作为编码器使用，其门限设置为 0，输出值设置为 1 和 0；解码端 Relay2 模块作为解码器使用， 其门限设置为 0.5 ，输出值分别为0.4 和-0.4 ；使用单位延时器Unit Delay作为预测滤波名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 15 页 - - - - - - - - - 沈阳理工大学通信系统课程设计报告8 器，初始状态均设置为零。使用 DPCM 编解码模块进行等价实现，DPCM 编码模块的设置是，预测器分子系数为 0,1 ，分母系数是 1，量化分割值为 0，码书为-0.4,0.4，解码器与编码器设置相同。仿真时间设置为0.02s ，即仿真前 20 个采样点。仿真结果如图 （7） 所示，采用 Simulink 基本模块实现的解码结果与编程法得到的波形相同。但是，由于初始值设置问题，采用DPCM 编解码模块得出的解码结果与采用Simulink 基本模块实现的解码结果在起始部分稍有不同，随着仿真时间的增加，两者输出结果相同。图 6 增量调制编码仿真测试模型其中 f(u)=sin(2*pi*50*u)+0.5*sin(2*pi*150*u) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 15 页 - - - - - - - - - 沈阳理工大学通信系统课程设计报告9 6 结果分析程序执行结果如图5 所示。从图中原信号和解码结果对比看，在输入信号变化平缓的部分， 编码器输出 1、0 交替码，相应的解码结果以正负阶距交替变化，形成颗粒噪声， 称空载失真； 在输入信号变化过快的部分，解码信号因不能跟踪上信号的变化而引起斜率过载失真。量化阶距越小， 则空载失真就越小， 但是容易发生过载失真； 反之，量化阶距增大， 则斜率过载失真减小， 但空载失真增大。如果量化阶距能根据信号的变化缓急自适应调整，则可以兼顾优化空载失真和过载失真，这就是自适应增量调制的意思。图 5 增量调制编码解码波形仿真结果（一）波形解析：第一个图形是原信号及离散样值第二个图形是编码输出二进制序列的波形名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 15 页 - - - - - - - - - 沈阳理工大学通信系统课程设计报告10 第三个图形解码结果和信号波形对比0.0040.006 为空载失真部分0.0090.012 为过载失真部分图 7 增量调制编码解码波形仿真结果（二） 6.1 增量调制存在的问题增量调制尽管有前面所述的不少优点，但它也有两个不足： 一个是一般量化噪声问题 ; 另一个是过载噪声问题。两者可统一称为量化噪声。观察图 1 可以发现，阶梯曲线（调制曲线）的最大上升和下降斜率是一个定值，只要增量 和时间间隔 t 给定，它们就不变。那么，如果原始模拟信号的变化率超过调制曲线的最大斜率， 则调制曲线就跟不上原始信号的变化，从而造成误差。我们把这种因调制曲线跟不上原始信号变化的现象叫做过载现象，由此产生的波形失真或者信号误差叫做过载噪声。另外，由于增量调制是利用调制曲线和原始信号的差值进行编码，也就是利用增量进行量化， 因此在调制曲线和原始信号之间存在误差，这种误差称为一般量化误差或一般量化噪声。两种噪声示意图如图8 所示。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 15 页 - - - - - - - - - 沈阳理工大学通信系统课程设计报告11 图 8 两种量化噪声示意图仔细分析两种噪声波形我们发现，两种噪声的大小与阶梯波的抽样间隔t和增量 有关。我们定义K 为阶梯波一个台阶的斜率式中, f s 是抽样频率。 该斜率被称为最大跟踪斜率。当信号斜率大于跟踪斜率时， 称为过载条件，此时就会出现过载现象；当信号斜率等于跟踪斜率时，称为临界条件；当信号斜率小于跟踪斜率时，称为不过载条件。可见，通过增大量化台阶（增量） 进而提高阶梯波形的最大跟踪斜率，就可以减小过载噪声； 而降低 则可减小一般量化噪声。 显然，通过改变量化台阶进行降噪出现了矛盾，因此， 值必须两头兼顾，适当选取。不过，利用增大抽样频率（即减小抽样时间间隔t ），却可以“左右逢源”，既能减小过载噪声，又可降低一般量化噪声。因此，实际应用中，M系统的抽样频率要比PCM 系统高得多（一般在两倍以上，对于话音信号典型值为16kHz和 32kHz）。 6.2 自适应增量调制增量调制中增量的幅值是固定的。若幅值选得过大， 粒状噪声过大； 若选得过小，超载噪声增加， 这给增量的幅度选择带来了一定的困难。为了解决这一问题， 可让增量的幅值在调制的过程中随着声音信号的变化自动地进行调制、变化，这就是自适应调制ADM(Adaptive Delta Modulation)。ADM 调制的基本原理是：在声音信号变化不大的情况下，取较小的增量幅值以抑制粒状噪声。 在声音信号变化较大的情况下， 预测信号跟不上声音信号的变化，应采取一定的算法增加增量的幅值，以此抑制超载噪声。调制过程中，增量的幅值随声音信号的变化自适应地变化。ttf (t)f (t)n(t)n(t)f (t)f (t)sKft名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 15 页 - - - - - - - - - 沈阳理工大学通信系统课程设计报告12 ADM 调制虽然能较好地克服超载噪声，解决粒状噪声和超载噪声的矛盾，但在声音信号从高速变化转向平坦处时，容易出现由于增量幅值过大而产生的噪声。6.2.1 自适应增量调制MATLAB 的仿真程序代码：Ts=1e-3; t=0Ts40Ts; x=sin(2pi50t)+0.5sin(2pi150t); x(2041)=0.2sin(2pi50t(2041); delta=0.4; D(1+length(t)=0; K=1.3; for k=1length(t) e(k)=x(k)-D(k); e_q(k)=delta(2(e(k)=0)-1); if k1 delta=delta(K.sign(e_q(k).e_q(k-1); end D(k+1)=e_q(k)+D(k); codeout(k)=(e_q(k)0); end 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 15 页 - - - - - - - - - 沈阳理工大学通信系统课程设计报告13 Dr(1+length(t)=0; delta=0.4; for k=1length(t) eq(k)=delta(2codeout(k)-1); if k1 delta=delta(K.sign(eq(k).eq(k-1); end xr(k)=eq(k)+Dr(k); Dr(k+1)=xr(k); end stairs(t,xr);hold on; plot(t,x); 自适应增量调制中，量化间距是自适应变化的：如果波形斜率陡峭，则连续输出的一串量化误差是同符号的，那么应使量化间距增大以减小斜率失真；如果波形平缓， 则连续输出的一串量化误差是正负符号交替的，这时减小量化间距就可以减小颗粒噪声。例如，一种较简单的自适应规则是)sgn(11nnnneeK（8）其中，自适应量化间距调整系数K 1。显然，当一串量化误差是同符号时，则sgn(1enne)0，于是1nn，即量化间距增加，反之，量化间距减少。程序运行结果如图 9 所示。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 15 页 - - - - - - - - - 沈阳理工大学通信系统课程设计报告14 图 9 自适应增量调制编码解码仿真结果波形解析：0.01 左右及 0.02 左右量化阶距增加0.015 左右及 0.025-0.035量化阶距减小，颗粒噪声随之减小6.3 增量调制应用的例子试建立 Simulink 模型，研究信道误码对增量调制的语音质量的影响。增量调制的采样率为 32kHz。仿真模型如图（ 10）所示，其中使用了PateTranstion模块将输入语音信号的采样率由 8000 次/s 升至 32000次/s ，然后进行增量调制。增量调制的预测器分子系数设置为 0,0.9以避免系统处于临界稳定状态。信道误码率可在01内任意设置。通过仿真聆听相应误码率下的恢复话音，主观感觉误码率在0.1时语音仍然具有相当的可懂度， 说明增量调制的抗噪声能力比PCM 强，但在无误码传输中，增量调制的解码音质不如PCM 。将信道误码率设置为0.5 ，则输出为纯噪声，相当于通信中断。有趣的是，名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 15 页 - - - - - - - - - 沈阳理工大学通信系统课程设计报告15 当将误码率设置为0.9 时，又可听到解码语音信号， 这说明增量调制能够抵抗信道传输中信号相位反转。 由于增量调制解码方法相当简单，实际中甚至直接使用积分器或低通滤波器对输入信号处理即可，也不需要时钟同步，抗干扰能力很强，所以在军用无线语音通信中得到了较为广泛的应用。心得体会由于基础不扎实， 以前 Matlab 并没怎么学， 这一次用上了， 还真是无从下手。在网上查看了很多资料， 本来以为随便找一个程序能运行就行了，但发现这样不行，还是必须自己会点， 因为根本不可能找到找到一个要求一摸一样的，最终还是要修改。参考了一些程序后，分析了一下，不会的就看书、百度，还是有点收获的。增量调制通信原理学过，基础还是没忘记，关键就是编程序了。由于对 Matlab 的陌生，好多函数找不到， 也不会用。 体现在这次课程设计中的就是不会求频谱。 问了好多同学， 还是有点模糊。有幸在看到一个求频谱的，原来就是使用的函数不会用。 虽然很简单， 但不会就是最难的事。 这次课程设计的确花了我很长的时间，但从中学到的东西还是令我欣慰的。7 参考文献1 郭文彬，桑林编著，通信原理-基于 Matlab 的计算机仿真，北京邮电大学出版社， 2006 2 曹志刚，钱亚生，现代通信原理，清华大学出版社，2003 年3 郭仕剑等， MATLAB7.x 数字信号处理，人民邮电出版社，2006 年4 张辉，曹丽娜著，通信原理学习指导，西安电子电子科技大学，2003 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 15 页，共 15 页 - - - - - - - - -