《MATLABQPSK在AWGN信道下的仿真》由会员分享，可在线阅读，更多相关《MATLABQPSK在AWGN信道下的仿真（10页珍藏版）》请在金锄头文库上搜索。
1、、基本原理1.1QPSK 简介在数字相位调制中,M个信号波形可表示为Sm(t )=Re二 g(t)cos 2 二 fctg t ej2”(mJl)/Mej2 二以2：.=g(t)cos m 1 cos2二 fct 一 g(t)sin M(m=1,2, , , M, 0 mt )(m -1)sin 2二 fct式中，g(t)是信号脉冲形状，入=2二(m-1)/M(m=1,2, , , M)是载波的M个可能的相位，用于传送发送信息。这些信号波形具有相等的能量，即T 21 T 21；=0 Sm(t)d- 0 g (t)dt =2而且这些信号波形可以表示为两个标准正交信号波形fl(t)和f2(t)的线
2、性星座图合，即Sm(tSm1f1(t)Sm2f2(t)，式中弾)二cif2(t)= -2 g(t)sin 2二 fct且二维向量 Sm二Sm1,Sm2为2兀_cosM(m1)-2g sinMg 亠2 (m1)(m=1,2,，M)其中当M=4时就是本文要讨论的4PSK( QPSK, QPSK勺载波相位有四种取值，每种取值代表两比特的信号。随着信号的改变,幅度恒定的载波信号的相位在四种取值间跳变。这四个相位的取值为间隔相等的值，比如,0，二/ 2,二,3二/ 2，每一个相位值对应于唯一的一对消息比特。有一种变形，称为 二/ 4 -QPSK是通过在每一个符号间隔的载波相位中引入附加的 /4相移来使符
3、号同步变得容易些。QPSK言号可以表示为：sqpsk(t)cos ct (i 0 _t _Ts , i=，2, 3, 4V Ts2式中Es为单位符号的信号能量，即0乞匕Ts时间内的信号能量；c为载波角频率,Ts为符号持续时间。QPSK言号可以看成是对两个正交的载波进行多电平双边带调制后所得信号的叠加，因此可以用正交调制的方法得到QPSK言号。QPSK言号的星座如图4.1.1所示：图1.1 QPSK信号星座图1.2 QPSK相位解调与检测从AWG信道中，在一个信号区间内接收到的带宽信号可以表示为r t 二 Um(t) n(t)二 Um(t) nc(t)cos(2 二 fct) -ns(t)sin
4、(2 二 fct)这里nc(t)和ns(t)是加性噪声的两个正交分量。可以将这个接收信号与 =(t)二 gT(t)cos(2二 fct),，2 二-gT(t)sin(2二 fct)给出的 (t)和 2(t)作相关，两个相关器的输出产生受噪声污损的信号分量，它们可表示为m “(仁込乜仁sin讥)mMM式中n,和ns定义为1% =2 ._：gT(t)nc(t)dt1足 J gT(t)ns(t)dt这两个正交噪声分量nc(t)和ns(t)是零均值，互不相关的高斯随机过程。这样,E(nc)二 E (ns) = 0和 E(ncns) =0。n和 ns 的方差是2E(nJ E(ns)_ 最佳检测器将接收信
5、号向量r投射到M个可能的传输信号向量 sm之一上去，并选取对应于最大投影的向量。据此，得到相关准则为C(r,Sm) =r *sm，m=0,1, ,，M-1由于全部信号都具有相等的能量，因此，对数字相位调制一种等效的检测器标准是计算接收 信号向量r=( rc, rs)的相位为片=arcta n并从信号集 sm中选取其相位最接近 齐的信号。在AWGN言道中，因为二相相位调制与二进制PAM是相同的，所以差错概率为F2=qI伍 1式中.是每比特的能量。四相相位调制可以看作两个正交载波上的二相No相位调制系统，所以1个比特的差错概率与二相相位调制是一样的。对于M - 4的符号差错概率不存在简单的闭式表达
6、式。对P.的一种好的近似式是式中k = log2 M比特/符号。1.3 QPSK调制解调原理四相相位键控(QPSK也称之为正交 PSK其调制原理如图 4.1.2所示。图1.2 QPSK调制原理图如果输入的二进制信息码流（假设+1V为逻辑1, -1V为逻辑0）串行进入比特分离器，产生2个码流以并行方式输出，分别被送入I （正交支路）通道及 Q （同相支路）通道，又各自经过一个平衡调制器，与一个和参考振荡器同频的正交的载波（sit和cost ）调制形成了四相相移键控信号即得到平衡器的输出信号后，经过一个带通滤波器， 然后再进入行信号叠加，可以得到已经调制的QPSK言号。QPSK勺4种（I , Q星
7、座图合为4种0 0，0 1 , 1 0 , 1 1）输出相位有相等的幅度，而且2个相邻的相位相差值为 90度，但是输出相位并不满足m（ m=0,1, ,M-1 ），信号相位移可以偏移 45度和-45度，接受M端仍可以得到正确的解码，实际中数字输入电压必须比峰值载波电压高出很多，以确保平衡器的正常工作。经过调制的信号通过信道传输到达用户端，需要进行解调，这样一过程是与调制相类似的逆过程。首先，QPSK信号经过功率分离器形成两路相同的信号，进入乘积检验波，用两个正交的载波信号（sin ，t和cos t）实现相干解调，然后各自通过一个低通滤波器滤波得到低频和直流的成分，再经过一个并行-串行变换器，得
8、到解调信号。QPSK的解调原理如图4.1.3所示。FPGA图1.3 QPSK解调原理图目前QPSK调制的实现主要是利用数字电路和专用芯片来完成，通常利用可编程数字电路对基带信号进行码元变换，差分编码，成型滤波等处理后得到同相分量和正交分量，然后将两路信号分量经过数模转换获得模拟信号送入一个正交相乘器与中频载波调制得到中频QPSK调制信号。该方法适合高码率数字信号的传输，但系统的开放性和灵活性较差。在解调过程中，若不考虑信道失真及噪声的影响，加到解调器输入端的接收信号在一个码元持续时间内可表示为：s(t)g(tkTs)cosJtk)二 g(t)cos( L)式中，g(t)为信号的包络；为码元中的
9、载波相位；c为载波角频率；该信号同时加到两个鉴相(相乘)器上，在上支路积分器输出电压为：TsUi = g(t)cos(，ct：k)cos ctdt0TsTs1 s1s=Jg(t)cos(2 t ZQdt 十 cos毋k Jg(t)dt2 o2o取样器在t= Ts时刻对U i进行取样，所得到的是两个电压的叠加，即前一积分在Ts时刻的积分值加上后一积分在 Ts时刻的积分值。当持续时间Ts内包含整数个载波周期时，前一积分在t= Ts时刻的积分值为0,这时测到取样值完全由后一个积分所决定。后一积分在t= Ts时刻的积分值是1/2cos ;与包络g(t)的面积的乘积。因此I支路取样器的输出电压与cos成
10、正比，即U cos k同理可得Q支路的输出电压与 sink成正比，即Uk 二 sin k若判决器按极性判决，正的取样值为“1”，负的取样值为“ 0”，则可将调相信号解调为相应的数字信号，再经并串变换即可恢复出与发送端完全相同的数字信号。二、仿真1、误码率、时域图、眼图仿真等程序：close allclcclear allSNR_DB=0:1:12;sum=1000000;data= ran dsrc(sum,2,0 1);a1,b1=find(data(:,1)=0&data(:,2)=0);message(a1)=-1-j;a2,b2=find(data(:,1)=0&data(:,2)=1
11、); message(a2)=-1+j;a3,b3=find(data(:,1)=1&data(:,2)=0); message(a3)=1-j;a4,b4=find(data(:,1)=1&data(:,2)=1);message(a4)=1+j;scatterplot(message)title(B 点信号的星座图 )A=1;Tb=1;Eb=A*A*Tb;P_signal=Eb/Tb;NO=Eb./(10.A(SNR_DB/10);P_noise=P_signal*NO;sigma=sqrt(P_noise);for Eb_NO_id=1:length(sigma) noise1=sigm
12、a(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/2&angle(receive)-pi&angle(receive)-pi/2&angle(receive)=0); remessage(1,m4)=1-j;redata
13、(m4,1)=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;end scatterplot(receive)title(C点信号的星座图)Pe=1-(1-1/2*erfc(sqrt(10.A(SNR_DB/10)/2).A2;Pbit=1/2*erfc(sqrt(10.A(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 onB点信号的星座團11-*1 1 11II08-I06-0.4-0.2-10-0.2-0.4-T-0.6-o.a-11-1-1ji|*-1-D.50
