基于Matlab的16QAM通信系统的设计与仿真报告

淮海工学院课程设计报告书课程名称：通信系统的计算机仿真设计题目：16QAM通信系统性能分析与MATLAB仿真系（院）：电子工程学院学期：2013-2014-2专业班级：姓名：学号：评语：成绩：签名：日期：课程设计报告专用纸第1页共17页基于Matlab的16QAM通信系统的设计与仿真1绪论1.1研究背景与研究意义应用MATLAB的编程方法和功能模块可以搭建各种仿真系统，还可以应用丰富的时间域、频率域、相位域的仿真测量仪器。许多新一代通信系统的系统级仿真程序出现在MATLAB软件的演示实例中，这使得学习的效率大为提高，对技术与系统的理解已经从概念深入到电路方案和选取层面。Simulink是Mathworks公司推出的基于Matlab平台的著名仿真环境。Simulink作为一种专业和功能强大且操作简单的仿真工具，目前已被越来越多的工程技术人员所青睐，它搭建积木式的建模仿真方式既简单又直观，而且已经在各个领域得到了广泛的应用。QAM（QuadratureAmplitudeModulation）：正交振幅调制。正交振幅调制，这是近年来被国际上移动通信技术专家十分重视的一种信号调制方式。QAM是数字信号的一种调制方式，在调制过程中，同时以载波信号的幅度和相位来代表不同的数字比特编码，把多进制与正交载波技术结合起来，进一步提高频带利用率。正交调幅是一种将两种调幅信号汇合到一个信道的方法，因此会双倍扩展有效带宽。正交调幅被用于脉冲调幅，特别是在无线网络应用。1.2课程设计的目的和任务随着现代通信技术的发展，特别是移动通信技术高速发展，频带利用率问题越来越被人们关注。在频谱资源非常有限的今天，传统通信系统的容量已经不能满足当前用户的要求。正交幅度调制QAM(QuadratureAmplitudeModulation)以其高频谱利用率、高功率谱密度等优势，成为宽带无线接入和无线视频通信的重要技术方案。首先介绍了QAM调制解调原理，提出了一种基于MATLAB的16QAM系统调制解调方案，包括串并转换，2-4电平转换，抽样判决，4-2电平转换和并串转换子系统的设计，对16QAM的星座图和调制解调进行了仿真，并对系统性能进行了分析，进而证明16QAM调制技术的优越性。课程设计的任务是：（1）掌握一般通信系统设计的过程，步骤，要求，工作内容及设计方法：掌握用计算机仿真通信系统的方法。（2）训练学生网络设计能力。（3）训练学生综合运用专业知识的能力，提高学生进行通信工程设计的能力。课程设计报告专用纸第2页共17页216QAM通信系统2.116QAM通信系统基本模型通信系统框图如图2.1所示:图2.116QAM通信系统基本模型信号源：模拟的正弦波语音信号4KHZ。Sample：抽样器，对模拟信号进行抽样，抽样频率8KHZ。A-law：量化器，A-law十三折线法。PCM：编码器，将量化后的信号进行PCM编码，变成一个传速为64Kbit/s的数字信号。信道编码：可以选择分组码，卷积码，汉明码。调制：QAM调制方法。信道：经过调制以后，通过信道。信道可以选择高斯加性白噪声信道，二进制对称信道，多径瑞利衰落信道，莱斯衰落信道等。设置不同的信道信噪比，对系统进行仿真，分析不同信噪比情况下的系统性能。解调：根据调制方式，选择对应解调方式。译码：根据信道编码方式，选择对应的信道解码方式。性能分析：信号经过调制，信道，解调过程。在接收端，将得到的数据与原始信号源数据比较，得到在特定信噪比下的误码率。改变系统信噪比，到系统的误码率曲线图，并给出各关键点信号图及星座图。2.216QAM通信系统的性能指标数字通信系统的有效性可用传输速率来衡量,可靠性可用差错率来衡量。2.2.1有效性码元传输速率RB简称传码率，又称符号速率等。它表示单位时间内传输码元的数目，单位是波特（Baud），记为B。RB=1/T(B)=64KB信号源抽样器量化器编码信道编码调制信道解调译码性能分析噪声课程设计报告专用纸第3页共17页信息传输速率Rb简称传信率，又称比特率等。它表示单位时间内传递的平均信息量或比特数，单位是比特/秒，可记为bit/s，或b/s，或bps。由于传输的是二进制码元，因此码元传输速率与信息传输速率相等，即：Rb=RB=64Kbit/s2.2.2可靠性误码率（码元差错率）Pe是指发生差错的码元数在传输总码元数中所占的比例，更确切地说，误码率是码元在传输系统中被传错的概率，即Pe=错误码元数/传输总码元数误信率（信息差错率）Pb是指发生差错的比特数在传输总比特数中所占的比例，即Pb=错误比特数/传输总比特数在二进制中有Pe，具体数值将在仿真中体现。316QAM通信系统主要模块3.1信源/信宿及其编译码3.1.1信源编码信号源是模拟的4KHz正弦语音信号，sample抽样器，对模拟信号进行抽样，抽样频率是8KHz。A-law量化编码器是采用十三折线法对其进行量化，然后将量化后的信号纪念性PCM编码，变成一个1个传输速率为64Kbit/s的数字信号。3.1.2信宿译码A-law量化译码器是采用十三折线法进行相应的量化译码，信宿可以是scope示波器，也可以是其他的信号输出。3.2基带信号处理本次试验我们采用卷积码和汉明码。汉明码属于线性分组码。汉明码的抗干扰能力较强。它能够纠正单个随机错误。设汉明码编码器的输入信号的长度为k，输出信号的长度为n,则产生的是一个（n,k）汉明码，其中n=2m-1,m3,并且满足n=k+m。MATLAB提供了一个函数”“gfrimfd(m,'min')”用来找到一个最小的本原多项式。卷积码是将发送的信息序列通过一个线性的、有状态的移位寄存器的输出完全取决于这段时间中的输入信息；而在卷积码中，任何一段规定时间内产生的n个码元不仅取决于这段时间内的k个信息位，而且取决于前N-1段时间内的信息位，这个N称为卷积码的约束长度。3.3调制/解调正交振幅键控是一种将两种调幅信号（2ask和2psk）汇合到一个信道的方课程设计报告专用纸第4页共17页法，因此会双倍扩展有效带宽。正交调幅被用于脉冲调幅，特别是在无线网络应用。如图3-3所示的是16QAM所谓调制与解调的示意图。正交调幅信号有两个相同频率的载波，但是相位相差90度（四分之一周期）。一个信号叫I信号，另一个信号叫Q信号。从数学角度将一个信号可以表示成正弦，另一个表示成余弦。两种被调制的载波在发射时已被混和。到达目的地后，载波被分离，数据被分别提取然后和原始调制信息相混和。QAM是用两路独立的基带信号对两个相互正交的同频载波进行抑制载波双边带调幅，利用这种已调信号的频谱在同一带宽内的正交性，实现两路并行的数字信息的传输。这里的正交是指两路载波信号之间的相互正交，即满足：220sincosTcTcttcc3.4信道信道可以选择高斯加性白噪声、二进制对称信道、多径瑞利（Rayleigh）衰落信道、赖斯（Racian）衰落信道等。本次课程设计采用的是加性高斯白噪声信道是信号传输中的最基本的一种信道，加性高斯白噪声（AWGN）信道用高斯分布的噪声信号叠加通过它的信号上，使通过该信道的信号产生与噪声均值相应的偏移，并且围绕平均值做随机波动。在该信道中，当噪声均值为零时，方差表现为噪声的功率大小。加性高斯白噪声模块其采样时间从输入信号中继承。当输入为实信号时，该模块将实高斯白噪声加入该实输入信号，并产生实输出信号；同理，当输入为复信号时，该模块将复高斯白噪声信号加入该复输入信号，并产生复输出信号。模块参数设置中，Initialseed（初始化种子）参数用于初始化噪声的产生，该参数可以是标量，也可以是与输入信号信道数量相匹配的向量形式。4MATLAB对16QAM通信系统的仿真4.1MATLAB主要模块及参数设置4.1.1信源/信宿及其编译码信源编码如图4.1所示:课程设计报告专用纸第5页共17页图4.1信源编码的MATLAB模块Pulsegenerator参数如表4-1所示：表4-1Pulsegenerator参数表Signalgenerator参数如表4-2所示：表4-2Signalgenerator参数表A-LawCompressor参数如表4-3所示：表4-3A-LawCompressor参数表设置项目参数设置PulsetypeTimebasedTimeUsesimulationtimeAmplitude1Period0.000125PulseWidth(%ofperiod)50Phasedelay0设置项目设置参数WaveformsineTime(t)UsesimulationtimeAmplitude1Frequency4000UnitsHertz设置项目设置参数Avalue87.6Peaksignalmagnitude1课程设计报告专用纸第6页共17页Quantizer参数如表4-4所示：表4-4Quantizer参数表设置项目设置参数Quantizationinterval1Sampletime(-1forinherited)0.00001信宿解码如图4.2所示:图4.2信宿解码的MATALAB模块A-LawExpander参数如表4-5所示：表4-5A-LawExpander参数表设置项目设置参数Avalue87.6Peaksignalmagnitude14.1.2基带信号处理汉明码编码译码模块如图4.3和图4.4所示：图4.3汉明码编码模块图4.4汉明码解码模块汉明码模块参数如表4-6所示：表4-6汉明码模块参数表设置项目设置参数CodewordlengthN7MessagelengthKGfprimfd(3,'min')课程设计报告专用纸第7页共17页卷积码编码译码模块如图4.5和图4.6所示：图4.5卷积码编码图4.6卷积码译码卷积码编码模块参数如表4-7所示：表4-7卷积码编码模块码参数表设置项目设置参数Trellisstructurepoly2trellis(7,[171133])ResetNoneViterbiDecoder参数如表4-8所示：表4-8ViterbiDecoder参数表设置项目设置参数Trellisstructurepoly2trellis(7,[171133])DecisiontypeUnquantizedNumberofsoftdecisionbits4Tracebackdepth16OperationmodeContinous4.1.3调制/解调调制解调模块如图4.7所示：图4.7调制/解调的MATLAB模块调制解调模块参数设置如表4-9所示：表4-9调制解调模块参数设置表设置项目设置参数M-arynumber16InputtypebitConstellationorderingBinaryNormalizationmethodAveragepower课程设计报告专用纸第8页共17页Averagepower10Phaseoffset(rad)piSamplespersymbol44.1.4信道信道模块参数如表4-10所示：表4-10信道模块参数表设置项目设置参数Initialseed67ModeSignaltonoiseratio(SNR)SNR(dB)10Inputsignalpower(watts)14.216QAM通信系统的仿真图和结果分析无信道编码总图4.8如下所示：图4.8无信道编码总图无信道编码仿真波形图4.9如下图所示：课程设计报告专用纸第9页共17页图4.9无信道编码总图仿真图卷积编码总图如图4.10所示：图4.10卷积码总图卷积码仿真波形图4.11如下所示：课程设计报告专用纸第10页共17页图4.11卷积码仿真波形图卷积码功率谱图4.12如下所示：图4.12卷积码功率谱图卷积码星座图4.13如下图所示：课程设计报告专用纸第11页共17页图4.13卷积码星座图卷积码眼图4.14如下所示：图4.14卷积码眼图汉明码总图4.15如下所示：课程设计报告专用纸第12页共17页图4.15汉明码总图汉明码仿真波形图4.16如下所示：图4.16汉明