基于某Matlab地2PSK,2DPSK仿真

1、实用标准文案摘要： Simulink 是 Mathworks 公司推出的基于Matlab 平台的著名仿真环境Simulin 作为一种专业和功能强大且操作简单的仿真工具，目前已被越来越多的工程技术人员所青睐，它搭建积木式的建模仿真方式既简单又直观，而且已经在各个领域得到了广泛的应用。本次课程设计是基于 MATLAB的 2PSK 和 2 DPSK 仿 真 ，通过系统分析，步骤来完成本次设计任务。通过课程设计从理论学习的轨道逐步引向实际应用，把理论上熟悉的定性分析、定量计算逐步和工程估算、实验调整等手段结合起来，掌握工程设计的步骤和方法，了解科学实验的程序和实施方法，为以后毕业设计和从事信息处理技术

2、的实际工作打下基础。关键词：MATLAB 2PSK, 2DPSK 仿真精彩文档实用标准文案目录第 1 章 MATLABU介 1第2章二进制相移键控 12.1 PSK调制原理 22.2 PSKW调原理 32.3 仿真结果及分析 4第3章二进制差分相移键控 63.1 DPSKM 制原理 63.2 DPSKB调原理 63.3 仿真结果及分析 8第3章总结 10附录 11参考文献 20致谢 错误！未定义书签。精彩文档实用标准文案第 1 章 Matlab 简介美国 Mathworks 公司于 1967年推出了矩阵实验室“Matrix Laboratory ”（缩写为Matlab ）这就是Matlab 最

3、早的雏形。开发的最早的目的是帮助学校的老师和学生更好的授课和学习。Matlab 是一种解释性执行语言，具有强大的计算、仿真、绘图等功能。Simulink是MATLA叶的一种可视化仿真工具，也是目前在动态系统的建模和仿真等方面应用最广泛的工具之一。 确切的说，Simulink 是一个用来对动态系统进行建模、仿真和分析的软件包，它支持线性和非线性系统，连续、离散时间模型，或者是两者的混合。系统还可以使多种采样频率的系统，而且系统可以是多进程的。在Simulink 环境中，它为用户提供了方框图进行建模的图形接口，采用这种结构画模型图就如同用手在纸上画模型一样自如、方便，故用户只需进行简单的点击和拖动

4、就能完成建模，并可直接进行系统的仿真，快速的得到仿真结果。但是 Simulink不能脱离MATLAB5独立工 作。Matlab 将高性能的数值计算和可视化集成在一起，并提供了大量的内置函数，从而被广泛地应用于科学计算、控制系统、信息处理等领域的分析、仿真和设计工作，而且利用 Matlab 产品的开放式结构，可以非常容易地对Matlab 的功能进行扩充，从而在不断深化对问题认识的同时，不断完善Matlab 产品以提高产品自身的竞争能力。利用M语言还开发了相应的Matlab专业工具箱函数供用户直接使用。这些工具箱应用的算法是开放的可扩展的，用户不仅可以查看其中的算法，还可以针对一些算法进行修改，甚

5、至允许开发自己的算法扩充工具箱的功能。目前Matlab 产品的工具箱有四十多个，分别涵盖了数据获取、科学计算、控制系统设计与分析、数字信号处理、数字图像处理、金融财务分析以及生物遗传工程等专业领域。第 2 章 二进制相移键控本次设计所用流程图如图2-1 所示。精彩文档实用标准文案s(t)cos ctn(t) cos ct图2-1 2PSK调制解调框图2.1 PSK调制原理在二进制数字调制中，当正弦载波的相位随二进制数字基带信号离散变化时，则产 生二进制移相键控(2PSK)信号。2PSK信号调制有两种方法，即模拟调制法和键控法。通 常用已调信号载波的0°和180 °分别表示二

6、进制数字基带信号的 1和0,模拟调制 法用两个反相的载波信号进行调制。2PSK以载波的相位变化作为参考基准的，当基带信 号为0时相位相对于初始相位为0° ,当基带信号为1时相对于初始相位为180°。键控法，是用载波的相位来携带二进制信息的调制方式。通常用 0°和180。来分 别代表0和1。其时域表达式为：e2PSK = '' ang(t - nT.) cos ct, n其中，2PSK的调制中an必须为双极性码。本次设计中采用模拟调制法。两种方法 原理图分别如图2-2和图2-3所示。双极性图2-2模拟调制法原理图精彩文档实用标准文案心电电美图2-3键

7、控法原理图2.2 PSK解调原理由于2PSK的幅度是恒定的，必须进行相干解调。经过带通滤波的信号在相乘器中 与本地载波相乘，然后用低通滤波器滤除高频分量，在进行抽样判决。判决器是按极性 来判决的。即正抽样值判为1,负抽样值判为00 2PSK信号的相干解调原理图如图2-4 所示，各点白波形如图2-5所示。由于2PSK信号的载波回复过程中存在着180。的相位模糊，即恢复的本地载波与所 需相干载波可能相同，也可能相反，这种相位关系的不确定性将会造成解调出的数字基 带信号与发送的基带信号正好相反，即“ 1”变成“0”吗"0”变成“1”，判决器输出 数字信号全部出错。这种现象称为 2PS3式的

8、“倒冗”现象或“反相工作”。但在本次 仿真中是直接给其同频同相的载波信号，所以不存在此问题。图2-4 2PSK的相干解调原理图精彩文档实用标准文案图2-5相干解调中各点波形图2.3 仿真结果及分析通过编写M文件程序（见附录），产生随机信号，按流程图2-1所示顺序对每一模 块编程后。程序中注有需注意语句及解释。运行程序，实现 2PSK的调制与解调过程。 本次设计采用模拟调制法和相干解调法。仿真后调制过程及解调过程的图形分别如图 2-6和图2-7所示。0-1I!I：1 I1w11111rl00.511-522 533.544.55基带信号反码统205001000150020002500300035

9、004000 45005000图2-6 2PSK模拟调制方法过程仿真图精彩文档实用标准文案图2-7 2PSK相干解调过程仿真图由图可以看出，产生的随机信号为“ 0011001100”，经过反相产生反码，并将原码跟 反码一起合成双极性码，与载波相乘后加入随机噪声。在解调时，与对应的载波相乘经 过低通滤波、抽样判决后，解调出基带信号与原基带信号一致，可以看出实现了本次设 计目的。精彩文档实用标准文案第3章二进制差分相移键控本次设计所用流程图如图3-1所示cos ctn(t) cos ct图3-1 2DPSK调制解调框图3.1 DPSK调制原理二进制差分相移键控。2DPSKT式是用前后相邻码元的载波

10、相对相位变化来表示数 字信息。假设前后相邻码元的载波相位差为 中，可定义一种数字信息与中之间的关系 为：m 0 （数字信息“0”）邛=g（数字信息”1?）邛为前一码元的相位。实现二进制差分相移键控的最常用的方法是：先对二进制数字基带信号进行差分编码，然后对变换出的差分码进行绝对调相即可。2DPS明制原理图如图3-2。图3-2 2DPSK调制原理框图3.2 DPSK解调原理2DPSKJ号解调有相干解调方式和差分相干解调。用差分相干解调这种方法解调时精彩文档实用标准文案不需要恢复本地载波，只要将 DPSK言号精确地延迟一个码元时间问隔，然后与 DPSK言 号相乘，相乘的结果就反映了前后码元的相对相

11、位关系，经低通滤波后直接抽样判决即 可恢复出原始的数字信息，而不需要在进行差分解码。相干解调码变换法及相干解调法的解调原理是，先对2DPSKJ号进行相干解调，恢复出相对码，再通过码反变换器变换为绝对码，从而恢复出发送的二进制数字信息。在解调过程中，若相干载波产生180。相位模糊，解调出的相对码将产生倒置现象， 但是经过码反变换器后，输出的绝对码不会发生任何倒置现象，从而解决了载波相位模 糊的问题。本次设计采用相干解调。两种解调方式的原理图如图3-3和图3-4所示。2DPSKS干解调各点波形图如图3-5所示。b带通 滤波器相乘器低通幄波器抽样判决器定时脉图3-3 2DPSK差分相干解调原理图2D

12、PSK信号带通底波器相乘器低通滤波器.抽样判决器码（反i变换器数据丽i本地截波图3-4 2DPSK相干解调原理图zVVaA/vWVxAAAAAAAy图3-5 2DPSK相干解调各点波形图精彩文档实用标准文案3.3 仿真结果及分析通过编写M文件程序（见附录），产生随机信号，按流程图3-1所示顺序对每一模 块编程后。程序中注有需注意语句及解释。本次设计采用相干解调法。运行程序，实现2DPSK勺调制与解调过程。仿真后调制过程及解调过程的图形分别如图3-6、图3-7和图3-8所示。由图可以看出，产生的绝对随机码为“ 0100111010',经码差分变化产生相对码，经反相产生相对码反码，分别与两

13、个载波相乘后加入随机噪声。在解调时，分 别与对应的载波相乘经过低通滤波、抽样判决后，解调出基带信号与原基带信号一致,可以看出实现了本次设计目的。精彩文档图3-6 2DPSK调制过程仿真图20*200 511.522.533.544 55加口桑声后信号IIILIIIEI00.511.522.533.5图3-7 2DPSK调制过程仿真图4.55实用标准文案0 1111-4-0 OS 11 522.S33.544 5 S距反娈推后波开负0 一 I III I IliditiiJ,J,1上10Q.511 522.533.544 55图3-8 2DPSK相干解调过程仿真图精彩文档实用标准文案第 3 章

14、总结通过在本次设计中的实践明白了自己知识上的误区，例如，在低通滤波的过程中，主要目的是滤去高频分量，滤去载波成分，所以对于低通滤波器的截止频率的设置较为关键。而在2PSK的调制与解调中所用信号为双极性的信号，因此要将在本次设计中产生的单极性信号经过码的变化形成双极性码来传输。本次设计只是按理论上的知识结构进行简单地系统构建，目的是明确数字基带传输的原理及过程，而对于具体问题，例如，实际中信道噪声一般为高斯白噪声，本次设计为简便并未采用而是用rand 函数产生了随机噪声信号。止匕外，还有在2PSK实际传输系统中，在恢复载波的过程中会出现“倒 冗现象”即相位模糊现象，但是在本次设计中直接在解调时给

15、其同频同相的载波，所以 不会出现此种现象，因为不必考虑。这是自己第一次利用Matlab 编程功能实现通信原理中基础知识系统地构建，以前都是利用Matlab 中 Simulink 模块搭建系统来实现，通过自己亲自去动手和调试我明白了实践的重要性，尤其是对程序的调试，更需要大量的时间反复上机运行，发现错误并改正， 这样也就加强了自己对程序分析的能力，更深刻地明白了通信原理中的知识内容，更进一步懂得了理论联系实际的含义，同时提高了自己的思考能力，使得自己对课本里的内容理解、记忆地更加透彻，这无论是在我以后的工作中或是学习中都是非常有用的。精彩文档实用标准文案附录2PS明制解调程序及注释clear a

16、llclose alli=10;j=5000;fc=4;% 载波频率fm=i/5;% 码元速率B=2*fm;t=linspace(0,5,j);a=round(rand(1,i); % 随机序列, 基带信号figure(3);stem(a);st1=t;for n=1:10if a(n)<1;for m=j/i*(n-1)+1:j/i*nst1(m)=0;endelsefor m=j/i*(n-1)+1:j/i*nst1(m)=1;endendendfigure(1);subplot(411);plot(t,st1);title(' 基带信号st1');axis(0,5,

17、-1,2);叫于psa的是双极性信号，因此对上面所求单极性信号取反来与之一起构成双极 性码精彩文档实用标准文案st2=t;for k=1:j;if st1(k)>=1;st2(k)=0;elsest2(k)=1;endend;subplot(412);plot(t,st2);title(' 基带信号反码st2');axis(0,5,-1,2);st3=st1-st2;subplot(413);plot(t,st3);title(' 双极性基带信号st3');axis(0,5,-2,2);s1=sin(2*pi*fc*t);subplot(414);plot

18、(s1);title(' 载波信号s1');e_psk=st3.*s1;figure(2);subplot(511);plot(t,e_psk);title('e_2psk');noise=rand(1,j);加入噪声psk=e_psk+noise; %subplot(512);plot(t,psk);title(' 加噪后波形');与载波相乘psk=psk.*s1;%subplot(513);plot(t,psk);title(' 与载波 s1 相乘后波形');精彩文档实用标准文案f,af = T2F(t,psk); %通过低通

19、滤波器t,psk = lpf(f,af,B);subplot(514);plot(t,psk);title(' 低通滤波后波形');for m=0:i-1;if psk(1,m*500+250)<0;for j=m*500+1:(m+1)*500;psk(1,j)=0;endelsefor j=m*500+1:(m+1)*500;psk(1,j)=1;endendendsubplot(515);plot(t,psk);axis(0,5,-1,2);title(' 抽样判决后波形')2DPSKI制解调程序及注释clear allclose alli=10;j

20、=5000;fc=4;%载波频率fm=i/5;%码元速率B=2*fm;t=linspace(0,5,j);a=round(rand(1,i);figure(4);stem(a);st1=t;for n=1:10精彩文档实用标准文案精彩文档全零矩阵if a(n)<1;for m=j/i*(n-1)+1:j/i*n st1(m)=0;endelsefor m=j/i*(n-1)+1:j/i*n st1(m)=1;endendendfigure(1);subplot(321);plot(t,st1);title(' 绝对码 ');axis(0,5,-1,2);b=zeros(1

21、,i);%b(1)=a(1);for n=2:10if a(n)>=1;if b(n-1)>=1b(n)=0;elseb(n)=1;endelseb(n)=b(n-1);endendst1=t;for n=1:10if b(n)<1;for m=j/i*(n-1)+1:j/i*n st1(m)=0;endelse实用标准文案for m=j/i*(n-1)+1:j/i*n st1(m)=1;endendendsubplot(323);plot(t,st1);title(' 相对码 st1');axis(0,5,-1,2);st2=t;for k=1:j;if s

22、t1(k)>=1;st2(k)=0;elsest2(k)=1;endend;subplot(324);plot(t,st2);title(' 相对码反码st2');axis(0,5,-1,2);s1=sin(2*pi*fc*t);subplot(325);plot(s1);title(' 载波信号s1');s2=sin(2*pi*fc*t+pi);subplot(326);plot(s2);title(' 载波信号s2');d1=st1.*s1;d2=st2.*s2;figure(2);subplot(411);plot(t,d1);tit

23、le('st1*s1');精彩文档实用标准文案精彩文档subplot(412);plot(t,d2);title('st2*s2');e_dpsk=d1+d2;subplot(413);plot(t,e_dpsk);title(' 调制后波形');noise=rand(1,j);dpsk=e_dpsk+noise;%subplot(414);plot(t,dpsk);title(' 加噪声后信号');dpsk=dpsk.*s1;%figure(3);subplot(411);plot(t,dpsk);title(' 与载波

24、相乘后波形');f,af=T2F(t,dpsk);%t,dpsk=lpf(f,af,B);subplot(412);plot(t,dpsk);title(' 低通滤波后波形');st=zeros(1,i);%for m=0:i-1;if dpsk(1,m*500+250)<0;st(m+1)=0;for j=m*500+1:(m+1)*500; dpsk(1,j)=0;endelsefor j=m*500+1:(m+1)*500;st(m+1)=1;dpsk(1,j)=1;endend加入噪声与载波si相乘通过低通滤波器全零矩阵实用标准文案endsubplot(4

25、13);plot(t,dpsk);axis(0,5,-1,2);title(' 抽样判决后波形')dt=zeros(1,i);%全零矩阵dt(1)=st(1);for n=2:10;if (st(n)-st(n-1)<=0&&(st(n)-st(n-1)>-1;dt(n)=0;elsedt(n)=1;endendst=t;for n=1:10if dt(n)<1;for m=j/i*(n-1)+1:j/i*nst(m)=0;endelsefor m=j/i*(n-1)+1:j/i*nst(m)=1;endendendsubplot(414);p

26、lot(t,st);axis(0,5,-1,2);title(' 码反变换后波形');冽用FFTU算信号的频谱并与信号的真实频谱的抽样比较。%却本文件T2F.m定义了函数T2F,计算信号的傅立叶变换。function f,sf= T2F(t,st)%This is a function using the FFT function to calculate a signal's Fourier精彩文档实用标准文案%Translation%Input is the time and the signal vectors,the length of time must greater%than 2%Output is the frequency and the signal spectrumdt = t(2)-t(1);T=t(end);df = 1/T;N = length(st);f=-N/2*df:df:N/2*df-df;sf = fft(st);sf = T/N*fftshift(s ) ；用到的低通滤波器函数function t,st=lpf(f,sf,B)%This function filter an input da