1、 QPSK 在 AWGN 信道下的仿真实验一、基本原理1.1QPSK 简介在数字相位调制中，M 个信号波形可表示为22(1)/Re()cos1cos2()sin(1)sin2cjftjmMmc ctgteft ftgmft （m=1,2,，M， ）0tT式中，g(t)是信号脉冲形状， （m=1,2,，M）是载波的 M 个可能的(1)/m相位，用于传送发送信息。这些信号波形具有相等的能量，即2200()()TTmgstdtd而且这些信号波形可以表示为两个标准正交信号波形 和 的线性星座图合，即1()ft2ft，式中12()()()mmstftsft1cos2cgfttft2()()incgftt
2、ft且二维向量 为1,2ms（m=1,2,，M）2cos(1),sin(1)ggmM 其中当 M=4 时就是本文要讨论的 4PSK（QPSK） ，QPSK 的载波相位有四种取值，每种取值代表两比特的信号。随着信号的改变，幅度恒定的载波信号的相位在四种取值间跳变。这四个相位的取值为间隔相等的值，比如，0， ，每一个相位值对应于唯一的一/2,3/对消息比特。有一种变形，称为 是通过在每一个符号间隔的载波相位中引入/4QPSK附加的 相移来使符号同步变得容易些。/4QPSK 信号可以表示为：，i=1，2，3，42()cos(1)2QPSKEttiT0StT式中 为单位符号的信号能量，即 时间内的信号
3、能量； 为载波角频率，S St c为符号持续时间。sTQPSK 信号可以看成是对两个正交的载波进行多电平双边带调制后所得信号的叠加，因此可以用正交调制的方法得到 QPSK 信号。QPSK 信号的星座如图 4.1.1 所示：图 1.1 QPSK 信号星座图1.2 QPSK 相位解调与检测从 AWGN 信道中，在一个信号区间内接收到的带宽信号可以表示为()cos2)(sin2)mccrtutnftft这里 和 是加性噪声的两个正交分量。()ctst可以将这个接收信号与 ， 给出的1()cos()Tctgtft2()sin2)Tcgtft和 作相关，两个相关器的输出产生受噪声污损的信号分量，它们可表
4、示为1()t2t 2(cossin)mcsmrsnM式中 和 定义为cs1()2cTcgtdt1()2sTsngtndt这两个正交噪声分量 和 是零均值，互不相关的高斯随机过程。这样，cs和 。 和 的方差是()0csEn()0sEcs220()csNn最佳检测器将接收信号向量 r 投射到 M 个可能的传输信号向量 之一上去，并选ms取对应于最大投影的向量。据此，得到相关准则为，m=0,1,，M-1(,)mCrs由于全部信号都具有相等的能量，因此，对数字相位调制一种等效的检测器标准是计算接收信号向量 r=( , )的相位为crsarctns并从信号集 中选取其相位最接近 的信号。msr在 AW
5、GN 信道中，因为二相相位调制与二进制 PAM 是相同的，所以差错概率为，式中 是每比特的能量。四相相位调制可以看作两个正交载波上的二相20bPQNb相位调制系统，所以 1 个比特的差错概率与二相相位调制是一样的。对于 M 4 的符号差错概率不存在简单的闭式表达式。对 的一种好的近似式是MP002sinibQNk式中 比特/符号。2logkM1.3 QPSK 调制解调原理四相相位键控（QPSK）也称之为正交 PSK，其调制原理如图 4.1.2 所示。图 1.2 QPSK 调制原理图如果输入的二进制信息码流（假设+1V 为逻辑 1，-1V 为逻辑 0）串行进入比特分离器，产生 2 个码流以并行方
6、式输出，分别被送入 I（正交支路）通道及 Q（同相支路）通道，又各自经过一个平衡调制器，与一个和参考振荡器同频的正交的载波（ 和 ）调sintcost制形成了四相相移键控信号即得到平衡器的输出信号后，经过一个带通滤波器，然后再进入行信号叠加，可以得到已经调制的 QPSK 信号。QPSK 的 4 种（I，Q 星座图合为 4 种0 0，0 1， 1 0，1 1）输出相位有相等的幅度，而且 2 个相邻的相位相差值为 90 度，但是输出相位并不满足 （m=0,1,M-1）,信号相位移可以偏移 45 度和-45 度，2mM接受端仍可以得到正确的解码，实际中数字输入电压必须比峰值载波电压高出很多，以确保平
7、衡器的正常工作。经过调制的信号通过信道传输到达用户端，需要进行解调，这样一过程是与调制相类似的逆过程。首先，QPSK 信号经过功率分离器形成两路相同的信号，进入乘积检验波，用两个正交的载波信号（ 和 ）实现相干解调，然后各自通过sintcost一个低通滤波器滤波得到低频和直流的成分，再经过一个并行-串行变换器，得到解调信号。QPSK 的解调原理如图 4.1.3 所示。图 1.3 QPSK 解调原理图目前 QPSK 调制的实现主要是利用数字电路和专用芯片来完成，通常利用可编程数字电路对基带信号进行码元变换，差分编码，成型滤波等处理后得到同相分量和正交分量，然后将两路信号分量经过数模转换获得模拟信
8、号送入一个正交相乘器与中频载波调制得到中频 QPSK 调制信号。该方法适合高码率数字信号的传输，但系统的开放性和灵活性较差。在解调过程中，若不考虑信道失真及噪声的影响,加到解调器输入端的接收信号在一个码元持续时间内可表示为： ()()cos()csSkkktgtTt式中，g(t)为信号的包络；为码元中的载波相位；k为载波角频率；c该信号同时加到两个鉴相（相乘）器上，在上支路积分器输出电压为： 100 0()cos)cs12o()2ss sTkTkkUgtttdtttgt取样器在 t= 时刻对 进行取样，所得到的是两个电压的叠加，即前一积分在 时刻的sTI sT积分值加上后一积分在 时刻的积分值
9、。当持续时间 内包含整数个载波周期时，前一积s sT分在 t= 时刻的积分值为 0，这时测到取样值完全由后一个积分所决定。后一积分在 t=s时刻的积分值是 与包络 g(t)的面积的乘积。因此 I 支路取样器的输出电压与sT1/2cosk成正比，即cok1coskU同理可得 Q 支路的输出电压与 成正比，即sinksik若判决器按极性判决，正的取样值为“1” ，负的取样值为“0” ，则可将调相信号解调为相应的数字信号，再经并串变换即可恢复出与发送端完全相同的数字信号。二、仿真1、误码率、时域图、眼图仿真等程序：close allclcclear allSNR_DB=0:1:12;sum=1000
10、000;data= randsrc(sum,2,0 1);a1,b1=find(data(:,1)=0message(a1)=-1-j;a2,b2=find(data(:,1)=0message(a2)=-1+j;a3,b3=find(data(:,1)=1message(a3)=1-j;a4,b4=find(data(:,1)=1message(a4)=1+j;scatterplot(message)title(B 点信号的星座图)A=1;Tb=1;Eb=A*A*Tb;P_signal=Eb/Tb;NO=Eb./(10.(SNR_DB/10);P_noise=P_signal*NO;sigm
11、a=sqrt(P_noise);for Eb_NO_id=1:length(sigma)noise1=sigma(Eb_NO_id)*randn(1,sum);noise2=sigma(Eb_NO_id)*randn(1,sum);receive=message+noise1+noise2*j;resum=0;total=0;m1=find(angle(receive)0);remessage(1,m1)=1+j;redata(m1,1)=1;redata(m1,2)=1;m2= find( angle(receive)pi/2remessage(1,m4)=1-j;redata(m4,1)=
12、1;redata(m4,2)=0;resum,ratio1=symerr(data,redata);pbit(Eb_NO_id)=resum/(sum*2);total,ratio2=symerr(message,remessage);pe(Eb_NO_id)=total/sum;endscatterplot(receive)title(C 点信号的星座图)Pe=1-(1-1/2*erfc(sqrt(10.(SNR_DB/10)/2).2;Pbit=1/2*erfc(sqrt(10.(SNR_DB/10)/2);figure(3)semilogy(SNR_DB,pe,:s,SNR_DB,Pe,-*,SNR_DB,pbit,-o,SNR_DB,Pbit,:+)legend(QPSK 仿真误码率,QPSK 理论误码率,QPSK 仿真误比特率,QPSK 理论误比特率,1)xlabel(信噪比/dB)ylabel(概率 P)gird on二、总结对实验结果的简单分析和说明由图可见 QPSK 仿真误码率曲线和理论误码率曲线重合在一起，QPSK 仿真误比特率曲线和理论误比特率曲线也重合在一起，误码率约是误比特率的两倍，说明实验方法是正确可行的。QPSK 信 号 的 误 码 率 ：QPSK 信号的误比特率：误码率是误比特率的两倍。2/21refcPe/21rfce
