基于matlab的通信系统的与仿真.doc

毕业论文基于MATLAB的通信系统的设计与仿真摘 要 通信系统是用以完成信息传输过程的技术系统的总称，作用是将信息从信源发送到一个或多个目的地。调制与解调在信息的传输过程中占据着重要的地位，是不可或缺的，因此研究系统的调制和解调过程就极为重要。MATLAB是集数值计算、图形绘制、图像处理及系统仿真等强大功能于一体的科学计算语言，它强大的矩阵运算和图形可视化的功能以及丰富的工具箱，为通信系统的调制和解调过程的分析提供了极大的方便。本论文首先介绍了通信系统的概念，进而引出调制和解调，然后介绍了我们常用的几种调制和解调的方法。由于MATLAB具有的强大功能所以详细介绍了MATLAB通信系统工具箱，并给出了基于MATLAB的通信系统的调制与解调的实现，运用MATLAB仿真软件进行仿真。【关键词】： 通信系统 调制与解调 MATLAB - 2 -ABSTRACTCommunication system is used to complete the information transmission process of the floorboard of the technology system, is used to send information from the source to one or more destinations. Modulation and demodulation in the process of information transmission occupies the important position, is indispensable, so the system of modulation and demodulation process is extremely important. MATLAB numerical computation, graphics rendering, image processing and system simulation, and powerful functions in the integration of scientific computing language, its powerful function of matrix and graph visualization and rich toolkit, for the analysis of the communication system of modulation and demodulation process provides a great convenience.This paper firstly introduces the concept of communication system, which lead to modulation and demodulation, and then introduces our commonly used several kinds of modulation and demodulation method.【Key words】：Communication system Modulation and demodulation MATLAB目 录前 言4第一章 通信系统的调制与解调5第一节 通信系统的概念5一、模拟通信系统5二、数字通信系统6三、QAM调制解调原理14第二节 多进制正交幅度调制及相干解调原理框图17一、正交调制原理框图17二、相干解调原理框图18第二章 MATLAB通信系统工具箱19第三章 MATLAB的模拟调制和解调实例21第一节 用MATLAB分析双边带幅度调制（DSM-AM）21第二节 用 MATLAB 分析相干解调过程24第三节 基于MATLAB的多进制正交幅度调制及相干解调设计与仿真26一、对系统进行分析与设计26二、随机信号的生成26三、星座图映射27四、波形成形（平方根升余弦滤波器）29五、调制30六、加入高斯白噪声之后解调31七、误码率曲线33八、16-QAM载波调制信号在AWGN信道下的性能33第四节 仿真结果及分析35结 论39致 谢40参考文献41前 言通信按照传统的理解就是信息的传输，在当今高度信息化得社会，信息和通信已经成为现代社会的“命脉”。信息作为一种资源，只有通过广泛的传播与交流，才能产生利用价值，促进社会成员之间的合作，推动社会生产力的发展，创造出巨大的经济效益。而通信作为传输信息的手段或方式，与传感技术、计算机技术相互融合，已经成为21世纪国际社会和世界经济发展的强大推动力，所以未来的通信对人们的生活方式和社会的发展将会产生更加重大和意义深远的影响。在信息传输过程中，要求天线的尺寸要和信号的波长相比拟，信号才能有效的被辐射。对于语音信号来说，相应的天线尺寸要在几十公里以上，实际上不可能实现，所以需要经过调制将信号频谱搬移到较高的频率范围，如果不进行调制就把信号直接辐射出去，那么各电台所发出信号的频率就会相同。调制作用的实质就是使相同频率范围的信号分别依托于不同频率的载波上，接收机就可以分离出所需的频率信号，不致互相干扰。有时信号过于复杂，人工计算其调制和解调过程较难实现，对其结果的分析又缺乏可视化的直观表现，影响了所得结果在实际生活中的应用，美国MathWorks公司开发的MATLAB解决了这一问题。它应用于自动控制、数学计算、信号分析、信号处理等诸多领域,也是国内高校和研究部门进行许多科学研究的重要工具。MATLAB的出现给通信系统的分析提供了极大的方便。42第一章 通信系统的调制与解调第一节 通信系统的概念通信是为了传输信息，通信系统就是将信息从信源发送到一个或多个目的地，对于电通信来说，首先要把消息转变成电信号，然后经过发送设备，将信号送入信道，在接收端利用接受设备对接收信号作相应的处理后，送给信宿再转换为原来的消息，这一过程可利用图1.1所示的通信系统一般模型来概括。接收设备信道发送设备 （发送端） （接收端） 信息源受信者噪声源图1.1 通信系统的一般模型通信(Communication)传输的消息是多种多样的，可以分成两大类：一类称为连续消息；另一类称为离散消息。消息的传递是通过它的物理载体电信号来实现的，按信号参量的取值不同，可以把信号分为两类：模拟信号和数字信号。通常按照信道中传输的是模拟信号还是数字信号，相应的可以把通信系统分为模拟通信系统和数字通信系统。一、模拟通信系统模拟通信系统是利用模拟信号来传递信息的通信系统，其模型如图1.2所示，其中包含两种重要变换。第一种变换是，在发送端把连续消息变换成原始电信号，在接收端进行相反的变换，这种变换由信源和信宿来完成，通常称为原始电信号为基带信号，基带的意思是指信号的频谱从零频附近开始。有些信道可以直接传输基带信号，而以自由空间作为信道的无线电传输却无法直接传输这些信号。因此，模拟通信系统中常常需要进行第二种变换：把基带信号变换成适合在信道中传输的信号，并在接收端进行反变换。完成这种变换和反变换的通常是调制器和解调器。噪声源调制器模拟信息源信道受信者解调器图1.2 模拟通信系统模型二、数字通信系统数字通信系统（Digital Communication System,DCS）是利用数字信号来传递信息的通信系统，如图1.3所示。数字通信所涉及的技术问题很多，其中主要有信源编码与译码、信道编码与译码、数字调制与解调、同步以及加密与解密等。数字调制信道编码加密信源编码信息源噪声源信道受信者信源译码信道译码解密数字解调图1.3 数字通信系统模型 与模拟通信相比，数字通信具有许多优良的特性：（1）、抗干扰能力强，且噪声不积累。数字通信系统中传输的是离散取值的数字波形，接收端的目标不是精确的还原被传输的波形，而是从收到噪声干扰的信号中判决出发送端所发送的是哪一个波形。（2）、传输差错可控。在数字同喜系统中，可以通过信道编码技术进行检错与纠错，降低误码率，提高传输质量。（3）、便于用现代数字信号处理技术对数字信息进行处理、变换、存储。这种数字处理的灵活性表现为可以将来自不同信源的信号综合到一起传输。（4）、易于集成，使通信设备微型化，重量轻。（5）、易于加密处理，且保密性好。 数字通信的缺点是，一般需要较大的传输带宽，设备复杂。近年来，随着大规模集成电路的出现，数字系统的设备复杂程度和技术难度大大降低了，数字传输方式日益受到欢迎。数字处理的灵活性使得数字传输系统中传输的数字信息既可以来自计算机、电传机等数据终端的各种数字代码，也可以来自模拟信号经过数字化处理后的脉冲编码（PCM）信号等。原理上，数字信息可以直接用数字代码序列表示和传输，但在世纪传输中，视系统的要求和信道情况，一般需要进行不同形式的编码，并且选用一组取值有限的离散波形来表示，这些取值离散的波形可以是未经调制的电信号，也可以是调制后的信号。未经调制的数字信号所占据的频谱是从零频或很低频率开始，称为数字基带信号。1、调制和解调的概念 调制：把信号转换成适合在信道中传输的形式的一种过程，广义的调制分为基带调制和带通调制（也称为载波调制）。在无线通信中和其他大多数场合，调制均指载波调制。载波调制，就是用调制信号去控制载波的参数的过程，使载波的某一个或某几个参数按照调制信号的规律而变化。调制信号是指来自信源的消息信号（基带信号），这些信号可以是模拟的，也可以是数字的。未受调制的周期性振荡信号称为载波，它可以是正弦波，也可以是非正弦波（如周期性脉冲序列）。载波调制后称为已调信号，它含有调制信号的全部特征。解调：将已调信号中的调制信号恢复出来，是调制的逆过程。调制方式不同，解调方法也不一样。解调可以分为正弦波解调（有时也称为连续波解调）。正弦波解调还可再分为幅度解调、频率解调和相位解调，此外还有一些变种如单边带信号解调、残留边带信号解调等。同样，脉冲波解调也可分为脉冲幅度解调、脉冲相位解调、脉冲宽度解调和脉冲编码解调等。调制方式有很多，根据调制信号是模拟信号还是数字信号，载波是连续波（通常是正弦波）还是脉冲序列，相应的调制方式有模拟连续波调制（简称模拟调制）、数字连续波调制（简称数字调制）、模拟脉冲调制和数字脉冲调制等。 最重要和最常用的模拟调制方式是用正弦波作为载波的幅度调制和角度调制。常见的调幅、双边带、单边带和残留边带等调制就是幅度调制的几个典型实例。解调的方法可以分为两类：相干解调和非相干解调（包络检波）。解调过程与采用何种解调方式有关，对于常规幅度调制，一般用包络检波进行解调，由于在这种解调方式中，接收机对载波频率和相应精度的了解是无关紧要的，所以解调过程相对简单。对于DSB-AM调制和SSB-AM调制，用相干解调的方法，它要求在接收机中有一个与载波同频同相的信号，接收机中产生所需要的正弦波振荡器，为本地振荡器。数字通信系统中，发送信号某一参数的离散取值（如正弦波的振幅、频率、相位等）于所承载消息的数字信号（即离散符号）之间通常是一种一一对应关系。如同模拟信号的频带传输时一样，传输数字信号时也有三种基本的调制方式：振幅键控（ASK）、频移键控（FSK）和相移键控（PSK）。它们分别对应于利用载波（正弦波）的振幅、频率和相位来承载数字基带信号，可以是模拟线性调制和角度调制的特殊情况。理论上数字调制与模拟调制在本质上没有什么不同，它们都属于正弦波调制。但是，数字调制是源信号为离散型的正弦波调制，而模拟调制则是源信号为连续型的正弦波调制，因而，数字调制具有由数字信号带来的一些特点。这些特点主要包括两个方面：第一，数字调制信号的产生，除把数字的调制信号当做模拟信号的特例而直接采用模拟调制方式产生数字调制信号外，还可以采用键控载波的方法。第二，对于数字调制信号的解调，为提高系统的抗噪声性能，通常采用与模拟调制系统中不同的解调方式。数字信息由二进制和多进制之分，因此，数字调制可分为二进制调制和多进制调制。在二进制调制中，信号参量只有两种可能取值；而在多进制调制中，信号参量可能有M（M2）中取值。2、幅度调制和相干解调的原理1）幅度调制（线性调制）的原理幅度调制是由调制信号去控制高频载波的幅度，使之随调制信号作线性变化的过程，设正弦型载波为 c(t)=Acos() (2-1)式中：A为载波幅度；为载波角频率；为载波初始相位（以后可以假定为为0，而不失讨论的一般性）。根据调制定义，幅度调制信号（已调信号）一般可表示成 Sm(t)=Am(t)cos () (2-2)式中：m(t)为基带调制信号。设调制信号m(t)的频谱为,则由上式（2.3-2）不难得到已调信号Sm(t)的频谱 = (2-3)由以上表示式可见，在波形上，幅度已调信号的幅度随基带信号的规律而呈正比地变化；在频谱结构上，它的频谱完全是基信号频谱在频域内的简单搬移（精确到常数因子）。由于这种搬移是线性的，因此，幅度调制通常又称为线性调制，这里的“线性”并不意味着已调信号与调制信号之间符合线性变换关系，实际上，任何调制过程都是一种非线性的变换过程。标准调幅就是常规双边带调制，简称调幅（AM）。假设调制信号m(t)的平均值为0，将其叠加一个直流偏量A0后与载波相乘（如图1.4所示），即可形成调幅信号。m(t)Sm(t) A0 图 1.4 AM调制模型其时域表示式为SAM(t)=A0+m(t) =A0 +m(t) (2-4)式中：A0为外加的直流分量；m(t)可以是确知信号，也可以是随机信号。若m(t)是确知信号，则AM信号的频谱为SAM()=A0+1/2(2-5)2) 相干解调相干解调也叫同步检波，解调与调制的实质一样，均是频谱搬移。调制是把基带信号的频谱搬到了载波位置，这一过程可以通过一个相乘器与载波相乘来实现。解调则是调制的反过程，即把在载频位置的已调信号的频谱搬回到原始基带位置，因此同样可以用相乘器与载波相乘来实现，相干解调的一般模型如图1.5所示。LPF Sm(t) Sp(t) Sd(t) c(t)= 图 1.5 相干解调器的一般模型相干解调时，为了无失真的恢复原基带信号，接收端必须提供一个与接收的已调载波严格同步（同频同相）的本地载波（称为相干载波），它与接收的已调信号相乘后，经低通滤波器取出低频分量，即可得到原始的基带调制信号。相干解调器适用于所有的线性调制信号的解调。送入解调器的已调信号的一般表达式为 Sm(t)=SI(t) +SQ(t) (2-6)与同频同相的相干载波c(t)相乘后，得 Sp(t)= Sm(t) =1/2 SI(t)+1/2 SI(t) +1/2 SQ(t) (2-7)经低通滤波器（LPF）后，得 Sd(t)=1/2 SI(t) （2-8）其中，SI(t)是m(t)通过一全通滤波器H1（）后的结果。因此，Sd(t)就是解调输出，即 Sd(t)=1/2 SI(t) m(t)由此可见，相干解调器适用于所有线性调制信号的解调，即对于AM、DSB、SSB和VSB都是适用的。只是AM信号的解调结果中含有直流成分A0，这时在解调后加上一个简单的隔直流电容即可。3、正交幅度调制解调原理1) 正交幅度调制技术正交振幅调制（Quadrature Amplitude Modulation,QAM）是一种振幅和相位联合键控。虽然MPSK和MDPSK等相移键控的带宽和功率方面都具有优势，即带宽占用小和比特噪声比要求低。但是由图1.6可见，在MPSK体制中，随着 图1.6 8PSK信号相位M的增大，相邻相位的距离逐渐减小，使噪声容限随之减小，误码率难于保证。为了改善在M大时的噪声容限，发展出了QAM体制。在QAM体制中，信号的振幅和相位作为两个独立的参量同时受到调制。这种信号的一个码元可以表示为 （21）式中：k=整数；和分别可以取多个离散值。 式（21）可以展开为 （22）令 Xk = Akcosqk， Yk = -Aksinqk则式（21）变为 （23）和也是可以取多个离散的变量。从式（23）看出，可以看作是两个正交的振幅键控信号之和。在式（21）中，若qk值仅可以取p/4和-p/4，Ak值仅可以取+A和-A，则此QAM信号就成为QPSK信号，如图1.7所示：图1.7 4QAM信号矢量图所以，QPSK信号就是一种最简单的QAM信号。有代表性的QAM信号是16进制的，记为16QAM，它的矢量图示于下图1.8中：Ak图1.8 16QAM信号矢量图图中用黑点表示每个码元的位置，并且示出它是由两个正交矢量合成的。类似地，有64QAM和256QAM等QAM信号，如图1.9、图1.10所示。它们总称为MQAM调制。由于从其矢量图看像是星座，故又称星座调制。 图1.9 64QAM信号矢量图 图1.10 256QAM信号矢量图16QAM信号的产生方法主要有两种。第一种是正交调幅法，即用两路独立的正交4ASK信号叠加，形成16QAM信号，如图1.11所示。AM图1.11 正交调幅法第二种方法是复合相移法，它用两路独立的QPSK信号叠加，形成16QAM信号，如图1.12所示。图中 AMAM图1.12 复合相移法虚线大圆上的4个大黑点表示一个QPSK信号矢量的位置。在这4个位置上可以叠加上第二个QPSK矢量，后者的位置用虚线小圆上的4个小黑点表示。三、QAM调制解调原理 1、QAM调制正交幅度调制QAM是数字通信中一种经常利用的数字调制技术，尤其是多进制QAM具有很高的频带利用率，在通信业务日益增多使得频带利用率成为主要矛盾的情况下，正交幅度调制方式是一种比较好的选择。正交幅度调制（QAM）信号采用了两个正交载波，每一个载波都被一个独立的信息比特序列所调制。发送信号波形如图1.13所示 图1.13 M=16QAM信号星座图式中和是电平集合，这些电平是通过将k比特序列映射为信号振幅而获得的。例如一个16位正交幅度调制信号的星座图如下图所示，该星座是通过用M4PAM信号对每个正交载波进行振幅调制得到的。利用PAM分别调制两个正交载波可得到矩形信号星座。QAM 可以看成是振幅调制和相位调制的结合。因此发送的QAM信号波形可表示为 如果那么QAM方法就可以达到以符号速率同时发送个二进制数据。图1.14给出了QAM调制器的框图。图1.14 QAM调制器框图2、QAM的解调和判决假设在信号传输中存在载波相位偏移和加性高斯噪声。因此r(t)可以表示为其中是载波相位偏移，且将接收信号与下述两个相移函数进行相关如图1.15所示，相关器的输出抽样后输入判决器。使用图1-15中所示的锁相环估算接收信号的载波相位偏移，相移和对该相位偏移进行补偿。图1.15 QAM信号的解调和判决假设图中1.15所示的时钟与接收信号同步，以使相关器的输出在适当的时刻及时被抽样。在这些条件下两个相关器的输出分别为其中噪声分量是均值为0，方差为的互不相关的高斯随机变量。最佳判决器计算距离量度 3、QAM的误码率性能矩形QAM信号星座最突出的优点就是容易产生PAM信号可直接加到两个正交载波相位上，此外它们还便于解调。对于下的矩形信号星座图（K为偶数），QAM信号星座图与正交载波上的两个PAM信号是等价的，这两个信号中的每一个上都有个信号点。因为相位正交分量上的信号能被相干判决极好的分离，所以易于通过PAM的误码率确定QAM的误码率。M进制QAM系统正确判决的概率是：。式中是进制PAM系统的误码率，该PAM系统具有等价QAM系统的每一个正交信号中的一半平均功率。通过适当调整M进制PAM系统的误码率，可得：式中是每个符号的平均信噪比。因此，因此M进制QAM的误码率为： 可以注意到，当K为偶数时，这个结果对情形时精确的，而当K为奇数时，就找不到等价的进制PAM系统。如果使用最佳距离量度进行判决的最佳判决器，可以求出任意K=1误码率的的严格上限。=其中，是每比特的平均信噪比。第二节 多进制正交幅度调制及相干解调原理框图一、正交调制原理框图图1.16 正交调制原理框图二、相干解调原理框图图1.17 相干解调原理框图第二章 MATLAB通信系统工具箱MATLAB通信系统工具箱主要应用于参数化模型，频谱分析和估计等。通信系统工具箱实际上是用MATLAB的基本语句变成的各种子程序集，用于解决某一方面的专门问题或实现某一类的新算法。MATLAB 的通信系统工具箱可以任意增减，不同的工具箱可以给不同领域的用户提供丰富强大的功能。MATLAB中的通信系统工具箱是目前比较成熟的通信系统仿真工具， MATLAB通信系统工具箱提供的函数主要用于处理通信及传输问题，为我们分析通信系统中每个过程提供了极大的帮助，使我们能够轻松直观的解决问题。常用的MATLAB函数主要有以下几项：1、figure功能：创建新的图形窗口(用于输出图形的窗口)。格式：figure说明：figure函数创建一个新的图形窗口, 并成为当前图形窗口, 所创建的图形窗口的序号是按同一MATLAB程序中创建的顺序号。2、plot功能：线型绘图函数。格式：plot(x)plot(x,y)说明：plot(x)是一种最简单的调用方式, x是长度为n的数值向量。plot(x)的作用是在坐标系中顺序地用直接连接顶点i,x(i),i=1,2,n生成一条折线。当向量元素充分多时, 即可生成一条光滑的曲线。3、subplot功能：多坐标设置与定位当前坐标系。格式：subplot(m,n,k)说明：subplot(m,n,k)将图形窗口分成m行n列m*n块子区域，按从上到下,从左到右的顺序,在第k块子区域定义一个坐标系, 使其成为当前坐标系, 随后的绘图函数将在该坐标系输出图形。4、ezplot功能：该命令用来绘制符号表达式的自变量和对应各函数值的二维曲线。格式： ezplot(F, xmin,xmax,fig)说明：其中F是要画的符号函数；xmin,xmax是绘图的自变量范围，fip是窗口。5、axis功能：该命令用来控制坐标轴的特性。格式： axis(xmin,xmax, ymin,ymax)说明：此为坐标范围,其中xminxmax ，yminymax 。7、angle功能：求复数的相角。格式：P=angle(Z)说明：当Z为复数矩阵时，用来求矩阵Z中每个元素的相角，相角位于-,。第三章 MATLAB的模拟调制和解调实例第一节 用MATLAB分析双边带幅度调制（DSM-AM）我们可以利用 MATLAB 强大的符号运算功能来进行运算,再根据MATLAB的可视化结果进行分析。在DSB-AM中，已调信号的时域表示为： u(t)=m(t)c(t)=Acm(t) 式中，m(t)是消息信号，c(t)=Ac 为载波，fc是载波的频率（单位：HZ）, 是初始相位。为了讨论方便取初相=0（以下类似）。随u(t)作傅里叶变换，即可得到信号的频域表示： U（f）=Ac/2 +Ac/2 传输带宽Bt是消息信号带宽W的两倍，即：Bt=2W。 0tt0/3 某消息信号 m(t)= -2 t0/3t2t0/3 其他用信号m(t)以DSB-AM方式调制载波c(t)=cos(2fct),所得到的已调制信号记为u(t).设t。=0.15s,fc=250Hz。试比较消息信号与已调信号，并绘制它们的频谱。运用如下 MATLAB程序: t0=0.15; 信号持续时间 ts=0.001; 采样时间间隔 Fc=250; 载波频率 FS=1/ts; 采样频率 df= 0.3; 频率分辨率 t=0:ts:t0; 时间矢量 m=ones(1,t0/(3*ts), -2*ones(1,t0/(3*ts),zeros(1,t0/(3*ts)+1); 定义信号序列 C=cos(2*pi*Fc.*t); 载波信号 u=m.*c; 调制信号 M,m,df1=fft_seq(m,ts,df); 傅里叶变换 M=M/FS; U,u,df1=fft_seq(m,ts,df); U=U/FS; C,c,df1=fft_seq(c,ts,df); f=0:df1:df1*(length(m)-1)-Fs/2; 频率矢量 subplot(2,2,1);plot(t,m(1:length(t); 未调制信号 title(未调制信号)； subplot(2,2,2);plot(t,u(1:length(t); 已调制信号 title(已调制信号)； subplot(2,2,3)；plot(f,abs(fftshift(M); 未调制信号频谱 title(未调制信号频谱)； subplot（2,2,4）；plot(f,abs(fftshift(U); 已调制信号频谱 title(已调制信号频谱)；傅里叶变换函数fft_seq,源代码如下： functionM,m,df=fft_seq(m,ts,df) M,m,df=fft_seq(m,ts,df) M,m,df=fft_seq(m,ts) M为输入序列m的傅里叶变换，ts为抽样间隔，输入df为频率分辨率 输出序列m按要求的频率分辨率df进行补零后的序列 输出df为最终的频率分辨率 fs=1/ts; if nargin=2,n1=0; else, n1=fs/df; end n2=length(m); n=2(max(nextpow2(n1),nexpow2(n2); M=fft(m,n); m=m,zeros(n-n2); df=fs/n;运行后得到的信号和调制信号如图3.1所示：图3.1 双边带幅度调制信号第二节 用 MATLAB 分析相干解调过程 0tt0/3例 某消息信号m (t)=-2t0/3t2t0/3用信号m(t)以DSB-AM方式调制载波c(t)=cos(2fct),所得到的已调制信号记为u(t).设t0=0.15s,fc=250Hz。得到DSM-AM的调制信号，现试对该调制信号进行相干解调，并绘出消息信号的时域频域曲线。MATLAB设计程序如下：t0=0.15; 信号持续时间 ts=1/1500; 采样时间间隔 Fc=250; 载波频率 FS=1/ts; 采样频率 df= 0.3; 频率分辨率 t=0:ts:t0; 时间矢量 m=ones(1,t0/(3*ts), -2*ones(1,t0/(3*ts),zeros(1,t0/(3*ts)+1); 定义信号序列 C=cos(2*pi*Fc.*t); 载波信号u=m.*c; 调制信号y=u.*c; 混频 M,m,df1=fft_seq(m,ts,df); 傅里叶变换 M=M/Fs; 缩放 U,u,df1=fft_seq(m,ts,df); U=U/FS; Y,y,df1=fft_seq(y,ts,df); Y=Y/Fs; f_cutoff=150; 滤波器的截止频率 n_cutoff=floor(150/df1); 设计滤波器 f=0:df1:df1*(length(m)-1)-Fs/2; 频率矢量 H=zero(size(f); H(1:n_cutoff)=2*ones(1,n_cutoff); H(length(f)-n_cutoff+1:length(f)=2*ones(1,n_cutoff); DEM=H.*Y; 滤波器输出的频谱 dem=real(ifft(DEM)*Fs; 滤波器的输出 subplot(2,2,1);plot(t,m(1:length(t); 未调制信号 title(未调制信号)； subplot(2,2,2);plot(t,dem(1:length(t); 解调信号 title(解调信号)； subplot(2,2,3)；plot(f,abs(fftshift(M); 未调制信号频谱title(未调制信号频谱)； subplot（2,2,4）；plot(f,abs(fftshift(DEM)；解调信号的频谱 title(已调制信号谱)；运行该程序后得到的信号和调制信号，以及信号调制前、后的频谱对比如图3.2所示：图3.2 相干解调信号第三节 基于MATLAB的多进制正交幅度调制及相干解调设计与仿真一、对系统进行分析与设计首先进行系统的分析的设计，整个设计分为如下几个部分：随机序列的产生，序列的串并和并串转换，16QAM调制，星座图的绘制，16QAM解调，加入噪声，误码率的测量及绘图。二、随机信号的生成利用Matlab中的random_binary函数来产生0、1等概分布的随机信号。源代码如下所示： random_binary.m%产生二进制信源随机序列function info=random_binary(N)if nargin = 0, %如果没有输入参数，则指定信息序列为10000个码元 N=10000;e