QPSK调制解调的仿真文档格式.docx
《QPSK调制解调的仿真文档格式.docx》由会员分享，可在线阅读，更多相关《QPSK调制解调的仿真文档格式.docx（34页珍藏版）》请在冰豆网上搜索。
因为在Matlab中应用程序对QPSK调制解调系统进行仿真，一方面降低了系统设计的复杂性，并且有效的克服了电子瓶颈的问题。
因为在电子电路中，总体电路对参数设置是非常敏感的，一旦参数设置出现小的偏差，将会影响到整个电路的结果。
另一方面程序的仿真过程中可以将主程序分为很多个子程序，将它们逐个进行仿真，这样不会影响到整体结果，可以逐步进行调试。
结合以上的比较论证以及自身的能力，最终我决定用程序法来完成此次对QPSK调制解调系统的仿真。
2仿真原理
四相相位调制解调是利用载波的四种不同相位差来表征输入的数字信息,是四进制移相键控。
QPSK是在M=4时的调相技术,它规定了四种载波相位,分别为45°
135°
225°
275°
调制器输入的数据是二进制数字序列,为了能和四进制的载波相位配合起来,则需要把二进制数据变换为四进制数据,这就是说需要把二进制数字序列中每两个比特分成一组,共有四种组合,即00,01,10,11,其中每一组称为双比特码元。
每一个双比特码元是由两位二进制信息比特组成,它们分别代表四进制四个符号中的一个符号。
QPSK中每次调制可传输2个信息比特,这些信息比特是通过载波的四种相位来传递的。
解调器根据星座图及接收到的载波信号的相位来判断发送端发送的信息比特。
数字调制用“星座图”来描述,星座图中定义了一种调制技术的两个基本参数:
（1）信号分布;
（2）与调制数字比特之间的映射关系。
星座图中规定了星座点与传输比特间的对应关系,这种关系称为“映射”,一种调制技术的特性可由信号分布和映射完全定义,即可由星座图来完全定义[3]。
在QPSK调制中,QPSK信号可以看作两个载波正交的2PSK调制器构成。
串/并变换器将输入的二进制序列分为速率减半的两个并行的双极性序列,然后分别对sinωct和cosωct调制,相加后得到QPSK调制信号。
QPSK同相支路和正交支路可分别采用相干解调方式解调,得到I（t）和Q（t）。
经抽样判决和并/串变换器,将上、下支路得到的并行数据恢复成串行数据。
2．1QPSK调制原理
在QPSK调制中,QPSK信号可以看作两个载波正交的2PSK调制器构成。
调制原理框图如图2.1所示
图2.1QPSK调制原理框图
原理分析：
基本原理及系统结构QPSK与二进制PSK一样，传输信号包含的信息都存在于相位中。
的别的载波相位取四个等间隔值之一，如л/4,3л/4,5л/4,和7л/4。
相应的，可将发射信号定义为
0≤t≤T
Si（t）＝
0（2.1.1）
其中，i＝1，2，2，4；
E为发射信号的每个符号的能量，T为符号持续时间，载波频率f等于nc/T，nc为固定整数。
每一个可能的相位值对应于一个特定的二位组。
例如，可用前述的一组相位值来表示格雷码的一组二位组：
10，00，01，11。
下面介绍QPSK信号的产生和检测。
如果a为典型的QPSK发射机框图。
输入的二进制数据序列首先被不归零（NRZ）电平编码转换器转换为极性形式，即负号1和0分别用
和－
表示。
接着，该二进制波形被分接器分成两个分别由输入序列的奇数位偶数位组成的彼此独立的二进制波形，这两个二进制波形分别用a1（t），和a2（t）表示。
容易注意到，在任何一信号时间间隔内a1（t），和a2（t）的幅度恰好分别等于Si1和Si2，即由发送的二位组决定。
这两个二进制波形a1（t），和a2（t）被用来调制一对正交载波或者说正交基本函数：
1（t）＝
，
2（t）＝
。
这样就得到一对二进制PSK信号。
1（t）和
2（t）的正交性使这两个信号可以被独立地检测。
最后，将这两个二进制PSK信号相加，从而得期望的QPSK。
2.2QPSK解调原理
在QPSK解调中,正交支路和同相支路分别设置两个相关器（或匹配滤波器）,得到I（t）和Q（t）,经电平判决和并/串变换后即可恢复原始信息从发射机发射的已调信号经过传输媒质传播到接收端,接收机接收到的己调信号为:
SQPSK（t）=I（t）cosωct+Q（t）sinωct（2.2.1）
I（t）、Q（t）分别为同相和正交支路,ωc为载波频率,那么相干解调后,同相支路相乘可得:
Ii（t）=SQPSK（t）cosωct=[I（t）cosωct+Q（t）sinωct]cosωct
=I（t）cos2wct+Q（t）sinwct/2
=I（t）/2-I（t）cos2ωct+Q（t）sin2ωct（2.2.2）正交支路相乘可得:
Qq（t）=SQPSK（t）sinωct=[I（t）cosωct+Q（t）sinωct]sinωct
=I（t）sinωct×
cosωct+Q（t）sin2ωct
=I（t）sin2ωct/2+Q（t）-Q（t）cos2ωct（2.2.3）
经过低通滤波器可得:
Ii（t）=0.5I（t）,Qq（t）=0.5Q（t）（2.2.4）
原理框图如图2.2所示：
1（t）同相信道门限＝0
2（t）正交信道门限＝0
图2.2QPSK解调原理框图
QPSK接收机由一对共输入地相关器组成。
这两个相关器分别提供本地产生地相干参考信号
2（t）。
相关器接收信号x（t），相关器输出地x1和x2被用来与门限值0进行比较。
如果x1>
0,则判决同相信道地输出为符号1；
如果x1<
0,则判决同相信道的输出为符号0。
；
类似地。
如果正交通道也是如此判决输出。
最后同相信道和正交信道输出这两个二进制数据序列被复加器合并，重新得到原始的二进制序列。
在AWGN信道中，判决结果具有最小的负号差错概率。
2.3调制解调系统总框图
图2.3QPSK调制与解调系统总框图
3仿真与分析
3.1程序仿真流程图
程序仿真流程图是程序仿真的步骤及大体框架，基于流程图，我们可以有效的避免程序次序的先后问题，将主程序分割成很多个子程序来逐一实现，严格按照流程图的次序进行编程仿真。
最终将各个子程序综合起来实现最终的程序仿真。
流程图如图3.1所示
图3.1程序仿真流程图
3.2QPSK调制方式的Matlab仿真
I路信号是用余弦载波，由2进制数据流的奇数序列组成；
Q路信号用正弦载波，由2进制数据流的偶数序列组成。
下面的a是Idata，b就是Qdata，它们分布与各自的载波相乘分别输出I路信号和Q路信号。
I路信号加上Q路信号就是QPSK输出信号。
当I路载波信号是0相位时为1，是180°
相位时为0；
当Q路载波信号是0相位时为1，是180°
相位时为0。
调制仿真图如图3.2所示：
图3.2matlab调制仿真图
3.3QPSK解调方式Matlab仿真
正交支路和同相支路分别设置两个相关器（或匹配滤波器）,得到I（t）和Q（t）,经电平判决和并/串变换后即可恢复原始信息。
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）;
图3.3matlab解调仿真图
根据图2.1和图2.2的流程框图设计仿真程序，得出结果并且分析如下：
3.4理想信道下的仿真
仿真结果出现的波形依次是基带信号，基带信号的功率谱密度，调制信号，调制信号的功率谱密度，解调输出，解调输出功率谱密度。
实验结果如图3.4所示
图3.4理想信道下的QPSK仿真结果
图3.4结果显示，完成了QPSK信号在理想信道上的调制，传输，解调的过程，由于调制过程中加进了载波，因此调制信号的功率谱密度会发生变化。
并且可以看出调制解调的结果没有误码。
3.5高斯信道下的仿真
在高斯信道下的仿真结果出现的波形依次是调制信号，调制信号功率谱密度，高斯噪声曲线，QPSK信号星座图（Awgn）。
结果如图3.4所示：
图3.5高斯信道下的QPSK仿真结果
由图3.5可以得到高斯信道下的调制信号，高斯噪声，调制输出功率谱密度曲线和QPSK信号的星座图。
在高斯噪声的影响下，调制信号的波形发生了明显的变化，其功率谱密度函数相对于图3.4中的调制信号的功率谱密度只发生了微小的变化，原因在于高斯噪声是一个均值为0的白噪声，在各个频率上其功率是均匀的，因此此结果是真确的。
星座图反映可接收信号早高斯噪声的影响下发生了误码，但是大部分还是保持了原来的特性。
3.6先通过瑞利衰落信道再通过高斯信道的仿真
信号先通过瑞丽衰落信道再通过高斯信道的仿真后出现的波形依次是调制信号（Ray+Awgn）,调试信号功率谱密度（Ray）,QPSK信号星座图（Awgn+Ray）,误码率曲线。
仿真结果如图3.6所示：
图3.6先通过瑞利衰落信道再通过高斯信道的仿真结果
由图3.6可以得到瑞利衰落信道前后的星座图，调制信号的曲线图及其功率谱密度。
最后显示的是高斯信道和瑞利衰落信道的误码率对比。
由图可知瑞利衰落信道下的误码率比高斯信道下的误码率高。
至此，仿真实验就全部完成。
4结束语
本次试验是基于matlab环境下的QPSK调制解调系统的仿真。
由此次仿真实验可知，由于高斯信道和瑞利信道的影响，波形发生了明显的变化，功率谱密度也发生变化，星座图可反映在噪声影响下产生的误码。
从仿真解调前与调制后的图对比可知仿真正确。
通过图3.3可知瑞利信道下的误码率和误比特率明显高于高斯信道，并与理论值基本符合。
随着通信事业的发展，通信系统的设计也会越来越复杂，通过计算机的仿真，可以大大地降低通信过程实验成本。
本文设计出一个QPSK仿真模型，以衡量QPSK在理想信道，高斯白噪声信道，以及先通过瑞利再通过高斯信道的性能，通过仿真，可以更好地了解QPSK系统的工作原理，而且为硬件的研制提供一定的参考作用。
这次的通信专业方向设计让我把通信原理的一些内容又重新复习了一遍，中间出了很多问题，比如用matlab7.0版本就不能把几个图整合到一起，也就是不能使用subplot函数，会出现Undefinedfunctionorvariable"
hanalysisparammenu"
的错误.而我后来又换成了matlab7.8就可以实现该功能。
我把程序发给老师后老师细心的帮我进行了排错，我错误的以为依靠高版本的优化功能把在低版本不能实现的功能强制执行，所以我还有很多地方需要学习，包括对误码率进行计算的两个函数我也是借鉴的网上的程序。
这次方向设计让我学会了很多，也认识到了自己还有很多方面的欠缺，希望在今后的学习和工作中好好努力，尽快的弥补自己的不足。
致谢
一个月的课程设计很快就过去了，在这期间我学习和掌握了使用MATLAB的工具箱simulink对数字通信系统的仿真。
回想起这期间的经历，感慨万千。
课程设计是对大学四年所学的知识的一次总结、一次检验。
同时也是锻炼和进一步提高自己自学能力的一次宝贵机会。
在做本次设计的初期阶段，由于自己在很的多知识点上有所遗忘，因此设计的进度相当缓慢，特别是在原理分析和simulink仿真上更是裹足不前。
我重点复习并新学习了相关章节的知识,并很快理解了要点，进度由此快速展开。
设计后期阶段，在设计思路上卡了壳。
此时，查阅了很多资料并多次请教老师，使我得以进行软件调试从而验证了设计的正确可行性并按时完成设计，在此特别感谢魏瑞老师的指导，此外还要感谢学校教研室给我们这次宝贵而难得的机会，也感谢帮助和支持我的同学。
在此表示衷心的感谢！
由于理论知识有限、文中难免出现不足之处，敬请各位老师批评指正。
参考文献
［1］张志勇等.《精通matlab6.5》[J].北京:
北京航空航天大学出版社,2003.3.
［2］JohnG.ProakisMasoudSalehiGerhardBauch.著刘树棠译comtemporarycommunicationsystemsusingmatlabandsimulink.北京:
电子工业出版社,2005.4.
［3］樊昌信,张甫詡、徐炳祥、吴成柯.《通信原理》[M].北京:
国防工业出版社,2001.5.
［4］王立宁,乐光新,詹菲等.《MATLAB与通信仿真》[M].北京:
人民邮电出版社,1999
［5］陈怀琛，吴大正，高西金《matlab及在信息课程中的应用》[J]北京：
电子工业出版社,2006
［6］楼天顺、刘晓东、李博涵《基于matlab7.x的系统分析与设计》[M].西安：
西安电子科技大学出版社,2005.5
［7］张肃文主编.《高频电子线路》[M].北京：
高等教育出版社,2004.11
［8］程佩青编著.《数字信号处理教程》[M].北京：
清华大学出版社,2006.3
［9］张辉、曹丽娜、王勇编著.《通信原理辅导》[M]西安.西安电子科技大学出版社,2004.11
[10]徐明远.《MATLAB仿真在通信工程中的应用》[J].西安.西安电子科技大学出版社,2005.6
[11]田丽华等.《编码理论》[J].西安西安电子科技大学出版社,2003.8
附录
%主文件
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%题目：
基于matlab环境下的QPSK仿真
%%%%%%%作者：
杨遇春
%%%%%%%完成日期：
2015-01-14
%%%%%%%邮箱：
710476235@
clear;
%%%%%%%%%%%%%%%%%%%%%
%%%初始化参数
T=1;
%基带信号宽度，也就是频率
fc=10/T;
%载波频率
ml=2;
%调制信号类型的一个标志位（选取2的原因见23行）
nb=100;
%传输的比特数
delta_T=T/200;
%采样间隔
fs=1/delta_T;
%采样频率
SNR=0;
%信噪比
t=0:
delta_T:
nb*T-delta_T;
%限定t的取值范围
N=length（t）;
%采样数
%%%%%%%%%%%%%%%%%%%调制部分
%基带信号的产生
data=randn（1,nb）>
0.5;
%调用一个随机函数（0or1），输出到一个1*100的矩阵
datanrz=data.*2-1;
%变成极性码
data1=zeros（1,nb/delta_T）;
%创建一个1*nb/delta_T的零矩阵
forq=1:
nb
data1（（q-1）/delta_T+1:
q/delta_T）=datanrz（q）;
%将极性码变成对应的波形信号
end
%将基带信号变换成对应波形信号
data0=zeros（1,nb/delta_T）;
data0（（q-1）/delta_T+1:
q/delta_T）=data（q）;
%发射的信号
data2=abs（fft（data1））;
%串并转换，将奇偶位数据分开
idata=datanrz（1:
ml:
（nb-1））;
%将奇偶位分开，因此间隔m1为2
qdata=datanrz（2:
nb）;
%QPSK信号的调制
ich=zeros（1,nb/delta_T/2）;
%创建一个1*nb/delta_T/2的零矩阵，以便后面存放奇偶位数据
fori=1:
nb/2
ich（（i-1）/delta_T+1:
i/delta_T）=idata（i）;
end
forii=1:
N/2
a（ii）=sqrt（2/T）*cos（2*pi*fc*t（ii））;
idata1=ich.*a;
%奇数位数据与余弦函数相乘，得到一路的调制信号
qch=zeros（1,nb/2/delta_T）;
forj1=1:
qch（（j1-1）/delta_T+1:
j1/delta_T）=qdata（j1）;
forjj=1:
b（jj）=sqrt（2/T）*sin（2*pi*fc*t（jj））;
qdata1=qch.*b;
%偶数位数据与余弦函数相乘，得到另一路的调制信号
s=idata1+qdata1;
%将奇偶位数据合并，s即为QPSK调制信号
ss=abs（fft（s））;
%快速傅里叶变换得到频谱
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%瑞利衰落信道和高斯信道
%瑞利衰落信道
ray_ich=raylrnd（0.8,1,nb/2/delta_T）;
ray_qch=raylrnd（0.8,1,nb/2/delta_T）;
Ray_idata=idata1.*ray_ich;
Ray_qdata=qdata1.*ray_qch;
Ray_s=Ray_idata+Ray_qdata;
%%%%%%%%%%%%%%%%%%%%%%%%%%
%高斯信道
s1=awgn（s,SNR）;
%通过高斯信道之后的信号
s11=abs（fft（s1））;
%快速傅里叶变换得到频谱
s111=s1-s;
%高斯噪声曲线
Awgn_s=awgn（Ray_s,SNR）;
%通过高斯信道再通过瑞利衰落信道
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%QPSK解调部分
%解调部分（高斯信道）
idata2=s1.*a;
%这里面其实隐藏了一个串并转换的过程
qdata2=s1.*b;
%对应的信号与正余弦信号相乘
idata3=zeros（1,nb/2）;
%建立1*nb数组，以存放解调之后的信号
qdata3=zeros（1,nb/2）;
%抽样判决的过程，与0作比较，data>
=0,则置1，否则置0
forn=1:
%A1（n）=sum（idata2（（n-1）/delta_T+1:
n/delta_T））;
ifsum（idata2（（n-1）/delta_T+1:
n/delta_T））>
=0
idata3（n）=1;
elseidata3（n）=0;
end
%A2（n）=sum（qdata2（（n-1）/delta_T+1:
ifsum（qdata
升级会员 