8PSK通信系统的蒙特卡罗仿真分析报告

1、设计要求:1、产生等概率且相互独立的二进制序列，画出波形;2、产生均值为0,方差为1的加性高斯随机噪声；3、进行8PSK调制，画出波形；4、进行蒙特卡罗分析；5、解调8PSK画出眼图。课程设计学生日志时间设计内容7.37.5翻阅Matlab、通信原理的有关资料7.6用Matlab编写并调试程序，实现调制、解调信号的时域波形7.7实行蒙特卡罗及眼图的仿真7.8写课程设计论文7.9修改论文7.10课程设计答辩课程设计考勤表周星期一星期二星期三星期四星期五2课程设计评语表指导教师评语:成绩：指导教师：年月日8PSK通信系统的蒙特卡罗仿真分析一、设计目的和意义（1）根据蒙特卡罗仿真方法基本思想，分析了

2、加性高斯噪声和单频干扰条件下的8PSK信号相位分布模型及其相关检测方法；（2）讨论了加性噪声和单频干扰模型，建立了相关检测系统蒙特卡罗仿真模型，模型通过动态系统仿真软件实现；（3）进行了仿真设计与应用系统性能分析。二、设计原理（1）蒙特卡罗概及原理蒙特卡罗法亦称为随机仿真方法，是一种与一般数值计算方法有本质区别的计算方法，它起源于早期的用几率近似概率的数学思想，它利用随机数学进行统计试验，以求得的统计特征值作为待解问题的数值解。蒙特卡罗法首先构造或描述各种概率过程，然后从这些概率过程中实现相应的概率分布，最后建立各种估计量，相当于对模拟实验的结果进行考察和登记，从中得到问题的解。对于数字通信系

3、统而言有效性用信息的传输速率表示而可靠性可用误码率来表示。信号在传送检波的过程中由于噪声的干扰可能会检码错误。这样如何统计分析信号传输过程中噪声干扰下的误码率问题就成了衡量该系统好坏的标志。蒙特卡罗方法正是基于随机数信号的产生原理，统计分析信号传输过程中由于噪声干扰而引起的误码率。其基本原理如下：由概率定义知，某事件的概率可以用大量试验中该事件发生的频率来估算，当样本容量足够大时，可以认为该事件的发生频率即为其概率。因此，可以先对影响其可靠度的随机变量进行大量的随机抽样，然后把这些抽样值一组一组地代入功能函数式，确定结构是否失效，最后从中求得结构的失效概率。蒙特卡罗法正是基于此思路进行分析的。

4、蒙特卡罗方法的基本思想是：当所求解问题是某种随机事件出现的概率，以这种事件出现的频率估计这一随机事件的概率，或者得到这个随机变量的某些数字特征,并将其作为问题的解.(OMSK)1QFSK)riftti图1通信系统仿真模型(2) 8PSK的调制8PSK是一种常用的多相键控，输入的二进制信息序列经串一并交换每次产生一个3位码组blb2b3，因此符号率为比特率的1/3。在b1b2b3控制下，同相路和正交路分别产生两个四电平基带信号l(t)和Q(t)。bl用于决定同相路的信号的极性，b2用于决定正交路信号的极性，b3则用于确定同相路和正交路信号的幅度。因此同相路和正交路的基带信号幅度是相互关联的，不能

5、独立选取。其调制与那里图如图1所示b、b、二进信息申一并执空换电平产生倒相靳Lg$电1F严生图28PSK调制（3）8PSK的解调8PSK言号可采用QPSKI号的相干解调器进行解调，区别在于判决电路二电平判决改为四电平判决，判决结果经逻辑运算后得到比特马组，再进行并/串变换，8PSK信号的另一种解调方案如图2所示，它采用两组正交相干解调器，其中一组参考载波相位为0和90，另一组参考载波相位为-45。和45,没个相干解调器后接一个二电平判决电路，对判决结果经逻辑运算后得到比特马组，在进行并/串变换，得到原始的串行二进制信息。8PSK调制是利用载波的8种不同相位（或相位差）来表征数字信息的调制方式，

6、和二进制调制一样，8PSK调制也分为绝对调相和相对调相。它把输入的二进制信号序列经过串并变换，每次把一个3位的码组映射为一个符号的相位，映射星座按自然码或者格雷码（文中采用自然码），符号速率是比特速率的1/3，在8PSK调制方式中，输入的串行二进制序列经串并转换，每次产生3位并行码组b1b2b3。bl和b3分别决定I路的极性和幅度，b2和b3分别决定Q路的极性和幅度，星座图如图4H100001图48PSK星座图4、眼图如图5所示，一个确定的数字序列在眼图空间对应固定的眼图结构。在不同干扰强度的情况下，基带信号对应的眼图结构会出现不同程度的变化，即干扰强度的不同会在眼图结构中得以体现，主要表现在

7、眼图的形状和亮、暗面积比会发生变化。基于眼图分析的干扰效果评估研究是一种数字信号基带传输受扰测度研究方法,主要思想既是通过眼图对基带信号的受扰程度进行表征。眼图的仿真模型如图6所示最佳取样时向最佳取样点零点位置的失真图5眼图模型图2askeyp-pstterndiagram(L9ILLLN13A1角1.11L7I.fiLS?lirtie/second图6仿真模型由图可知“眼睛”的张开程度可以作为基带传输系统性能的一种度量，它不但反映串扰的大小，而且也可以反映信道噪声的影响。眼图张开部分的宽度决定了接受波形可以不受串扰影响而抽样、再生的时间间隔。显然，抽样的最佳时刻是“眼睛”张开最大的时刻；“眼

8、睛”在特定抽样时刻的张开高度决定了系统的噪声容限；“眼睛”的闭合斜率决定了系统对抽样定时误差的敏感程度，斜率愈大则对定时误差愈敏感。三、详细设计步骤1、MATLA库中的高斯随机数发生函数randn是应用给出的最低标准随机数发生器以及极坐标法，将均匀分布的随机数映射成高斯分布的随机数。2、产生二进制序列后当进行8PS调制时伴随着高斯随机噪声，这时的信号便会有不同概率误差，我们就用蒙特卡罗来估计分析，存在的误差是怎么的，有多少,基于蒙特卡罗的思想原则和应用特征，结合通信中调制解调理论，利用蒙特卡罗思想在通信系统建模的应用，并给出仿真结果。蒙特卡罗法的步骤是：（1）构造实际问题的目标概率模型；（2）

9、根据目标概率模型的特点，确定概率密度函数和累积概率分布函数；（3）统计试验结果，给出问题的解和精度估计.仿真平台为matlab7.0,验证蒙特卡罗方法的好坏主要是能否把仿真结果确切显示出来.在计算信号误差函数时为了提高实验效果，在伪随机采样时就对概率密度函数和累积概率分布函数作好统计，这样得到的采样信号统计特征量，才能防止随着伪随机数数目增加而处理时间以指数级增加仿真结果见图10蒙特卡罗法的主要思路是：一、针对所要求解或所要验证的问题建立一个简单且便于实现的概率统计模型；二、对模型中的随机变量建立抽样方法，在计算机上进行随机模拟，抽取足够的随机数,并对有关的事件进行统计；三、对模拟试验结果加以

10、分析说明3、数字调相或相移键控PSK）最简单的形式为，利用=进制数字信号对两个同频反相正弦渡进行控制不断切换合成调相渡，即BPSK8BPSK是一种常用的八相相移键控，有八种可能输出的相位，它不是恒幅度调制，不仅在相位中而且在幅度中亦包台信息，一个符号携带三个比特的信息8PSK的调制信号可以看成是对两个正交载波进行四电平双边带调制后所得的两路四进制幅度键控4ASK信号的叠加，其实现如图3所示。4、在传输信道中引入噪声，解调出的信号以眼图的形式输出，改变干扰信号的强度,录取不同干扰强度下的眼图，利用眼图来观察其噪声带来的误码率。计算对应眼图的眼皮厚度与眼图开启度。改变信噪比的大小，相应的给出眼皮厚

11、度和眼图开启度的值及相应的受扰后的眼图。结合眼皮厚度与眼图开启度两个参数，建立二者的比值（眼图厚度比）与干扰强度的对应关系。通过大量的仿真数据获得拟合曲线。四、设计结果及分析图7二进制序列及8PSK波形FiyiLire1-inikJlriserr10居Qekioojindowtieti鼻1Q冈呂b|區富|口困o-S2一dFU.5Q6TimeQ.3OSTint：5)QQ9-n|x10.500.6-1-1.5-20d98G09S50JS90J99505Timefs)mpnadErVfl*T4*-44-*4-rl-l-PHVIVI1VI1i*tiI-BinifrI*IWWIRvSVg.w. I344

12、-曇*AAAI!41*VIV1R.aVI11vl11MiEIIIIRUi孚VVVV9VjVwrrvtvvHIX1I：!0;Transmitted_data_bits=data_2(1:(fs)/N:end)figure(1)subplot(2,1,1)plot(td,data2)axis(01-22);gridonxlabel(Time)ylabel(Amplitude)title(输入数据)twopi_fc_t=(1:fs/2)*2*pi*fcarr/fs;a=1;phi=0;cs_t=a*cos(twopi_fc_t+phi);sn_t=a*sin(twopi_fc_t+phi);cs_t

13、=cs_t;sn_t=sn_t;si=cs_t.*lsymbols;sq=sn_t.*Qsymbols;sumiq=si+sq;gridonsubplot(2,1,2)plot(tiq,sumiq)axis(.498.502-22);gridonxlabel(Time(s)ylabel(Amplitude)title(8PSK波形);M=2;seed=7864554321;SNRpBit=60;SNR=SNRpBit/log2(M);numPlot=25;source=randsrc(numSymb,1,0:M-1);n=randn(1,100)*sqrt(1);figure(2)stem(n

14、)title(高斯随机噪声)figure(4)subplot(1,1,1)plot(tiq,sumiq1)axis(.498.502-22);gridonxlabel(Time(s)ylabel(Amplitude)title(叠加噪声后的8PSK波形)figure(5)xlabel(Time)ylabel(Voltage)out=;fori=1:1000:10000out=outbitout(i:i+999)bitout1(i:i+999);endout1=;fori=1:1000:20000ifsum(out(i:i+999)0outa(i:i+999)=1;out1=out11;else

15、outa(i:i+999)=-1;out1=out1-1;endendplot(td,outa)title(解调后的8PSK信号);gridon;axis(01-22);N=fs/fcarr;functionpb,ps=cm_sm32(snr_in_dB)N=10000;E=1;snr=10A(snr_in_dB/10);sgma=sqrt(E/(3*2*snr);s000=10;s001=cos(pi/4)sin(pi/4);s01仁01;s010=cos(3*pi/4)sin(3*pi/4);s110=-10;s11仁cos(5*pi/4)sin(5*pi/4);s101=0-1;s100

16、=cos(7*pi/4)sin(7*pi/4);fori=1:N,temp=rand;if(temp0.125),dsource1(i)=0;dsource2(i)=0;dsource3(i)=0;elseif(temp0.25),dsource1(i)=0;dsource2(i)=0;dsource3(i)=1;elseif(temp0.375),dsource1(i)=0;dsource2(i)=1;dsource3(i)=0;elseif(temp0.5),dsource1(i)=0;dsource2(i)=1;dsource3(i)=1;elseif(temp0.625),dsourc

17、e1(i)=1;dsource2(i)=0;dsource3(i)=0;elseif(temp0.75),dsource1(i)=1;dsource2(i)=0;dsource3(i)=1;elseif(temp0.875).dsource1(i)=1;dsource2(i)=1;dsource3(i)=0;elsedsource1(i)=1;dsource2(i)=1;dsource3(i)=1;end;end;numofsymbolerror=0;numofbiterror=0;fori=1:N,n(1)=gngauss(sgma);n(2)=gngauss(sgma);if(dsourc

18、e1(i)=0)&(dsource2(i)=0)&(dsource3(i)=0),r=s000+n;elseif(dsource1(i)=0)&(dsource2(i)=0)&(dsource3(i)=1),r=s001+n;elseif(dsource1(i)=0)&(dsource2(i)=1)&(dsource3(i)=0),r=s010+n;elseif(dsource1(i)=0)&(dsource2(i)=1)&(dsource3(i)=1),r=s011+n;elseif(dsource1(i)=1)&(dsource2(i)=0)&(dsource3(i)=0),r=s100+

19、n;elseif(dsource1(i)=1)&(dsource2(i)=0)&(dsource3(i)=1),r=s101+n;elseif(dsource1(i)=1)&(dsource2(i)=1)&(dsource3(i)=0).r=s110+n;elser=s111+n;end;cOOO=dot(r,sOOO);cOO仁dot(r,s001);cO1O=dot(r,sO1O);cO11=dot(r,sO11);c1OO=dot(r,s1OO);c1O1=dot(r,s1O1);c11O=dot(r,s11O);c111=dot(r,s111);c_max=max(cOOOcOO1cO

20、1OcO11c1OOc1O1c11Oc111);if(cOOO=c_max),decis1=O;decis2=O;decis3=O;elseif(cOO1=c_max),decis1=O;decis2=O;decis3=1;elseif(cO1O=c_max),decis1=O;decis2=1;decis3=O;elseif(cO11=c_max),decis1=O;decis2=1;decis3=1;elseif(c1OO=c_max),decis1=1;decis2=O;decis3=O;elseif(c1O1=c_max),decis1=1;decis2=O;decis3=1;else

21、if(c11O=c_max),decis1=1;decis2=1;decis3=O;elsedecis1=1;decis2=1;decis3=1;end;symbolerror=0;if(decis1=dsource1(i),numofbiterror=numofbiterro叶1;symbolerror=1;end;if(decis2=dsource2(i),numofbiterror=numofbiterro叶1;symbolerror=1;end;if(decis3=dsource3(i),numofbiterror=numofbiterro叶1;symbolerror=1;end;if

22、(symbolerror=1),numofsymbolerror=numofsymbolerror+1;end;end;ps=numofsymbolerror/N;pb=numofbiterror/(3*N);functionchannelout=pskmoto(input,No)u=input;sgma=sqrt(No);%经过高斯信道fori=1:length(u)n=gngauss(sgma);if(dsource1(i)=0)&(dsource2(i)=0)&(dsource3(i)=0)r=s000+n;elseif(dsource1(i)=0)&(dsource2(i)=0)&(d

23、source3(i)=1)r=sOO1+n;elseif(dsource1(i)=0)&(dsource2(i)=1)&(dsource3(i)=0)r=sO1O+n;elseif(dsource1(i)=0)&(dsource2(i)=1)&(dsource3(i)=1)r=s011+n;elseif(dsource1(i)=1)&(dsource2(i)=0)&(dsource3(i)=0)r=s100+n;elseif(dsource1(i)=1)&(dsource2(i)=0)&(dsource3(i)=1)r=s101+n;elseif(dsource1(i)=1)&(dsource

24、2(i)=1)&(dsource3(i)=0)r=s110+n;elser=s111+n;end%计算互相关度c000=dot(r,s000);c001=dot(r,s001);c010=dot(r,s010);c011=dot(r,s011);c100=dot(r,s100);c101=dot(r,s101);c110=dot(r,s110);c111=dot(r,s111);%对第i个符号进行判决c_max=max(c000c001c010c011c100c101c110c111);if(c000=c_max)decis1(i)=0;decis2(i)=0;decis3(i)=0;elseif(c001=c_max)decis1(i)=0;decis2(i)=0;decis3(i)=1