基于某Matlab的QPSK调制解调仿真设计与研究.doc

word基于Matlab/Simulink 的QPSK调制解调仿真设计与研究目录摘要2第一章前言21.1 专业设计任务与要求21.2 Matlab简介21.3 Matlab下的simulink简介31.4 通信系统模型3第二章QPSK调制42.1 QPSK介绍42.2 QPSK调制原理42.2.1 相乘法42.2.2 选择法52.3 QPSK调制原理框图62.4 QPSK调制方式的Matlab仿真62.5 QPSK调制方式Matlab-simulink仿真72.5.1 simulink调制建模72.5.2 simulink调制仿真结果8第三章QPSK解调13QPSK解调原理133.2 QPSK解调原理框图13QPSK解调方式Matlab仿真13QPSK解调方式的Matlab-simulink仿真143.4.1 QPSK解调建模143.4.2 传输信道16仿真结果163.5 仿真结果分析18第四章QPSK通信系统性能分析19第五章结论19参考文献20附录20摘要正交相移键控QPSK，是一种数字调制方式。QPSK技术具有抗干扰能力好、误码率低、频谱利用效率高等一系列优点。论文主要介绍了正交相移键控(QPSK)的概况，以与正交相移键控QPSK的调制解调概念和原理，利用Matlab中M文件和Simulink模块对QPSK的调制解调系统进展了仿真，对QPSK在高斯白噪声信道中的性能进展了，分析了解Simulink中涉与到QPSK的各种模块的功能。【关键词】Matlab QPSK Simulnk 仿真第一章 前 言专业设计任务与要求1了解并掌握QPSK调制与解调的根本原理；2在通信原理课程的根底上设计与分析简单的通信系统；编写程序进展仿真，根据实验结果能分析所设计系统的性能。4学习MATLAB的根本知识，熟悉MATLAB集成环境下的Simulink的仿真平台。5利用通信原理相关知识在仿真平台中设计QPSK调制与解调仿真系统并用示波器观察解调后的波形6在指导教师的指导下，独立完成课程设计的全部容，能正确的阐述和分析设计和实验结果。1.2 Matlab简介MATLAB是MATrix LABoratory的缩写，是一款由美国Math Works公司出品的商业数学软件。MATLAB 是一种用于算法开发、数据可视化、数据分析以与数值计算的高级技术计算语言和交互式环境。除了矩阵运算、绘制函数/数据图像等常用功能外，MATLAB还可以用来创建用户界面与与调用其它语言(包括C，C+和FORTRAN)编写的程序。尽管MATLAB主要用于数值计算,但是因为大量的额外的工具箱它也适合于不同领域的应用,如控制系统设计与分析、图像处理和信号处理和通信、金融建模和分析等。除了一个完整的Simulink包,提供了一个可视化的开发环境,通常用于系统仿真、动态/嵌入式系统开发等。1.3 Matlab下的simulink简介Simulink是MATLAB最重要的组件之一，它提供一个动态系统建模、仿真和综合分析的集成环境。在该环境中只要通过简单的鼠标操作，就可以构造出复杂的系统。Simulink提供了一个建立模型方块图的图形用户接口，这个创建过程只需单击和拖动鼠标操作就能完成，它提供了一种更快捷、直接明了的方式，而且用户可以立即看到系统的仿真结果。Simulink具有适应面广、结构和流程清晰与仿真精细、效率高、贴近实际、等优点，基于以上优点Simulink已被广泛应用于控制理论和数字信号处理的复杂仿真和设计。同时有大量的第三方软件和硬件应用于Simulink。1.4 通信系统模型通信系统就是传递信息所需要的一切技术设备和传输媒质的总和，包括信息源、发送设备、信道、接收设备和信宿(受信者) ，它的一般模型如图1.4.1所示。图1.4.1 通信系统一般模型模拟通信系统是利用模拟信号来传递消息的通信系统，其模型如图1.4.2所示。图1.4.2 模拟通信系统模型第二章 QPSK调制2.1 QPSK介绍Quadrature Phase Shift Keying通过使用载波的四个各不一样的相位差来表示输入的信息，是具有四进制的相移键控。QPSK是在M=4时的数字的调相技术，它通过约定的四种载波相位，分别为45°，°，225°，275°，输入数据为二进制的数字序列，因为载波相位是四进制的，所有我们需要把二进制的数据变为四进制的，即把二进制序列中每两个比特分成一组，四种排列组合，即00，01，10，11，双比特码元即为一组。每两位二进制信息比特构成每一组，它们分别表示着着四个符号中的某一个符号。2.2 QPSK调制原理QPSK的调制有两种产生方法相乘电路法和选择法。相乘法输入信号是二进制不归零的双极性码元，它通过“串并变换电路变成了两路码元。变成并行码元后，每个码元的持续时间是输入码元的两倍。用两路正交载波去调制并行码元。选择法QPSK的调制中，QPSK信号可以看成是两个载波正交的2PSK信号调制器构成。原理分析如下：根本原理和系统结构QPSK与二进制PSK一样，传输信号包含的信息都存在于相位中。个别的载波相位取四个等间隔值之一，如/4、3/4、5/4、7/4。相应的，可将发射信号定义为：其中，i1，2，3，4；E是发射信号的每个符号的能量，T为符号的持续时间，载波频率f等于nc/T，nc为固定整数。每一个可能的相位值对应于一个特定的二位组。下面介绍QPSK信号的产生和检测。如图为典型的QPSK发射机框图。输入的二进制数据序列首先被不归零(NRZ)电平编码转换器转换为极性形式，即负号1和0分别用和-表示。该二进制波形被分接器分成两个分别由输入序列的奇数位偶数位组成的彼此独立的二进制波形，这两个二进制波形分别用a1(t)和a2(t)表示。此时，在任何一信号时间间隔a1(t)，和a2(t)的幅度恰好分别等于Si1和 Si2，即由发送的二位组决定。这两个二进制波形a1(t)和a2(t)被用来调制一对正交载波：,。这样就得到一对二进制PSK信号。和的正交性使这两个信号可以被独立地检测。最后，将这两个二进制PSK信号相加，从而得期望的QPSK。2.2.2 选择法输入基带信号经过串并变换后用于控制一个相位选择电路，按照当时的输入双比特ab，决定选择哪个相位的载波输出。QPSK调制原理框图图2.3 调制原理框图2.4 QPSK调制方式的Matlab仿真I 路信号是用余弦载波，由2进制数据流的奇数序列组成；Q路信号用正弦载波，由2进制数据流的偶数序列组成。下面的a是Idata，b就是Qdata，它们分布与各自的载波相乘分别输出 I 路信号和 Q 路信号。I 路信号加上Q路信号就是QPSK输出信号。当 I 路载波信号是0相位时为1，是180°相位时为0；当Q路载波信号是0相位时为1，是180°相位时为0。2.4 matlab调制仿真图2.5 QPSK调制方式Matlab-simulink仿真2.5.1 simulink调制建模调制框图1)产生需要的信号源在搭建QPSK调制解调系统中使用伯努力信号发生器产生随机的01比特序列，每两比特代表就一个符号。Bernoulli Binary Generator模块利用伯努利分布的原理，相应得到参数为p的伯努利分布。伯努利分布的均值1 - p和方差p(1 p)的。一个零概率参数指定p。本次实验中的p设置为，即0和1等概。采样时间可根据需要进展设置，例如测误码率时采样时间设为。信号源参数设置2串并变换我们先通过使用buffer 这个模块来实现将信号源信号转变为两路信号。Buffer 模块可以重新分配缓冲区块的输入样本，用到了Demux，可以将一个复合输入转化为多个单一输出，即可以输出多个采样率较低的帧信号。但会产生与缓冲区容量一样的时延。所以，我们可以设置buffer的参数容量为2。图2.5.3 Buffer的参数设置3单极性信号转化为双极性信号因为QPSK的调制信号要求的是双极性信号，所以用伯努利随机生成二进制Generator模块产生的信号必须经过转化才能够被使用。利用加法模块和常数产生模块将1和0的序列各自减去1/2，再利用比例运算模块乘以2，就得到了1和-1 的双极性序列。4调制模块分别将两路信号乘以相位相差 / 2的载波，然后相加。载波由正弦信号发生器产生。正弦波模块的参数设置为可基于时间的模式，时间设为使用仿真时间，我们设载波信号的幅度为1，载波频率可根据需要来进展设置，两路载波同频正交，相位相差 / 2。我们设上支路的相位为0，下支路的相位为 / 2。2.5.2 simulink调制仿真结果图2.5.8 下支路载波第三章 QPSK解调3.1 QPSK解调原理QPSK接收机由一对共输入地相关器组成。这两个相关器分别提供本地产生地相干参考信号和。相关器接收信号x(t),相关器输出地x1和x2被用来与门限值0进展比拟。如果x1>0,如此判决同相信道地输出为符号1；如果x1<0 ,如此判决同相信道的输出为符号0。如果正交通道也是如此判决输出。最后同相信道和正交信道输出这两个二进制数据序列被复加器合并，重新得到原始的二进制序列。在AWGN信道中，判决结果具有最小的负号过失概率。用两路具有相互正交特性的载波来解调信号，可以别离这两路正交的2PSK信号。相干解调后，并行码元经过并/串变换后，最终得到串行的数据流。3.2 QPSK解调原理框图3.3 QPSK解调方式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.3 matlab解调仿真图3.4 QPSK解调方式的Matlab-simulink仿真3.4.1 QPSK解调建模首先将从高斯信道送过来的信号分别乘以与调制时的载波同频的载波，且相位相差为 / 2的载波。解调可以使用相关器或者匹配滤波器进展解调，本次实验使用的是相关器，这时信号需要通过设置的积分器。因为积分器设置为使用积分器时，需要在时间t=T 时使得积分器复位，所以需要设置积分模块续设置在时钟下降沿时复位。并需要设置参数为使用外部信号，此时，时钟设置为与该支路码元时间一样，即是发送信元时间的两倍输入。然后积分后的信号经过采样并保持模块，即sample and hold 模块，此时，设置这个模块为触发上升沿，同样使用时钟设置为与该支路码元时间一样，即是发送信元时间的两倍输入。此时各路传输的信号相应地会使一个单位的时延产生。然后使用autothreshold 模块，该模块根据输入的信号数据自动设置出阀值，由此，可对输入信号做出判定，再输出相应的二进制比特序列，并可输出阀值。最后使用N-sample switch 模块来实现并串转换，因为我们最后需要的是最原始的信号，将两路信号合二为一，在第一路信号发出一个样本时间后，样本时间设置为发送信元时间，开关会自动转换到第二路信号，此时换做第二路信号输入，一个码元时间后模块重置，如此循环。同样的，此模块也需要两倍的发送信元时间输入。积分器设置采样和保持设置传输信道本次实验使用的是高斯信道和理想信道。实验所需的高斯噪声我们可以由高斯信道模块来提供，用到了Zero-Order Hold，和子模块，即SubSystem，通过子模块建立新的封装(Mask)功能模块其中参数设置噪比为Es/No， Es/No为信号能量比噪声功率谱密度。AWGN信道模块可以将加性高斯白噪声加到一个实数的或复数的输入信号。现在输入信号是实数，这个模块增加了实数的高斯噪声，产生一个实数的输出信号。此块继承它的输入信号的采样时间。模块使用信号处理模块随机产生的噪声。初始种子可以是一个标量或矢量的长度相匹配的输入信号通道数。种子的详细资料初次，查看随机源模块库文件参考页面中设置的信号处理。该端口的数据类型都继承自该驱动器的信号块。注意权力的所有值假设一个1欧姆的标称阻抗。 仿真结果信号源和转变后的双极性信号经过高斯信道后的调制信号上支路积分和采样后的信号上支路判决后的信号下支路积分和采样后的信号下支路判决后的信号3.5 仿真结果分析从上述图中可以看出，因为整个系统模块有引进噪声，以与电路使用了积分、采样保持模块，复原后的信号幅度差异较大，通过判决门限后，得到原来的二进制信号。我们从仿真结果图中看出，信道噪声功率谱密度越大，信号信噪比越小，误码率越高，这也符合实际情况。仿真的各种条件都是理想化的，除了噪声之外不会发生任何错误，和实际情况相比，在一样的信噪比之下，比特错误率理应要小的多，但是仿真所得结果的误码率偏大。第四章 QPSK通信系统性能分析1从仿真结果中看出，信号信噪比越小，错误率越高，跟实际情况比拟符合。2由于仿真中各种条件都是理想化的，包括数据在传输过程中，除了噪声影响以外不会发生任何错误，所以相对实际情况来说，在一样信噪比下，比特错误率要小的多，但是仿真所得结果与事实规律并不违背。总结：第五章 结论在搭建QPSK调制解调系统中首先遇到的问题就是串并转换问题，simulink中没有可以直接解决这个问题的模块，经过网上的搜索与查阅相关资料，最终选取buffer模块，完成可串并转换的功能，但是带人了两个单位的延迟。之后便是示波器的显示问题，此次产生的信号采样速率设置为1秒，载波频率1HZ，对仿真过程中的波形分析带来很多的方便。3在完成的QPSK系统的simulink仿真中，由于方便起见，载波大小设为1HZ，与实际的调制解调中所用载波相差甚远，但假如将载波设置变大，会对实验仿真和分析带来极大的不便。参考文献1 苗长云，保锁.现代通信原理与应用(第2版)M：电子工业，2009.2 樊昌信，曹丽娜通信原理(第6版)M国防工业，2009.13J电子设计工程,2011(9).4J.科学技术与工程，2010(5)5J.科技广场,2006(1).附 录1QPSK调制M文件程序clear allclose all % x1是类似1 1 -1 -1 -1 -1 1 1的分布,作用是控制相位的180°反转。%由于仿真中载波的频率是f=1Hz，所以1s的间隔有一个完整周期的正弦波。t=-1:0.01:7-0.01; % t共800个数据，-17st1=0:0.01:8-0.01; %t1也是800个数据点 ，0 8stt=length(t); % tt=800x1=ones(1,800);for i=1:tt if (t(i)>=-1 & t(i)<=1) | (t(i)>=5& t(i)<=7); x1(i)=1; else x1(i)=-1; endendt2 = 0:0.01:7-0.01; %t2是700个数据点，是QPSK_rc绘图的下标t3 = -1:0.01:7.1-0.01; %t3有810个数据点，是i_rc的时间变量t4 = 0:0.01:8.1-0.01; %t4有810个数据点，是q_rc的时间变量tt1=length(t1);x2=ones(1,800); %x2是类似于1 1 -1 -1 1 1 1 1的分布,作用是控制相位的180°反转for i=1:tt1 if (t1(i)>=0 & t1(i)<=2) | (t1(i)>=4& t1(i)<=8); x2(i)=1; else x2(i)=-1; endendf=0:0.1:1;xrc=0.5+0.5*cos(pi*f); %xrc是一个低通特性的传输函数y1=conv(x1,xrc)/5.5; %y1和x1 实际上没什么区别，仅仅是上升沿、下降沿有点过渡带y2=conv(x2,xrc)/5.5; % y2和x2 实际上没什么区别，仅仅是上升沿、下降沿有点过渡带 n0=randn(size(t2); f1=1;i=x1.*cos(2*pi*f1*t); % x1就是I dataq=x2.*sin(2*pi*f1*t1); %x2就是Q dataI=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); % y1就是I data，i_rc可能是贴近实际的波形,i如此是理想波形q_rc=y2.*sin(2*pi*f1*t4); %y2就是Q data，q_rc可能是贴近实际的波形,q如此是理想波形I_rc=i_rc(101:800);Q_rc=q_rc(1:700);QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc);QPSK_rc_n1=QPSK_rc+n1;figure(1)subplot(3,1,1);plot(t3,i_rc);axis(-1 8 -1 1);ylabel('a序列');subplot(3,1,2);plot(t4,q_rc);axis(-1 8 -1 1);ylabel('b序列');subplot(3,1,3);plot(t2,QPSK_rc);axis(-1 8 -1 1);ylabel('合成序列');2QPSK解调M文件程序clear allclose allbit_in = randint(1e3, 1, 0 1);bit_I = bit_in(1:2:1e3); %bit_I为奇数序列，奇数序列是同相分量，以cos为载波bit_Q = bit_in(2:2:1e3); %bit_Q是bit_in的所有偶数下标组成的偶数序列，以sin为载波data_I = -2*bit_I+1; % 将bit_I中的1变成-1，0变成1; 注意data_I是500点data_Q = -2*bit_Q+1; %将bit_Q中的1变成-1，0变成1data_I1=repmat(data_I',20,1); %将500行的列向量data_I的共轭转置data_I复制为20*500的矩阵，20行数据是一样的。data_Q1=repmat(data_Q',20,1);for i=1:1e4 %data_I2是将data_I1这个20*500的矩阵拉长为1*10000的行向量data_I2(i)=data_I1(i); data_Q2(i)=data_Q1(i);end;f=0:0.1:1;xrc=0.5+0.5*cos(pi*f); data_I2_rc=conv(data_I2,xrc)/5.5; % data_I2_rc就是Idatadata_Q2_rc=conv(data_Q2,xrc)/5.5; % data_Q2_rc就是Qdataf1=1;t1=0:0.1:1e3+0.9; % 10010个数据，长度和data_I2_rc以与data_Q2_rc一样n0=rand(size(t1); %n0是1*10010的均匀分布的噪声I_rc=data_I2_rc.*cos(2*pi*f1*t1);Q_rc=data_Q2_rc.*sin(2*pi*f1*t1);QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc);QPSK_rc_n0=QPSK_rc+n0;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:10010);t2=0:0.05:1e3-0.05; %t2有20000个点t3=0:0.1:1e3-0.1; %t3有10000个点data_recover=; % 抽样判决for i=1:20:10000 %data_recover是待判决的20000点data_recover=data_recover I(i:1:i+19) Q(i:1:i+19); end;bit_recover=;for i=1:20:20000 if sum(data_recover(i:i+19)>0 %20点为同一个值，20点数据叠加后与阈值0比拟 data_recover_a(i:i+19)=1; %data_recover_a是并/串转换后的20000点 bit_recover=bit_recover 1; %bit_recover是1000点数据 else data_recover_a(i:i+19)=-1; bit_recover=bit_recover -1; endenderror=0;dd = -2*bit_in+1; % 将bit_in中的1变成-1，0变成1ddd=dd' %ddd是1表示0，-1表示1的原始序列，1000个点ddd1=repmat(ddd,20,1); %ddd1是20*1000的矩阵for i=1:2e4 ddd2(i)=ddd1(i); %将ddd1拉直为1*20000的行向量ddd2endfor i=1:1e3 if bit_recover(i)=ddd(i) error=error+1; endendp=error/1000;figure(1)subplot(4,1,1);plot(t2,ddd2);axis(0 100 -2 2);title('原序列');subplot(4,1,2);plot(t1,I_demo);axis(0 100 -2 2);title('I 支路解调');subplot(4,1,3);plot(t1,Q_demo);axis(0 100 -2 2);title('Q 支路解调');subplot(4,1,4);plot(t2,data_recover_a);axis(0 100 -2 2);title('解调后序列');23 / 24