在matlab上的的QPSK调制与解调仿真概要Word文档下载推荐.docx
《在matlab上的的QPSK调制与解调仿真概要Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《在matlab上的的QPSK调制与解调仿真概要Word文档下载推荐.docx(20页珍藏版)》请在冰豆网上搜索。
3.3、先通过瑞利衰落信道再通过高斯信道的仿真8
4致谢9
参考文献10
附录11
1.前言
1.1QPSK系统的应用背景简介
QPSK是英文QuadraturePhaseShiftKeying的缩略语简称,意为正交相移键控,是一种数字调制方式。
在19世纪80年代初期,人们选用恒定包络数字调制。
这类数字调制技术的优点是已调信号具有相对窄的功率谱和对放大设备没有线性要求,不足之处是其频谱利用率低于线性调制技术。
19世纪80年代中期以后,四相绝对移相键控(QPSK)技术以其抗干扰性能强、误码性能好、频谱利用率高等优点,广泛应用于数字微波通信系统、数字卫星通信系统、宽带接入、移动通信及有线电视系统之中。
1.2QPSK实验仿真的意义
通过完成设计内容,复习QPSK调制解调的基本原理,同时也要复习通信系统的主要组成部分,了解调制解调方式中最基础的方法。
了解QPSK的实现方法及数学原理。
并对“通信”这个概念有个整体的理解,学习数字调制中误码率测试的标准及计算方法。
同时还要复习随机信号中时域用自相关函数,频域用功率谱密度来描述平稳随机过程的特性等基础知识,来理解高斯信道中噪声的表示方法,以便在编程中使用。
理解QPSK调制解调的基本原理,并使用MATLAB编程实现QPSK信号在高斯信道和瑞利衰落信道下传输,以及该方式的误码率测试。
复习MATLAB编程的基础知识和编程的常用算法以及使用MATLAB仿真系统的注意事项,并锻炼自己的编程能力,通过编程完成QPSK调制解调系统的仿真,以及误码率测试,并得出响应波形。
在完成要求任务的条件下,尝试优化程序。
通过本次实验,除了和队友培养了默契学到了知识之外,还可以将次实验作为一种推广,让更多的学生来深入一层的了解QPSK以至其他调制方式的原理和实现方法。
可以方便学生进行测试和对比。
足不出户便可以做实验。
1.3实验平台和实验内容
1.3.1实验平台
本实验是基于Matlab的软件仿真,只需PC机上安装MATLAB6.0或者以上版本即可。
(本实验附带基于MatlabSimulink(模块化)仿真,如需使用必须安装simulink模块)
1.3.2实验内容
1.构建一个理想信道基本QPSK仿真系统,要求仿真结果有
a.基带输入波形及其功率谱
b.QPSK信号及其功率谱
c.QPSK信号星座图
2.构建一个在AWGN(高斯白噪声)信道条件下的QPSK仿真系统,要求仿真结果有
a.QPSK信号及其功率谱
b.QPSK信号星座图
c.高斯白噪声信道条件下的误码性能以及高斯白噪声的理论曲线,要求所有误码性能曲线在同一坐标比例下绘制
3验可选做扩展内容要求:
构建一个先经过Rayleigh(瑞利衰落信道),再通过AWGN(高斯白噪声)信道条件下的条件下的QPSK仿真系统,要求仿真结果有
b.通过瑞利衰落信道之前和之后的信号星座图,前后进行比较
c.在瑞利衰落信道和在高斯白噪声条件下的误码性能曲线,并和二.2.c中所要求的误码性能曲线在同一坐标比例下绘制
2系统实现框图和分析
2.1、QPSK调制部分
原理框图如图1所示
1(t)=
2(t)=
图2-1
原理分析:
基本原理及系统结构
QPSK与二进制PSK一样,传输信号包含的信息都存在于相位中。
的别的载波相位取四个等间隔值之一,如л/4,3л/4,5л/4,和7л/4。
相应的,可将发射信号定义为
0≤t≤T
Si(t)=
0。
,其他
其中,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.2、QPSK解调部分
原理框图如图2所示:
1(t)同相信道门限=0
2(t)正交信道门限=0
图2-2
QPSK接收机由一对共输入地相关器组成。
这两个相关器分别提供本地产生地相干参考信号
2(t)。
相关器接收信号x(t),相关器输出地x1和x2被用来与门限值0进行比较。
如果x1>
0,则判决同相信道地输出为符号1;
如果x1<
0,则判决同相信道的输出为符号0。
;
类似地。
如果正交通道也是如此判决输出。
最后同相信道和正交信道输出这两个二进制数据序列被复加器合并,重新得到原始的二进制序列。
在AWGN信道中,判决结果具有最小的负号差错概率。
3实验结果及分析
根据图1和图2的流程框图设计仿真程序,得出结果并且分析如下:
3.1、理想信道下的仿真
实验结果如图3所示
图3-1
实验结果分析:
如图上结果显示,完成了QPSK信号在理想信道上的调制,传输,解调的过程,由于调制过程中加进了载波,因此调制信号的功率谱密度会发生变化。
并且可以看出调制解调的结果没有误码。
3.2、高斯信道下的仿真
结果如图4所示:
图3-2
由图4可以得到高斯信道下的调制信号,高斯噪声,调制输出功率谱密度曲线和QPSK信号的星座图。
在高斯噪声的影响下,调制信号的波形发生了明显的变化,其功率谱密度函数相对于图1中的调制信号的功率谱密度只发生了微小的变化,原因在于高斯噪声是一个均值为0的白噪声,在各个频率上其功率是均匀的,因此此结果是真确的。
星座图反映可接收信号早高斯噪声的影响下发生了误码,但是大部分还是保持了原来的特性。
3.3、先通过瑞利衰落信道再通过高斯信道的仿真。
实验结果如图5所示:
图3-3
由图5可以得到瑞利衰落信道前后的星座图,调制信号的曲线图及其功率谱密度。
最后显示的是高斯信道和瑞利衰落信道的误码率对比。
由图可知瑞利衰落信道下的误码率比高斯信道下的误码率高。
至此,仿真实验就全部完成。
4致谢
感谢指导老师赵老师对我们的指导,帮我们解决了不少的问题。
也感谢队友之间的相互合作。
希望以后再接再厉,努力学习。
参考文献
[1]曹志刚,钱亚生.QPSK及其解调.北京:
清华大学出版社,1992
[2]黄正.BPSK,QPSK及其解调.电光系统,2003,103
(1):
43-47
[3]黄文梅,熊桂林,杨勇.信号分析与处理—MATLAB语言及应用.长沙:
国防科技大学出版社,2000.
[4]邓华..MATLAB通信仿真及应用实例详解:
北京:
人民邮电出版社,2003
[5]李贺冰.Simulink通信仿真教程.北京:
国防工业出版社,2006.5
附录
%调相法
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
end
t1=[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);
I=i(101:
800);
Q=q(1:
700);
QPSK=sqrt(1/2).*I+sqrt(1/2).*Q;
QPSK_n=(sqrt(1/2).*I+sqrt(1/2).*Q)+n0;
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*b