4ASK数字频带系统仿真.docx

1、4ASK数字频带系统仿真*实践教学*兰州理工大学计算机与通信学院2015年春季学期通信系统仿真训练课程设计报告题 目 4ASK数字调制系统的设计与仿真 班级： 通信工程12级（1）班 姓名： 朱芮碧 设计质量（30分）： 学号： 12250108 说明书质量（10分）： 同组成员： 汪永前高泽雷张佩君李慧慧梁汉聪刘昆鑫樊丰收摘要本文在分析通信系统的原理上，利 用MATLAB来仿真一个基于4ASK（四进制振幅键控） 调制解调的数字通信系统。 该系统的仿真设计为数字通信系统的优化及设计方法提供了科学 理论依据。该仿真基本能够满足题目的任务要求，完成了数字通信系统的仿真与测试。通过 实际方案的分析比

2、较、参数计算、仿真测试等环节，初步巩固了课程所学的有关理论知识， 加深了对数字通信的理解和掌握；掌握了 4ASK的原理和应用，学习了 MATLAB件的基本使 用；学会运用MATLAB件进行一些仿真和设计。最后，在此基础上进行了满足题目要求的优 化与调试，完成了题目的设计任务。该数字通信系统利用MATLAB件对基带信号通过四进制振幅键控来调制和解调，经过产 生信息序列、变为基带信形、调制、加噪、解调、接收信息序列、误码等七个步骤完成通信 系统的设计。本文着重介绍了 数字通信系统 的设计思路、方案论证、信号分析，通过 MATLAB 的软件调试，实现了设计之初的目标，论证了设计方案的可行性。通过数字

3、通信系统的仿真 与测试，加深对基本原理的了解，增强了我们的实践能力。摘要 1一、前言 3二、系统概述 42.1课题研究背景 42.2设计概要 42.3设计任务与设计要求 4三、系统设计原理 53.1数字通信系统 53.24ASK 信号的原理 63.34ASK 调制解调原理 7四、 4ASK调制解调系统仿真 94.1调制系统仿真 94.2解调系统仿真 104.3误码率分析 11五、 总结 13六、参考文献 14附录 15源代码 15一、前言当今社会已经步入信息时代 . 在各种信息技术中， 信息的传输及通信起着支撑作用。 而对 于信息的传输，开发更加完善、快捷的通信系统已经变得越来越重要，数字通信

4、就成为重要 的手段。因此，数字信号的调制就显得非常重要。实际通信中的许多信道都不能直接传送基 带信号，必须用基带信号对载波波形的某些参量进行控制，使得载波的这些参量随基带信号 的变化而变化，即正弦载波调制。在数字通信系统中，有二进制数字调制和多进制调制。多 进制数字调制与二进制数字调制相比又具有如下两个特点，在相同的码元传输率下，多进制 系统的信息传输率比二进制系统的高，在相同的信息速率下，多进制信号码元的持续时间要 比二进制的长，因此会增加码元的能量，减小信号特性引起的码间干扰的影响，利用层次化 和模块化的设计方法，通过 MATLAB软件平台设计并实现了多进制幅移键控( M-aryAmpli

5、tude-Shift Keying ， MAS)K 中的四电平调制( 4-ary Amplitude Shift Keying ， 4ASK) 的调制系统和解调系统。二、系统概述2.1课题研究背景远古时代，远距离的传递消息是以书信的形式来完成的 , 这种通信方式明显具有传递时间 长的缺点。为了在尽量短的时间内传递尽量多的消息 , 人们不断地尝试所能找到的各种最新技 术手段。 1837 年发明的莫尔斯电磁式电报机标志着电通信的开始 ,之后，利用电进行通信的 研究取得了长足的进步。 1866年利用海底电缆实现了跨大西洋的越洋电报通信。 1876 年贝尔 发明了电话 , 利用电信号实现了语音信号的有

6、线传递使信息的传递变的既迅速又准确， 这标志 着模拟通信的开始， 由于它比电报更便于交流使用， 所以直到 20 世纪前半叶这种采用模拟技 术的电话通信技术比电报的到了更为迅速和广泛的发展。 1937 年瑞威斯发明的脉冲编码调制 标志数字通信的开始。 20 世纪 60 年代以后集成电路、电子计算机的出现，使得数字通信迅 速发展。在 70 年代末在全球发展起来的模拟移动电话在 90 年代中期被数字移动电话所代替 现有的模拟电视也正在被数字电视所代替。数字通信的高速率和大容量等各方面的优越性也 使人们看到了它的发展前途。2.2设计概要本次课设主要通过研究4ASK信号的调制解调，首先通过对二进制 2A

7、SK的分析来研究出 四进制4ASK的变化，对2ASK的基带信号和传输的载波信号，以及其波形图进行分析，从而 掌握多进制的振幅键控（MASK调制解调的原理及其实现方法，然后利用 MATLAB仿真实现4ASK的调制与解调，并仿真4ASK载波信号在高斯白噪声下的误码率和误比特率的性能，同 时给出调制信号、载波信号及已调信号的波形图和频谱图。 最后根据仿真的波形图来分析 4ASK 的性能特点，以及对以后信道的传输有更重要的意义和频带利用率，资源有效充分利用全方 面的来考虑4ASK的用途。2.3设计任务与设计要求1.设计一个数字通信仿真系统，调制为 4ASK；2.能显示信息序列、调制信号等信号波形；3.

8、解调部分，相关解调和非相关解调都可的，相关解调和非相关解调都做；4.对在不同信噪比情况下的误码指标。过程：产生信息序列变为基带波形调制加噪解调接受信息序列误码（注意其中信号的变化）三、系统设计原理3.1数字通信系统数字通信(digital telecommunications )是用数字信号作为载体来传输消息，或用数字信号对载波进行数字调制后再传输的通信方式。它可传输电报、数字数据等数字信号，也 可传输经过数字化处理的语声和图像等模拟信号。数字通信系统的模型如图 3.1所示：图3.1数字通信系统模型图中信源输出的是模拟信号，经过数字终端的信源编码器成为数字信号，终端输出的数 字信号，经过信道编

9、码器后变成适合于信道传输的数字信号，然后由解调器把数字信号调制 到系统所使用的数字信道上，再传输到接收端，经过相反的转换后最终送到信宿。数字通信系统各部分作用：1、 信源：把原始信息变换成原始电信号。2、 信源编码：1实现模拟信号的数字化传输即完成 A/D变化。2提高信号传输的有效性。即在保证一定传输质量的情况下，用竟可能少的数字脉冲来 表示信源产生的信息。信源编码也称作频带压缩编码或数据压缩编码。3、 信道编码：1信源编码的目的：信道编码主要解决数字通信的可靠性问题。2信道编码的原理：对传输的信息码元按一定的规则加入一些冗余码(监督码)，形成新 的码字，接收端按照约定好的规律进行检错甚至纠错

10、。3信道编码又称为差错控制编码、抗干扰编码、纠错编码 。4、 数字调制：1数字调制技术的概念：把数字基带信号的频谱搬移到高频处，形成适合在信道中传输 的频带信号。2数字调制的主要作用：提高信号在信道上传输的效率，达到信号远距离传输的目的。3基本的数字调制方式：振幅键控 ASK频移键控FSK相移键控PSK5、 同步：1同步的概念：指通信系统的收、发双方具有统一的时间标准，使它们的工作“步调一 致”。2同步的作用：对于数字通信时是至关重要的。如果同步存在误差或失去同步，通信过 程中就会出现大量的误码，导致整个通信系统失效。6、 信道： 信道是信号传输媒介的总称，传输信道的类型有无线信道(如电缆、光

11、纤)和有线信道(如自由空间)两种。7、 噪声源： 通信系统中各种设备以及信道中所固有的，为了分析方便，把噪声源视为各处噪声的集中表现而抽象加入到信道。数字通信系统较模拟通信系统而言，具有抗干扰能力强、差错可控、便于加密、易于实 现集成化、便于与计算机连接等优点。因而，数字通信更能适应对通信技术的越来越高的要 求。近二十年来，数字通信发展十分迅速，在整个通信领域中所占比重日益增长，在大多数 通信系统中已代替模拟通信，成为当代通信系统的主流。3.2 4ASK 信号的原理实际通信中的许多信道都不能直接传送基带信号，必须用基带信号对载波波形的某些参 量进行控制，使得载波的这些参量随基带信号的变化而变化

12、，即正弦载波调制。在数字通信 系统中，有二进制数字调制和多进制调制。多进制数字调制与二进制数字调制相比又具有如 下两个特点，在相同的码元传输率下，多进制系统的信息传输率比二进制系统的高，在相同 的信息速率下，多进制信号码元的持续时间要比二进制的长，因此会增加码元的能量，减小 信号特性引起的码间干扰的影响，利用层次化和模块化的设计方法，通过MATLAB件平台设 计并实现了多进制幅移键控 (M-ary Amplitude-Shift Keying，MASK中的四电平调制(4-ary Amplitude Shift Keying ，4ASK的调制系统和解调系统。本文首先介绍了四电平调制和解调的原理，

13、随后介绍载波产生、振幅调制、振幅判别等功能模块的设计，最后给出了整体调 制解调的模块图和仿真波形及测试结果。多进制数字幅度调制( 4ASK又称为四电平调制，它是二进制数字幅度调制方式的推广。四进制幅度调制信号的载波振幅有四种取值，在一个 码元期间内，发送其中的一种幅度的载波信号。 MAS已调信号的表示式为：st = Dan ? g(t - nTb )s(t)这里为四进制数字基带信号，式中，g(t)是高度为1、宽度为Tb的门函数；an有4种取值 0,1,2,3 ，出现的概率分别为 P0,P1,P2,P3,且P0+P1+P2+P3=1.图3.2和图3.3分别为四进制 数字基带信号和已调信号的波形图

14、。图3.3已调信号波形图3.34ASK调制解调原理图3.4数字调制信号基本结构数字调制与模拟调制原理是相同的，一般可以采用模拟调制的方法实现数字调制。 但是， 数字基带信号具有与模拟基带信号不同的特点，其取值是有限的离散状态。这样，可以用载 波的某些离散状态来表示数字基带信号的离散状态。基本的三种数字调制方式是 ：振幅键控(ASK)、移频键控(FSK)和移相键控(PSK或DPSK)多进制数字振幅调制又称多电平振幅调制， 它用高频载波的多种振幅去代表数字信息。MAS碉制方式就是使载波的幅值随着基带信号的 变化而变化。MASK勺表达式如下所示:其中Ao为基带信号的电平，为载波顺率。由上式可以看出，

15、如果其中电平是0的多进 制信号,只要让载波信号与多进制信号通过乘法器即可调制完成。如果两个电平都不是 0,只要让载波信号的振幅固定，通过乘法器与多进制信号相乘就行。四电平振辐调制，高频载波有四种：振幅为0、1A、2A和3A,分别代表数字信息0、1、 2、3或者双比特二进制输入信息00、01、10、11进行振幅调制。而解调的两种方法包络检波和相干解调的原理图如图 3.5和图3.6所示。 图3.5相干法解调示意图图3.6包络检波法示意图4ASK的基带信号只有“ 0”、“ 1”、“2”、“3”，四个电平值，它与载波相乘的结果相当于 将载波关断，或者接通放大。它的实际意义是当调制的数字信号为“ 3”时

16、，假设传输振幅为126个量化单位的载波，则当调制的数字信号分别为“ 2”、1、0时，传输振幅分别为 84、42，0个量化单位的载波。其典型波形如图1所示。4ASK的键控调制原理如图2所示载波通过基带信号的控制选择不同的开关，当基带信号是“ 0”时，调制信号的幅度为0个量化单位；当基带信号分别是“ 1”、“2”、“3”，选择开关接通相应的乘法器，则调制信号的幅 度就分别为载波信号幅度的1,2, 3 倍。对于4ASK的解调，我们采用最高幅值判别的方法，在调制信号中检测出最大幅值，随后 根据最大幅值与基带信号的对应关系就可以解调出来。四、4ASK调制解调系统仿真4.1调制系统仿真在软件设计过程中，先

17、设计出二进制基带，信号使用语句: subplot(3,1,1);plot(t,a(ceil(100*t+0.1)/5)然后在此基础上生成四进制基带信号，使用语句 ：sym( n+1)=a(2* n+1)*2+a(2* n+2); subplot(3,1,2);plot(t,sym(ceil(10*t+0.01)具体结果如图4.1所示：图4.1二进制基带信号变为四进制基带信号将生成的四进制基带信号进行4ASK调制，然后对4ASK信号进行功率谱密度分析，加入 高斯白噪声噪经过信道，到达解调部分。运行相关MATLAB?序，得到4ASK调制波形如图4.2 所示，将前后波形进行对比、分析。图4.2 4A

18、SK调制结果4.2解调系统仿真1.米用相关解调法（或同步检测法）进行解调，信号经过相关流程后，得到4ASK相关解调波形如图4.3所示图4.3 4ASK相关解调结果2.采用非相关解调法（或包络检波法）进行解调，信号经过相关流程后，得到4ASK非相关 解调波形如图4.4所示图4.4 4ASK非相关解调结果输入基带信号的直流分量最大，经 4ASK调制后载波频率分量占主要成分，故 4ASK信号 的功率谱密度在f= 2Hz出现峰值。将加噪后的4ASK信号进行相干解调后通过低通滤波后得 到的信号发生明显失真，通过对失真的信号抽样判决后很好的恢复出原基带信号。4.3误码率分析误码率统计中，直接调用 MATL

19、A的函数randn产生均值为0,方差为1的加性高斯随机 噪声，对每一次加入噪声后解调输出的四进制序列与输入的四进制序列进行对比，计算解调 后的误码数及其比特率，结果如图 4.5所示。图4.5误码率分析五、总结通过前面的设计分析和原理叙述，我们最终可以总结以下：首先，在为期三个星期的课程设计中，我们终于完成了课程设计的基本要求，完成了题 目的设计任务。本次课设主要通过研究 4ASK信号的调制解调，首先通过对二进制 2ASK的分 析来研究出四进制4ASK的变化，对2ASK的基带信号和传输的载波信号，以及其波形图进行 分析，从而掌握多进制的振幅键控（MASK调制解调的原理及其实现方法，然后利用 MA

20、TLAB 仿真实现4ASK的调制与解调，并仿真4ASK载波信号在高斯白噪声下的误码率和误比特率的 性能，同时给出调制信号、载波信号及已调信号的波形图和频谱图。最后根据仿真的波形图 来分析4ASK的性能特点，以及对以后信道的传输有更重要的意义和频带利用率， 资源有效充分利用全方面的来考虑4ASK的用途。其次，通过这次课程设计，我懂得了理论与实际相结合的重要性，只有理论知识是远远 不够的。只有把所学到的理论知识与实践相结合起来，从理论中得出结论，将结论辅助与理 论，才能真正学到知识。所以提高自己的实际动手能力与独立思考的能力是十分重要的，在 设计的过程中遇到很多的问题，由于对以前所学过的知识理解不

21、深刻，掌握不牢靠，导致在 设计过程中出现一些错误。还有一些小的问题，这是由于个人的不认真思考造成的，这也让 我认识到整体设计一定要仔细，要踏踏实实的多查资料，多问，多看。还要善于同别人交流 与合作，善于获取各种有用的资源。最后，在课程设计与论文整理期间，老师给了我很大的支持和帮助，才使得课程设计得 以顺利完成，在此谨向导师表示忠心的感谢和崇高的敬意。同时，对我而言，知识上的收获 重要，精神上的丰收更加可喜。挫折是一份财富，经历是一份拥有。这次实习必将成为我人 生旅途上一个非常美好的回忆！六、参考文献1. 李建新现代通信系统分析与仿真一MATLAB信工具箱.西安：西安电子科技大学出版社,2000

22、 .2. 樊昌信 . 通信原理 . 北京： 国防工业出版社 ,2002.3. 刘敏. MATLAB 通信仿真与应用 . 北京： 国防工业出版社 .4. 曹志刚等著 . 现代通信原理 . 北京：清华大学出版社 .2001. 55. 吴伟陵等著 . 移动通信原理 . 北京：电子工业出版社， 2005.附录：源代码：clear all;close all;A=1;fc = 2; % 设定载波频率 2HzN_sample = 100;% 插入点数为 2*200-1N = 100;% 考察的码元个数Ts = 1;% 码元周期dt = Ts/fc/N_sample; % 波形采样间隔t = 0:dt:N*

23、Ts-dt;Lt = length(t);%二进制基带信号变为四进制基带信号M=4;d=1;t=0:1/1e3:0.999;a=randint(1,20,2);for n=0:9sym(n+1)=a(2*n+1)*2+a(2*n+2);ends=sym(ceil(10*t+0.01).*cos(2*pi*100*t);subplot(2,1,1);plot(t,a(ceil(100*t+0.1)/5);axis(0,1,-0.2,1.2);subplot(2,1,2);plot(t,sym(ceil(10*t+0.01);%产生四进制信源d =rand(1,N);d1=floor(d*4);d

24、d = sigexpand(d1,fc*N_sample);gt = ones(1,fc*N_sample); %NRZ 波形d_NRZ = conv(dd,gt);% 基带信号 d_NRZ dd*gt(t-nTs)ht = A*cos(2*pi*fc*t);% 载波s_4ask = d_NRZ(1:Lt)*ht; % 进行 4ASK调希9 %4AS信号 s_4ask=刀 dd*gt(t-nTs)*Acos(2*pi*fc*t) figure(1)subplot(2,2,1);plot(t,d_NRZ(1:length(t);axis(0 50 0 4); xlabel( 输入信号);subp

25、lot(2,2,2);plot(t,s_4ask);axis(0 50 -4 4); xlabel(4ASK);f, s_4askf=T2F(t,s_4ask);subplot(2,2,3);plot(f,10*log10(abs(s_4askf).A2/Ts);axis(-10 10 -20 40); xlabel(4ASK 功率谱密度 );%加噪t4ask=awgn(s_4ask,15);% 在信号 s_4ask 中加入高斯白噪声 , 信噪比 15dB。 subplot(2,2,4);plot(t,t4ask)axis(0 50 -4 4); xlabel( 加噪 4ASK);%相干解调后

26、通过低通滤波器tt4ask=2*t4ask.*cos(2*pi*fc*t);% 相干解调f, tt4askf=T2F(t,tt4ask);B=2/Ts;%B=2*fs,fs=N/T,T=N*Tst st=lpf(f,tt4askf,B);% 低通figure(2)subplot(2,2,1)plot(t,tt4ask);axis(0 50 0 8); xlabel( 相干解调输出信号 ); subplot(2,2,2) plot(t,st);axis(0 50 0 4); xlabel( 通过低通滤波后信号 ); %抽样判决for k=1:floor(length(st)/(2*N_sampl

27、e)TT4ask(k)=floor(st(k+(k-1)*(2*N_sample)+N_sample)+0.5);% 在每个码元中间抽一次样，得到的值四舍五入endtt= sigexpand(TT4ask,fc*N_sample);gtt = ones(1,fc*N_sample); %NRZ 波形 tt_NRZ = conv(tt,gtt);subplot(2,2,3) plot(t,tt_NRZ(1:length(t);axis(0 50 0 4); xlabel( 抽样判决后信号 );f, stf=T2F(t,tt_NRZ);subplot(2,2,4)plot(f,10*log10(a

28、bs(stf)A2/Ts);axis(-10 10 -20 40); xlabel( 输入信号功率谱密度 );%非相干解调后通过低通滤波器tt4ask=2*t4ask.*cos(2*pi*fc*t);% 非相干解调f, tt4askf=T2F(t,tt4ask);B=2/Ts;%B=2*fs,fs=N/T,T=N*Tst st=lpf(f,tt4askf,B);% 低通figure(3)subplot(2,2,1)s=abs(tt4ask);plot(t,s);axis(0 50 0 8); xlabel( 非相干解调输出信号 ); subplot(2,2,2)plot(t,st);axis(

29、0 50 0 4); xlabel( 通过低通滤波后信号 );%抽样判决for k=1:floor(length(st)/(2*N_sample) TT4ask(k)=floor(st(k+(k-1)*(2*N_sample)+N_sample)+0.5);% 在每个码元中间抽一次样，得到的值四舍五入endtt= sigexpand(TT4ask,fc*N_sample);gtt = ones(1,fc*N_sample); %NRZ 波形tt_NRZ = conv(tt,gtt); subplot(2,2,3)plot(t,tt_NRZ(1:length(t);axis(0 50 0 4);

30、 xlabel( 抽样判决后信号 );f, stf=T2F(t,tt_NRZ); subplot(2,2,4)plot(f,10*log10(abs(stf)A2/Ts);axis(-10 10 -20 40); xlabel( 输入信号功率谱密度 );%T2函数function f, sf=T2F(t, st)% This is a function using the FFT function to calculate a signals Fourier% Translation% Input is the time and the signal vectors, the length of time must greater% than 2% Output is the frequency and the signal spectrumdt=t(2)-t(1);T=t(end); df=1/T;N=length(st);f=-N/2*df:df:N/2*df-df; sf=fft(st);sf=T/N*fftshift(sf);end%F2函数function t, st=F2T(f, sf)% This function calculate the time signal using ifft function for the input % signals