增量调制仿真设计

1、增量调制仿真设计1. 课程设计目的(1) 加深理解增量调制编译码的基本原理。(2) 培养独立开展科研的能力和编程能力。(3) 了解系统的过载特性，动态范围以及最大虽化信噪比等三大指标的测试方法。2. 课程设计要求(1) 掌握课程设计的相关知识、概念清晰。(2) 程序设计合理、能够正确运行。相关知识3.1增量调制简介增虽调制简称或增虽脉码调制方式(DM,它是继PCM后出现的又一种模拟信号数字化的方法。1946年由法国工程师DeLoraine提出，目的在于简化模拟信号的数字化方法。主要在军事通信和卫星通信中广泛使用，有时也作为高速大规模集成电路中的A/D转换器使用。对模拟信号采样，并用每个样值与它

2、的预测值的差值对周期脉冲序列进行调制，简称增M或DM已调脉冲序列以脉冲的有、无来表征差值的正负号，也就是差值只编成一位二进制码。增虽调制的基本原理是于1946年提出的，它是一种最简单的差值脉冲编码。早期的语言增虽调制编码器是由分立元件组成的。随着模拟集成电路技术的发展，70年代末出现了音节压扩增虽调制集成单片，80年代出现了瞬时压扩集成单片，单片内包括了开关电容滤波器与开关电容积分器，集成度不断提高，使增虽调制的编码器的体积减小，功耗降低。3.2基本概念在PCM系统中，为了得到二进制数字序列，要对虽化后的数字信号进行编码，每个抽样虽化值用一个码组(码字)表示其大小。码长一般为7位或8位，码长越

3、大，可表示的虽化级数越多，但编、解码设备就越复杂。那么能否找到其它更为简单的方法完成信号的模/数转换呢？我们看一下图1。图中在模拟信号f(t)的曲线附近，有一条阶梯状的变化曲线f'(t),f'(t)与f(t)的形状相似。显然，只要阶梯“台阶和时间间隔t足够小，则f(t)与f(t)的相似程度就会提高。对f'(t)进行滤波处理，去掉高频波动，所得到的曲线将会很好地与原曲线里合，这意味着f'(t)可以携带f(t)的全部信息(这一点很S要)。因此，f'(t)可以看成是用一个给定的“台阶对f(t)进行抽样与虽化后的曲线。我们把“台阶”的高度C称为增虽，用“T表示正

4、增虽，代表向上增加一个C；用“0”表示负增虽，代表向下减少一个Co则这种阶梯状曲线就可用一个“0”、“T数字序列来表示(如图(1)所示)，也就是说，对f(t)的编码只用一位二进制码即可。此时的二进制码序列不是代表某一时刻的抽样值，每一位码值反映的是曲线向上或向下的变化趋势。这种只用一位二进制编码将模拟信号变为数字序列的方法(过程)就称为增虽调制(DeltaModulation)，缩写为DMI或AM调制。增虽调制最早由法国人DeLoraine于1946年提出，目的是简化模拟信号的数字化方法。其主要特点是：在比特率较低的场合，虽化信噪比高于PCM抗误码性能好。能工作在误比特率为102?103的信道

5、中，而PCM则要求信道的误比特率为104?106。设备简单、制造容易它与PCM的本质区另U是只用一位二进制码进彳亍编码，但这一位码不表示信号抽样值的大小，而是表示抽样时刻信号曲线的变化趋向。t0Tl010111111000二进制码序列编码后的数字信号图1增量调制波形示意示4课程设计分析4.1M的调制原理如何在发送端形成f'(t)信号并编制成相应的二元码序列呢？仔细分析上图(1),比较在每个抽样时刻t处的f(t)和f'(t)的值可以发现,当f(it)>f'(it_)时，上升一个c,发“1”码；当f(it)<f'(it_)时，下降一个c,发“0”码。f&

6、#39;(it_)是第i个抽样时刻前一瞬间的虽化值。根据上述分析，我们给出增虽调制器框图如图2所示。图2增量调制原理框图t=2At时,e(2At)=f(2At)-fz(2At_)<0，则Po(2At)=0;(3)t=3At时，e(3At)=f(3At)-f'(3At_)>0，则Po(3At)=1;(4)t=4At时,e(4At)=f(4At)-fz(4At_)<0，则Po(4At)=0;(5)t=5At时,e(5At)=f(5At)-fz(5At_)>0，则Po(5At)=1;(6)0旦0R二二Cf'(it_)可以由编码输出的二进制序列反馈到一个理想的积

7、分器以(b)RC积分器后得到。由于该积分器又具有解码功能，因此又称为本地解码器(译码器)。f(it)和f'(it_)的差值，可以用一个比较电路(减法器)来完成。员化编码可以用一个双稳判决器来执行，并生成双极性二进制码序列。具体调制过程描述如下：设f'(0-)=0(即t=0时刻前一瞬间的虽化值为零)，因此有t=0时，e(O)=f(O)-f'(0-)>0,则Po(0)=1t=At时，e(t)=f(t)-f'(t_)>0，贝UPo(t)=1t=6At时，e(6At)=f(6At)-f'以此类推，即可得到如图(6At_)>0，则Po(6At)=

8、1;(7)3所示的波形。会发现图3中的f'(t)和图1的波形不一样。其实，图1的阶梯波只是为了形象地说明增虽调制原理，而实际积分器的输出波形如图3d所示。(a)抽样脉冲(b)样值信号1T(t)i.:?:4闩i5t了:e;宕':矣t：2di:.t14-|at1112t%t|_1122tiSt14t0tJ3t15t6t78t：9t1"n(C)tp(t)110101111110002XL35XL67匚9XLTo2XL34f(t)(d)积分器输出信号图3增量调制过程示意图4.2M的解调原理为了完成整个通信过程，发送端调制出的信号必须在接收端通过解调恢复出原始模拟信号?M信号的

9、解调比较简单，用一个和本地解码器一样的积分器即可。在接收端和发送端的积分器一般都是一个RC积分器。解调过程就是图43中的积分过程。当积分器输入“T码时，积分器输出产生一个正斜变的电压并上升一个虽化台阶c;而当输入“0”码时，积分器输出电压就下降一个虽化台为了保证解调质虽，对解码器有两个要求：每次上升或下降的大小要一致，即正负斜率大小一样(2)解码器应具有“记忆”功能，即输入为连续“T或“0”码时，输出能连续上升或下降。对积分器的输出信号进行低通滤波,滤除波形中的高频成分，即可得到与原始模拟信号十分近似的解调信号，如图4'-IPo(tf(t)低通滤沌(a)增虽解调器(译码器)框图图4增量

10、调制译码(解调)示意示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;fork=1:length(t)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);endsubplot(3,1,1);plot(t,x,'-o');axis(020*Ts,-22);holdon;subplot(3,1

11、,2);stairs(t,codeout);axis(020*Ts,-22);Dr(1+length(t)=0;fork=1:length(t)eq(k)=delta*(2*codeout(k)-1);xr(k)=eq(k)+Dr(k);Dr(k+1)=xr(k);endsubplot(3,1,3);stairs(t,xr);holdon;subplot(3,1,3);plot(t,x);5.2增量调制Simulink仿真实现采用Simulink基本模块实现和采用DPCM?解码模块实现。仿真测试模型如图(6)所示。仿真步进设置为0.001s,模型中所有需要设置采样时间的地方均设置采样时间为0.

12、001s。在增虽调制部分，Relay模块作为虽化器适应，其门限设置为0,输出值分别设置为0.4和-0.4;Relay作为编码器使用，其门限设置为0,输出值设置为1和0;解码端Relay2模块作为解码器使用，其门限设置为0.5,输出值分别为0.4和-0.4;使用单位延时器UnitDelay作为预测滤波器，初始状态均设置为零使用DPCM?解码模块进行等价实现，DPCMS码模块的设置是，预测器分子系数为0,1，分母系数是1,i化分割值为0,码书为卜0.4,0.4，解码器与编码器设置相同。仿真时间设置为0.02s,即仿真前20个采样点。仿真结果如图（7）所示，采用Simulink基本模块实现的解码结果

13、与编程法得到的波形相同。但是，由于初始值设置问题，采用DPCM编解码模块得出的解码结果与采用Simulink基本模块实现的解码结果在起始部分稍有不同，随着仿真时间的增加，两者输出结果相同。Wu-rtitledfieviewSimulationFn?itH?p|s;口(三)*0IOIR显址10T%OtJ期图6增量调制编码仿真测试模型其中f(u)=sin(2*pi*50*u)+0.5*sin(2*pi*150*u)6结果分析程序执行结果如图5所示。从图中原信号和解码结果对比看，在输入信号变化平缓的部分，编码器输出1、0交替码，相应的解码结果以正负阶距交替变化，形成颗粒噪声，称空载失真；在输入信号变

14、化过快的部分，解码信号因不能跟踪上信号的变化而引起斜率过载失真。虽化阶距越小，则空载失真就越小，但是容易发生过载失真；反之，H化阶距增大，则斜率过载失真减小，但空载失真增大。如果虽化阶距能根据信号的变化缓急自适应调整，则可以兼顾优化空载失真和过载失真，这就是自适应增H调制的意思。图5增量调制编码解码波形仿真结果（一）波形解析:第一个图形是原信号及离散样值第二个图形是编码输出二进制序列的波形第三个图形解码结果和信号波形对比0.0040.006为空载失真部分0.0090.012为过载失真部分图7增量调制编码解码波形仿真结果（二）6.1增量调制存在的问题增虽调制尽管有前面所述的不少优点，但它也有两个

15、不足：一个是一般虽化噪声问题；另一个是过载噪声问题。两者可统一称为虽化噪声。观察图1可以发现，阶梯曲线（调制曲线）的最大上升和下降斜率是一个定值，只要增虽c和时间间隔t给定，它们就不变。那么，如果原始模拟信号的变化率超过调制曲线的最大斜率，则调制曲线就跟不上原始信号的变化，从而造成误差。我们把这种因调制曲线跟不上原始信号变化的现象叫做过载现象，由此产生的波形失真或者信号误差叫做过载噪声。另外，由于增虽调制是利用调制曲线和原始信号的差值进行编码，也就是利用增虽进行虽化，因此在调制曲线和原始信号之间存在误差，这种误差称为一般虽化误差或一般虽化噪声。两种噪声示意图如图8所示。n(t)f(t)f(t)

16、图8两种量化噪声示意图仔细分析两种噪声波形我们发现，两种噪声的大小与阶梯波的抽样间隔t和增虽c有关。我们定义K为阶梯波一个台阶的斜率一fs式中,fs是抽样频率。该斜率被称为最大跟踪斜率。当信号斜率大于跟踪斜率时，称为过载条件，此时就会出现过载现象；当信号斜率等于跟踪斜率时，称为临界条件；当信号斜率小丁跟踪斜率时，称为不过载条件。可见，通过增大虽化台阶（增虽）C进而提高阶梯波形的最大跟踪斜率，就可以减小过载噪声；而降低C则可减小一般虽化噪声。显然，通过改变虽化台阶进行降噪出现了矛盾，因此，C值必须两头兼顾，适当选取。不过，利用增大抽样频率（即减小抽样时间间隔At）,却可以“左右逢源”，既能减小过

17、载噪声，又可降低一般虽化噪声。因此，实际应用中，M系统的抽样频率要比PCM系统高得多（一般在两倍以上，对于话音信号典型值为16kHz和32kHz）。6.2自适应增量调制增虽调制中增虽的幅值是固定的。若幅值选得过大，粒状噪声过大；若选得过小，超载噪声增加，这给增虽的幅度选择带来了一定的困难。为了解决这一问题，可让增虽的幅值在调制的过程中随着声音信号的变化自动地进行调制、变化,这就是自适应调制ADM（AdaptiveDeltaModulation）。ADMM制的基本原理是：在声音信号变化不大的情况下，取较小的增虽幅值以抑制粒状噪声。在声音信号变化较大的情况下，预测信号跟不上声音信号的变化，应采取一

18、定的算法增加增虽的幅值，以此抑制超载噪声。调制过程中，增虽的幅值随声音信号的变化自适应地变化。ADM调制虽然能较好地克服超载噪声，解决粒状噪声和超载噪声的矛盾，但在声音信号从高速变化转向平坦处时，容易出现由于增虽幅值过大而产生的噪声。621自适应增量调制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;fork=1length(t)e(k)=x(k)-D(k);e_q(k)=delta(2(e(k)=0

19、)-1);ifk1delta=delta(K.Asign(e_q(k).e_q(k-1);endD(k+1)=e_q(k)+D(k);codeout(k)=(e_q(k)0);endDr(1+length(t)=O;delta=0.4;fork=1length(t)eq(k)=delta(2codeout(k)-1);ifk1delta=delta(K.Asign(eq(k).eq(k-1);endxr(k)=eq(k)+Dr(k);Dr(k+1)=xr(k);endstairs(t,xr);holdon;plot(t,x);自适应增虽调制中，虽化间距是自适应变化的：如果波形斜率陡峭，则连续输

20、出的一串虽化误差是同符号的，那么应使虽化间距增大以减小斜率失真；如果波形平缓，则连续输出的一串虽化误差是正负符号交替的，这时减小虽化间距就可以减小颗粒噪声。例如，一种较简单的自适应规则是nn1KSgn(enen1)(8)其中，自适应虽化间距调整系数K1。显然，当一串虽化误差是同符号时，则sgn(nn1)>0,于是nn1,即虽化间距增加，反之，虽化间距减少。程序运行结果如图9所示。图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内任意设置。通过仿真聆听相应误码率下的恢复话音