qpsk调制与解调Word文件下载.docx
《qpsk调制与解调Word文件下载.docx》由会员分享，可在线阅读，更多相关《qpsk调制与解调Word文件下载.docx（23页珍藏版）》请在冰豆网上搜索。
基于MATLAB的QPSK仿真设计与实现
1.概述
1.1QPSI系统的应用背景简介
QPSK是英文QuadraturePhaseShiftKeying的缩略语简称，意为正交相移键控，是一种数字调制方式。
在19世纪80年代初期，人们选用恒定包络数字调制。
这类数字调制技术的优点是已调信号具有相对窄的功率谱和对放大设备没有线性要求,不足之处是其频谱利用率低于线性调制技术。
19世纪80年代中期以
后,四相绝对移相键控（QPSK技术以其抗干扰性能强、误码性能好、频谱利用率高等优点，广泛应用于数字微波通信系统、数字卫星通信系统、宽带接入、移动通信及有线电视系统之中。
1.2QPSK实验仿真的意义
通过完成设计内容，复习QPS调制解调的基本原理，同时也要复习通信系统的主要组成部分，了解调制解调方式中最基础的方法。
了解QPS的实现方法
及数学原理。
并对“通信”这个概念有个整体的理解，学习数字调制中误码率测试的标准及计算方法。
同时还要复习随机信号中时域用自相关函数，频域用功率谱密度来描述平稳随机过程的特性等基础知识，来理解高斯信道中噪声的表示方法，以便在编程中使用。
理解QPS调制解调的基本原理，并使用MATLA编程实现QPS信号在高斯信道和瑞利衰落信道下传输，以及该方式的误码率测试。
复习MATLA编程的基础
知识和编程的常用算法以及使用MATLA仿真系统的注意事项，并锻炼自己的编程能力，通过编程完成QPS调制解调系统的仿真，以及误码率测试，并得出响应波形。
在完成要求任务的条件下，尝试优化程序。
通过本次实验，除了和队友培养了默契学到了知识之外，还可以将次实验作为一种推广，让更多的学生来深入一层的了解QPS以至其他调制方式的原理和实现方法。
可以方便学生进行测试和对比。
足不出户便可以做实验。
1.3实验平台和实验内容
1.3.1实验平台
本实验是基于Matlab的软件仿真，只需PC机上安装MATLAB6.0或者以上版本即可。
（本实验附带基于MatlabSimulink（模块化）仿真，如需使用必须安装
simulink模块）
1.3.2实验内容
1.构建一个理想信道基本QPSI仿真系统,要求仿真结果有
a.基带输入波形及其功率谱
b.QPSK信号及其功率谱
c.QPSK信号星座图
2.构建一个在AWGN高斯白噪声）信道条件下的QPSI仿真系统，要求仿真结果有
a.QPSK信号及其功率谱
b.QPSK信号星座图
c.高斯白噪声信道条件下的误码性能以及高斯白噪声的理论曲线，要求所有
误码性能曲线在同一坐标比例下绘制
3验可选做扩展内容要求：
构建一个先经过Rayleigh（瑞利衰落信道），再通过AWGN高斯白噪声）信道条件下的条件下的QPSI仿真系统，要求仿真结果有
b.通过瑞利衰落信道之前和之后的信号星座图，前后进行比较
c.在瑞利衰落信道和在高斯白噪声条件下的误码性能曲线，并和二.2.c中
所要求的误码性能曲线在同一坐标比例下绘制
2.系统实现框图和分析
2.1QPSK调制部分
原理框图如图1所示
二进制
数据序列
2（t）=2tsin（2何）
图2-1QPSK调制
原理分析：
基本原理及系统结构
QPSK与二进制PSK-样，传输信号包含的信息都存在于相位中。
的别的载波相位取四个等间隔值之一，如ji/4,3畀/4,5ji/4,和7畀/4。
相应的，可将发射信号定义为
{
v2E/tcos[2ft（2i1）/4]0<
t<
T
。
，其他
图2-2载波信号
其中，i二1,2,2,4;
E为发射信号的每个符号的能量，T为符号持续时间，载波频率f等于nc/T，nc为固定整数。
每一个可能的相位值对应于一个特定的二位组。
例如，可用前述的一组相位值来表示格雷码的一组二位组：
10，00，01,11。
下面介绍QPSI信号的产生和检测。
如果a为典型的QPSKg射机框图。
输入的二进制数据序列首先被不归零（NRZ电平编码转换器转换为极性形式，即负号1和0分别用.Eb和—.Eb表示。
接着，该二进制波形被分接器分成两个分别由输入序列的奇数位偶数位组成的彼此独立的二进制波形，这两个二进制波形分别用al（t），和a2（t）表示。
容易注意到，在任何一信号时间间隔内al（t）,
和a2（t）的幅度恰好分别等于Si1和Si2，即由发送的二位组决定。
这两个二进制波形al（t），和a2（t）被用来调制一对正交载波或者说正交基本函数：
1（t）=.'
2tcos（2何），2（t）=2tsin（2何）。
这样就得到一对二进制PSK信号。
1（t）和2（t）的正交性使这两个信号可以被独立地检测。
最后，将这两个二进制PSK信号相加，从而得期望的QPSK
2.2QPSK解调部分
，原理框图如图2所示:
2（t）正交信道门限=0
图2-3QPSK的解调
QPSK接收机由一对共输入地相关器组成。
这两个相关器分别提供本地产生地相干参考信号1（t）和2（t）。
相关器接收信号x（t），相关器输出地x1和x2被用来与门限值0进行比较。
如果x1>
0,则判决同相信道地输出为符号1;
如
果x1<
0，则判决同相信道的输出为符号0。
类似地。
如果正交通道也是如此判决输出。
最后同相信道和正交信道输出这两个二进制数据序列被复加器合并，重
新得到原始的二进制序列。
在AWG信道中，判决结果具有最小的负号差错概率。
3、实验结果及分析
根据图2-1和图2-2的流程框图设计仿真程序，得出结果并且分析如下:
3.1、理想信道下的仿真，实验结果如图3所示
O■60
甚帝信号功率谱密空
100
图3-1理想信道下的仿真
实验结果分析:
如图上结果显示，完成了QPSK信号在理想信道上的调制，传输，解调的过程，由于调制过程中加进了载波，因此调制信号的功率谱密度会发生变化。
并且可以看出调制解调的结果没有误码。
3.2、高斯信道下的仿真，结果如图4所示:
00
5
IJj
GPSK恬号星座图<
.^n>
辰斯操岁曲裁
5.，
n
-it-
T—
十理主愷〔址射端）*实尿｛接脱端、
-1n
图3-2高斯信道下的仿真
由图3-2可以得到高斯信道下的调制信号，高斯噪声，调制输出功率谱密度曲线和QPSK信号的星座图。
在高斯噪声的影响下，调制信号的波形发生了明显的变化，其功率谱密度函
数相对于图1中的调制信号的功率谱密度只发生了微小的变化，原因在于高斯噪
声是一个均值为0的白噪声，在各个频率上其功率是均匀的，因此此结果是真确
的。
星座图反映可接收信号早高斯噪声的影响下发生了误码，但是大部分还是保
持了原来的特性。
3.3、先通过瑞利衰落信道再通过高斯信道的仿真。
实验结果如
图5所示：
010D70G3004CI0600
匚
—TI
L'
J
020切6DBD100
润制信号功率谱港度（R钟）
QPSKf^号星座图（期£
汁沖昕》r,,.
图3-3先通过瑞利衰落信道再通过高斯信道的仿真实验结果分析：
由图3-3可以得到瑞利衰落信道前后的星座图，调制信号的曲线图及其功率
谱密度。
最后显示的是高斯信道和瑞利衰落信道的误码率对比。
由图可知瑞利衰
落信道下的误码率比高斯信道下的误码率高。
至此，仿真实验就全部完成。
参考文献:
[1]陈杰等.MATLAB宝典.电子工业出版社
[2]刘波，文忠，曾涯编.MATLAB信号处理•北京电子工业出版社
[3]万永革•数字信号处理的MATLAB实现•北京科学出版社
[4]网上资料
附录
MATLAB程序
%调相法
clearall
closeall
t=[-1:
0.01:
7-0.01];
tt=length（t）;
x1=ones（1,800）;
fori=1:
tt
if（t（i）>
=-1&
t（i）<
=1）|（t（i）>
=5&
=7）;
x1（i）=1;
elsex1（i）=-1;
end
t仁[0:
8-0.01];
t2=0:
7-0.01;
t3=-1:
7.1-0.01;
t4=0:
8.1-0.01;
tt1=length（t1）;
x2=ones（1,800）;
tt1
if（t1（i）>
=0&
t1（i）<
=2）|（t1（i）>
=4&
=8）;
x2（i）=1;
elsex2（i）=-1;
f=0:
0.1:
1;
xrc=0.5+0.5*cos（pi*f）;
y1=conv（x1,xrc）/5.5;
y2=conv（x2,xrc）/5.5;
n0=randn（size（t2））;
f1=1;
i=x1.*cos（2*pi*f1*t）;
q=x2.*sin（2*pi*f1*t1）;
l=i（101:
800）;
Q=q（1:
700）;
QPSK=sqrt（1/2）.*l+sqrt（1/2）.*Q;
QPSK_n=（sqrt（1/2）.*l+sqrt（1/2）.*Q）+nO;
n1=randn（size（t2））;
i_rc=y1.*cos（2*pi*f1*t3）;
q_rc=y2.*sin（2*pi*f1*t4）;
I_rc=i_rc（101:
Q_rc=q_rc（1:
QPSK_rc=（sqrt（1/2）.*I_rc+sqrt（1/2）.*Q_rc）;
QPSK_rc_n1=QPSK_rc+n1;
figure
（1）subplot（4,1,1）;
plot（t3,i_rc）;
axis（[-18-11]）;
ylabel（'
a序列'
）;
subplot（4,1,2）;
plot（t4,q_rc）;
b序列'
subplot（4,1,3）;
plot（t2,QPSK_rc）;
合成序列'
subplot（4,1,4）;
plot（t2,QPSK_rc_n1）;
加入噪声'
%设定T=1,加入高斯噪声
%调制
bit_in=randint（1e3,1,[01]）;
bit_I=bit_in（1:
2:
1e3）;
bit_Q=bit_in（2:
data_I=-2*bit_I+1;
data_Q=-2*bit_Q+1;
data_I1=repmat（data_I'
20,1）;
data_Q1=repmat（data_Q'
1e4
data_I2（i）=data_I1（i）;
data_Q2（i）=data_Q1（i）;
end;
data」2_rc=conv（data」2,xrc）/5.5;
data_Q2_rc=conv（data_Q2,xrc）/5.5;
t1=0:
1e3+0.9;
nO=rand（size（t1））;
I_rc=data_l2_rc.*cos（2*pi*f1*t1）;
Q_rc=data_Q2_rc.*sin（2*pi*f1*t1）;
QPSK_rc_nO=QPSK_rc+nO;
%解调
I_demo=QPSK_rc_n0.*cos（2*pi*f1*t1）;
Q_demo=QPSK_rc_n0.*sin（2*pi*f1*t1）;
%低通滤波
I_recover=conv（I_demo,xrc）;
Q_recover=conv（Q_demo,xrc）;
I=I_recover（11:
10010）;
Q=Q_recover（11:
0.05:
1e3-0.05;
t3=0:
1e3-0.1;
%抽样判决
data_recover=[];
20:
10000
data_recover=[data_recoverI（i:
1:
i+19）Q（i:
i+19）];
bit_recover=[];
20000
ifsum（data_recover（i:
i+19））>
data_recover_a（i:
i+19）=1;
bit_recover=[bit_recover1];
else
i+19）=-1;
bit_recover=[bit_recover-1];
error=0;
dd=-2*bit_in+1;
ddd=[dd'
];
ddd仁repmat（ddd,20,1）;
2e4
ddd2（i）=ddd1（i）;
1e3
ifbit_recover（i）~=ddd（i）
error=error+1;
p=error/1000;
figure
（1）
subplot（2,1,1）;
plot（t2,ddd2）;
axis（[0100-22]）;
title（'
原序列'
解调后序列’）;
subplot（2,1,2）;
plot（t2,data_recover_a）;
%设定T=1,不加噪声
t=0:
I_demo=QPSK_rc.*cos（2*pi*f1*t1）;
Q_demo=QPSK_rc.*sin（2*pi*f1*t1）;
ddd=-2*bit_in+1;
ddd1=repmat（ddd'
10,1）;
subplot（4,1,1）;
plot（t3,l）;
axis（[020-66]）;
subplot（4,1,2）;
plot（t3,Q）;
subplot（4,1,3）;
plot（t2,data_recover）;
subplot（4,1,4）;
plot（t,ddd2）;
%QPSK误码率分析
SNRindB仁0:
10;
SNRindB2=0:
length（SNRindB1）
[pb,ps]=cm_sm32（SNRindB1（i））;
smld_bit_err_prb（i）=pb;
smld_symbol_err_prb（i）=ps;
length（SNRindB2）
SNR=exp（SNRindB2（i）*log（10）/10）;
theo_err_prb（i）=Qfunct（sqrt（2*SNR））;
end;
QPSK误码率分析'
semilogy（SNRindB1,smld_bit_err_prb,'
*'
axis（[01010e-81]）;
holdon;
%semilogy（SNRindB1,smld_symbol_err_prb,'
o'
semilogy（SNRindB2,theo_err_prb）;
legend（'
仿真比特误码率'
'
理论比特误码率'
holdoff;
function[y]=Qfunct（x）
y=（1/2）*erfc（x/sqrt
（2））;
function[pb,ps]=cm_sm32（SNRindB）N=10000;
E=1;
SNR=10A（SNRindB/10）;
sgma=sqrt（E/SNR）/2;
s00=[10];
s0仁[01];
si仁[-10];
s10=[0-1];
N
dsource1（i）=[1011000101101011];
numofsymbolerror=0;
numofbiterror=0;
n=sgma*randn（size（sOO））;
if（（dsource1（i）==0）&
（dsource2（i）==0））
r=sOO+n;
elseif（（dsource1（i）==0）&
（dsource2（i）==1））
r=s01+n;
elseif（（dsource1（i）==1）&
r=s10+n;
r=s11+n;
cOO=dot（r,sOO）;
cO仁dot（r,s01）;
c10=dot（r,s10）;
c11=dot（r,s11）;
c_max=max（[cOOc01c10c11]）;
if（c00==c_max）
decis1=0;
decis2=0;
elseif（c01==c_max）
decis2=1;
elseif（c10==c_max）
decis1=1;
decis1=1;
symbolerror=0;
if（decis1~=dsource1（i））
numofbiterror=numofbiterror+1;
symbolerror=1;
end;
if（decis2~=dsource2（i））
symbolerror=1;
if（symbolerror==1）
numofsymbolerror=numofsymbolerror+1;
ps=numofsymbolerror/N;
pb=numofbiterror/（2*N）;
升级会员 