《信息与通信系统仿真》课件8 通信系统仿真方法

1、别志松 通信系统仿真方法2主要内容通信系统仿真的多学科性通信系统仿真的层次性蒙特卡洛仿真原理二进制反极性基带通信系统的蒙特卡洛仿真蒙特卡洛仿真精度分析与仿真量确定方法方差缩减方法：重要性采样通信系统仿真的多学科性其他系统的仿真一样，通信系统仿真同样包含建立系统模型、开发仿真模型、实现仿真模型、执行仿真和仿真验证等步骤，在通信系统仿真的这一系列过程中，需要用到多个领域、多个科目的知识多学科性通信系统仿真数论概率论与随机过程数理统计数值分析信号 与系统数字信号处理排队论估计理论通信原理程序设计算法与数据结构分布式计算多学科性信号与系统中关于线性系统的理论给出了一整套时域和频域的分析方法，为建立通信

2、链路的系统模型打下了基础；通信原理的重要性显而易见，理解调制、均衡等各种通信技术是建立适当系统模型的关键；数字信号处理技术常用来开发构成通信系统仿真模型的算法，而且现代数字通信系统中的各种处理与数字信号处理技术本身就密切相关；多学科性与底层不同，通信系统的上层的主要建模与分析工具是排队论；估计理论和数理统计能够用于对特定仿真结果的有效性和可信度进行评估；通信系统仿真大部分情况下是随机仿真，其中的很多量都是随机变量，仿真要处理的信号和噪声等都是随机过程的样本函数，因此概率论和随机过程是通信系统仿真中的一种常用数学工具；多学科性通信系统仿真需要产生满足某些条件的伪随机变量，这些伪随机变量的产生方法

3、通常是以数论为基础的；与其他仿真一样，通信系统仿真中经常会涉及差分方程求解、曲线拟合和插值等操作，这些方法都源于数值分析。为了实现仿真模型，不可避免的会用到计算机科学方面的一些理论，程序设计理论、算法与数据结构理论的利用都是解决仿真模型实现问题的保证；现代通信系统仿真的规模和仿真量越来越大，分布式计算方法需要计算机科学中的分布式计算理论作指导。8主要内容通信系统仿真的多学科性通信系统仿真的层次性蒙特卡洛仿真原理二进制反极性基带通信系统的蒙特卡洛仿真蒙特卡洛仿真精度分析与仿真量确定方法方差缩减方法：重要性采样通信系统仿真的层次性特点通信系统的复杂性广义上讲，通信系统指的是遍布全球的通信网络，包括

4、有线通信网、移动通信网、卫星通信网、广播网和数据通信网等。每个网络都包括很多节点，节点之间通过底层链路和上层协议互相联系，构成一个非常复杂的系统。每个链路的情况又各不相同，具有不同的信道、不同的调制编码技术、不同的干扰情况、不同的器件等。通信系统仿真的层次性特点对如此复杂的通信系统用单个模型进行建模仿真显然是不现实的。通信系统仿真通常采用分层抽象仿真的方法。所谓分层抽象仿真是整个复杂系统的仿真分为若干个层次，分别评估系统的不同性能指标，在进行上层仿真时，下层通过相应的接口为上层提供一个抽象模型。由于各个层次仿真具有不同的特点和仿真目的，因此采用的仿真机制往往也各不相同。通信系统仿真的常见层次举

5、例网络（系统级）仿真链路级仿真同步单元仿真滤波器仿真射频电路仿真参数估计仿真网络（系统级）仿真网络（系统级）仿真用于数据包或信息在网络节点之间的流动和处理过程。主要评估指标包括网络的吞吐量、时延、时延抖动、丢包率和资源利用效率等。网络仿真常用语评估处理器配置、协议参数和缓冲器（Buffer）大小等资源配置的合理性。网络（系统级）仿真所采用的仿真机制应当是离散事件驱动的动态仿真机制。网络仿真需要链路级仿真为其提供一个抽象的接口模型。这个模型通常表征的是一些系统参数与误码率（或误码块率）之间的关系，而不关心链路级以下滤波、调制、编码、均衡等细节。链路级仿真链路级仿真 是与具体通信技术紧密相关的一个

6、仿真层次，用于仿真承载信息的波形在各种信道条件下的传输情况。链路级仿真模型通常需要包含调制解调、编译码、均衡等收发两端的各种信号处理过程。对数字通信系统而言，链路级仿真的评估指标通常是误码率、误帧率和误码块率等。链路级仿真通常采用的是同步数据流驱动的仿真机制。链路级仿真为系统仿真提供抽象接口。亚链路级仿真链路级以下的仿真层次变种比较多，也可以归属于链路级的子类。常见的有同步单元仿真、滤波器仿真、射频电路仿真和参数估计部分仿真等。这些仿真将一些相对独立的模块从链路级仿真中独立出来，用于评估一些模块独有的性能指标。亚链路级仿真同步单元仿真的评估指标包括时频跟踪精度、捕获时间、锁相环的失锁概率等。滤

7、波器仿真的评估指标主要是滤波器特性与设计指标的吻合程度。射频电路仿真的评估指标是电路的动态范围和线性范围等。参数估计部分仿真的性能指标则主要是估计参数的均方误差。亚链路级仿真该层次仿真任务的多样性使得它们所采用的仿真机制也不尽相同。如滤波器的仿真属于时间驱动的确定性仿真。同步部分的仿真采用时间驱动的连续系统级仿真机制或数据流驱动的仿真机制。参数估计部分则采用与链路级仿真类似的数据流驱动仿真机制。亚链路级仿真这些仿真能够为链路级仿真提供一些简化的抽象模型。如同步单元仿真可以为链路级仿真提供同步误差模型。滤波器仿真为链路级仿真提供传递函数模型。参数估计部分为链路级仿真提供估计误差模型。射频电路仿真

8、可以为链路级仿真提供放大器的非线性特性。这些模型的抽象能够在简化链路级仿真实现的同时保证一定的仿真准确性。19主要内容通信系统仿真的多学科性通信系统仿真的层次性蒙特卡洛仿真原理二进制反极性基带通信系统的蒙特卡洛仿真蒙特卡洛仿真精度分析与仿真量确定方法方差缩减方法：重要性采样通信系统仿真的共性特征通信系统仿真虽然种类繁多，方法各异，但绝大多数都是随机仿真，通常都采用蒙特卡洛仿真方法。蒙特卡洛方法概述蒙特卡洛(Monte Carlo)方法，又称随机抽样或统计模拟方法，泛指所有基于统计采样进行数值计算的方法。二战期间，这种方法首先被用于解决原子弹研制中的一个关键问题。后来N.Metropolis用驰

9、名世界的赌城摩纳哥的Monte Carlo来命名这种方法，为它蒙上了一层神秘色彩。蒙特卡洛方法概述随着现代计算机技术的飞速发展，蒙特卡洛方法已经在统计物理、经济学、社会学甚至气象学等方面的科学研究中发挥了极其重要的作用。将蒙特卡洛方法用于仿真即为蒙特卡洛仿真。蒙特卡洛方法概述蒙特卡洛方法适用于两类问题： -本身就具有随机性的问题； -能够转化为概率模型进行求解的确定性问题。蒙特卡洛方法求解问题的三个步骤： -构造或描述概率过程 -从已知概率分布抽样 -建立估计量蒙特卡洛方法求解问题步骤构造或描述概率过程实际上就是建立随机试验模型.构造概率过程是对确定性问题而言的.描述概率过程是对随机性问题而言

10、的.不同的问题所需要建立的随机试验模型各不相同。蒙特卡洛方法求解问题步骤所谓的从已知概率分布抽样指的是随机试验过程。随机模型中必须要包含某些已知概率分布的随机变量或随机过程作为输入。进行随机试验的过程就是对这些随机变量的样本或随机过程的样本函数作为输入产生相应输出的过程，因此通常被称为对已知概率分布的抽样。如何产生已知分布的随机变量或随机过程是蒙特卡洛方法中的一个关键问题。蒙特卡洛方法求解问题步骤最后一个步骤是获得估计量。蒙特卡洛方法所得到的问题的解总是对真实解的一个估计，本身也是一个随机变量，这个随机变量是由随机试验模型输出通过统计处理得到的。例：用蒙特卡洛方法估计 值1、构造随机试验模型可

11、以将 值估计问题归结为对四分之一圆面积的估计问题。随机试验模型：在正方形区域中随机撒点，统计落在1/4圆区域内的点数例：用蒙特卡洛方法估计 值例：用蒙特卡洛方法估计 值2、从已知概率分布抽样根据实验模型，需要在正方形区域内均匀撒点。这实际上就是要对一个二维均匀分布进行采样。这个二维均匀分布可以通过对两个独立均匀分布采样来实现。在计算机上实现，每次采样就是产生两个0,1区间上均匀分布的随机数。3、获得统计量例：用蒙特卡洛方法估计 值m=input(enter M, the number of experiments:);n=input(enter N, the number of trial:)

12、;z=zeros(1,m); %innitializationdata=zeros(n,m);for j=1:m x=rand(1,n); y=rand(1,n); k=0; for i=1:n if x(i)2+y(i)2Eb1-Eb判决+输出该模型中包括对两个随机过程的采样，一个是模拟二进制信源的贝努里分布，另一个是模拟最佳接收机输出端所包含高斯噪声。由于假设信源和信道均无记忆，且信源和信道相互独立，因此我们对这两个随机过程进行采样，只需要产生一系列相互独立、满足概率分布的随机变量即可。我们所关心的事件是判决结果与实际传送比特相反，建立的估计量使用错误事件的出现相对频率近似系统的误码率，即

13、：模型分析42源代码%Monte Carlo Simulation to estimate bit error rate performance of a binary%antipodal communication systemclcclearN = 100, 1000, 10000, 100000,1000000;EbN0_indB = -2:2:8;%Ebn0 = 10.(EbN0_indB/10); E=1;times = 50;err_rate = zeros(length(Ebn0),length(N),times); 43源代码for iiii=1:length(N) for i

14、ii=1:length(Ebn0) sigma = E/sqrt(2*Ebn0(iii); for ii = 1:10 r = rand(1,N(iiii); source = double(r=0.5); x=1-2*source; noise = randn(1,N(iiii)*sigma; y=E*x+noise; result = double(y1.e-6); err_rate(iii,iiii,ii) = error_num/N(iiii); end endend44相对误差的方差45现象可以看出，在同一信噪比条件下，仿真量越大相对误差的偏差越小；如果固定仿真量，相对误差的方差则随

15、着信噪比的提高而迅速增大。对于这种问题至少要做多少次仿真才能得到接近于真实解的估计值呢？4647主要内容通信系统仿真的多学科性通信系统仿真的层次性蒙特卡洛仿真原理二进制反极性基带通信系统的蒙特卡洛仿真蒙特卡洛仿真精度分析与仿真量确定方法方差缩减方法：重要性采样仿真量与精度之间关系的定性分析蒙特卡洛仿真方法的本质是在计算机上进行的随机试验和结果统计分析过程，试验次数越多，得到的数据样本就越多，根据这些样本所得到的统计结果精度和可信程度就越高，这一点非常容易理解。但是仿真量到底应该多大呢？48仿真量问题的关键：仿真精度的定量分析仿真量是我们在进行蒙特卡洛仿真方案设计时不能回避的一个问题。仿真量太小

16、，会造成仿真精度不能满足要求。仿真量过大则需要消耗大量的计算资源。因此往往需要根据仿真精度要求确定仿真量，即需要解决仿真精度的定量分析问题。49仿真精度的衡量指标仿真精度的两种衡量指标 -绝对精度和相对精度设数据的精确值是 ，通过仿真得到的估计值为 ， 是一个服从某种分布的随机变量。如果估计值 有概率1-a落在某一区间 上，称区间 为置信度1-a的置信区间，将置信区间长度的一半，即 称为绝对精度，绝对精度 与真值 的比值称为相对精度。50由置信度和绝对精度确定仿真次数在很多仿真场景中每次蒙特卡洛采样可以看做一次独立的贝努里试验。 -例如通信中传输一个数据符号，传输可能是正确的也可能是错误的；

17、-每次电话呼叫有可能接通也有可能呼叫阻塞； -通过随机试验法求圆周率时每次投下的点可能在四分之一圆内也可能在四分之一圆外等等。设一次独立的贝努里试验中事件E发生的概率为p，那么n次独立的贝努里试验中E发生的次数k服从二项分布，事件E发生的次数恰好为k次的概率是：51由置信度和绝对精度确定仿真次数如果以相对频率作为对概率p的估计，绝对精度要求小于 ，即k落在这个区间范围内的概率即置信度：给定置信度和绝对精度，可计算出所需要仿真的最少次数n。52由置信度和绝对精度确定仿真次数但是，直接按上式计算比较复杂，特别是当需要的次数n比较大时，式中的组合数计算难以实现，这种情况下可通过近似方法进行计算。53

18、由置信度和绝对精度确定仿真次数近似为正态分布根据大数定理，当试验次数 ，试验中事件发生次数k服从均值为np方差为np(1-p)的正态分布，即其中54由置信度和绝对精度确定仿真次数这样，给定置信度1-a和绝对精度以及事件的概率值p，就可以求解方程得出最小仿真次数n。如果概率值p未知，可用频率估计代替。55例子已知某通信系统的设计传输错误概率为0.001，为了至少有95%的把握使仿真得到的错误概率与真值之间的误差在2* 之内，问至少需要多少次仿真（即传输多少个独立符号）？delta=2e-4;p=1e-3;alpha=0.05;%显著性水平n=ceil(2*p*(1-p)/delta2*(erfi

19、nv(1-alpha)2)n=9594156泊松分布近似除了利用正态分布来近似之外，还可以用泊松分布来近似。泊松定理指出，在随机试验中事件的发生概率很小，随机试验次数很多的情况下，试验中事件的发生次数k近似服从参数为 的泊松分布，即F(x)是参数为np的泊松分布的分布函数。57例子在前一例子的仿真系统中，设计传输错误率为 ，置信区间为 ，若总试验次数（独立传输符号数）为95941次，求仿真结果的置信度（分别用泊松分布和正态分布近似）。delta=2e-4;%绝对误差p=1e-3;n=95941；P_delta_poiss=poisscdf(n*p+n*delta,n*p)-poisscdf(n

20、*p-n*delta,n*p);P_delta_norm=normcdf(n*p+n*delta,n*p)-normcdf(n*p-n*delta,n*p);58P_delta_poiss=0.9538P_delta_norm=0.9500由置信度和相对精度计算仿真次数在很多情况下，相对精度比绝对精度更有意义，因此在实际仿真中往往需要通过置信度和相对精度确定最小仿真次数。给定仿真的相对精度要求 ，则 ，将之代入到 并整理得到已知相对精度条件下的最小仿真次数为59由置信度和相对精度计算仿真次数反过来，如果给定仿真次数和置信度，则仿真结果的相对精度也可以计算出来 其中pn的物理意义是n次试验中事件

21、出现的平均次数。可以看出，在统计误码率时，平均误码数越多，统计结果的相对精度越高。60例子画出置信度为90%、95%和99%条件下试验中事件发生次数np与相对精度r之间的关系曲线。clear;alpha=0.1,0.05,0.01;pn=1 10 100 1000 10000 100000;for i=1:3 r(:,i)=sqrt(2./pn).*erfinv(1-alpha(i);end61例子f1=loglog(pn,r);set(f1(1),Marker,o);set(f1(1),Marker,square);set(f1(1),Marker,);legend(alpha=0.1, a

22、lpha=0.05, alpha=0.01);xlabel(多次试验中事件发生的次数np);xlabel(相对精度r);62事件发生次数与相对精度的关系图63分析如果要求试验结果的相对精度提高，那么就要使试验中观察得到事件发生次数成平方数量级增加。在事件发生概率较小的情况下，将导致总试验次数大大增多。这就是蒙特卡洛通常都需要较大数量仿真的原因。64例子一个通信系统，设传输错误概率很小，如果在仿真中每观察到10个、100个和1000个误码就进行一次误码率的统计，问得到的结果在95%置信度条件下的相对精度是多少？以统计次数代替平均出错次数alpha=0.05;err_num=10,100,1000;r=sqrt(2./err_num)*erfinv(1-alpha)65结论可见，对于误码率仿真统计而言，只有平均出错个数大于100，才能将相对误差超过20%的可能性控制在5%以内。因此，在实际通信系统仿真中，通常要求出错个数超过100.比如对传输误码率为 的通信系统进行仿真时，要求仿真量大于 。66蒙特卡洛仿真的优缺点蒙特卡洛方法能应用于一切系统中而无需考虑其结构和复杂度。但是，传统的