基于Matlab调制与解调系统仿真Word格式文档下载.docx

1、 基于Matlab的QPSK调制与解调系统仿真基本要求：分析QPSK调制、解调原理，画出系统框图,推导星座图（设载波相位为0）；产生随机二进制序列，作为待传输的数据,数据的速率是5Mbit/s；载波频率为50MHz,完成QPSK调制，分别画出两路信号经过乘法器后的时域波形，并给出分析、解释；产生加性高斯白噪声，画出时域波形及频谱,将加性高斯白噪声作为传输信道的干扰,画出接收端信号的时域图；完成QPSK解调,输出解码序列，画出信号经过乘法器，低通滤波器后的波形.计算系统信噪比、传输误码率扩展要求：计算在不同系统信噪比下，传输误码率;画出系统信噪比，传输误码率的曲线。课 程 设 计 学 生 日 志

2、时间设计内容6-21-6-23熟悉MATLAB软件与数字信号处理得有关资料6-246-28用MATLAB编写并调试程序，实现QPSK调制与解调6-297-1书写课程设计的论文 7-2完善论文，准备答辩7-4答辩周星期一星期二星期三星期四星期五课 程 设 计 考 勤 表课 程 设 计 评 语 表指导教师评语： 成绩: 年 月 日一、设计目的和意义 目的:用MATLAB设计QPSK调制与解调系统. 意义:掌握用MATLAB工具来实现数字基带信号的调制与解调，包括QPSK的调制与解调的原理,工作流程,实现方法等。二、设计原理带通二进制键控系统中,每个码元只能传输1b信息，其频带利用率不高.为了提高频

3、带利用率，最有效的办法是使一个码元传输多个比特的信息。这就是多进制键控体制。多进制数字调制是利用多进制数字基带信号去调制载波的振幅，频率或相位.因此，相应地有多进制数字振幅调制（MASK）,多进制数字频率调制（MPSK）以及多进制数字相位调制（MPSK）等. QPSK（4PSK，正交相移键控）又叫四相绝对相移调制，是最常用的MPSK,分为p/2系统和p/4系统两种。它是利用载波的四种不同相位来表征数字信息。由于每一种载波相位代表两个比特信息，故每个四进制码元又被称为双比特码元。QPSK系统可采用多种调制方式：正交调制法，相位选择法，脉冲插入法。在这篇课程设计中采用了正交调制法进行设计。QPSK

4、系统在用正交调制部分需要进行串/并变换，其中串/并变换电路将QPSK调制的两位编码按比特分开,走上下两路,分成的两路序列速率减半，电平发生器分别产生双极性二电平信号I（t）和Q（t）,然后各自去调制相互正交的正弦波,再进行矢量合成，即得到QPSK信号。 图2-1 QPSK调制电路在解调部分可以用两个正交的载波信号实现相干解调。正交路和同相路分别设置两个相关器，得到I（t）和Q（t）,经电平判决和并/串变换即可恢复原始信息。 图22 QPSK解调电路数字调制用“星座图”来描述，星座图中定义了一种调制技术的两个基本参数：（1）信号分布；（2）与调制数字比特之间的映射关系。星座图中规定了星座点与传输

5、比特间的对应关系,这种关系称为“映射，一种调制技术的特性可由信号分布和映射完全定义,即可由星座图来完全定义。 QPSK的星座图相位有0，p/2，p，p/4四种，它们分别对应于信息11,01,00和10。为了能和四进制的载波相位配合起来,则需要把二进制数据变换为四进制数据，这就是说需要把二进制数字序列中每两个比特分成一组,共有四种组合，即00，01,10，11,其中每一组称为双比特码元.每一个双比特码元是由两位二进制信息比特组成，它们分别代表四进制四个符号中的一个符号.QPSK中每次调制可传输2个信息比特，这些信息比特是通过载波的四种相位来传递的。解调器根据星座图及接收到的载波信号的相位来判断发

6、送端发送的信息比特。 图23 QPSK星座图三、详细设计步骤 1 根据QPSK的调制及解调原理及原理框图，MATLAB程序来仿真这个系统应遵循以下几个步骤： 调制部分:（1）串/并变换 （2）极性转换 （3）电平产生 （4）两分路分别与载波相乘 （5）合并两路信号 信道部分：加入高斯白噪声 解调部分:（1）接收到的信号分别乘以正弦信号及余弦信号 （2）两路信号分别进行抽样判决 （3）并/串变换 2 用MATLAB程序来实现QPSK调制与解调系统仿真nb=8; %生成二进制序列的长度ml =2; %分奇偶位SNR=10；data=randint（1,nb，2）； % 调用一个随机函数（0 or

7、1），产生长度为nb的二进制序列bit=； %存放产生的二进制波形for n=1:length（data） if data（n）=0； %这位为0，产生100个0（产生波形） bita=zeros（1,100）； else data（n）=1; %为1，产生100个1 bita=ones（1,100）； end bit=bit bita；endfigure（1）plot（bit） ； title（产生的二进制波形）;grid on;axis（0 100*length（data） -2 2）；% 串并转换,将奇偶位数据分开idata=data（1:ml:（nb1）； 将奇偶位分开，因此间隔m1为

8、2 qdata=data（2：nb）; %偶数位 偶数位正弦调制f1=20; 单路信号的速率为载波的1/20s=qdata；t1=0：2*pi/99:2*pi； %采样频率为基带信号的100倍 cp1=；mod1=；mod11=;bit1=；for n1=1：length（s） if s（n1）=0； cp11=-ones（1,100）； bit11=zeros（1,100）; else s（n1）=1; cp11=ones（1，100）; bit11=ones（1,100）; c1=sin（f1*t1）； cp1=cp1 cp11; 基带信号的抽样波形 mod1=mod1 c1； 载频信号的

9、波形 bit1=bit1 bit11； 正弦分路的波形 bpsksin=cp1。mod1; QPSK调制信号figure（3）subplot（2，1,1）;plot（bit1，LineWidth，1。5）; title（二进制信息grid on；axis（0 100length（s） -2.5 2.5）;subplot（2，1，2）；plot（bpsksin，LineWidth,1.5）； title（正弦分路信号）;axis（0 100*length（s） 2.5 2.5）;基数位余弦调制f2=20;s=idata;t2=0:2pi/99：2*pi; %采样频率为基带信号的100倍cp2=;

10、mod2=；mod21=；bit2=;for n2=1： if s（n2）=0; cp21=-ones（1，100）； bit21=zeros（1，100）; else s（n2）=1； cp21=ones（1,100）; bit21=ones（1，100）; c2=cos（f2t2）； cp2=cp2 cp21； mod2=mod2 c2; bit2=bit2 bit21； %正弦分路的波形 bpskcos=cp2.mod2； %QPSK调制信号figure（2）subplot（2,1，1）；plot（bit2,LineWidth,1.5）；二进制信息axis（0 100length（s）

11、2。5 2.5）；subplot（2，1,2）；plot（bpskcos，LineWidth，1.5）;余弦分路信号axis（0 100length（s） -2.5 2.5）；s=bpsksin+bpskcos; %两路信号合并高斯白噪声s1=awgn（s,SNR）； 通过高斯信道之后的信号s11=abs（fft（s1）； % 快速傅里叶变换得到频谱 s111=s1s; 高斯噪声曲线ss=abs（fft（s111）; 快速傅里叶变换得到频谱h = spectrum.welch; % 类似于C语言的宏定义，方便以下的调用figure（4）subplot（2,1,1） ；plot（s111）,ti

12、tle（加性高斯白噪声时域波形subplot（2,1,2） ;plot（abs（ss）,title（加性高斯白噪声频谱figure（6）psd（h，s111,fs，100）,title（加性高斯白噪声功率谱密度）; 解调部分（高斯信道）idata2=s1。*mod2; %余弦分路相干解调qdata2=s1.mod1； 正弦分路相干解调figure（5） subplot（3,1，1） ;plot（s1）; title（接收端信号的时域图）；subplot（3，1,3）;plot（qdata2）；接收端正弦分路信号subplot（3,1,2） ；plot（idata2）;接收端余弦分路信号）；id

13、ata3=zeros（1，nb/2）； 建立1*nb数组，以存放解调之后的信号qdata3=zeros（1，nb/2）； 抽样判决的过程，与0作比较，data=0,则置1，否则置0nb/2 if sum（idata2（100（n-1）+1） ： （100n）=0 %100个抽样点相加与0比较 idata3（n）=1; 大于等于0为1 else idata3（n）=0； %小于 0为1 if sum（qdata2（100*（n1）+1） ： （100n）=0 qdata3（n）=1; else qdata3（n）=0；end 将判决之后的数据存放进数组demodata=zeros（1,nb）；d

14、emodata（1：（nb1）=idata3; 存放奇数位demodata（2：nb）=qdata3； % 存放偶数位bit3=;for n3=1：length（demodata） if demodata（n3）=0； bit31=zeros（1,100）; else demodata（n3）=1； bit31=ones（1,100）; bit3=bit3 bit31；figure（7）plot（bit3） ; title（解调出的波形axis（0 100*length（demodata） -2 2）;penumber,perat=symerr（data,demodata） %计算误码率解误码

15、个数四、设计结果及分析 图41 产生的随机二进制序列 图4-2 基数位的二进制序列及调制后的波形 图43偶数位的二进制序列及调制后的波形 图4-2和图43的两路信号分别是产生的二进制序列的偶数位和基数位，把信号（10000101）分为偶数位（0011）和基数位（1000）后基带信号的速率减半，即为2。5Mbit/s。基数位的基带信号用余弦信号调制，而偶数位的基带信号用正弦信号调制.当基带信号从1变为0或从0变为1时,正弦或余弦信号的相位翻转180度。图4-4 加性高斯白噪声的时域波形及频谱 图45接收端信号的时域图及正余弦两路信号的波形 图4-6 加性高斯白噪声的功率谱密度 图47 接收端最终

16、解调后的二进制序列QPSK信号的误码率的计算公式: 图4-8 误比特数及误码率 图49 传输误码率的曲线五、体会这次的课程设计中我收获了很多的MATLAB实际操作的技巧，我学会了如何逐步调试MATLAB程序以及如何使用MATLAB的help命令来学会使用它自带的很多有用的函数。通过一次次的调试以及排除错误，我对QPSK的调制与解调原理以及如何实现它有了更加深入的认识.因为我对MATLAB工具箱中的有关通信系统的函数了解很少，在这次的课程设计中，我到了图书馆查阅了大量的关于MATLAB以及通信系统的资料，这教会了我如何在短时间内找到我所需要的资料以及如何去消化这些资料并把有用的内容用到自己的课程设计中。六、参考文献赵静 张瑾 高新科基于MATLAB的通信系统仿真北京航空航天大学出版社程佩青 数字信号处理教程 清华大学出版社邓薇 MATLAB函数速查手册人民邮电出版社张德丰 MATLAB通信工程仿真 机械工业出版社李妍 MATLAB通信仿真开发手册 国防工业出版社