cha3蒙特卡罗-课件课件.pptVIP

郑宏云njtu1第三章蒙特卡罗仿真方法

郑宏云NJTU 7/19/2025

郑宏云njtu2蒙特卡罗方法的基本概念蒙特卡罗方法的应用蒙特卡罗积分蒙特卡罗方法仿真通信系统

郑宏云njtu3一、蒙特卡罗方法基本概念蒙特卡罗（MonteCarlo，MC）方法，又称随机抽样或统计试验方法利用随机数进行数值模拟的方法蒙特卡罗仿真是指利用蒙特卡罗方法估计系统参数的仿真蒙特卡罗估计则是指通过内在随机试验来估计参数值的过程

郑宏云njtu41、基本思想例1：产品合格率。抽查N个产品，n个合格品，合格率pN数目越大，n/N作为合格率估计值越准确

郑宏云njtu5相对频率随机实验和事件A重复N次随机实验，记录事件A发生次数NA计算事件A发生的相对频率事件A发生的概率P（A）=是P（A）的估计值,因为Ninf是随机变量，决定了估计的质量

郑宏云njtu6MC基本思想 当要求解的问题是某种事件出现的频率，或者是某个随机变量的期望值时，可以通过某种“试验”的方法，得到这种事件出现的概率，或这个随机数的平均值，再进行统计平均求其近似值，作为问题的解。抓住事物运动的几何数量和几何特征，利用数学方法来模拟；以一种概率模型为基础，按照该模型描绘的过程，通过模拟实验的结果，作为问题的近似解

郑宏云njtu72、无偏和一致性估计无偏估计平均意义下,估计值是正确的一致性估计估计值具有小方差如果是无偏的,估计值聚集在真值的周围,且具有较小的散布范围,估计值有用的具有无偏和一致性性质的估计误差零均值，且方差在N趋于inf时收敛于0

郑宏云njtu83、蒙特卡罗估计例：估计Pir=1正方形面积Ab=1饼区面积Apie=(1/4)*[pi*r^2]=pi/4落点均匀分布，落入正方形和饼区的数量之比，近似等于两个区的面积之比，即

郑宏云njtu9MC估计的重要特征测试条件2个计数器

郑宏云njtu10二、MC的应用--仿真通信系统AWGN信道,独立丢失事件“N个字符中出现Ne个错误”服从二项分布误码率估计Q:无偏?一致?PE:单次发送时的差错概率

郑宏云njtu11无偏，一致注意：假设条件

郑宏云njtu12数字通信系统MC仿真假设条件发送端没有脉冲成形AWGN信道数据序列独立同分布没有滤波——无码间干扰

郑宏云njtu13通信系统仿真模型数据信源发送机比较高斯白噪声接收机延时高斯白噪声同相分量和正交分量表示法:描述的是信号的信号空间分量,而不是时域波形的采样点;好处是信号空间分量只要用每个发送符号的一个样点就可给出，运行快

郑宏云njtu14通信系统仿真模型(续)为确定以Eb/N0为函数的BER,固定Eb,改变N0归一化fs和Eb复包络用于确定仿真中的噪声标准偏差

郑宏云njtu15例1BPSKd[n]=0d[n]=10数据信源发送机比较高斯白噪声接收机高斯白噪声Ac=1,km=pi

郑宏云njtu16BPSK

郑宏云njtu17%file:c9_MCBPSK.msnrdB_min=-3;snrdB_max=8;snrdB=snrdB_min:1:snrdB_max;Nsymbols=input(Enternumberofsymbols);snr=10.^(snrdB/10);h=waitbar(0,SNRIteration);len_snr=length(snrdB);forj=1:len_snr%waitbar(j/len_snr);sigma=sqrt(1/(2*snr(j)));error_count=0;fork=1:Nsymbolsd=round(rand(1));%datax_d=2*d-1;%tansmitteroutputn_d=sigma*randn(1);%noisey_d=x_d+n_d;%receiverinputify_d0d_est=1;elsed_est=0;end

郑宏云njtu18if(d_est~=d)error_count=error_count+1;%errorcounterendend%simulationloopendserrors(j)=error_cou