《实验名称 4PSK通信系统的蒙特卡洛仿真.docx》由会员分享，可在线阅读，更多相关《实验名称 4PSK通信系统的蒙特卡洛仿真.docx（10页珍藏版）》请在三一办公上搜索。
1、通信系统的数字调制实验名称：4PSK通信系统的蒙特卡洛仿真（包含在无线通信的数字调制技术实验里面）一. 实验目的1）学会运用MATLAB分析通信系统。2）了解蒙特卡洛仿真原理。3）掌握调相原理。4）学会分析系统的误码率及误比特率。5）理论联系实际，能运用所学知识分析具体问题。二. 实验仪器：1）PC 机。2）采用matlab或其它软件工具编程实现三. 实验意义在数字通信中，有些场合可以不经过载波调制和解调过程而让基带信号直 接进行传输。称为基带传输系统。与之对应，把包括了载波调制和解调过程的 传输系统称为频带传输系统。无论是基带传输还是频带传输，基带信号处理是 必须的组成部分。因此掌握数字基带
2、传输的基本理论十分重要，它在数字通信 系统中具有普遍意义。PSK是利用载波的不同相位表示相应的数字信息。对于M相相位调制来说 M-2K，这里K是每个传输符号的信息比特数4PSK是M=4的载波相位调制。 蒙特卡洛法的实质是利用服从某种分布的随机数来模拟现实系统中可能出现的 随机现象。其基本思想是：为了求解数学、物理、工程技术以及生产管理等方 面的问题，首先建立一个概率模型或随机过程，使它的参数等于问题的解；然 后通过对模型或过程的观察或抽样试验来计算所求随机参数的统计特征，最后 给出所求解的近似值，解的精度可用估计值的标准误差来表示。四. 实验原理1、移相键控移相键控即受键控的载波相位调制是按基
3、带脉冲改变的一种数字调制方 式。其中，四相移相键控制（4PSK）的应用广泛，它是用4种不同相位代表4种 不同相位的信息，因此对于输入的二进制数字序列应该先分组，将每两个比特 编为一组；然后用4种不同的相位对其表征。例如，若输入的二进制数字信息 序列为，则可将他们分成10, 11，00, 10,，然后用4种不同的相位对其 表征。采用相位选择法产生4PSK信号，以实现4PSK调制器的设计。2、MPSK信号所谓多进制数字调制，就是利用多进制数字基带信号去调制高频载波的某 个参量，如幅度、频率或相位的过程。根据被调参量的不同，多进制数字调制 可分为多进制幅度键控（MASK）、多进制频移键控（MFSK）
4、以及多进制相移 键控（MPSK或MDPSK）。也可以把载波的两个参量组合起来进行调制，如把 幅度和相位组合起来得到多进制幅相键控（MAPK）或它的特殊形式多进制正 交幅度调制（MQAM）等。由于多进制数字已调信号的被调参数在一个码元间隔内有多个取值，因 此，与二进制数字调制相比，多进制数字调制有以下几个特点：（1）在码元速率（传码率）相同条件下，可以提高信息速率（传信率）， 使系统频带利用率增大。码元速率相同时，血进制数传系统的信息速率是二进 制的汨倍。在实际应用中，通常取仙=2七*为大于1的正整数。（2）在信息速率相同条件下，可以降低码元速率，以提高传输的可靠性。 信息速率相同时，M进制的码
5、元宽度是二进制的1吧财倍，这样可以增加每个 码元的能量，并能减小码间串扰影响等。正是基于这些特点，使多进制数字调制方式得到了广泛的使用。不过， 获得以上几点好处所付出的代价是，信号功率需求增加和实现复杂度加大。MPSK信号将M进制符号与M个载波相位对应，可以写成如下形式：队（Q = Re血/&=g（t）cos 2” 十常= g（/）cos |（wy）cos 2”盘g（E）sin 咨（小-）成n 2fct五、设计步骤1、4PSK调相波的产生当传双比特码元00时，发送起始相位为0的载波；当传双比特码元10时，发送起始相位为n/2的载波；当传双比特码元11时，发送起始相位为n的载波；当传双比特码元0
6、1时，发送起始相位为3n/2的载波。按这种定义方式定义的4PSK系统称n/2系统（也称A方式）。对4PSK信号还有另外一种定义方式：称n/4系统或B方式。具体定义如下：当传双比特码元00时，发送起始相位为45的载波；当传双比特码元10时，发送起始相位为135的载波；当传双比特码元11时，发送起始相位为225的载波；当传双比特码元01时，发送起始相位为315的载波。（1）调相法调相法产生QPSK,由两路相互正交的2PSK相加构成。串/并变换器将输入的二进制序列变成两路并行的双极性序列，将这两路信 号分别进行2PSK调制后再相加，这样就得到的四相移相信号。每一个支路信 号的速率是原始基带信号速率的
7、一半。二进制数码一串/# 变换（1）（b）图1调相法（2）相位选择法用相位选择法产生QPSK信号。四相载波发生器分别送出调相所需的四种 不同相位的载波。按照串/并变换器输出的双比特码元的不同，逻辑选相电路输 出响应相位的载波。图2相位调相法（3）QPSK信号的解调四相调相信号可以看作是两个正交2PSK信号的合成，所以它可以采用与 2PSK信号相似的解调方法进行解调，即由两个2PSK信号相干解调器构成。图3相干解调2、蒙特卡洛仿真套胃 srfc sin 式因而得到广4PSK信号采用相干解调时系统的误码率为 们中，产为信噪比。多相制是一种频带利用率较高的高效率传输方式。再加之有较好的抗噪声性能，
8、泛的应用，而MDPSK比MPSK用得更广泛一些。for i=1:length (SNRindBl),pb , ps=cm_sm32 (SNRindBl (i);smld_bit_err_prb (i)=pb;smld_symbol_err_prb (i)=ps;end;for i=1:length (SNRindB2),SNR=exp(SNRindB2(i)*log(10)/10);theo_err_prb(i)=Qfunct(sqrt(2火SNR);仿真比特和符号比特率，信噪比，理论比特误码率分别画出曲线。六、设计结果及分析1、4PSK调相波产生图4 4PSK的调相波2、蒙特卡洛仿真图5 4
9、PSK蒙特卡洛仿真附：程序代码%对M=4的PSK通信系统进行蒙特卡罗仿真%clear all echo off SNRindB1=0:2:10; SNRindB2=0:0.1:10;for i=1:length (SNRindB1),pb , ps=cm_sm32 (SNRindB1 (i);仿真比特和符号误码率smld_bit_err_prb (i)=pb; smld_symbol_err_prb (i)=ps;end;for i=1:length (SNRindB2), SNR=exp(SNRindB2(i)*log(10)/10);信%噪上匕theo_err_prb(i)=Qfunct(
10、sqrt(2*SNR);理论比特误码率end;%随后绘图曲线%semilogy (SNRindB1,smld_bit_err_prb,*r); holdsemilogy (SNRindB1,smld_symbol_err_prb,Og); semilogy (SNRindB2,theo_err_prb); grid xlabel(Eb/No in dB) ylabel(error probability) title(4PSK通信系统的蒙特卡洛仿真)%gtext(注：一”理论误码率；“大”误比特率；“o”误符号率)functionpb,ps=cm_sm32(snr_in_dB) N=10000
11、;E=1;snr=10八(snr_in_dB/10); sgma=sqrt(E/snr)/2; s00=10;s01=01;s11=-1 0; s10=0 -1; for i=1:N, temp=rand; if(temp0.25), dsource1(i)=0; dsource2(i)=0; elseif(temp0.5), dsource1(i)=0; dsource2(i)=1; elseif(temp0.75),dsource1(i)=1;dsource2(i)=0;elsedsource1(i)=1;dsource2(i)=1;end;end;numofsymbolerror=0;n
12、umofbiterror=0;for i=1:N,n(1)=gngauss(sgma);n(2)=gngauss(sgma);if(dsource1(i)=0) & (dsource2(i)=0), r=s00+n;elseif(dsource1(i)=0) & (dsource2(i)=1), r=s01+n;elseif(dsource1(i)=1) & (dsource2(i)=0), r=s10+n;elser=s11+n;end;c00=dot(r , s00);c01=dot(r , s01);c10=dot(r , s10);c11=dot(r , s11);c_max=max(
13、c00 c01 c10 c11);if(c00=c_max)decis1=0;decis2=0;elseif(c01=c_max)decis1=0;decis2=1;elseif(c10=c_max)decis1=1;decis2=0;elsedecis1=1;decis2=1;end;symbolerror=0;if(decis1=dsource1(i),numofbiterror=numofbiterror+1; symbolerror=1;end;if(decis2=dsource2(i),numofbiterror=numofbiterror+1; symbolerror=1;end;
14、if(symbolerror=1),numofsymbolerror=numofsymbolerror+1;end;end;ps=numofsymbolerror/N;pb=numofbiterror/(2*N);%functionpb,ps=cm_sm32(snr_in_dB)clear allsnr_in_dB=1N=10000;E=1;snr=10A(snr_in_dB/10);sgma=sqrt(E/snr)/2;s00=1 0;s01=0 1;s11=-1 0;s10=0 -1;for i=1:N,temp=rand;if(temp0.25),dsource1(i)=0;dsourc
15、e2(i)=0;elseif(temp0.5),dsource1(i)=0;dsource2(i)=1;elseif(temp0.75),dsource1(i)=1;dsource2(i)=0;elsedsource1(i)=1;dsource2(i)=1;end;end;numofsymbolerror=0;numofbiterror=0;for i=1:N,n(1)=gngauss(sgma);n(2)=gngauss(sgma);if(dsource1(i)=0) & (dsource2(i)=0), r=s00+n;elseif(dsource1(i)=0) & (dsource2(i
16、)=1), r=s01+n;elseif(dsource1(i)=1) & (dsource2(i)=0), r=s10+n;elser=s11+n;end;c00=dot(r , s00);c01=dot(r , s01);c10=dot(r , s10);c11=dot(r , s11);c_max=max(c00 c01 c10 c11);if(c00=c_max)decis1=0;decis2=0;elseif(c01=c_max)decis1=0;decis2=1;elseif(c10=c_max)decis1=1;decis2=0;elsedecis1=1;decis2=1;end
17、;symbolerror=0;if(decis1=dsource1(i),numofbiterror=numofbiterror+1;symbolerror=1;end;if(decis2=dsource2(i),numofbiterror=numofbiterror+1;symbolerror=1;end;if(symbolerror=1),numofsymbolerror=numofsymbolerror+1;end;end;ps=numofsymbolerror/N;pb=numofbiterror/(2*N);function y=Qfunct(x)y=(1/2)*erfc(x/sqrt(2);functiongsrv1,gsrv2=gngauss(m,sgma)if nargin=0,m=0;sgma=1;elseif nargin=1,sgma=m;m=0;end;u=rand;z=sgma*(sqrt(2*log(1/(1-u);u=rand;gsrv1=m+z*cos(2*pi*u);gsrv2=m+z*sin(2*pi*u);
