通信系统仿真技术 第4章蒙特卡洛仿真与随机数产生.ppt

第4章 蒙特卡洛仿真与随机数产生 4.1 蒙特卡洛仿真原理 4.2 随机数产生 4.3 独立随机序列的产生 4.4 相关随机序列的产生 4.5 随机数产生器的测试 4.1 蒙特卡洛仿真原理 4.2 随机数产生 4.3 独立随机序列的产生 4.4 相关随机序列的产生 4.5 随机数产生器的测试 1、m序列的产生 定义：由带线性反馈的移位寄存器产生的周期最长的一种序列。 由任何4级反馈移位寄存器产生的序列的周期最长为15。 一个n级反馈移位寄存器可能产生的最长周期等于(2n-1)。反馈电路如何连接才能使移位寄存器产生的序列最长，这就是本小节后面将要讨论的主题。 线性反馈以为寄存器的结构如下： 递推方程：对于任意一状态ak 特征方程： 母函数： 产生这个最大周期（p=2n-1）序列的充要条件，就是线性反馈移位寄存器的特征多项式f(x)为本原多项式。 f(x)多项式所需要满足的条件如下： ①f(x)为既约的多项式，即不能分解因子的多项式； ②f(x)能够被 整除，其中p=2n-1； ③但f(x)不能被 整除，其中 。 2、m序列的性质 ①均衡特性 ②游程分布 ③移位相加特性 ④自相关函数 ⑤功率谱密度 ⑥伪噪声特性 ⑦m序列具有n比特所有可能的组合 4.3.3 M进制伪随机序列 如果要仿真MASK、MFSK、MPSK以及QAM等系统，就需要利用M进制PN序列。 M进制反馈移位寄存器的特征多项式为 逆多项式为 系数ci可在共有q个元素的有限区域内取值（q进制），这个有限域定义为GF(q)。 以GF(q)域上的本原多项式构建的M进制反馈移位寄存器，能够产生一个最大长度为的q进制序列。其递推描述关系式为 如果要仿真16-QAM系统，可以采用上述类似的方法，在GF(4)中选择系数，确定本原多项式，构建两个8进制反馈移位寄存器电路，来表示I和Q信道。 在许多应用中，经常需要产生具有某种特定自相关函数，或者产生具有某种特定功率谱密度形式的输出序列，例如，当需要仿真随机时变的移动通信信道，其功率谱密度函数 当然还可以利用相关函数形式来表示相关序列。 除了单一采样值相关性研制外，还存在矢量形式。例如，多经效应的仿真。 矢量随机过程的相关性主要反映在两方面： ①时间轴上的相关特性，这种相关特性是瞬时的，并引出了功率谱密度这个概念； ②各个元素间的相关性，它代表了不同方向上的相关，通常出现较多的是空间维上的相关等。 4.4.1 相关高斯标量序列 不相关的高斯序列适当的处理，就可以变成具有相关特性的高斯序列。 ①经过线性变换得到，基于相关函数的时域法； ②经过滤波器得到，基于功率谱密度的频域法。 1、时域法 时域法通常采用ARMA过程，构建离散的时间模型。问题是确定a b 。 X(n)为输入模型的不相关高斯序列，通常假设其均值为零，方差为 ；Y(n)为希望产生的相关高斯序列。Z变换后得： 如果ARMA(N，M)模型就退化成为AR(N)模型，其Y(n)随机序列产生模型为 计算式（4.4-5）的自相关函数可得 两个等式写成矩阵形式，得 计算步骤： ①已知Y(n)的自相关函数和 ； ②求出ak ； ③将高斯序列X(n)转换成Y(n)序列。 2、频域法 如果已知输出序列Y(n)的功率谱密度，就可以设计一个频域滤波器，把不相关的高斯序列X(n)变成相关的序列Y(n)，这种方法通常被称为滤波法。 这时 选择H(s)作为滤波器的传输函数，为了保证该滤波器因果并且稳定，则H(s)的所有极点都应当位于以虚轴为中心的s左半平面。 得到了滤波器传输函数，用IIR或FIR设计实现。 在许多情况下，输出谱密度可以凭借经验形式给出，在这时可以采用下面方法进行计算： ①凭经验，对谱密度函数进行因式分解； ②利用滤波器进行确定。 4.4.2 相关高斯矢量序列 l个标量高斯随机过程的采样值就可以认为是矢量随机过程Y(k)的元素，即