基于Matlab的脉冲编码调制PCM系统设计与仿真.docx
- 文档编号:13232983
- 上传时间:2023-06-12
- 格式:DOCX
- 页数:31
- 大小:159.67KB
基于Matlab的脉冲编码调制PCM系统设计与仿真.docx
《基于Matlab的脉冲编码调制PCM系统设计与仿真.docx》由会员分享,可在线阅读,更多相关《基于Matlab的脉冲编码调制PCM系统设计与仿真.docx(31页珍藏版)》请在冰点文库上搜索。
基于Matlab的脉冲编码调制PCM系统设计与仿真
课程设计(论文)任务书
欧阳光明(2021.03.07)
信息工程学 院 通信工程专 业142 班
一、课程设计(论文)题目 脉冲编码调制(PCM)系统设计与仿真
二、课程设计(论文)工作自1月3日起至1月13日止。
三、课程设计(论文)地址:
图书馆、寝室、通信实验室(4410)。
四、课程设计(论文)内容要求:
1.本课程设计的目的
(1)使学生掌握通信系统各功能模块的基本工作原理;
(2)培养学生采取Matlab与Simulink相结合对各种编码与解码进行仿真的办法;
(3)培养学生对PCM的理解能力;
(4)能提高和挖掘学生对所学知识的实际应用能力即立异能力;
(5)提高学生的科技论文写作能力。
2.课程设计的任务及要求
1)基本要求:
(1)学习Matlab与Simulink仿真软件的使用;
(2)对PCM,DPCM,ΔM编码与解码各功能模块的工作原理进行阐发;
(3)提出各种编码与解码电路的设计计划,选用合适的模块;
(4)对所设计系统进行仿真;
(5)并对仿真结果进行阐发。
a.采样定理的原理仿真
b.PCM编码与解码
c.DPCM编码与解码;增量调制(至少选做一种)
2)立异要求:
3)课程设计论文编写要求
(1)要依照书稿的规格打印誊写结业论文
(2)论文包含目录、绪论、正文、小结、参考文献、谢辞、附录等
(3)结业论文装订按学校的统一要求完成
4)辩论标准:
(1)完成原理阐发(20分)
(2)系统计划选择(30分)
(3)仿真结果阐发(30分)
(4)论文写作 (20分)
5)参考文献:
(1)王俊峰.《通信原理MATLAB仿真教程》人民邮电出版社第1版..11.1
(2)赵静.《基于MATLAB的通信系统仿真》北京航空航天年夜学出版社
6)课程设计进度安插
内容天数 地址
构思及收集资料2 图书馆
仿真5 实验室
撰写论文3 实验室
学生签名:
1月3日
课程设计(论文)评审意见
(1)完成原理阐发(20分):
优( )、良( )、中( )、一般( )、差( );
(2)系统计划选择(30分):
优( )、良( )、中( )、一般( )、差( );
(3)仿真结果阐发(30分):
优( )、良( )、中( )、一般( )、差( );
(4)论文写作 (20分):
优( )、良( )、中( )、一般( )、差( );
(5)格式规范性及考勤是否降品级:
是( )、否( )
评阅人:
职称:
副教授
1月13日
摘要
本设计结合PCM的抽样、量化、编码原理,利用MATLAB软件编程和绘图功能,完成了对脉冲编码调制(PCM)系统的建模与仿真阐发。
课题中主要分为三部分对脉冲编码调制(PCM)系统原理进行建模与仿真阐发,辨别为采样、量化和编码原理的建模仿真。
同时仿真阐发了采样与欠采样的波形、均匀量化与A律13折线非均匀量化的量化性能及其差别。
通过对脉冲编码调制(PCM)系统原理的仿真阐发,设计者对PCM原理及性能有了更深刻的认识,并进一步掌握MATLAB软件的使用。
关键词:
脉冲编码调制(PCM)均匀与非均匀量化MATLAB仿真
Abstract
Inthisdesign,combinationtheSimulinkemulatationfunctionandtheSfunction’sspreadfunctionofMATLABsoftware,havecompletedthesystematicemulatationandmodelingforpulsecodemodulation(PCM).Inthisdesign,divideinto3partsmainly,emulatetobuildmouldandemulateanalysisfortheprincipleofpulsecodemodulation(PCM)systematic.Theyaremodelingandemulatationofsampling,quantizingandecoding.Atthesametime,emulatetoanalysethewaveformofsamplingandowesampling,thequantizingerrorofuniformquantizingandnonuniformquantizing.Throughthisdesign,thedesignerhasamoreprofoundunderstandingofPCMprinciplesandperformance,andfurthermastertheuseofMATLABsoftware.
Keywords:
Pulsecodingmodulation(PCM)uniformandnonuniformquantitativeMATLABsimulation
1绪论
数字通信作为一种新型的通信手段,早在20世纪30年代就已经提出。
在1937年,英国人里费(A.H.Reeves)提出了脉冲编码调制(PCM)方法。
从此揭开了近代数字传输的序幕。
PCM系统的优点是:
抗干扰性强;失真小;传输特性稳定,远距离再生中继时噪声不累积,并且可以采取有效编码、纠错编码和保密编码来提高通信系统的有效性、可靠性和保密性。
另外,由于PCM可以把各种消息(声音、图像、数据等等)都变换成数字信号进行传输,因此可以实现传输和交换一体化的综合通信方法,并且还可以实现数据传输与数据处理一体化的综合信息处理。
故它能较好地适应信息化社会对通信的要求。
PCM的缺点是传输带宽宽、系统较庞杂。
可是,随着数字技术的奔腾成长这些缺点也不重要。
因此,PCM是一种极有成长前途的通信方法。
2PCM脉冲编码原理
2.1模拟信号的抽样及频谱阐发
2.1.1信号的采样
离散时间信号通常是有连续时间信号经周期采样获得的。
完成采样功能的器件称为采样器,下图所示为采样器的示意图。
图中Xa(t)暗示模拟信号,Xa(nt)暗示采样信号,T为采样周期,n=0,1,2,…。
一般可以把采样器视为一个每隔T秒闭合一次的电子开关S。
在理想情况下,开关闭合时间τ满足τ< 实际采样过程可视为脉冲调幅过程,Xa(t)为调制信号,被调脉冲载波p(t)是周期为T、脉宽为τ的周期脉冲串。 当τ→0时的理想采样情况是实际采样的一种科学的、实质的笼统,同时可使数学推导获得简化。 下面主要讨论理想采样。 图3.1采样器示意图及波形图 2.1.2抽样定理 抽样也称取样、采样,是把时间连续的模拟信号变换为时间离散信号的过程。 抽样定理是指: 一个频带限制在(0,fH)内的时间连续信号m(t),如果以T≤1/2fH秒的间隔对它进行等间隔抽样,则m(t)将被所获得的抽样值完全确定。 这意味着,若m(t)的频谱在某一角频率ωH上为零,则m(t)中的全部信息完全包含在其间隔不年夜于1/2fH秒的均匀抽样序列里。 换句话说,在信号最高频率分量的每一个周期内起码应抽样两次。 根据抽样脉冲的特性,抽样分为理想抽样、自然抽样(亦称曲顶取样)、瞬时抽样(亦称平顶抽样);根据被抽样信号的性质,抽样又分为低通抽样和带通抽样。 虽然抽样种类很多,可是间隔一按时间,抽样连续信号的样值,把信号从时间上离散,这是各种抽样共同的作用,抽样是模拟信号数字化及时分多路的理论基础。 我们考察一个频带限制在(0,fH)赫的信号m(t)。 假定将信号m(t)和周期性冲击函数δ(t)相乘,如图所示,乘积函数即是均匀间隔为T秒的冲激序列,这些冲激的强度即是相应瞬时上的m(t)值,它暗示对函数m(t)的抽样。 我们用ms(t)暗示此已抽样的函数,即有 ms(t)=m(t)δ(t) 上述关系如下图所示。 图3.2抽样示意图 2.1.3采样信号的频谱阐发 频谱阐发自然要使用快速傅里叶变换FFT了,对应的命令即fft,简单使用办法为: Y=fft(b,N),其中b即是采样数据,N为fft数据采样个数。 一般不指定N,即简化为Y=fft(b)。 Y即为FFT变换后获得的结果,与b的元素数相等,为单数。 以频率为横坐标,Y数组每个元素的幅值为纵坐标,画图即得数据b的幅频特性;以频率为横坐标,Y数组每个元素的角度为纵坐标,画图即得数据b的相频特性。 对现实中的情况,采样频率fs一般都是由采样仪器决定的,即fs为一个给定的常数;另一方面,为了获得一定精度的频谱,对频率辩白率F有一个人为的规定,一般要求F<0.01,即采样时间ts>100秒;由采样时间ts和采样频率fs即可决定采样数据量,即采样总点数N=fs*ts。 这就从理论上对采样时间ts和采样总点数N提出了要求,以包管频谱阐发的精准度。 2.2量化 2.2.1量化的界说 模拟信号进行抽样以后,其抽样值还是随信号幅度连续变更的,即抽样值m(kT)可以取无穷多个可能值,如果用N个二进制数值信号来代表该样值的年夜小,以便当用数字传输系统来传输该样值的信息,那么N个二进制信号只能同M=2^N个电平样值相对应,而不克不及同无穷多个电平值相对应。 这样一来,抽样值必须被划分红M个离散电平,此电平被称作量化电平。 或者说,采取量化抽样值的办法才干够利用数字传输系统来实现抽样值信息的传输。 利用预先规定的有限个电平来暗示模拟抽样值的过程称为量化。 抽样是把一个时间连续信号变换成时间离散的信号,而量化则是将取值连续的抽样变换成取值离散的抽样。 通常,量化器的输入是随机模拟信号。 可以用适当速率对此随机信号m(t)进行抽样,并依照预先规定,将抽样值m(kT)变换成M个电平q1,q2,…,qM之一,有 mq(kTs)=qi,若mi1≤m(kTs) 2.2.2量化的分类 (1)依照量化级的划分方法分,有均匀量化和非均匀量化。 均匀量化: 把输入信号的取值域按等距离联系的量化称为均匀量化。 在均匀量化中,每个量化区间的量化电平在各区间的中点。 其量化间隔Δv取决于输入信号的变更规模和量化电平数。 当信号的变更规模和量化电平数确定后,量化间隔也被确定。 上述均匀量化的主要缺点是,无论抽样值的年夜小如何,量化噪声的均方根都固定不变。 因此,当信号较小时,则信号量化噪声功率比也就很小,这样,对弱信号时的信号量噪比就很难达到给定的要求。 通常,把满足信噪比要求的输入信号取值规模界说为静态规模。 可见,均匀量化是的信号静态规模将受到较年夜的限制。 为了克服这一个缺点,实际中往往采取非均匀量化。 非均匀量化: 非均匀量化是根据信号的不合区间来确定量化间隔的。 对信号取值小的区间,其量化间隔也小;反之,量化间隔就年夜。 它与均匀量化相比,有两个突出的优点。 首先,当输入量化器的信号具有非均匀散布的几率密度时,非均匀量化器的输出端可以获得较高的平均信号量化噪声功率比;其次,非均匀量化时,量化噪声功率的均方根基本上与信号抽样值成比例。 因此量化噪声对年夜、小信号的影响年夜致相同,即改良了小信号时的信号量噪比。 罕见的非均匀量化有A律和μ率等,它们的区别在于量化曲线不合。 μ压缩律: 所谓μ压缩律就是压缩器的压缩特性具有如下关系的压缩律: 式中y为归一化的压缩器输出电压,x为归一化的压缩器输入电压,μ为压扩参数,暗示压缩的水平。 由于上式暗示的是一个近似对数关系,因此这种特性也称为近似对数压扩律,其压缩特性曲线如下图所示。 由图可知,当μ=0时,压缩特性是通过原点的一条直线,故没有压缩效果;当μ值增年夜时,压缩作用明显,对改良小信号的性能也有利。 一般当μ=100时,压缩器的效果就比较理想了。 另外,需指出,μ律压缩特性曲线是以原点奇对称的,图中只画出了正向部分。 图3.3μ压缩律特性 A压缩律: 所谓A压缩律也就是压缩器具有如下特性的压缩律: 其中,A为压缩系数;y为归一化的压缩器输出电压;x为归一化的压缩器输入电压。 图画出了A为某一取值的归一化压缩特性。 A律压缩特性是以原点奇对称的,为了简便,图中只给出了正半轴部分。 图3.4A压缩律特性 上图中,x和y都在1和+1之间,取量化级数为N(在y标的目的上从1到+1被均匀划分为N个量化级),则量化间隔为 当N很年夜时,在每一量化级中压缩特性曲线可看作是直线,因此有 式中,xi为第i个量化级间隔的中间值。 因此 (3.1) 为了使量化信噪比不随信号x变更,也就是说在小信号时的量化信噪比不因x的减小而变小,即应使各量化级间隔与x成线性关系,即 则式3.1可写成 (3.2) 即 其中k为比例常数。 当量化级数很年夜时,可以将它看成连续曲线,因而式(3.2)成为线性微分方程 解此微分方程 (3.3) 其中c为常数。 为了满足归一化要求,当x=1时,y=1,代入式(3.3)可得 故所得结果为 即 (3.4) 如果压缩特性满足上式,就可获得理想的压缩效果,其量化信噪比和信号幅度无关。 满足上式的曲线如下图所示,由于其没有通过坐标原点,所以还需要对它作一定的修改。 图3.5理想压缩特性曲线 A律压缩特性就是对式(3.4)修改后的函数。 在上图中,通过原点作理想压缩特性曲线的切线oc,将oc、cd作为实际的压缩特性。 修改以后,必须用两个不合的方程来描述这段曲线,以切点c为分界点, 线段oc的方程: 设切点c的坐标为(x1,y1)斜率为 则由式(3.4)可得 (3.5) 所以线段oc的方程为 所以当x=x1时,y1=1/k时,有 因此有 所以,切点坐标为(exp[(k1)],1/k),令 则 将它代入式(3.5),就可获得以切点c为鸿沟的 段的方程为 (3.6) 因cd段的方程,满足式(3.4),所以由该式可得 (3.7) 由以上阐发可见,经过修改以后的理想压缩特性与图5中所示的曲线近似,而式(3.6)式(3.7)和式(3.4)完全一样。 13折线: 实际中,A压缩律通常采取13折线来近似,13折线法如图747所示,图中先把 轴的[0,1]区间分为8个不均匀段。 图3.613折线示意图 其具体分法如下: a.将区间[0,1]一分为二,其中点为1/2,取区间[1/2,1]作为第八段; b.将剩下的区间[0,1/2]再一分为二,其中点为1/4,取区间[1/4,1/2]作为第七段; c.将剩下的区间[0,1/4]再一分为二,其中点为1/8,取区间[1/8,1/4]作为第六段; d.将剩下的区间[0,1/8]再一分为二,其中点为1/16,取区间[1/16,1/8]作为第五段; e.将剩下的区间[0,1/16]再一分为二,其中点为1/32,取区间[1/32,1/16]作为第四段; f.将剩下的区间[0,1/32]再一分为二,其中点为1/64,取区间[1/64,1/32]作为第三段; g.将剩下的区间[0,1/64]再一分为二,其中点为1/128,取区间[1/128,1/64]作为第二段; h.最后剩下的区间[0,1/128]作为第一段。 然后将y轴的[0,1]区间均匀地分红八段,从第一段到第八段辨别为[0,1/8],(1/8,2/8],(2/8,3/8],(3/8,4/8],(4/8,5/8],(5/8,6/8],(6/8,7/8],(7/8,1]。 辨别与x轴的八段一一对应。 采取上述的办法就可以作出由八段直线构成的一条折线,该折线和A压缩律近似,图3.6中的八段线段的斜率辨别为: 表1各段落的斜率 段落 1 2 3 4 5 6 7 8 斜率 16 16 8 4 2 1 1/2 1/4 从上表中可以看出,除一、二段外,其他各段折线的斜率都不相同。 图748中只画出了第一象限的压缩特性,第三象限的压缩特性的形状与第一象限的压缩特性的形状相同,且它们以原点为奇对称,所以负标的目的也有八段直线,总共有16个线段。 但由于正向一、二两段和负向一、二两段的斜率相同,所以这四段实际上为一条直线,因此,正、负双向的折线总共由13条直线段构成,这就是13折线的由来。 从A律压缩特性中可以看出,取A=87.6主要基于下述两个原因: 1使压缩特性曲线在原点邻近的斜率为16; 2当用13折线迫近时, 的八段量化分界点近似为1/2^n(n=0,1,2,…,7)。 从表1可以看出,当要求满足x=1/2^n时,相应有y=1n/8代入式中,有 因此有 将上式代入式(7.416),就可以获得对应A=94.4时的压缩特性 (3.8) 此压缩特性如果用13折线迫近,除第一段落起始点外,其余各段落的分界点的x、y都应满足式(3.8)。 在13折线中,第一段落起始点要求的x、y都应该为零,而若依照式(3.8)计算时,当x=0时,y→∞;而当y=0,x=1/2^8。 因此,需要对式(3.8)的压缩特性曲线作适当的修正,我们可以在原点和点(1/2^7,1/8)之间用一段直线取代原来的曲线,这段直线的斜率是1/8÷1/2^7=16。 为了找到一个能够暗示修正后的整个压缩特性曲线的方程,将式(3.8)酿成 (3.9) 从上式中可以看出,它满足x=0时,y=0;x=1时,y=1。 虽然式(3.9)在其他点上会有误差,但x在区间(1/128,1]内,1+255x都能和原来的256x比较接近。 所以,在绝年夜部分规模内的压缩特性仍和A律压缩特性很是接近,只有在x→0的小信号部分和A律压缩特性有些不同。 若在式(3.9)中,令μ=255,则式(3.9)可写成 (3.10) 式(3.10)的压缩特性与μ律压缩特性完全一致。 (2)依照量化的维数分,量化分为标量量化和矢量量化。 标量量化是一维的量化,一个幅度对应一个量化结果。 而矢量量化是二维甚至多维的量化,两个或两个以上的幅度决定一个量化结果。 以二维情况为例,两个幅度决定了平面上的一点。 而这个平面事先依照几率已经划分为N个小区域,每个区域对应着一个输出结果(码数,codebook)。 由输入确定的那一点落在了哪个区域内,矢量量化器就会输出那个区域对应的码字(codeword)。 矢量量化的好处是引入了多个决定输出的因素,并且使用了几率的办法,一般会比标量量化效率更高。 2.2.3MATLAB的A律13折线量化 在MATLAB中编写法度实现A律对数量化,并输出13折线对数量化特性曲线如图所示,法度见第4章设计内容。 图3.7A律13折线量化特性曲线 2.3PCM编码 2.3.1编码的界说 量化后的抽样信号在一定的取值规模内仅有有限个可取的样值,且信号正、负幅度散布的对称性使正、负样值的个数相等,正、负向的量化级对称散布。 若将有限个量化样值的绝对值从小到年夜依次排列,并对应地依次付与一个十进制数字代码(例如,付与样值0的十进制数字代码为0),在码前以“+”、“-”号为前缀,来区分样值的正、负,则量化后的抽样信号就转化为按抽样时序排列的一串十进制数字码流,即十进制数字信号。 简单高效的数据系统是二进制码系统,因此,应将十进制数字代码变换成二进制编码。 根据十进制数字代码的总个数,可以确定所需二进制编码的位数,即字长。 这种把量化的抽样信号变换成给定字长的二进制码流的过程称为编码。 话音PCM的抽样频率为8kHz,每个量化样值对应一个8位二进制码,故话音数字编码信号的速率为8bits×8kHz=64kb/s。 量化噪声随量化级数的增多和级差的缩小而减小。 量化级数增多即样值个数增多,就要求更长的二进制编码。 因此,量化噪声随二进制编码的位数增多而减小,即随数字编码信号的速率提高而减小。 自然界中的声音很是庞杂,波形极其庞杂,通常我们采取的是脉冲代码调制编码,即PCM编码。 PCM通过抽样、量化、编码三个步调将连续变更的模拟信号转换为数字编码。 2.3.2码型的选择 经常使用的二进制码型有自然二进制码和折叠二进制码两种。 折叠码优点: 只需对单极性信号进行,再增加最高位来暗示信号的极性;小信号的抗噪性能强,年夜信号的抗噪性能弱。 2.3.3PCM脉冲编码的原理 若信源输出的是模拟信号,如德律风机传送的话音信号,模拟摄象机输出的图像信号等,要使其在数字信道中传输,必须在发送端将模拟信号转换成数字信号,即进行A/D变换,在接收端则要进行D/A。 对语音信号最典范的数字编码就是脉冲编码调制(PCM)。 所谓脉冲编码调制: 就是将模拟信号的抽样量化值转换成二进制码组的过程。 下图给出了脉冲编码调制的一个示意图。 图3.8脉冲编码调制示意图 假设模拟信号m(t)的求值规模为[4V,+4V],将其抽样值按8个量化级进行均匀量化,其量化间隔为1s,因此各个量化区间的端点依次为4、3、2、1、0、1、2、3、4V,8个量化级的电平辨别为3.5、2.5、1.5、0.5、0.5、1.5、2.5和3.5V。 PCM系统的原理方框图如下图所示。 图中,输入的模拟信号m(t)经抽样、量化、编码后变换成数字信号,经信道传送到接收真个译码器,由译码器还原出抽样值,再经低通滤波器滤出模拟信号m^(t)。 其中,量化与编码的组合通常称为A/D变换器;而译码与低通滤波的组合称为D/A变换。 图3.9PCM通信系统方框图 3PCM的MATLAB实现 3.1PCM抽样的MATLAB实现 PCM抽样的MATLAB法度设计按如下步调进行: (1)确定输入的模拟信号为sa(200t); (2)根据输入的模拟信号,确定抽样频率,对输入信号进行抽样,并将正常抽样和会产生失真的抽样进行比较,对抽样定理加以验证; (3)编写法度,画出满足采样定理和不满足的时、频域图形。 PCM抽样的MATLAB实现源法度如下: functionsample() t0=10;%界说时间长度 ts=0.001;fs=1/ts; t=[t0/2: ts: t0/2];%界说时间序列 df=0.5;%界说频率辩白率 x=sin(200*t);m=x./(200*t+eps); w=t0/(2*ts)+1;%确定t=0的点 m(w)=1;%修正t=0点的信号值 m=m.*m; [M,mn,dfy]=fft_seq(m,ts,df);%傅立叶变换 M=M/fs; f=[0: dfy: dfy*length(mn)dfy]fs/2;%界说频率序列 figure (1) subplot(2,1,1);plot(t,m); xlabel('时间');ylabel('幅值');title('原始信号(fh=200/2piHz)的波形'); axis([0.15,0.15,0,1.5]); subplot(2,1,2); plot(f,abs(fftshift(M))); xlabel('频率');ylabel('幅值'); axis([500,500,0,0.03]);title('原始信号的频谱'); t0=10;%信号继续的时间 ts1=0.005;%满足抽样条件的抽样间隔 fs1=1/ts1; t1=[t0/2: ts1: t0/2];%界说满足抽样条件的时间序列 x1=sin(200*t1);m1=x1./(200*t1+eps);w1=t0/(2*ts1)+1; m1
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 Matlab 脉冲 编码 调制 PCM 系统 设计 仿真