基于-matlab的QPSK和BPSK信号性能比较仿真.docx
《基于-matlab的QPSK和BPSK信号性能比较仿真.docx》由会员分享，可在线阅读，更多相关《基于-matlab的QPSK和BPSK信号性能比较仿真.docx（19页珍藏版）》请在冰点文库上搜索。
目 录
第一章 概述 2
第二章 QPSK通信系统原理与仿真 2
2.1QPSK系统框图介绍 2
2.2QPSK信号的调制原理 3
2.2.1QPSK信号产生方法 3
2.2.2QPSK星座图 4
2.3QPSK解调原理及误码率分析 4
2.3.1QPSK解调方法 4
2.3.2QPSK系统误码率 5
2.4QPSK信号在AWGN信道下仿真 5
第三章 BPSK通信系统原理与仿真 6
3.1BPSK信号的调制原理 6
3.2BPSK解调原理及误码率分析 7
第四章 QPSK与BPSK性能比较 8
4.1QPSK与BPSK在多信道下比较仿真 8
4.1.1纵向比较分析 8
4.1.2横向比较分析 10
4.2仿真结果分析 10
4.2.1误码率分析 10
4.2.2频带利用率比较 10
附 录 11
代码1 11
代码2 11
代码3 14
代码4 16
第一章 概述
QPSK是英文QuadraturePhaseShiftKeying的缩略语简称，意为正交相移键控，是一种数字调制方式。
它以其抗干扰性能强、误码性能好、频谱利用率高等优点，广泛应用于数字微波通信系统、数字卫星通信系统、宽带接人、移动通信及有线电视系统之中。
BPSK是英文BinaryPhaseShiftKeying的缩略语简称，意为二相相移键控，是利用偏离相位的复数波浪组合来表现信息键控移相方式的一种。
它使用了基准的正弦波和相位反转的波浪，使一方为0，另一方为1，从而可以同时传送接受2值（1比特）的信息。
本文所研究的QPSK系统与二进制的BPSK系统相比，具有以下特点:
1.在传码率相同的情况下，四进制数字调制系统的信息速率是二进制系统的2倍。
2.在相同信息速率条件下，四进制数字调制系统的传码率是二进制系统的1/4倍，这一特点使得四进制码元宽度是二进制码元宽度的2倍，码元宽度的加大，可增加每个码元的能量，也可减小码间串扰的影响。
3.由于四进制码元速率比二进制的降低，所需信道带宽减小。
4.在接收系统输入信噪比相同的条件下，四进制数字调制系统的误码率要高于二进制系统。
5.四进制数字调制系统较二进制系统复杂，常在信息速率要求较高的场合。
基于以上优点，在数字信号的调制方式中QPSK（QuadraturePhaseShiftKeying）四相移键控是目前最常用的一种卫星数字信号调制方式,它具有较高的频谱利用率、较强的抗干扰性，在电路上实现也较为简单，因而被WCDMA和CDMA2000等第三代移动通信系统采用。
第二章 QPSK通信系统原理与仿真
2.1QPSK系统框图介绍
在图2.1的系统中，发送方，QPSK数据源采用随机生成，信源编码采用差分编码，编码后的信号经QPSK调制器，经由发送滤波器进入传输信道。
接收方，信号首先经过相位旋转，再经匹配滤波器解调，经阈值比较得到未解码的接收信号，差分译码后得到接收信号，与信源发送信号相比较，由此得到系统误码率，同时计算系统误码率的理论值，将系统值与理论值进行比较。
对于信道，这里选取的是加性高斯白噪声（Additive White Gaussian Noise）以及多径
噪声
衰减
发送滤波器
QPSK调
制器
差分编码器
Rayleigh衰落信道（MultipathRayleighFadingChannel）。
计算BER
QPSK
数据源
符号同步器
相位同步
相位旋转
积分与清除
阈值比较
差分译码器
图2.1QPSK系统框图
在实验中，选用的是差分码。
差分码又称为相对码，在差分码中利用电平跳变来分别表示1或0，分为传号差分码和空号差分码。
传号差分码：
当输入数据为“1”时，编码波型相对于前一码电平产生跳变；输入为“0”时，波型不产生跳变。
空号差分码：
当输入数据为“0”时，编码波型相对于前一码电平产生跳变；输入为“1”时，波型不产生跳变。
2.2QPSK信号的调制原理
2.2.1QPSK信号产生方法
QPSK信号通常由图2.2所示的调相法产生，输入的二进制序列经过一个串/并转换器后分为两路二进制序列，这里假设两路序列为a，b，并通过平衡调制器分别对同相载波及正交载波进行二相调制，得到如图2.3中的虚线矢量。
QPSK信号每个码元含有2比特的信息，如图2.3所示现用ab代表这两个比特。
两个比特有4种组合，即00、01、10和11。
它们和相位之间的关系通常都按格雷码的规律安排，如表2-
1所示。
图2.2调相法产生QPSK信号 图2.3信号矢量
2.2.2QPSK星座图
在本次系统仿真中使用仿真代码产生QPSK信号星座图，具体代码如附录代码1。
运行得到的信号星座图如图2.4，2.5：
表2-1QPSK编码规则
a
b
qk
A方式
B方式
0
0
0o
225o
1
0
90o
315o
1
1
180o
45o
0
1
270o
135o
表中qk－一组间隔均匀的受调制相位
2.3QPSK解调原理及误码率分析
2.3.1QPSK解调方法
由于QPSK信号可以看作两个正交2PSK信号的叠加，解调框图如图2.6，用相干解调方法，即用两路正交的相干载波，可以很容易的分离出这两路正交的2PSK信号。
解调后的两路基带信号码元a和b，经过并串变换后，成为串行数据输出。
图2.4A方式信号星座图 图2.5B方式信号星座图
2.3.2QPSK系统误码率
在QPSK体制中，由其矢量图（图2.7）可以看出，因噪声的影响使接收端解调时发生错误判决，是由于信号矢量的相位发生偏离造成的。
例如，设发送矢量的相位为45o，它代表基带信号码元“11”，若因噪声的影响使接收矢量的相位变成135o，则将误判为“01”。
当各个发送矢量以等概率出现时，合理的判决门限应该设在和相邻矢量等距离的位置。
在图中对于矢量
“11”来说，判决门限应该设在0o和90o。
当发送“11”时，接收信号矢量的相位若超出这一
9
范围（图
0o
图2.6QPSK信号解调原理框图 图2.7QPSK噪声容限
中阴影区），则将发生错判。
设f（q）为接收矢量（包括信号和噪声）相位的概率密度，则发生错误的概率为：
p/2
Pe=1-ò0
f（q）dq
省略计算f（q）和Pe的繁琐过程，直接给出计算结果：
Pe=1-[1
2
erfc
r/2]2
上式计算出的是QPSK信号的误码率。
若考虑其误比特率，正交的两路相干解调方法和
2PSK中采用的解调方法一样。
所以其误比特率的计算公式也和2PSK的误码率公式一样。
2.4QPSK信号在AWGN信道下仿真
仿真代码见附录代码2。
运行结果如图2.8：
图2.8AWGN环境下QPSK系统的仿真和理论结果
由运行结果可以看出，仿真值与理论值基本吻合，误码率值随信噪比增大而下降。
第三章 BPSK通信系统原理与仿真
3.1BPSK信号的调制原理
s（t）
双极性不归零
乘法器
e2PSK（t）
码型变换
coswct
图3.1模拟调制方法
0
p
e2PSK（t）
s（t）
1800移相
coswct
开关电路
图3.2键控调制方法
BSPK信号通常有2种调制方式，分别如图3.1、3.2所示。
在2PSK中，通常用初始相位0和p分别表示二进制“1”和“0”。
因此，2PSK信号的时
域表达式为
e2PSK（t）=Acos（wct+jn）
式中，jn表示第n个符号的绝对相位：
j=ì0，发送“0”时
n í
îp,
发送“1”时
因此，上式可以改写为
e （t）=ì
Acoswct,
概率为P
2PSK
í-Acoswt,
概率为1-P
î c
3.2BPSK解调原理及误码率分析
BPSK信号的解调器原理方框图3.3和波形图3.4，这里给出的解调框图采用相干解调的接收方法。
e2PSK（t）
a
相乘器
带通滤波器
coswct b
c d e
抽样判决器
低通滤波器
输出
定时脉冲
图3.3BPSK信号相干解调框图
波形图中，假设相干载波的基准相位与2PSK信号的调制载波的基准相位一致（通常默认为
0相位）。
但是，由于在2PSK信号的载波恢复过程中存在着的相位模糊，即恢复的本地载波与
所需的相干载波可能同相，也可能反相，这种相位关系的不确定性将会造成解调出的数字基带信号与发送的数字基带信号正好相反，即“1”变为“0”，“0”变为“1”，判决器输出数字信
t
t
t
t
1
0
0
1
1
Ts
1 0 0 1 1
a
b
c
d
e t
图3.4相干解调波形示意图
号全部出错。
这种现象称为2PSK方式的“倒π”现象或“反相工作”。
这也是2PSK方式在实际中很少采用的主要原因。
另外，在随机信号码元序列中，信号波形有可能出现长时间连续的正弦波形，致使在接收端无法辨认信号码元的起止时刻。
第四章 QPSK与BPSK性能比较
4.1QPSK与BPSK在多信道下比较仿真
QPSK仿真程序见代码3。
BPSK仿真程序见代码4。
下面对QPSK与BPSK在多信道下的BER-SNR曲线进行横向及纵向的比较，如表4-1所示。
4.1.1纵向比较分析
图中依次给出LOS信道、莱斯平坦衰落信道、莱斯频率选择性衰落、瑞利平坦衰落信道、瑞利频率选择性衰落。
比较频率选择性衰落和其相应的平坦衰落，可发现频率选择性衰落比相应的平坦衰落的性能差。
比较莱斯信道和瑞利信道，可发现不存在视距分量的瑞利信道比存在视距分量的任何莱斯信道的性能差。
表4-1QPSK与BPSK在多信道下性能比较
信道
QPSK
BPSK
LOS
莱斯平坦衰落
莱斯频率选择性衰落
信道
QPSK
BPSK
瑞利平坦衰落
瑞利频率选择性衰落
4.1.2横向比较分析
相同信道下，BPSK调制的系统误码率小于QPSK调制。
本次仿真过程中，莱斯平坦衰落信道下两者误码率差别显著。
4.2仿真结果分析
4.2.1误码率分析
由前面介绍的误码率内容可以看出，QPSK判决门限为90o，BPSK的判决门限为180o。
因此相同系统情况下的误码率BPSK优于QPSK。
4.2.2频带利用率比较
在传码率相同的情况下，四进制数字调制系统的信息速率是二进制系统的2倍。
频带利用
率公式
h=Rb
B
在相同信号速率的情况下，QPSK和BPSK系统的带宽是相同的，但是由于QPSK每个信
号都是四进制的，QPSK每个信号包含2bit信息，所以比特率就是BPSK的两倍，因而其频带利用率即为BPSK的两倍。
BPSK系统理论的频带利用率最大为1，但是在实际的实现中不能达到1，而在QPSK系统中，频带利用率可以超过1。
附 录
代码1
functionqpskconstellationA（M）M=4;
x=[0:
M-1];
scatterplot（pskmod（x,M））; % A方式QPSK信号gridon;
return
functionqpskconstellationB（M）M=4;
x=[0:
M-1];
scatterplot（pskmod（x,M,pi/4））; % B方式QPSK信号gridon;
return
代码2
function[x,bits]=random_binary（nbits,nsamples）
x=zeros（1,nbits*nsamples）;bits=round（rand（1,nbits））;form=1:
nbits
forn=1:
nsamples
index=（m-1）*nsamples+n;x（1,index）=（-1）^bits（m）;
end
end
return
function[c,lags]=vxcorr（a,b）
%计算ab的互相关系数，返回C是列向量，长度为ab的长度和减1a=a（:
）;
b=b（:
）;
M=length（a）;maxlag=M-1;
lags=[-maxlag:
maxlag]';A=fft（a,2^nextpow2（2*M-1））;%快速傅立叶变换B=fft（b,2^nextpow2（2*M-1））;
c=ifft（A.*conj（B））;%快速傅立叶反变换
c=[c（end-maxlag+1:
end,1）;c（1:
maxlag+1,1）];[nrnc]=size（a）;
if（nr>nc）
c=c';lags=lags.';
endreturn
functionBER_MC=MCQPSKrun（N,Eb,N0,ChanAtt,TimingBias,TimingJitter,PhaseBias,PhaseJitter）
fs=1e+6; % 采样速率
SymRate=1e+5; % 信号速率
Ts=1/fs; % 采样周期
TSym=1/SymRate; % 信号周期
SymToSend=N; % 发送信号数
ChanBW=4.99e+5; % 带宽
MeanCarrierPhaseError=PhaseBias; % 载波相位均值
StdCarrierPhaseError=PhaseJitter; % 相位误差标准差MeanSymbolSyncError=TimingBias; % 符号同步误差均值StdSymbolSyncError=TimingJitter; % 符号同步误差标准差ChanGain=10^（-ChanAtt/20）; % 信道增益
TxBitClock=Ts/2; % 发送机时钟
RxBitClock=Ts/2; % 接收机时钟
%接收机输入端噪声标准差和信号幅度
RxNoiseStd=sqrt（（10^（（N0-30）/10））*（fs/2））;TxSigAmp=sqrt（10^（（Eb-30）/10）*SymRate）;
%分配缓存SampPerSym=fs/SymRate;
probe1=zeros（（SymToSend+1）*SampPerSym,1）;probe1counter=1;probe2=zeros（（SymToSend+1）*SampPerSym,1）;probe2counter=1;
%计已传输信号的个数
TxSymSent=1;
RxSymDemod=0;
%发送和接收数据缓冲区[unused,SourceBitsI]=random_binary（SymToSend,1）;[unused,SourceBitsQ]=random_binary（SymToSend,1）;
%差分编码TxBitsI=SourceBitsI*0;TxBitsQ=SourceBitsQ*0;fork=2:
length（TxBitsI）
TxBitsI（k）=or（and（not（xor（SourceBitsI（k）,SourceBitsQ（k）））,xor（SourceBitsI（k）,TxBitsI（k-1）））,and（xor（SourceBitsI（k）,SourceBitsQ（k））,xor（SourceBitsQ（k）,TxBitsQ（k-1））））;
TxBitsQ（k）=or（and（not（xor（SourceBitsI（k）,SourceBitsQ（k）））,xor（SourceBitsQ（k）,TxBitsQ（k-1）））,and（xor（SourceBitsI（k）,SourceBitsQ（k））,xor（SourceBitsI（k）,TxBitsI（k-1））））;
end
%产生复信号
TxBits=（（TxBitsI*2）-1）+（sqrt（-1）*（（TxBitsQ*2）-1））;RxIntegrator=0; %初始化接收机积分器
TxBitClock=2*TSym; %初始化发送机
%设计信道滤波器，产生滤波器参数序列[b,a]=butter（2,ChanBW/（fs/2））;
b=[1];
a=[1];
[junk,FilterState]=filter（b,a,0）;
%开始仿真循环
whileTxSymSentTxBitClock=TxBitClock+Ts;ifTxBitClock>TSym
TxSymSent=TxSymSent+1;TxBitClock=mod（TxBitClock,TSym）;TxOutput=TxBits（TxSymSent）*TxSigAmp;
end
%信号经过信道滤波器[Rx,FilterState]=filter（b,a,TxOutput,FilterState）;
%加高斯白噪声Rx=（ChanGain*Rx）+（RxNoiseStd*（randn（1,1）+sqrt（-1）*randn（1,1）））;
%基于接收机载波同步误差的相位旋转
PhaseRotation=exp（sqrt（-1）*2*pi*（MeanCarrierPhaseError+（randn（1,1）*StdCarrierPhaseError））/360）;Rx=Rx*PhaseRotation;
probe1（probe1counter）=Rx;probe1counter=probe1counter+1;
%更新接收机积分清除器RxIntegrator=RxIntegrator+Rx;probe2（probe2counter）=RxIntegrator;probe2counter=probe2counter+1;
%更新接收机时钟，判断是不是适合采样RxBitClock=RxBitClock+Ts;
RxTSym=TSym*（1+MeanSymbolSyncError+（StdSymbolSyncError*randn（1,1）））;ifRxBitClock>RxTSym%解调信号
RxSymDemod=RxSymDemod+1;RxBitsI（RxSymDemod）=round（sign（real（RxIntegrator））+1）/2;RxBitsQ（RxSymDemod）=round（sign（imag（RxIntegrator））+1）/2;RxBitClock=RxBitClock-TSym;
RxIntegrator=0;
end
end
%差分解码SinkBitsI=SourceBitsI*0;SinkBitsQ=SourceBitsQ*0;
fork=2:
RxSymDemodSinkBitsI（k）=or（and（not（xor（RxBitsI（k）,RxBitsQ（k）））,xor（RxBitsI（k）,RxBitsI（k-
1）））,and（xor（RxBitsI（k）,RxBitsQ（k））,xor（RxBitsQ（k）,RxBitsQ（k-1））））;SinkBitsQ（k）=or（and（not（xor（RxBitsI（k）,RxBitsQ（k）））,xor（RxBitsQ（k）,RxBitsQ（k-
1）））,and（xor（RxBitsI（k）,RxBitsQ（k））,xor（RxBitsI（k）,RxBitsI（k-1））））;end
%在输入和输出100字节中寻找最佳时延[C,Lags]=vxcorr（SourceBitsI（10:
110）,SinkBitsI（10:
110））;[MaxC,LocMaxC]=max（C）;
BestLag=Lags（LocMaxC）;
%调整时延
ifBestLag>0
SourceBitsI=SourceBitsI（BestLag+1:
length（SourceBitsI））;SourceBitsQ=SourceBitsQ（BestLag+1:
length（SourceBitsQ））;
elseifBestLag<0
SinkBitsI=SinkBitsI（-BestLag+1:
length（SinkBitsI））;SinkBitsQ=SinkBitsQ（-BestLag+1:
length（SinkBitsQ））;
end
%将序列调整成相同长度TotalBits=min（length（SourceBitsI）,length（SinkBitsI））;TotalBits=TotalBits-20;SourceBitsI=SourceBitsI（10:
TotalBits）;SourceBitsQ=SourceBitsQ（10:
TotalBits）;SinkBitsI=SinkBitsI（10:
TotalBits）;SinkBitsQ=SinkBitsQ（10:
TotalBits）;
Errors=sum（SourceBitsI~=SinkBitsI）+sum（SourceBitsQ~=SinkBitsQ）;BER_MC=Errors/（2*length（SourceBitsI））;
return
functionMCQPSKBER
Eb=22:
0.5:
26;
N0=-50;
ChannelAttenuation=70;
EbN0dB=（Eb-ChannelAttenuation）-N0;EbN0=10.^（EbN0dB./10）;
BER_T=0.5*erfc（sqrt（EbN0））;N=round（100./BER_T）;BER_MC=zeros（size（Eb））;fork=1:
len
升级会员 