QPSK在AWGN信道下的仿真

QPSK在AWGN信道下的仿真实验

一、 基本原理 1.1QPSK简介

在数字相位调制中，M个信号波形可表示为

j2(m1)/Mj2fctsmtRegtee2g(t)cos2fctm1M22g(t)cosm1cos2fctg(t)sin(m1)sin2fctMM

（m=1,2,…，M，0tT）

式中，g(t)是信号脉冲形状，m2(m1)/M （m=1,2,…，M）是载波的M个可能的相位，用于传送发送信息。

这些信号波形具有相等的能量，即

1T21g(t)dtg022

2sm(t)dt0T而且这些信号波形可以表示为两个标准正交信号波形f1(t)和f2(t)的线性星座图合，即

sm(t)sm1f1(t)sm2f2(t)，式中

f1(t)2gg(t)cos2fct

且二维向量

f2(t)2gg(t)sin2fct

smsm1,sm2为

gg22smcos(m1),sin(m1)2M2M （m=1,2,…，M）

其中当M=4时就是本文要讨论的4PSK（QPSK），QPSK的载波相位有四种取值，每种取值代表两比特的信号。随着信号的改变，幅度恒定的载波信号的相位在四种取值间跳变。这四个相位的取值为间隔相等的值，比如，0，/2,,3/2，每一个相位值对应于唯一的一对消息比特。有一种变形，称为/4QPSK是通过在每一个符号间隔的载波相位中引入附加的/4相移来使符号同步变得容易些。

QPSK信号可以表示为：

SQPSK(t)2Escosct(i1)Ts2 0tTS，i=1，2，3，4

式中ES为单位符号的信号能量，即0tTS时间内的信号能量；c为载波角频率，Ts为符号持续时间。

QPSK信号可以看成是对两个正交的载波进行多电平双边带调制后所得信号的叠加，因此可以用正交调制的方法得到QPSK信号。

QPSK信号的星座如图4.1.1所示：

图1.1 QPSK信号星座图

1.2 QPSK相位解调与检测

从AWGN信道中，在一个信号区间内接收到的带宽信号可以表示为

rtum(t)n(t)um(t)nc(t)cos(2fct)ns(t)sin(2fct)

这里nc(t)和ns(t)是加性噪声的两个正交分量。

可以将这个接收信号与1(t)gT(t)cos(2fct)，2gT(t)sin(2fct)给出的1(t)和2(t)作相关，两个相关器的输出产生受噪声污损的信号分量，它们可表示为

rsmn(scos2m2mncssinns)MM

式中nc和ns定义为

1ncgT(t)nc(t)dt2

1gT(t)ns(t)dt2

ns这两个正交噪声分量nc(t)和ns(t)是零均值，互不相关的高斯随机过程。这样，

E(nc)E(ns)0和E(ncns)0。nc和ns的方差是

N02

E(nc2)E(ns2)最佳检测器将接收信号向量r投射到M个可能的传输信号向量{sm}之一上去，并选取对应于最大投影的向量。据此，得到相关准则为

C(r,sm)r•sm，m=0,1,…，M-1

由于全部信号都具有相等的能量，因此，对数字相位调制一种等效的检测器标准是计算接收信号向量r=(rc,rs)的相位为

rsrc

rarctan并从信号集{sm}中选取其相位最接近r的信号。

在AWGN信道中，因为二相相位调制与二进制PAM是相同的，所以差错概率为

2bP2QN0，式中b是每比特的能量。四相相位调制可以看作两个正交载波上的二相相位

调制系统，所以1个比特的差错概率与二相相位调制是一样的。对于M4的符号差错概率不存在简单的闭式表达式。对PM的一种好的近似式是

2sPM2QsinNM02kb2QsinNM0

式中klog2M比特/符号。

1.3 QPSK调制解调原理

四相相位键控（QPSK）也称之为正交PSK，其调制原理如图4.1.2所示。

图1.2 QPSK调制原理图

如果输入的二进制信息码流（假设+1V为逻辑1，-1V为逻辑0）串行进入比特分离器，产生2个码流以并行方式输出，分别被送入I（正交支路）通道及Q（同相支路）通道，又各自经过一个平衡调制器，与一个和参考振荡器同频的正交的载波（sint和cost）调制形成了四相相移键控信号即得到平衡器的输出信号后，经过一个带通滤波器，然后再进入行信号叠加，可以得到已经调制的QPSK信号。QPSK的4种（I，Q星座图合为4种

[0 0]，[0 1]，[1 0]，[1 1]）输出相位有相等的幅度，而且2个相邻的相位相差值为90度，但是输出相位并不满足

m2mM（m=0,1,…,M-1）,信号相位移可以偏移45度和-45度，

接受端仍可以得到正确的解码，实际中数字输入电压必须比峰值载波电压高出很多，以确保平衡器的正常工作。经过调制的信号通过信道传输到达用户端，需要进行解调，这样一过程是与调制相类似的逆过程。首先，QPSK信号经过功率分离器形成两路相同的信号，进入乘积检验波，用两个正交的载波信号（sint和cost）实现相干解调，然后各自通过一个低通滤波器滤波得到低频和直流的成分，再经过一个并行-串行变换器，得到解调信号。QPSK的解调原理如图4.1.3所示。

图1.3 QPSK解调原理图

目前QPSK调制的实现主要是利用数字电路和专用芯片来完成，通常利用可编程数字电路对基带信号进行码元变换，差分编码，成型滤波等处理后得到同相分量和正交分量，然后将两路信号分量经过数模转换获得模拟信号送入一个正交相乘器与中频载波调制得到中频QPSK调制信号。该方法适合高码率数字信号的传输，但系统的开放性和灵活性较差。

在解调过程中，若不考虑信道失真及噪声的影响,加到解调器输入端的接收信号在一个

s(t)g(tkTS)cos(ctk)k码元持续时间内可表示为：g(t)cos(ctk)

式中，g(t)为信号的包络； k为码元中的载波相位；

c为载波角频率；

该信号同时加到两个鉴相（相乘）器上，在上支路积分器输出电压为：

TsU1g(t)cos(ctk)cosctdt0s1s1g(t)cos(2tk)dtcoskg(t)dt2020

TT取样器在t=Ts时刻对UI进行取样，所得到的是两个电压的叠加，即前一积分在Ts时刻的积分值加上后一积分在Ts时刻的积分值。当持续时间Ts内包含整数个载波周期时，前一积分在t=Ts时刻的积分值为0，这时测到取样值完全由后一个积分所决定。后一积分在t=Ts时刻的积分值是1/2cosk与包络g(t)的面积的乘积。因此I支路取样器的输出电压与cosk成正比，即

U1cosk 同理可得Q支路的输出电压与sink成正比，即

Uksink

若判决器按极性判决，正的取样值为“1”，负的取样值为“0”，则可将调相信号解调为相应的数字信号，再经并串变换即可恢复出与发送端完全相同的数字信号。

二、仿真

1、误码率、时域图、眼图仿真等

%QPSK误码率在瑞利高斯信道下的比较

clear;

close all;

N=10000;

Fd=1;

Fs=1*Fd;

M=4;

SNR_db=0:2:14;

for n=1:length(SNR_db)

Eb_N0=10^(SNR_db(n)/10);

sgma=sqrt(1/(2*Eb_N0));

x=randint(N,1,M);

y=dmodce(x,Fd,Fs,'psk',M);

a11=randn(1,N);

b11=randn(1,N);

%产生低通滤波器

b=-ones(1,Fs);b=b/Fs;a=1;

fad1=filter(b,a,a11);

fad2=filter(b,a,b11);

ynoise=y.*(abs(fad1+j*fad2))'+sqrt(Fs/Fd)*sgma*(randn(length(y),1)+j*randn(length(y),1));

ynoise1=y+sqrt(Fs/Fd)*sgma*(randn(length(y),1)+j*randn(length(y),1))

z=ddemodce(ynoise,Fd,Fs,'psk',M);

z1=ddemodce(ynoise1,Fd,Fs,'psk',M);

[numbers,p1m(n)]=symerr(x,z);

[numbers,p2m(n)]=symerr(x,z1)

end;

semilogy(SNR_db,p1m,'*-',SNR_db,p2m,'rx-');

xlabel('信噪比（dB）');

ylabel('误码率');

legend('瑞利信道仿真得到的误码率','高斯信道仿真得到的误码率');

axis([0 15 10^(-4) 1]);

%高斯信道下QPSK星座图程序

clear;

close all;

N=10000;

Fd=1;

Fs=1*Fd;

M=4;

for SNR_db=0:10:20

Eb_N0=10^(SNR_db/10);

sgma=sqrt(1/(8*Eb_N0));

x=randint(N,1,M);

y=dmodce(x,Fd,Fs,'psk',M);

ynoise=y+sqrt(Fs/Fd)*sgma*(randn(length(y),1)+j*randn(length(y),1));

figure(SNR_db+1);

axis([-1.2 1.2 -1.2 1.2]);

hold on;

for i=0:M-1

plot(cos(2*pi*i/M),sin(2*pi*i/M),'.','MarkerSize',20);

end

plot(ynoise,'+');

hold off;

end;

三、总结