1、完整word版8PSK通信系统的蒙特卡罗仿真分析(完整word版)8PSK通信系统的蒙特卡罗仿真分析8PSK通信系统的蒙特卡罗仿真分析一、设计原理 1. 设计整体方案根据题目要求和对设计原理的分析，设计的整体方案为：首先，产生一个等概率且相互独立的二进制序列，将二进制序列通过串-并转换为一个三位码组。通过调用Matlab系统函数，对输入信号进行调制和滤波，再送入信道传输，在在信道中对信号采用加性高斯白噪声进行干扰，然后进行解调，同时计数误比特和误符号数，并根据蒙特卡罗方法统计分析信号传输过程中由于噪声干扰作用下的误比特率和误符号率。最后画出各种波形。 图1 8PSK通信系统的蒙特卡罗仿真分析2
2、. 蒙特卡罗法的原理蒙特卡罗法又称随机抽样或统计试验方法，属于计算数学的一个分支，它是在本世纪四十年代中期为了适应当时原子能事业的发展而发展起来的。传统的经验方法由于不能逼近真实的物理过程，很难得到满意的结果，而蒙特卡罗方法 由于能够真实地模拟实际物理过程，故解决问题与实际非常符合，可以得到很圆满的结果。这也是我们采用该方法的原因。 蒙特卡罗方法的基本原理及思想如下： 当所要求解的问题是某种事件出现的概率，或者是某个随机变量的期望值时，它们可以通过某种“试验”的方法，得到这种事件出现的频率，或者这个随机变数的平均值，并用它们作为问题的解。这就是蒙特卡罗方法的基本思想。蒙特卡罗方法通过抓住事物运
3、动的几何数量和几何特征，利用数学方法来加以模拟，即进行一种数字模拟实验。它是 以一个概率模型为基础，按照这个模型所描绘的过程，通过模拟实验的结果，作为问题 的近似解。二、详细设计步骤根据整体设计方案，对各个设计模块进行具体设计。设计步骤如下：1 二进制序列的产生、串/并转换及二/十进制转换二进制的产生即是本系统的信息源的产生，串/并转换和二/十进制转换以备在调制的时候数据调用。在本系统中，首先，用Matlab库函数来直接产生一个串行二进制序列，并将其保存在一个向量中，以备函数调用；其次，将产生的串行二进制码转换成3行并行码，即是一个符号；再次，将二进制序列转换成十进制数（0、1、2、7）；最后
4、，调用Matlab绘图函数绘制等概率且相互独立的二进制序列波形。模块程序如下： %数据源的产生x=round(rand(1,3*N); %产生独立串行随机二进制for i=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)=bin2dec( num2str(D) );%二十进制转换 end%#以下是画出独立随机二进制序列程序for i=1:length(x), %计算码元的值 if(x(i)=1), %如果信息为1 for j0=1
5、:grid0, %该码元对应的点值取1 y(i-1)*grid0+j0)=1; end else for j0=1:grid0, %反之，信息元为0，码元对应点值取0 y(i-1)*grid0+j0)=0; end;end;endfigure(2); %图2plot(y,b); %画出独立随机二进制序列title(信号源);xlabel(X序列);ylabel(Y序列);axis(0 20000 -0.5 2) ;28psk调制、滤波、添加加性高斯白噪声及相应波形产生本系统中，首先，8psk调制、滤波通过调用Matlab库函数y=dmodce(x, Fd, Fs, method, )，从而在输
6、出端输出8psk波形（眼图）；其次，通过调用子函数产生高斯白噪声；最后，在加入高斯白噪声后，绘制受噪声干扰的星座图和眼图。模块程序如下：for SNR_db=0:10:20N3=0.1*SNR_db;Eb_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
7、,M); %8PSK解调figure(N3+3); % 循环画出不同加加性高斯随机噪声后的星座图axis(-1.2 1.2 -1.2 1.2);hold on;for i=0:M-1 plot(cos(2*pi*i/M),sin(2*pi*i/M), .,MarkerSize ,20); %产生星座图endplot(ynoise, +); %画出加噪声后的星座图title(加噪声后不同星座图);hold onfor j1=1:Nn1(N*N3+j1)=y0(j1); %存调制后的数据于n1n2(N*N3+j1)=ynoise(j1);%存加噪声后的数据于n2end %画出受到噪声干扰的信号向量
8、；end;hold ondelay = 3; % 提升余弦滤波器的延时 rcv = rcosflt(n1,Fd,Fs,fir/normal,.5,delay); rcvy = rcosflt(n2,Fd,Fs,fir/normal,.5,delay); propdelay = delay .* Fs/Fd + 1; % 繁殖延时 rcv1 = rcv(propdelay:end-(propdelay-1),:); rcvy1 = rcvy(propdelay:end-(propdelay-1),:); N0 = Fs/Fd; offset1 = 0;% 偏移值为2 h1 = eyediagra
9、m(rcv1,N0,1/Fd,offset1); set(h1,Name,Eye Diagram Displayed with No Offset); title(加噪前的眼图);h2 = eyediagram(rcvy1,N0,1/Fd,offset1); set(h2,Name,Eye Diagram Displayed with No Offset); title(加噪后的眼图); %# 以下是画出加加性高斯随机噪声的波形程序 B=10*N;for i=1:B n00(i)=10*gngauss(sgma);endfigure(8);hold on;for i=1:B plot(n00)
10、;end;title(加性高斯随机噪声);xlabel(时间t);ylabel(幅度y);axis(0 B -2.5 2.5);3加性高斯白噪声的产生及其波形的产生和绘制 本模块，调用一个子函数，产生均值为0，方差为1的高斯白噪声，并且绘出相应的高斯白噪声波形。模块程序如下：%# 以下是画出加加性高斯随机噪声的波形程序 B=10*Nfor i=1:B n00(i)=10*gngauss(sgma);endfigure(8);hold on;for i=1:B plot(n00);end;title(加性高斯随机噪声);xlabel(时间t);ylabel(幅度y);axis(0 B -2.5
11、2.5);48PSK的解调 本模块是将加入加性高斯随机噪声的信号与可能发送的8种发射信号相位向量am（m=0，1，2，3，4，5，6，7）进行相关，再通过判决器取出相关值最大的那组am，作为检测器的输出信号，完成相干检测，将该信号进行并串变换，其结果与原信号d进行比较。模块程序如下：%信号映射a000=1 0;a001=cos(pi/4) sin(pi/4);a011=0 1;a010=cos(3*pi/4) sin(3*pi/4);a110=-1 0;a111=cos(5*pi/4) sin(5*pi/4);a101=0 -1;a100=cos(7*pi/4) sin(7*pi/4); fo
12、r i=1:N, n(1)=gngauss(sgma); n(2)=gngauss(sgma); %存高斯白噪声数据于n中 if(d1(i)=0)&(d2(i)=0)&(d3(i)=0), %加噪声 r=a000+n; elseif(d1(i)=0)&(d2(i)=0)&(d3(i)=1), r=a001+n; elseif(d1(i)=0)&(d2(i)=1)&(d3(i)=0), r=a010+n; elseif(d1(i)=0)&(d2(i)=1)&(d3(i)=1), r=a011+n; elseif(d1(i)=1)&(d2(i)=0)&(d3(i)=0), r=a100+n; el
13、seif(d1(i)=1)&(d2(i)=0)&(d3(i)=1), r=a101+n; elseif(d1(i)=1)&(d2(i)=1)&(d3(i)=0), r=a110+n; else r=a111+n; end; %对应相关函数 c000=dot(r,a000); c001=dot(r,a001); c010=dot(r,a010); c011=dot(r,a011); c100=dot(r,a100); c101=dot(r,a101); c110=dot(r,a110); c111=dot(r,a111); c_max=max(c000 c001 c010 c011 c100 c
14、101 c110 c111); %判决器取出相关值最大，作为检测器的输出信号 if(c000=c_max), de1=0;de2=0;de3=0; %将该信号进行并串变换 elseif(c001=c_max), de1=0;de2=0;de3=1; elseif(c010=c_max), de1=0;de2=1;de3=0; elseif(c011=c_max), de1=0;de2=1;de3=1; elseif(c100=c_max), de1=1;de2=0;de3=0; elseif(c101=c_max), de1=1;de2=0;de3=1; elseif(c110=c_max),
15、 de1=1;de2=1;de3=0; else de1=1;de2=1;de3=1; end;5蒙特卡罗分析本模块完成对加高斯白噪声后的误码和误符号的概率计算，以及绘制出相应的蒙特卡罗分析图。首先，是检测加高斯白噪声后的对应码元是否错误，如果出错，则误码和误符号计数加1，算出总共的误码数和误符号数；然后，由误码数和误符号数计算出误码率和误符号率；最后，绘出8PSK的蒙特卡洛图。模块程序如下：for k=1:length(SNRindB1),snr=10(SNRindB1(k)/10); %信噪比sgma=sqrt(E/(3*2*snr); % %噪声标准均方差numofsymbolerror
16、=0; %无符号率计数初值numofbiterror=0; %无码率计数初值symbolerror=0; %误码标志 if(de1=d1(i), numofbiterror=numofbiterror+1; %有误码则加一 symbolerror=1; %置误码标志 end; if(de2=d2(i), numofbiterror=numofbiterror+1; symbolerror=1; end; if(de3=d3(i), numofbiterror=numofbiterror+1; symbolerror=1; end; if(symbolerror=1), numofsymbole
17、rror=numofsymbolerror+1; end;end;ps=numofsymbolerror/N; %计算误码率pb=numofbiterror/(3*N); %计算误符号率 smld_symbol_err_prb(k)=ps; %将误码率值存入向量smld_symbol_err_prb中 smld_bit_err_prb(k)=pb; %将误符号率值存入向量 smld_bit_err_prb中end;%#以下是画出理论误码率，误比特率和误符号率的程序for i=1:length(SNRindB2), SNR=exp(SNRindB2(i)*log(10)/10); %对数转换 %
18、信噪比 theo_err_prb(i)=2*Qfunct(sqrt(6*SNR)*sin(pi/8)/3;%理论比特误码率 end;figure(6);semilogy(SNRindB1,smld_bit_err_prb,*); %“*”画出误b比特率对应曲线hold onsemilogy(SNRindB1,smld_symbol_err_prb,O); %“o”画出误符号率对应曲线semilogy(SNRindB2,theo_err_prb); %“”画出理论比误码率曲线 axis(0,10,10(-5),100); %设置横纵轴间距大小gridxlabel(Eb/No in dB) %x轴
19、名称ylabel(error probability) %y轴名称title(8PSK的蒙特卡洛图)三、设计结果及分析1 通过对二进制序列产生、串/并转换及二/十进制转换模块程序的调试绘出波形如图2： 图2 等概率且相互独立的二进制波形2通过调试绘制出调制输出的8psk波形（如图3）、加高斯白噪声前的眼图（如图4）、加高斯白噪声后的眼图（如图5）以及加高斯白噪声后的星座图（如图6）： 图3 调制后输出8PSK波形 图4 加高斯白噪声前的眼图 图5 加高斯白噪声后的眼图 图6 加高斯白噪声后的星座图3通过对加性高斯白噪声的产生及其波形的产生和绘制模块程序的调试，绘出高斯随机噪声的波形如图6. 图
20、7 加性高斯随机噪声波形4. 通过对8psk蒙特卡罗分析模块程序的调试得到波形如图7. 图8 8PSK的蒙特卡罗图 经过对系统设计的分析，程序设计和调试，得到了正确的结果。本设计以产生一个等概率且相互独立的二进制序列，经串并变换形成三路并行的二进制符号，三路二进制符号经过8PSK调制，映射为相应的二维相位向量，与干扰噪声向量叠加后，形成受污染的接收信号，输出的接收信号以向量的形式表达，在检测器中与可能发送的8种发射信号相位向量an（n=0，1，2，3，4，5，6，7）进行相关，判决器取出相关值最大的那组an，作为检测器的输出信号，完成相干检测，将该信号进行并串变换，其结果与原信号d进行比较。如果进行N次试验，符号不一致的次数计数为M，则符号错误概率Ps的估值为M/N。对8PSK系统，如果发送N个符号，相当于发送了3N个比特。总的比特错误次数计数为M，则比特误码率Pb的估值为M/2N。最后通过图形进行形象的分析误比特率和误符号率。
