ch通信系统的蒙卡仿真.ppt

也就是 在AWGN环境下，概率通常用高斯Q函数来表示。于是有 其中 是高斯Q函数的定义式。注意，在信号能量相等的情况下，阈值T为零，所以v2＝-v1。因此确定了v1和?n，就完全确定了系统的误比特率。为了确定差错概率，我们只要开发仿真程序来估计v1和?n即可，不需要使用蒙特卡罗方法对差错发生的次数进行计数。 可以通过一个无噪声的仿真来确定vl值。如果去掉信道中的噪声，图10-11(a)中的两个条件pdfs退化成冲激函数(?n=0)，如图10-11(b)所示。每个冲激函数都有单位面积，冲激函数的位置分别确定了v1和v2。 只要简单地仿真噪声通过的那部分系统，就可以确定?n的值，这里考虑的系统中包含了用积分--清除符号检测器来建模的接收机。假设这部分系统具有传递函数H(f)，假设输入到匹配滤波器接收机中的白噪声（?-相关）具有双边功率谱密度No/2，随机变量Vk的方差为 正如在第7章所见，等价噪声带宽定义为 这也是接收机的等价噪声带宽。于是有差错概率为 注意到，尽管系统是加性高斯白噪声的，非线性放大器仍可能会影响系统的性能，因为非线性放大器会影响发送信号的形状，而这又会影响v1值。 可以很容易地把结果扩展到记忆长度为M个符号的情况。 现在考虑有发送滤波器的情况。这个滤波器的影响是，在时间上将所发送符号的能量扩展到发送符号周期以外，从而引起码间干扰。如果滤波器的记忆长度是两个符号周期，则与发送符号相对应的差错概率不仅依赖于当前的符号，而且和前面发送的一个符号有关。因此在计算差错概率时牵涉到四个概率密度函数，而不是如图10-11所示的两个概率密度函数，这种情况如图10-12(a)所示。和以前相同，执行无噪声的仿真可得出v1、v2、v3和v4。系统的差错概率为 (二)等效噪声源 在应用半解析方法时，我们采用了等效噪声源的思想。我们已经看到判决统计Vk是三个分量的函数，也就是说，Vk=f(Sk,Dk,Nk)，其中Sk是由于信号产生的，Dk来自于由系统因素引起的失真（如ISI），而Nk是由噪声引起的。 Sk和Dk的影响可由蒙特卡罗仿真确定，而如前所见，由Nk表示的噪声的影响则用解析的方法处理。 如果执行无噪声的仿真，则所得的充分统计量Vk,nf只是Sk和Dk的函数。该统计量加上一个方差由式（10-11)给出的随机变量Nk，可得 可将随机变量Nk视作来自于如图10-13所示等价噪声源ne（t）的一个样本。该等价噪声源包含了反映在积分--清除检测器的积分器输出中的，由热噪声、干扰和其他信道损伤所产生的总影响。如果信道噪声是白噪声，则可用冲激响应，或等价地用式(10-11)所定义的传递函数，来把信道噪声变换到积分器的输出中。 （三） PSK系统误比特率的半解析估计方法 现在我们简要地考虑开发一个算法，以利用半解析仿真来确定BPSK系统中的误比特率，我们采用一种很容易就可扩展到QPSK系统的方法进行处理。考虑如图10-14所示的信号星座，发送信号的点用S1和S2表示，相应的判决区域用D1和D2表示。如果发送的是Si而接收信号落在区域Di中，则接收机作出了正确的判决；不然，就发生了差错。在图10-14中，我们假设发送的是S1而接收到的是S1。正如在上节中所讨论的那样，由于存在符号间干扰、非线性或者损害信号的其他因素，S1跟S1会不同。用dx表示S1和S1之差。在发送S1的条件下，条件差错概率为 使用高斯Q函数来表示，上式变为 因此，已知由蒙特卡罗仿真求得的S1以及?n，就可以确定条件误比特率。在确定?n时，可由仿真的冲激响应h[n]求出BN值。 * * 一、两个蒙特卡罗实例 二、半解析方法 （一）基本考虑 （二）等效噪声源 （三）PSK系统误比特率的半解析估计方法 （四）QPSK系统误比特率的半解析估计方法 （五）数据序列的选择 引 言 本章将对上一章探讨的基本蒙特卡罗（MC）方法作进一步的扩展。 一、将考虑两个通信系统的仿真实例。第一个是相移键控（PSK）数字通信系统，它是本书后面部分开发的仿真系统的基本构建模块。接下来讨论的是比较复杂的差分QPSK系统的仿真，其中考虑了相位和符号同步误差的影响。 二、将重点讨论半解析（SA）方法，即把蒙特卡罗仿真和解析的方法结合在一起。 蒙特卡罗方法几乎不需要任何的数学分析，并且只要系统方框图中每个功能模块的信号处理算法是已知的就能应用。因此，蒙特卡罗仿真是一个很通用的工具，但它的缺点是仿真运行时间长，因此需要在仿真精度和仿真运行时间之间作个基本的折中。 半解析方法需要更高水平的分析，但得到的回报是大大地降低了仿真运行时间。另外，运行一次蒙特卡罗仿真得到的是单个Eb/No误比特率估计值，而半解析仿真得到的是以Eb/No为函数的完整的误比特率曲线。然而，