8psk通信系统的蒙特卡罗仿真分析

《8psk通信系统的蒙特卡罗仿真分析》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、8PSK通信系统的蒙特卡罗仿真分析一、设计原理1.设计整体方案根据题目要求和对设计原理的分析，设计的整体方案为：首先，产生一个等概率且相互独立的二进制序列，将二进制序列通过串-并转换为一个三位码组。通过调用Matlab系统函数，对输入信号进行调制和滤波，再送入信道传输，在在信道中对信号采用加性高斯白噪声进行干扰，然后进行解调，同时计数误比特和误符号数，并根据蒙特卡罗方法统计分析信号传输过程中由于噪声干扰作用下的误比特率和误符号率。最后画出各种波形。图18PSK通信系统的蒙特卡罗仿真分析2.蒙特卡罗法的原理蒙

2、特卡罗法又称随机抽样或统计试验方法，属于计算数学的一个分支，它是在本世纪四十年代中期为了适应当时原子能事业的发展而发展起来的。传统的经验方法由于不能逼近真实的物理过程，很难得到满意的结果，而蒙特卡罗方法由于能够真实地模拟实际物理过程，故解决问题与实际非常符合，可以得到很圆满的结果。这也是我们采用该方法的原因。蒙特卡罗方法的基本原理及思想如下：当所要求解的问题是某种事件出现的概率，或者是某个随机变量的期望值时，它们可以通过某种“试验”的方法，得到这种事件出现的频率，或者这个随机变数的平均值，并用它们作为问题的

3、解。这就是蒙特卡罗方法的基本思想。蒙特卡罗方法通过抓住事物运动的几何数量和几何特征，利用数学方法来加以模拟，即进行一种数字模拟实验。它是以一个概率模型为基础，按照这个模型所描绘的过程，通过模拟实验的结果，作为问题的近似解。一、详细设计步骤根据整体设计方案，对各个设计模块进行具体设计。设计步骤如下：1．二进制序列的产生、串/并转换及二/十进制转换二进制的产生即是本系统的信息源的产生，串/并转换和二/十进制转换以备在调制的时候数据调用。在本系统中，首先，用Matlab库函数来直接产生一个串行二进制序列，并将其保

4、存在一个向量中，以备函数调用；其次，将产生的串行二进制码转换成3行并行码，即是一个符号；再次，将二进制序列转换成十进制数（0、1、2、…、7）；最后，调用Matlab绘图函数绘制等概率且相互独立的二进制序列波形。模块程序如下：%数据源的产生x=round(rand(1,3*N));%产生独立串行随机二进制fori=1:N,D(1)=x(3*i-2);%D暂存二进制数以待转换D(2)=x(3*i-1);D(3)=x(3*i);d1(i)=D(1);d2(i)=D(2);d3(i)=D(3);Bin2Dec(i

5、)=bin2dec(num2str(D));%二十进制转换end%##################以下是画出独立随机二进制序列程序fori=1:length(x),%计算码元的值if(x(i)==1),%如果信息为1forj0=1:grid0,%该码元对应的点值取1y((i-1)*grid0+j0)=1;endelseforj0=1:grid0,%反之，信息元为0，码元对应点值取0y((i-1)*grid0+j0)=0;end;end;endfigure(2);%图2plot(y,'b');%画出独立随

6、机二进制序列title('信号源');xlabel('X序列');ylabel('Y序列');axis([020000-0.52]);2．8psk调制、滤波、添加加性高斯白噪声及相应波形产生本系统中，首先，8psk调制、滤波通过调用Matlab库函数y=dmodce(x,Fd,Fs,method,…)，从而在输出端输出8psk波形（眼图）；其次，通过调用子函数产生高斯白噪声；最后，在加入高斯白噪声后，绘制受噪声干扰的星座图和眼图。模块程序如下：forSNR_db=0:10:20N3=0.1*SNR_db;E

7、b_N0=10^(SNR_db/10);sgma=sqrt(1/(8*Eb_N0));%nc、ns的均方差y0=dmodce(Bin2Dec,Fd,Fs,'psk',M);%8PSK调制ynoise=sqrt(Fs/Fd)*sgma*(randn(length(y0),1)+j*randn(length(y0),1));%加加性高斯随机噪声%z=ddemodce(ynoise,Fd,Fs,'psk',M);%8PSK解调figure(N3+3);%循环画出不同加加性高斯随机噪声后的星座图axis([-1.2

8、1.2-1.21.2]);holdon;fori=0:M-1plot(cos(2*pi*i/M),sin(2*pi*i/M),'.','MarkerSize',20);%产生星座图endplot(ynoise,'+');%画出加噪声后的星座图title('加噪声后不同星座图');holdonforj1=1:Nn1(N*N3+j1)=y0(j1);%存调制后的数据于n1n2(N*N3+j1)=ynoise(j1);