基于QAM调制的无线衰落信道的能分析与仿真课程设计

1、实践教学*兰州理工大学计算机与通信学院通信系统仿真训练题 目：基于QAM调制的无线衰落信道的性能分析与仿真专业班级：通信工程（1）班摘 要本课程设计主要是为了研究无线衰落信道的特性，通过深入分析了无线衰落信道的衰落特征以及衰落信道对信号传输产生的影响，设计了一种基于QAM调制的无线衰落信道的仿真，主要是针对16QAM调制下的瑞利衰落信道、高斯衰落信道的模型在MATLAB环境下进行模拟仿真，然后观察数字信号经过这两类信道后的时域波形、频谱以及误码率，通过仿真，对比分析无线衰落信道对信号的影响,有助于对无线信道的进一步理解和掌握以及工程实践提供一定的指导。关键词：QAM调制；瑞利衰落；高斯噪声；目

2、录摘要2目录3前言1一、QAM调制解调2调制简介21.2 QAM调制21.3 16QAM调制11.3.1 16QAM调制解调原理11.3.2 16QAM星座图映射31.3.3 16QAM抽样判决31.3.4 16QAM误码率曲线41.3.5 16QAM信号与其它调制信号的性能比较4二、衰落信道62.1 通信系统信道模型及其分类62.1.1 恒参信道62.1.2 随参信道7瑞利衰落信道的统计模型82.3 通信系统中噪声概述及高斯白噪声92.3.1 噪声概述92.3.2 高斯白噪声9通信系统的主要性能11三、无线衰落信道仿真程序设计12仿真设备12系统整体框图123.3 QAM调制解调13瑞利信道

3、仿真模块13四、设计与仿真144.1 16QAM调制模块的建立与仿真144.1.1 信号源14串并转换14 2-4电平转换154.1.4 增加载波164.1.5 调制信号形成174.2 16QAM调制信号的噪声叠加184.3 16QAM解调模块的建立与仿真194.3.1 滤波器194.3.2 4-2电平转换204.3.3 并串转换204.4 QAM性能分析224.4.1 16QAM抗噪声性能仿真22总结23参考文献24附录25致谢40前言无线通信是当今社会最重要的通信方式之一。在进行无线通信系统的设计时，首先需要考虑的是信道的传输特性，因此无线衰落信道的建模与仿真研究对于无线通信有着重要的意义

4、。针对无线衰落信道已经提出了许多的仿真模型，其中针对宽带短波信道ITS模型和陆地移动卫星信道的多状态模型在两种通信方式下得到了广泛的应用。影响无线通信性能的主要因素有:多径效应、多普勒效应和阴影效应，无线衰落信道进行建模时主要考虑这几方面的影响，并通过现有的数学理论与实测数据尽可能精确的进行曲线拟合，得到想要的信道模型。在卫星移动通信系统、陆地移动通信系统中其电波传播方式主要以视距传播为主。由于多径和接收端运动等因素的影响，使得无线信道对接收信号在时间、频率和角度上造成了色散，这种色散表现在接收信号幅度上就是所谓的信号衰落。因此，多径效应对通信质量有着至关重要的影响，根据不同的无线环境，接收信

5、号包络一般服从几种典型分布，如瑞利分布、莱斯分布等。在本文中，专门针对接收信号包络服从瑞利分布的信道进行建模仿真，为实际的通信系统设计提供理论参考和支持。一、QAM调制解调1.1调制简介调制在通信系统中的作用至关重要。所谓调制，就是把信号转换成适合在信道中传输的形式的一种过程。广义的调制分为基带调制和带通调制（也称载波调制）。载波调制，就是用调制信号去控制载波的参数的过程，即使载波的某一个或某几个参数暗中啊调制信号的规律而变化。调制信号是指来自信源的消息信号（基带信号），这些信号可以是模拟的，也可以是数字的。未受调制的周期性震荡信号称为载波，它可以是正弦波，也可以使非正弦波（如周期性脉冲序列）

6、。载波调制后称为已调信号，它含有调制信号的全部特征。解调（也称检波）则是调制的逆过程，其作用是将已调信号中的调制信号恢复出来。解调的方法可分为两类：相干解调和非相干解调（包络检波）。相干解调时，为了无失真地恢复原基带信号，接收端必须提供一个与接收的已调载波严格同步（同频同相）的本地载波。本课题采用的是相干解调。数字调制具有3种基本方式：数字振幅调制、数字频率调制、数字相位调制，这3种数字调制方式都存在不足之处，如：频谱利用率低、抗多径抗衰弱能力差、功率谱衰减慢、带外辐射严重等。为了改善这些不足，在恒参信道中，提出了正交振幅调制（QAM）方式，它具有高的频谱利用率，因此在卫星通信和有线电视网络高

7、速数据传输等领域得到广泛应用。 QAM调制正交振幅调制（QAM）是一种矢量调制，它是将输入比特先映射（一般采用格雷码）到一个复平面（星座）上，形成复数调制符号。正交调幅信号有两个相同频率的载波，但是相位相差90度（四分之一周期，来自积分术语）。一个信号叫I信号，另一个信号叫Q信号。从数学角度将一个信号可以表示成正弦，另一个表示成余弦。两种被调制的载波在发射时已被混和。到达目的地后，载波被分离，数据被分别提取然后和原始调制信息相混和。 这样与之作幅度调制（AM）相比，其频谱利用率高出一倍。QAM是用两路独立的基带信号对两个相互正交的同频载波进行抑制载波双边带调幅，利用这种已调信号的频谱在同一带宽

8、内的正交性，实现两路并行的数字信息的传输。该调制方式通常有二进制QAM（4QAM）、四进制QAM（l6QAM）、八进制QAM（64QAM）、，对应的空间信号矢量端点分布图称为星座图，分别有4、16、64、个矢量端点。目前QAM最高已达到1024QAM。样点数目越多，其传输效率越高。但并不是样点数目越多越好，随着样点数目的增加，QAM系统的误码率会逐渐增大，所以在对可靠性要求较高的环境，不能使用较多样点数目的QAM。对于4QAM，当两路信号幅度相等时，其产生、解调、性能及相位矢量均与4PSK相同。 a 4QAM星座图b 16QAM星座图 图1-1 QAM星座图QAM采用格雷编码，采用格雷码的好处

9、在于相邻相位所代表的两个比特只有一位不同，由于因相位误差造成错判至相邻相位上的概率最大，故这样编码使之仅造成一个比特误码的概率最大。下图以16QAM为例，显示了编码图1-2 16QAM编16QAM调制1.3.1 16QAM调制解调原理16QAM是两路4ASK信号的叠加，其演变方式可以有以下两种：（1）正交调幅法，由两路独立的正交4ASK信号叠加而成； 图1-3 正交调幅（2）复合相移法，由两路独立的QPSK信号叠加而成。图中虚线大圆上的4个大黑点表示第一个QPSK信号矢量的位置，在这4个位置上可以叠加上第二个QPSK矢量，后者的位置用虚线小圆上的4个小黑点表示。图1-4 复合相移法在QAM体制

10、中，信号的振幅和相位作为两个独立的参量同时受到调制。这种信号的一个码元可以表示为Sk(t)=Akcos(0t+k) kT<t(k+1)T （1-1）式中，k取整数；Ak和qk分别可以取多个离散值。上式可以展开为Sk(t)=Akcoskcos0tAksinksin0t （1-2）令Xk=Akcosk Yk=-Aksink则信号表示式变为Sk(t)= Xkcos0t+Yk sin0t（1-3）Xk和Yk也是可以取多个离散值的变量。从上式看出，k(t)可以看作是两个正交的振幅键控信号之和。本课题采用了正交调幅法。在发送端调制器中串/并变换使得信息速率为Rb的输入二进制信号分成两个速率为Rb/2

11、的二进制信号，2/4电平转换将每个速率为Rb/2的二进制信号变为速率为Rb/8的电平信号，然后分别与两个正交载波相乘，再相加后即得16QAM信号下图是16QAM的调制框图。串并变换2-4电平变换2-4电平变换LPFLPF载波发生器90°相移Icoswct-sinwct已调信号Y基带信号XQ图1-5 16QAM调制框图解调是调制的逆过程，在接收端解调器中可以采用正交的相干解调方法。接受到的信号分两路进入两个正交的载波的相干解调器，再分别进入判决器形成L进制信号并输出二进制信号，最后经并/串变换后得到基带信号。下图为16QAM解调框图：并串变换4-2电平变换4-2电平变换LPFLPF载波

12、恢复90°相移QIcoswct-sinwctEPF抽样判决抽样判决解调输出定时脉冲接收信号图1-6 16 QAM解调框图1.3.2 16QAM星座图映射将等概分布的0、1信号映射到16QAM星座图上。每四个bit构成一个码子，具体实现的方法是，将输入的信号进行串并转换分成两路，分别叫做I路和Q路。再把每一路的信号分别按照两位格雷码的规则进行映射，这样实际上最终得到了四位格雷码。为了清楚说明，参看表1-1表1-1两位格雷码的映射规律两位0、1码映射后（按格雷码）0 0-30 1-11 111 03 16QAM抽样判决经过前边的匹配滤波器解调或者称为相关解调产生了一组向量，在这里就是一个

13、一维的向量，根据最大后验概率（MAP）准则（由于各个信号的先验概率相等，所以页可以认为是最大似然准则），得到了最小距离检测。具体在本仿真系统中，判断为各个信号的门限如表2所示。判决后得到的数据再按照格雷码的规则还原成0、1信号，最终将两路0、1信号合成一路0、1信号，用来同最初的信号一起决定误码率。表1-2判决电平对应表判决前的信号的幅度对应的判决后的幅度31131.3.4 16QAM误码率曲线对于16QAM信号星座图等效为在两个正交载波上的两个PAM信号，其中每一个具有4个信号点。因为在解调器中可以将相位正交的两个信号分量完全分开，所以QAM的错误概率可以由PAM的错误概率求得。16QAM系

14、统的正确判决概率是（1-4）式中，是4元PAM的错误概率，在等效QAM系统的每一个正交信号中，4元PAM具有一半的平均功率，通过适当的修改4元PAM的错误概率，可以得到（1-5）其中是平均符号SNR。因此，16QAM的错误概率是（1-6） 16QAM信号与其它调制信号的性能比较16QAM和16PSK星座图中相邻点欧氏距离直接代表这噪声容限的大小。按最大振幅相等，画出16QAM信号和16PSK信号的星座图。设其最大振幅为AM，则16PSK信号的相邻矢量端点的欧氏距离等于（1-7） 而16QAM信号的相邻点欧氏距离等于（1-8）d2和d1的比值就代表这两种体制的噪声容限之比。AM d2AM d1图

15、1-7 欧氏距离按上两式计算，d2超过d1约1.57 dB。但是，这时是在最大功率（振幅）相等的条件下比较的，没有考虑这两种体制的平均功率差别。16PSK信号的平均功率（振幅）就等于其最大功率（振幅）。而16QAM信号，在等概率出现条件下，可以计算出其最大功率和平均功率之比等于倍，即2.55 dB。因此，在平均功率相等条件下，16QAM比16PSK信号的噪声容限大4.12 dB。二、衰落信道2.1 通信系统信道模型及其分类各种发送信息传送到既定的信宿，可选用适于传输的物理媒体，完成通信功能。连接发信号与收信号设备、适用于不同类型通信业务的各种物理媒体通称为信道。信道可分为有界与无界两大类，即通

16、常所说的有线信道与无线信道。前者如双绞线、电缆、光纤、波导等，后者为自由空间提供的各种频段或波长的电磁波传播信道。根据各种信道不同的特征和参量及其变化情况，又将它们分为恒参信道和随参信道。前者如有线信道、微波与卫星信道等，后者如无线系统的短波和超短波散射信道。收信者一般地，如单指传输媒体而言称为狭义信道。在具体的通信系统构成中，往往把信源发出的模拟信号和数字编码基带信号视为信息部分，从调制器到接收端解调器这一中间变换历程中。经过了包括物理媒体在内的线路设备(如交换、放大、中继等中间部件)传输路径，因此将图3-1所表示的调制信道和编码信道称为广义信道。编码器调制器线路设备译码器信源解调器调制信道

17、编码信道图2-1 信道结构图2 恒参信道恒参信道是指由架空明线、电缆、中长波地波传输，超短波及微波视距传输，人造卫星中继，光导纤维以及光波视距传输等传输媒体构成的信道。恒参信道以有线信道为最典型，其特征参数主要是频率特征，如幅度频率特征与相位频率特征及频率漂移等。反映在时域，如信道时延、抖动，尚有电平波动和非线性等。其中，幅度频率特性，就理想而言，可表示为理想传输函数。2 随参信道随参信道包括短波电离层反射、超短波流星余迹散射，超短波及微波对流层散射，超短波电离层散射以及超短波超视距绕射等传输媒质所分别构成的信道。随参信道的特性比恒参信道要复杂的多。由于地面以上不同高度大气的电离层浓度不同，并

18、随机流性变化，对短波传输具有反射作用，对超短波具有对流层散射作用。乘性干扰的现象表现为各种类型的衰落，在多径信道中，发送端发出的信号通过多个反射之后沿多条路径到达接收端，这些路径具有不同的时延和不同的接收强度，它们之间的相互作用就形成了衰落。衰落可分为以下三种情况：(1) 慢衰落。它是由电离层随机变化引起的衰落。(2) 快衰落。它是由于短波信道的多径引起的衰落，对信号的影响更为严重。由于发射电磁波束传播到不同的电离层，反射后到地面构成不同的信号路径而至接收天线，也可能是地面反射回电离层经两次反射后到达接收天线。(3) 选择性衰落。频率选择性衰落是由于多径衰落导致的幅度随机性起伏衰落和相位随机性

19、变化所引起的。下面我们来介绍两个非常重要的衰落信道。(1) 瑞利衰落信道。瑞利衰落信道是移动通信中相当重要的衰落信道，它在很大程度上影响着移动通信系统的质童。在移动通信中，发送端和接收端都可能处在不停的运动状态之中，发送端和接收端之间的相对运动将产生多普勒频移。 多普勒频移与运动速度和方向有关，它的计算公式为式(3-1)。（2-1）其中，v是发送端和接收端的相对运动速度，是运动方向与发送端和接收端连线之间的夹角，是载波的波长。(2) 伦琴衰落信道。在移动通信系统中，如果发送端和接收端存在一条占优势的视距传播路径，这种信道就可以模拟成伦琴衰落信道。当发送端和接收端既存在视距传播路径，又有多条反射

20、路径时，它们之间的信道可以利用Simulink中的伦琴衰落信道模块和多径瑞利衰落信道模块的组合来进行仿真。瑞利衰落信道的统计模型衰落信道的统计特征的模型可以用几种概率分布表示。当信道中传送到接收机的信号的散射分量数目很大时，如电离层和对流层中的信号传播，应用中心极限定理可得到信道冲激响应的高斯过程模型，如果该过程是零均值的，那么任何时刻信道响应的包络都具有Rayleigh概率分布，而相位在（0,2）区间内是均匀分布的。瑞利分布的概率密度函数为：（2-2）是包络检波之前所接收电压信号的均方根值，是包络检波之前所接收信号包络的时间平均功率。瑞利分布的均值为： （2-3）瑞利分布的方差为：（2-4）

21、瑞利信号的均方根为：（2-5）瑞利分布的概率密度函数为：图2-2瑞利分布的概率密度函数2.3 通信系统中噪声概述及高斯白噪声2.3.1 噪声概述信号在信道中传输，要受到信道特性及噪声的影响，致使信号接收造成某些差错。此时受到干扰及噪声影响的信号可称为受扰信号。加性干扰的噪声其来源可分为人为噪声、自然噪声和内部噪声三个方面。通常将宇宙噪声、散弹噪声和热噪声归为起伏噪声，它们的统计特性基本上是高斯分布。2.3.2 高斯白噪声在研究通信系统时，为了分析方便，把噪声假想成一种理想化的形式，认为它通信的频段不受限于实际通信系统的频段，并包括电磁辐射全部可见频率，这种噪声称为白噪声。由于高斯过程的普遍存在

22、和高斯噪声过程在通信中的重要意义，我们归纳它们的统计特征的特点。高斯随机过程的一维统计特性只取决于均值和方差，二维统计特性主要取决于自协方差或自相关函数。高斯过程若广义平稳，则同时也严格平稳。窄带高斯噪声及其同相、正交分量均值皆为0，方差均等于窄带高斯噪声本身的方差，同相分量与正交分量不相关且统计独立。窄带噪声包络为瑞利分布，相位为均匀分布。载波信号加窄带噪声一般为赖斯分布，当信号幅度很大时，为高斯分布：幅度很小时，接近于瑞利分布；信号包络与噪声包络相差不大时，为赖斯分布。本设计中采用了加性高斯白噪声，它的幅度分布服从高斯分布，而它的功率谱密度又是均匀分布的，所以r（t）与s（t）不会有任何失

23、真。在信道中加入高斯白噪声，可以观察信号经过信道的变化情况。产生高斯白噪声序列的方法：假设信道以高斯白噪声相加来恶化信号，如下图所示。信道高斯白噪声n(t)发送信号s(t)接收信号r(t)=s(t)+n(t)图2-3通过AWGN信道的接收信号模型在间隔内，接收信号可以表示为：（） （2-6）其中n(t)表示具有功率密度谱（W/Hz）的加性高斯白噪声的样本函数。噪声性能指标：信噪比SNR： (dB)（2-7）噪声平均功率：（2-8）(其中为功率谱密度，为带宽)通信系统的主要性能通信系统的性能指标涉及其有效性、可靠性、适应性、经济性、标准性、可维护性等。尽管不同的通信业务对系统性能的要求不尽相同，

24、但从研究信息传输的角度来说，通信的有效性和可靠性是主要的矛盾所在。所谓有效性是指传输一定信息量时所占用的信道资源（频带宽度和时间间隔），或者说是传输的“速度”问题；而可靠性则是指接收信息的准确程度，也就是传输的“质量”问题。这两个问题相互矛盾而又相对统一，并且还可以进行互换。由于模拟通信系统和数字通信系统之间的区别，两者对有效性和可靠性的要求及度量的方法不尽相同。模拟通信系统的有效性可用有效传输频带来度量，同样的消息用不同的调制方式，则需要不同的频带宽度。可靠性通常用接收端解调器输出信噪比来度量。输出信噪比越高，通信质量就越好。不同调制方式在同样信道信噪比下所得到的解调后的输出信噪比是不同的。

25、数字通信系统的有效性可用传输速率和频带利用率来衡量，对于其可靠性可用差错率来衡量。差错率常用误码率和误信率表示。其中，误码率的计算公式如下： （2-9）三、无线衰落信道仿真程序设计3.1仿真设备装有MATLAB的PC机一台信号源16QAM调制瑞利信道高斯信道16QAM解调误码率、性能对比分析3.2系统整体框图图3-1系统整体框图串并变换2-4电平转换载波调制载波解调抽样判决4-2电平转换并串变换发送序列信道输出序列3.3 QAM调制解调图3-2QAM调制解调流程图瑞利信道仿真模块利用窄带高斯过程的特性，其振幅服从瑞利分布。首先产生独立的复高斯噪声的样本，并经过FFT后形成频域的样本，然后与S（

26、f）开方后的值相乘，以获得满足多普勒频谱特性要求的信号，经IFFT后变换成时域波形，再经过平方，将两路的信号相加并进行开方运算后，形成瑞利衰落的信号r(t)。图3-3瑞利信道仿真方法四、设计与仿真4.1 16QAM调制模块的建立与仿真16QAM调制解调原理进行了分析，这一章我将对系统仿真框图中的各个模块进行简单的介绍： 信号源本程序中，信号源为8位二进制代码-1 1 -1 1 1 1 -1 -1，sigexpand函数的作用是将代码扩展为码元宽度为1的双极性波形，如下图所示：图4-1星座图 串并转换信号源通过串并变换，将原来的一路信源信号变成两路信号，分别为上支路信号和下支路信号，独立地进行调

27、制和解调。串并变换的规则是根据序列编号的奇偶行，将编号为奇的码元编成一路信号，将编号为偶的码元编成一路信号。经过串并转换后，并行输出的每一路码元传输速率降为原来的一半即Rb/2. 输入d：-1 1 -1 1 1 1 -1 -1 上支路d_NRZ1：-1 -1 1 -1 下支路d_NRZ2：1 1 1 -1图4-2 串并转换后上下支路信号时域波形图4.1.3 2-4电平转换2-4电平转换就是将输入信号的2电平信号状态经过转换后变成相应的4电平信号。这里选择的映射关系如下所示：表4-1 电平映射关系映射前数据 双极性 电平/V 00 -1 -1 -3 01 -1 1 -1 10 1 -1 1 11

28、 1 1 3根据以上的映射关系，可得到上下支路分别为上支路d_NRZ1：-1 -1 1 -1 下支路d_NRZ2：1 1 1 -12-4电平转换信号： 上支路d_NRZ1：-1 -1 1 -1 2-4电平转换后：-3 1 下支路d_NRZ2：1 1 1 -12-4电平转换后：3 1图4-3 2-4电平转换后上下支路信号时域波形图 这里4电平信号的码元传输速率已降为Rb/4。 增加载波在本课题中，选用的载波是载波幅度A=1，载波频率fc=2Hz,上支路分量的载波是h1t=A*cos(2*pi*fc*t),正交分量的载波是h2t=A*sin(2*pi*fc*t)。上下支路信号在加载波之前还经过平滑

29、处理，以滤除较高频率的信号，使实验结果更加理想。上下支路信号加载波后的图形为：图4-4上下支路调制信号时域波形图 调制信号形成上下支路调制信号形成后，将两个分量相加，既可得到16QAM调制信号，如下图所示：图4-5 已调信号波形图4.2 16QAM调制信号的噪声叠加（1）高斯噪声本次仿真采用的噪声是高斯白噪声，这是一种最常见的噪声，白噪声的功率谱密度在所有频率上均为一常数，且仅在t=0时才相关，而在任意两个时刻的随机变量都是不相关的。对已调制信号可采用awgn函数添加加性高斯噪声。y = awgn(m,n,p) 产生一个m行n列的高斯白噪声的矩阵，p以dBW为单位指定输出噪声的强度。为使解调效

30、果较好，采用噪声的强度较小，设置Pn=-10dB。（2）瑞利信道仿真 另外产生一个瑞利信道，并使已调信号通过此信道，统计此路解调信号的误码率，并与高斯信道下的进行对比。Matlab仿真步骤：（1）首先确定接收机速率，载波频率，多普勒频移和的N点；（2）根据书上多普勒平移公式，计算出的正频率部分频域值，对其取反，合成得到全频域的（3）在matlab生成两个独立的复数高斯噪声：对于噪声一，在时域信号产生N点高斯噪声信号，对该号作快速傅里叶变换（FFT）后得到复数高斯信号频谱；同样的方法再合成高斯噪声二；（4）将多普勒功率谱分别与复数高斯噪声信号相乘，再求快速傅里叶反变换（IFFT），得到时域信号；

31、（5）对IFFT后的信号求模后分别平方相加，最后开方后即得到瑞利衰落信号。chan = rayleighchan(1/10000,100); %产生瑞利信道fadedSig = filter(chan,s_16qam); %使信号通过瑞利信道4.3 16QAM解调模块的建立与仿真系统先前所得的16QAM调制信号通过高斯白噪声信道以后便可以解调了。本文所采用的解调器原理为相干解调法，即已调信号与载波相乘，送入到低通滤波器，其对应原理图中信号输入并与载波相乘后通过LPF的部分，输出送入到判决器判决，再经4-2电平转换和并串转换即可得到解调信号。 滤波器IIR滤波器采用的巴特沃斯低通滤波器有现成的模

32、型，我们可以加以利用，因此在本文涉及的仿真中滤波器均选择巴特沃斯低通滤波器。巴特沃斯低通滤波器的标准形式为：b,a=butter(N,W0); y=filter(b,a,x);N表示要选取的低通滤波器的阶数，W0表示滤波器的截止频率，b,a为滤波器返回的特性参数。在第二行程序中，x表示输入序列，y表示输出序列。整个函数表示信号通过滤波器的过程。图4-6上下支路通过低通滤波器信号时域波形图4.3.2 4-2电平转换抽样判决是在每个码元中间抽样，并用几个判决语句进行判决，例如当码元幅度大于1时，判为3。4-2电平转换是2-4电平转换的逆过程，其映射关系如下图所示： 表4-2电平映射关系映射前数据电

33、平/V双极性 -3 00-1 -1 -1 01-1 1 1 101 -1 3 111 1图4-7上下支路抽样判决及4-2转换后信号时域波形图并串转换经过并串转换即可得到解调信号，将两路并联信号经过转换成为一路信号，下图为基带信号和解调信号的对比及它们的频率谱密度图图4-8高斯噪声下的基带信号与解调信号图4-9瑞利信道下的基带信号与解调信号可以看出，加高斯白噪声时并没有错码，经统计得此时误码率为0。而当已调信号通过瑞利信道时，从4-9中可以看出出现了明显的误码，统计得误码率为32%。因此无线通信中的瑞利信道对信号传输造成的影响更大。4.4 QAM性能分析4.4.1 16QAM抗噪声性能仿真对于Q

34、AM，可以看成是由两个相互正交且独立的多电平ASK信号叠加而成。因此，利用多电平误码率的分析方法，可得到M进制QAM的误码率为1： (4.1) 式中，Eb为每码元能量，n0为噪声单边功率谱密度。通过调整高斯白噪声信道的信噪比snr（Eb/No），可以得到如图5-2所示的误码率图：图4-10 QAM信号误码率分析可见16QAM信号的误码率随着信噪比的增大而逐渐减小，这与理论分析是完全一致的。总结本课程设计主要是在MATLAB环境下，设计仿真了基于QAM调制的无线衰落信道模型，首先对信号源产生的二进制序列进行QAM调制，然后使已调信号分别通过瑞利，高斯信道，在进行QAM解调，对解调后的信号进行误码

35、率分析。通过观察本次课程设计的仿真结果，对比分析了16QAM调制信号通过瑞利信道，高斯信道后的输出波形，及误码率可明显看出，瑞利信道使传输的信号产生更高的误码率。通过观察16QAM调制解调信号的频谱以及对查阅的资料进行研究发现16QAM的抗噪声性能相当优越。此次课程设计完成了对基于QAM调制的无线衰落信道模型的仿真，实现了设计要求的目标。通过本次课程设计使我对MATLAB语言有了更进一步的掌握，对无线衰落信道的性能有了更清楚的了解。参考文献 1BELLO. P. A. Characterization of randomly time variant lin2ear channelsJ . I

36、EEE Trans. Communications Systems , 1963 ,11 (4) :360 - 393.2CLARKE ,R. H.A statistical theory of mobile radio receptionJ . Bell Systems Technical Joumal ,1968 ,47 (6) : 957 - 1000.3樊昌信，曹丽娜编著，通信原理，北京：国防工业出版社，2010：73-844樊昌信，通信原理，北京：国防工业出版社，2002:70-80:67-905曹志刚等著，现代通信原理，北京：清华大学出版社，2001.5。：65-786吴伟陵等著，

37、移动通信原理，北京：电子工业出版社，2005：30-577.李建新，现代通信系统分析与仿真-MATLAB 通信工具箱，西安：西安电子科技大学出版社，2000：13-458潘子宇，Matlab通信仿真设计指导书，南京工程学院，2011：24-369刘敏，MATLAB通信仿真与应用，北京：国防工业出版社：56-60附录源程序代码如下1、将输入的序列扩展成间隔为N-1个0的序列function out=sigexpand(d,M)N=length(d);out=zeros(M,N);out(1,:)=d;out=reshape(out,1,M*N);2、计算信号的傅里叶变化functionf,sf=

38、T2F(t,st);dt=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); 主函数代码fc=2; % 载波频率N_sample=8; % 基带码元抽样点数N=8; % 码元数Ts=1; % 码元宽度dt=Ts/fc/N_sample; % 抽样时间间隔T=N*Ts; % 信号持续时间长度t=0:dt:T-dt; % 时间向量Lt=length(t); % 时间向量长度tx1=0; % 时域波形图横坐标起点tx2=8; %时域波形图横坐标终点ty1=-4.5;

39、 %时域波形图纵坐标起点ty2=4.5; %时域波形图纵坐标终点fx1=-10; %功率谱图横坐标起点fx2=10; %功率谱图横坐标终点fy1=-40; %功率谱图纵坐标起点fy2=25; %功率谱图纵坐标终点%产生二进制信源d=-1,1,-1,1,1,1,-1,-1dd=sigexpand(d,fc*N_sample); % 双极性gt=ones(1,fc*N_sample); % NRZ波形d_NRZ=conv(dd,gt); % 基带信号figure(1);subplot(2,3,1);plot(t,d_NRZ(1:Lt);axis(tx1,tx2,ty1,ty2);xlabel(&#

40、39;时间(S)');ylabel('幅度');title('基带信号时域波形图');grid;figure(5);subplot(2,2,1);plot(t,d_NRZ(1:Lt);axis(tx1,tx2,ty1,ty2);xlabel('时间(S)');ylabel('幅度');title('基带信号时域波形图');grid;f1,d_NRZf=T2F(t,d_NRZ(1:Lt);figure(5);subplot(2,2,2);plot(f1,10*log10(abs(d_NRZf).2/T);ax

41、is(fx1,fx2,fy1,fy2);xlabel('频率(Hz)');ylabel('功率谱密度(dB/Hz)');title('基带信号功率谱图');grid;% 串并转换d1=;d2=;for i=1:N/2 d1(i)=d(2*(i-1)+1); d2(i)=d(2*(i-1)+2);end% 上支路dd1=sigexpand(d1,2*fc*N_sample); gt1=ones(1,2*fc*N_sample); d_NRZ1=conv(dd1,gt1); figure(1);subplot(2,3,2);plot(t,d_NRZ1

42、(1:Lt);axis(tx1,tx2,ty1,ty2);xlabel('时间(S)');ylabel('幅度');title('串并转换后上支路信号时域波形图');grid;% 下支路dd2=sigexpand(d2,2*fc*N_sample); d_NRZ2=conv(dd2,gt1); figure(1);subplot(2,3,5);plot(t,d_NRZ2(1:Lt);axis(tx1,tx2,ty1,ty2);xlabel('时间(S)');ylabel('幅度');title('串并转换后

43、下支路信号时域波形图');grid;% 载波h1t=cos(2*pi*fc*t);h2t=sin(2*pi*fc*t);figure(1);subplot(2,3,4);plot(t,h1t);axis(tx1,tx2,ty1,ty2);xlabel('时间(S)');ylabel('幅度');title('载波信号时域波形图');grid;% 上下支路2-4电平转换d11=; % 上支路d22=; % 下支路·for m=1:N/4; d11(m)=2*d1(2*m-1)+d1(2*m); d22(m)=2*d2(2*m-1)

44、+d2(2*m);enddd11=sigexpand(d11,4*fc*N_sample); % 上支路gt2=ones(1,4*fc*N_sample); d_NRZ11=conv(dd11,gt2); figure(1);subplot(2,3,3);plot(t,d_NRZ11(1:Lt);axis(tx1,tx2,ty1,ty2);xlabel('时间(S)');ylabel('幅度');title('2-4电平转换后上支路信号时域波形图');grid;figure(3);subplot(2,2,1);plot(t,d_NRZ11(1:L

45、t);axis(tx1,tx2,ty1,ty2);xlabel('时间(S)');ylabel('幅度');title('2-4转换后上支路信号时域波形图');grid;dd22=sigexpand(d22,4*fc*N_sample); % 下支路d_NRZ22=conv(dd22,gt2); figure(1);subplot(2,3,6);plot(t,d_NRZ22(1:Lt);axis(tx1,tx2,ty1,ty2);xlabel('时间(S)');ylabel('幅度');title('2-4

46、转换后下支路信号时域波形图');grid;figure(3);subplot(2,2,3);plot(t,d_NRZ22(1:Lt);axis(tx1,tx2,ty1,ty2);xlabel('时间(S)');ylabel('幅度');title('2-4转换后下支路信号时域波形图');grid;% 对上下支路信号进行平滑pd_NRZ11=d_NRZ11(1:Lt);pd_NRZ22=d_NRZ22(1:Lt);tao=3/16;tm=0:dt:8*Ts*tao-dt;for w=2:N/4 phxh1=d11(w-1)+(d11(w)-

47、d11(w-1)*0.5*(1+sin(pi*tm/8*tao*Ts); phxh2=d22(w-1)+(d22(w)-d22(w-1)*0.5*(1+sin(pi*tm/8*tao*Ts); for k=1:8*tao*fc*N_sample pd_NRZ11(w-1)*64+k)=phxh1(k); pd_NRZ22(w-1)*64+k)=phxh2(k); endend% 生成16QAM信号s_16qam1=pd_NRZ11 .* h1t; %生成上支路频带信号s_16qam2=pd_NRZ22 .* h2t; %生成下支路频带信号figure(2);subplot(2,4,1);plo

48、t(t,s_16qam1);axis(tx1,tx2,ty1,ty2);xlabel('时间(S)');ylabel('幅度');title('上支路已调信号时域波形图');grid;f2,s_16qam1f=T2F(t,s_16qam1);figure(2);subplot(2,4,5);plot(f2,10*log10(abs(s_16qam1f).2/T);axis(fx1,fx2,fy1,fy2);xlabel('频率(Hz)');ylabel('功率谱密度(dB/Hz)');title('上支路已

49、调信号功率谱图');grid;figure(2);subplot(2,4,2);plot(t,s_16qam2);axis(tx1,tx2,ty1,ty2);xlabel('时间(S)');ylabel('幅度');title('下支路已调信号时域波形图');grid;f3,s_16qam2f=T2F(t,s_16qam2);figure(2);subplot(2,4,6);plot(f3,10*log10(abs(s_16qam2f).2/T);axis(fx1,fx2,fy1,fy2);xlabel('频率(Hz)')

50、;ylabel('功率谱密度(dB/Hz)');title('下支路已调信号功率谱图');grid;s_16qam=s_16qam1+s_16qam2;figure(2);subplot(2,4,3);plot(t,s_16qam);axis(tx1,tx2,ty1,ty2);xlabel('时间(S)');ylabel('幅度');title('已调信号时域波形图');grid;f4,s_16qamf=T2F(t,s_16qam);figure(2);subplot(2,4,7);plot(f4,10*log10

51、(abs(s_16qamf).2/T);axis(fx1,fx2,fy1,fy2);xlabel('频率(Hz)');ylabel('功率谱密度(dB/Hz)');title('已调信号信号功率谱图');grid; % 信道加入高斯白噪声进行接受解调% 产生高斯白噪声 m1=1;p1=-10;noise =wgn(m1,Lt,p1);%产生瑞利信道chan = rayleighchan(1/10000,100);%使信号通过瑞利信道fadedSig = filter(chan,s_16qam); % 接收信号y_16qam = s_16qam +

52、 noise;figure(2);subplot(2,4,4);plot(t,y_16QAM);axis(tx1,tx2,ty1,ty2);xlabel('时间(S)');ylabel('幅度');title('接收信号时域波形图');grid;f5,y_16qamf=T2F(t,y_16qam);figure(2);subplot(2,4,8);plot(f5,10*log10(abs(y_16qamf).2/T);axis(fx1,fx2,fy1,fy2);xlabel('频率(Hz)');ylabel('功率谱密度(

53、dB/Hz)');title('接收信号功率谱图');grid;% 相干解调% 通过乘法器1r_16QAM11 = y_16qam .* h1t;%通过低通滤波器b,a=butter(3,0.1);r_16QAM11=filter(b,a,r_16QAM11);figure(3);subplot(2,2,2);plot(t,r_16QAM11)axis(tx1,tx2,ty1,ty2);xlabel('时间(S)');ylabel('幅度');title('上支路通过低通滤波器信号时域波形图');grid;%抽样判决dd111=r_16QAM11(2*fc*N_sample:4*fc*N_sample:end); %在每个码元中间抽样 dd1111=;for n=1:N/4 if dd111(n)>1 dd1111(n)=3; elseif 0<dd111(n) & dd111(n)<1 dd1111(n)=1; elseif -1<dd111(n) & dd111(n)<0 dd1111(n)=-1; else dd1111(n)=-3;