物理层仿真 实验.docx

1、物理层仿真物理层仿真 实验实验物理层仿真实验系别：通信工程系专业：通信工程系 x 级学号：姓名：吕 XX实验时间：2014 年 6 月 30日撰写日期：2014 年 7 月 3日 实验一：随机信号的产生1.1 实验目的掌握各种伪随机序列的产生方法。1.2 实验原理1.2.1 Wichmann-Hill算法产生均匀分布随机变量该算法是通过将 3 个周期相近的随机数发生器产生的数据序列进行相加，进而得到更大的周期的数据序列。定义三个随机数发生器：以上三式中均需要设定一初始值，这三个初始值一般称为种子。产生的三个序列的周期分别是：30269、30307、30323。将这三个序列组合相加即可得到一个周

2、期更大的均匀分布随机序列：1.2.2 逆变换发产生 Rayleigh 分布随机变量逆变换法的基本思想如图 1.1所示，条件是产生的随机变量的分布函数具有闭合表达式。图 1.1 将一个不相关均匀分布的随机序列 U 映射到一个具有概率分布函数 Fx(x)的不相关序列随机序列 XRayleigh 分布的分布函数：根据上面的逆变换方法有：因此，根据上式即可将均匀分布的随机变量映射为 Rayleigh 分布的随机变量。1.2.3 根据 Rayleigh 分布随机变量产生 Gaussian分布随机变量基于 Rayleigh 随机变量，可以方便的产生 Gaussian分布随机变量。关系如下：或者其中 u1

3、和 u2 分别是两个（01）之间均匀分布的随机变量，产生的 X和 Y 均为高斯随机变量。1.3 实验内容根据实验原理部分，完成以下实验任务。任务一：采用 Willmann-Hill算法产生 10000 个均匀分布的随机变量，根据两组种子做出随机序列的直方图（命令 hist），两组种子自行设定。任务 2：使用逆变换法产生其它分布的随机变量（Rayleigh）设定 sigma2为 0.5、1 和 2，分别做出 Rayleigh 随机变量序列的直方图，观察图形。任务 3：通过 Rayleigh 分布随机变量产生 Gaussian分布随机变量，设定 sigma2为1 产生标准正态分布的随机变量，计算其

4、方差并做直方图。改变 sigma2 值观察直方图的变化。1.4 实验结果与分析任务一：均匀分布两组直方图如下 由上图可看出，产生的随机序列满足（01）均匀分布。任务二：Rayleigh 分布直方图 由这三张图可看出，图像大致满足瑞利分布，且方差不同，瑞利分布的横坐标不同。方差越大，图像在横坐标轴上的跨度越大任务三：通过瑞利分布随机变量产生高斯随机分布变量，方差为 1计算其方差得：varx=1.0033；vary=1.0079与所设定的 sigma2=1的值相接近。实验二：QPSK 系统仿真2.1 实验目的掌握简单调制方法的基带仿真实现，以及 AWGN 信道和 Rayleigh 平坦衰落信道的建

5、模，并完成在这两种信道下的误码率仿真。2.2 实验原理QPSK 发射信号的等效基带格式和星座图Sk=A(cosk+jsink)其中 A 为信号幅度，k 为调制相位，取值范围为(/4,3/4,5/4,7/4),比特与相位之间的映射关系如图 2.1 所示。图 1.2 QPSK、8PSK、16QAM信号星座图2.3 实验内容任务 1：根据参考教材第一章的 BPSK 示例，采用实验一产生随机数的方法重新改写程序，并仿真的 BER 结果，将结果与示例仿真结果进行比较。（提示：主要修改的内容有比特的产生、AWGN 噪声的产生，即只要涉及需要随机数的地方均改用实验一中的方法得到）。任务 2：参考 BPSK

6、程序，搭建 QPSK、8PSK、16QAM 的基带仿真程序，仿真在AWGN 信道和 Rayleigh 平坦衰落信道下的误码率性能，比较它们的误比特率和带宽效率，解释它们误码率性能差别的原因。2.4 实验结果与分析任务一：用实验一的方法产生的随机数和用 matlab 自带函数产生的随机数经过BPSK 通信系统仿真的 BER结果的比较。如图，对这两种方法产生的随机数进行 BPSK 仿真，其误码率几乎一致。任务二：BPSK、QPSK、8PSK、16QAM 在 AWGN信道下的仿真：以上四图分别是 BPSK、QPSK、8PSK、16QAM 在伪随机数下 AWGN信道中的仿真曲线，可见其结果与理论曲线几

7、乎重叠，这是在满足大数定理下的仿真结果。BPSK、QPSK、8PSK、16QAM 在 Rayleigh信道下的仿真曲线：可见各个调制在 AWGN 和 Rayleigh信道下的仿真结果几乎一致。接下来为节省仿真时间，在比较他们之间的性能差异将直接利用理论值曲线进行比较。BPSK、QPSK、8PSK、16QAM 仿真在 AWGN信道下的误码率性能比较。如上图所示，在相同信噪比下，BPSK 和 QPSK 的误码率性能最好，8PSK 次之，QAM 最差，8PSK 和 QAM 的误码率性能大致在同一数量级。但是从带宽效率看，QAM 的带宽效率最高为 4,8PSK 的带宽效率为 3，QPSK 为 2。说明

8、在相同信噪比下，带宽效率和误码率是一对矛盾互补的，它们分别对应着通信系统的有效性和可靠性。QPSK、8PSK、6QAM 在 Rayleigh 衰落信道下的误码率性能，并比较它们的比特率。如上图可看出，在 Rayleigh 平坦衰落信道下 BPSK、QPSK和 8PSK的性能较在AWGN 信道下有较大的减少，要得到相同的误码率，Rayleigh 信道下要使用比 AWGN信道下更大的信噪比。在相同信噪比下，QPSK 的误码率会比 8PSK 小，但小的不是很多，而 8PSK 的带宽效率比 QPSK 高，所以在 Rayleigh 平坦衰落信道下，采用 8PSK 调制比 QPSK 更好。实验三：直接序列

9、扩频抗 ISI 研究3.1 实验目的掌握直接序列扩频技术(DSSS)原理及基带仿真实现，理解扩频码设计对直接序列扩频抗多径干扰的影响。3.2 实验原理定义：扩展频谱技术一般是指用比信号带宽宽得多的频带宽度来传输信息的技术。直接序列扩频就是用比信息速率高很多倍的伪随机噪声码（PN）与信号相乘来达到扩展信号的带宽。上图是直接序列扩频的系统框图，从图可以看到，与传统窄带系统相比主要引入了扩频和解扩过程(通过扩频码实现)，广义来说可以把扩频/解扩也堪称一种调制/解调。通过扩频码的作用窄带信号变成了宽带信号。因此，扩频系统具有良好的抗窄带干扰和抗多径(ISI)干扰的能力。下图展示了 DSSS 系统抗窄带

10、干扰的原理，其中 S(f)是调制信号的频谱，经过扩频进入接收端的信号频谱为 S(f)*Sc(f)，由于扩频码具有宽频特性，因此，信号带宽被展宽，如图 1.3(b)所示，I(f)表示窄带干扰的频谱，其幅度超过有用信号，但经过解扩后，如图 1.3(c)所示，由于扩频码的相关性，调制信号的功率得到恢复，而干扰信号的频谱被扩散，达到了抗窄带干扰的效果。类似的由于扩频码的自相关特性，DSSS 系统还有有效的抗 ISI 干扰，其频域原理如图 3.3 所示，其中 H(f)=+e-j2f 表示 ISI信道的频域传递函数。在直接序列扩频系统中，其扩频码的设计是整个系统的关键，扩频码的选择要求具有好的互相关性和自

11、相关性，即要能有效的抗多用户干扰和抗 ISI 干扰。为了抗多用户干扰，需要尽量选择互相正交的扩频码，为了抗 ISI 干扰，需要扩频码的自相关特性尽量满足一定的特性。常用的扩频码序列有 m序列、Golden 序列等等。接收机可采用相干接收机、RAKE 接收机。3.3 实验内容根据 DSSS 技术原理搭建仿真系统，并完成下面的任务。任务一：参考后面提供的程序，采用 m序列替代程序中的随机序列作为扩频码，仿真其性能，完成与随机序列扩频性能之间的比较；任务二：在 ISI 信道下，研究 m序列与 golden序列抗 ISI 的性能。3.4 实验结果与分析任务一：比较 m 序列与随机序列的扩频性能高斯信道

12、下采用 m 序列扩频和采用随机序列扩频，信噪比相同时，它们的误码率是大致也是相同的。说明在无码间串扰的信道下 m序列扩频与一般随机序列的性能没有显著差别。任务二：在 ISI 信道下，m序列与 golden序列抗 ISI 性能比较。如上图所示，在 ISI 信道下，采用 m序列和 golden序列扩频都能够减少通信系统的误码率。但是 golden 序列扩频的性能明显优于 m序列的性能。m 序列需要在足够大信噪比下才能得到满意的误码率。4 实验小结4.1 实验一实验小结在通信系统中的各种信号以及噪声、信道衰落都可以看作相应的随机过程，通信系统仿真就是要通过计算机产生各种随机信号，并依据通信原理对这些

13、随机信号进行相应的处理，获取期望的结果。这次的实验中，我掌握了 Wichmann-Hill算法，并利用此算法产生的均匀分布的随机变量变换得到了 Rayleigh 分布的随机变量，以及 Gussian分布的随机变量，为后续实验奠定了基础。4.2 实验二实验小结在本次实验中，我利用上一个实验一产生的伪随机数，对 BPSK、QPSK、8PSK和 16QAM 通信系统进行了仿真，并且得到了和理论曲线一致的仿真曲线。在实验过程中的第一个问题：我一开始采用 10000个随机数作为信号源，发现仿真的结果与理论相距甚远。后来通过分析，明白在误码率达到 10(-3)10(-5)时，仿真的结果在满足大数定理的情况

14、下即错误个数尽可能的大，实验数据至少要 10(7)数量级。因此后来的实验中，我都采用另一种仿真模式，即设定错误限，在每次达到错误限之前一直产生伪随机源和伪随机噪声，直到达到大数定理，才完成仿真。通过设定错误限来控制仿真的曲线稳定性。最后得到几乎与理论曲线一致的仿真曲线。在实验过程中的第二个问题：在进行仿真的过程中，我发现 16QAM 的仿真曲线与理论曲线不一致，虽然仿真曲线和理论曲线在数量级上一致，但仿真曲线的误码率总是高于理论曲线，并与其几乎平行。这与其他的仿真中仿真曲线几乎与理论曲线重合的现象不一致，让我很费解。后来分析，觉得这种现象产生的原因是，16QAM 的星座图不是最好的星座图，在译

15、码过程中，产生错误时，译码结果没有最小限度的减小错误的比特，也就是星座图上相邻两个码点距离过大。通过修改星座图，我发现仿真曲线与理论曲线又一致了，下图是我修改后的星座图。所有相邻的码字距离皆为 1。这样能最小化错误的比特数目。4.3 实验三实验小结这个实验中，我又重新上网搜索了 m序列和 golden序列的概念，并通过网上的一些描述原理，产生了两组 m 序列和 golden序列作为扩频码，完成了实验的仿真。通过修改实验书上的代码，我对一些 matlab 矩阵实验平台上工具又有了进一步的了解。比如 ISI 信道的建模，可以直接通过卷积实现，而且实验书上的例程是利用错误限来控制仿真数据尽量的满足大

16、数定理。这也给了我很多的启发。虽然实验中的一些原理不是很清晰，但是让我有了总体的把握，并加深了我对通信系统的理解。附录：代码实验一：lab1.mclc;clear all;%lab1 的实验一%种子设定x(1)=22;y(1)=21;z(1)=23;u(1)=mod(x(1)/30269+y(1)/30307+z(1)/30323),1);%随机数发生器for i=2:10000 x(i)=mod(171*x(i-1),30269);y(i)=mod(170*y(i-1),30307);z(i)=mod(172*z(i-1),30323);u(i)=mod(x(i)/30269+y(i)/30307+z(i)/30323),1);end%种子设定x(1)=9;y(1)=8;z(1)=10;v(1)=mod(x(1)/30269+y(1)/30307+z(1)/30323),1);%随机数发生器for i=2:10000 x(i)=mod(171*x(i-1),30269);y(i)=mod(170*y(i-1),30307);z(i)=mod(172*z(i-1),30323);v(i