QAM调制与解调的MATLAB实现及调制性能分析
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
通信原理课程设计报告书
课题名称 16QAM 调制与解调 的MATLAB 实现及调制性能分析
姓 名
学 号 学 院 通信与电子工程学院
专 业 通信工程 指导教师
李梦醒
2012年 01 月 01日
※※※※※※※※※ ※※ ※
※ ※
※
2009级通信工程专业
通信原理课程设计
16QAM调制与解调的MATLAB实现及调制性能分
析
1 设计目的
(1)掌握16QAM调制与解调的原理。
(2)掌握星座图的原理并能熟悉星座图的应用。
(3)熟悉并掌握MATLAB的使用方法。
(4)通过对16QAM调制性能的分析了解16QAM调制相对于其它调制方式的优缺点。
2 设计原理
正交振幅调制(Quadrature Amplitude Modulation,QAM)是一种振幅和相位联合键控。虽然MPSK和MDPSK等相移键控的带宽和功率方面都具有优势,即带宽占用小和比特噪声比要求低。但是由图1可见,在MPSK体制中,随着
8/ 5π
8/
3π
8/
π
8/
7π
8/
9π
8/
11π
8/
13π
8/15π
图 1 8PSK 信号相位
M 的增大,相邻相位的距离逐渐减小,使噪声容限随之减小,误码率难于保证。为了改善在M 大时的噪声容限,发展出了QAM 体制。在QAM 体制中,信号的振幅和相位作为两个独立的参量同时受到调制。这种信号的一个码元可以表示为
0()cos() (1)k k k s t A t kT t k T ωθ=+<≤+ (2—1)
式中:k=整数;k A 和k θ分别可以取多个离散值。 式(2—1)可以展开为
00()cos cos sin sin k k k k k s t A t A t θωθω=- (2—2)
令 X k = A k cos k , Y k = -A k sin
k
则式(2—1)变为
00()cos sin k k k s t X t Y t ωω=+ (2—3)
k X 和k Y 也是可以取多个离散的变量。从式(2—3)看出,()k s t 可以看作是两个
正交的振幅键控信号之和。
在式(2—1)中,若k 值仅可以取
/4和-/4,A k 值仅可以取+A 和-A ,
则此QAM 信号就成为QPSK 信号,如图2所示:
图2 4QAM 信号矢量图
所以,QPSK 信号就是一种最简单的QAM 信号。有代表性的QAM 信号是16进制的,记为16QAM ,它的矢量图示于下图中:
图3 16QAM 信号矢量图
图中用黑点表示每个码元的位置,并且示出它是由两个正交矢量合成的。类似地,有64QAM 和256QAM 等QAM 信号,如图4、图5所示。它们总称为MQAM
调制。由于从其矢量图看像是星座,故又称星座调制。
16QAM 信号的产生方法主要有两种。第一种是正交调幅法,即用两路独立的正交4ASK 信号叠加,形成16QAM 信号,如图6所示。第二种方法是复合相
图4 64QAM 信号矢量图
图5 256QAM 信号矢量图
图6 正交调幅法
移法,它用两路独立的QPSK信号叠加,形成16QAM信号,如图7所示。图中
图7 复合相移法
虚线大圆上的4个大黑点表示一个QPSK信号矢量的位置。在这4个位置上可以叠加上第二个QPSK矢量,后者的位置用虚线小圆上的4个小黑点表示。
3设计过程
设计思路
由设计原理中可知MQAM调制又称为星座调制,故我们在设计16QAM调制系统时就可以星座图来进行编程。下面我们就借用如图8所示的星座图设计一个16QAM调制系统。
图8 16QAM星座
在图中共有16个点,每个点用4个比特表示,代表调制以后的一个矢量位置(这个点拥有唯一的振幅与相位)。因此我们可以把横轴看作是实轴,纵轴看作虚轴。由于每个点与跟它相邻的四个点是等距,且设为2,则每个点都可用一个虚数进行表示。例如点0000可用-1-j表示,这个虚数的模就是相当于16QAM 信号的振幅,相角就相当于16QAM信号的相位。所以16QAM的调制过程就可以用如下语句进行描述:
if A(1,b:c)==[0 0 0 0]
B(k)=-1-1i;
elseif A(1,b:c)==[0 0 0 1]
B(k)=-3-1i;
elseif A(1,b:c)==[0 0 1 0]
B(k)=-1-3i;
elseif A(1,b:c)==[0 0 1 1]
B(k)=-3-3i;
elseif A(1,b:c)==[0 1 0 0]
B(k)=1-1i;
elseif A(1,b:c)==[0 1 0 1]
B(k)=1-3i;
elseif A(1,b:c)==[0 1 1 0]
B(k)=3-1i;
elseif A(1,b:c)==[0 1 1 1]
B(k)=3-3i;
elseif A(1,b:c)==[1 0 0 0]
B(k)=-1+1i;
elseif A(1,b:c)==[1 0 0 1]
B(k)=-1+3i;
elseif A(1,b:c)==[1 0 1 0]
B(k)=-3+1i;
elseif A(1,b:c)==[1 0 1 1]
B(k)=-3+3i;
elseif A(1,b:c)==[1 1 0 0]
B(k)=1+1i;
elseif A(1,b:c)==[1 1 0 1]
B(k)=3+1i;
elseif A(1,b:c)==[1 1 1 0]
B(k)=1+3i;
elseif A(1,b:c)==[1 1 1 1]
B(k)=3+3i;
end
当调制以后的信号经过信道加噪以后,我们必须对其进行解调。由于加噪了的缘故,调制以后的信号不再是原来的信号,而应该有不同。因此在解调时不能简单的将上述过程逆转,即不能由-1-j就判断为0000。而应该对虚数的实部和虚部设定一个范围后再进行判断。这个范围边界的选取原理我们可以借用量化的概念,取相邻虚数的实部的平均数和虚部的平均数。以下为16QAM解调过程的程