QPSK通信系统的_MonteCarlo_仿真.doc

QPSK通信系统的 MonteCarlo 仿真 实验报告 山东大学 荣林林 20071001089 一、实验目的 1、提高独立学习的能力； 2、培养发现问题、解决问题和分析问题的能力； 3、学习Matlab 的使用； 4、掌握4PSK通信系统的Monte Carlo仿真方法； 5、掌握4PSK通信系统的的组成原理； 6、比较编码信号和为编码信号在随机信道中的传输，加深对纠错编码的理解； 二、系统框图及编程原理 实验原理 PSK是利用载波的不同相位表示相应的数字信息。对于二进制相位调制（M=2)来说，两个载波相位是0和π。对于M相相位调制来说M=2k，这里k是每个传输符号的信息比特数。4PSK是M=4的载波相位调制。 这里，将理论差错概率与仿真的差错概率比较，进一步观察仿真与理论值之间的差别。同时，用不同的判决准则对接受信号进行判决。并比较两种判别方法的差别。 一．QPSK调制原理 1．信号能量分析 一组M载波相位调制信号波形的一般表示式为 m=[0，M-1] 式中是发送滤波器的脉冲形状，它决定了传输信号的频谱特性，A是信号的幅度。注意到，PSK信号对所有m都具有相等的能量，即 代表每个传输符号的能量。 2．噪声分析 传输信号的信道假设被加性噪声n(t)所污损，这样信号在接收端将产生误码。因为n(t)是功率谱为的白高斯过程的一个样本函数，所以噪声分量就是零均值高斯型的，即 3．信号判决分析 最佳检测器将接收信号向量r投射到M个可能的传输信号向量{}之一上去，并选取对应于最大投影的向量。据此，得到相关准则为 m=[0,M-1] 检测器观察到接收信号向量，并计算r在4种可能的信号向量上的投影。根据选取对应于最大投影的信号点作为判决，从而判决出信号。 同时，检测器的判决准则也可采用最小距离法，即利用星座图上符号间的距离进行判决，从而得到判决结果。 二．Monte Carlo 仿真过程 仿真框图如图（一） 图（一） 用于 Monte Carlo仿真的4PSK系统的方框图 如图所示，利用一个随机数发生器，产生（0，1）范围内的随机数。再将这个范围分成四个相等的区间（0，0.25），（0.25，0.5），（0.5，0.75），（0.75，1.0），这些子区间分别对应于00，01，11，10信息比特对，再用这些比特对来选择信号相位向量。加性噪声的同相分量和正交分量 ，在上面讨论过，即为零均值，方差为 的统计独立的高斯随机变量。在检测器观察到的接收信号向量 ，利用上面讨论的两种检测方法，得到判决结果，并与传输符号作比较，最后对符号差错和比特差错计数 三、实验内容及程序分析（以下程序皆以N=1000为例） %映射比较子函数 %函数分为四步第一步产生随机序列，进行4PSK映射。 %第二步：调用高斯高斯噪声子函数，产生正交两路高斯噪声，与输出符号序列相加 %第三步：检测接受信号 %第四步：计算误码率和误比特率 N=1000; %符号长度 E=1; %计算噪声方差 sgma=input(方差=); sgma=sqrt(sgma); % 4PSK比特映射 s00=[1 0]; s01=[0 1]; s11=[-1 0]; s10=[0 -1]; %第一步产生随机序列，进行4PSK映射。生成随机信源 for i=1:N, %生成随机信源 temp=rand; if (temp0.25), % With probability 1/4, source output is 00. dsource1(i)=0; dsource2(i)=0; elseif (temp0.5), % With probability 1/4, source output is 01. dsource1(i)=0; dsource2(i)=1; elseif (temp0.75), % With probability 1/4, source output is 10. dsource1(i)=1; dsource2(i)=0; else % With probability 1/4, source output is 11. dsource1(i)=1; dsource2(i)=1; end; end; numofsymbolerror=0; %检测错误并计算错误率 numofbiterror