2PSK数字传输系统设计与仿真

通信系统设计与仿真实践课程设计实验报告课题名称：2PSK数字传输系统设计与仿真专业班级：姓名：学号：起止时间：浙江科技学院信息与电子工程学院目录一、课题内容………………………………………..….……...1二、设计目的………………………………………..….….…..1三、设计要求…………………………………………………..1四、实验条件…………………………………………..…..…..2五、系统设计………………………………………….…...…..21.通信系统的基本原理…………………………………….22.所设计子系统的原理……………………….…...….……4六、详细设计与编码…………………………….…………….51.设计方案……………………………….…….…..……..52.编程工具的选择…………………………………….…..73.编码与测试……………………………………....……..84.编码与调试过程………………………………..….…..135.运行结果及分析……………………………..…….…..14七、设计心得………………………………………..………..21八、参考文献……………………………….………..……….211一、课题内容使用Matlab进行2PSK的调制解调系统设计与仿真，能输出调制前的基带信号、调制后的2PSK信号和叠加噪声后的2PSK信号波形、解调器在接收到信号后解调的各点的信号波形以及眼图和星座图，并对仿真结果进行分析。二、设计目的1、综合应用《Matlab原理及应用》、《信号与系统》、《通信原理》等多门课程知识，使学生建立通信系统的整体概念；2、培养学生系统设计与系统开发的思想；3、培养学生利用软件进行通信仿真的能力。三、设计要求1、每2人一组，组内成员进行各自分工，分别完成不同子系统的详细功能；2、对通信系统有整体的较深入的理解，深入理解自己仿真部分的原理的基础，画出对应的通信子系统的原理框图；3、提出仿真方案；4、完成仿真软件的编制；5、仿真软件的演示；6、提交详细的设计报告。2四、实验条件计算机、Matlab软件五、系统设计1．通信系统的原理通信系统的一般模型图1通信系统的一般模型信息源：消息的生成者或来源；发送设备：将信源输出的信号变为适合信道传输的发射信号，且发送信号包含了原始信号的一切信息；信道：传输信号的通道，可以是有线的，也可以是无线的；噪声源：在信道中传输，噪声是绝不可避免的，噪声又可为加性噪声（线性的噪声）和乘性噪声（非线性的噪声），一般我们只考虑加性噪声；接收设备：从接收信号中提取我们所希望的信号，并将其转换成适合输出传感器的形式；受信者：消息接收者。在通信系统中，按信号参量的取值方式不同可把信号分为两类，3即模拟信号和数字信号，再按照信道中传输信号的特征，来分为模拟通信系统和数字通信系统。下面分别来介绍模拟通信系统与数字通信系统：模拟通信系统模型图2模拟通信系统模型调制器:将原始电信号变换成其频带适合信道传输的信号；解调器:在接收端将信道中传输的信号还原成原始的电信号。数字通信系统模型图3数字通信系统模型信源编码/译码：主要完成A/D变换、数据压缩、加密等三部分功能；信道编码/译码：数据扩张，增加保镖，提高抗干扰能力等；信息源调制器信道解调器受信者噪声源模拟通信系统信道编码器基带脉冲生成器数字调制器信道数字解调器采样判决器信道译码器信源译码器受信者数字通信系统噪声源信源编码器信息源4基带脉冲生成器/抽样判决器：数字基带传输系统的重要部分；数字调制器/解调器：(1)将数字基带信号的频谱搬移到信道的频带之内，以便于在信道中传输；(2)便于对信道进行频分复用，从而有效地利用信道的频率资源。数字解调是数字调制的反过程。模拟通信系统与数字通信系统的比较模拟通信系统和数字通信系统各有自己的一些特点，但相对于模拟通信系统，数字通信系统有许多独特的优点，例如：数字通信系统抗噪能力强、便于加密、易复用也易压缩、易于集成化、便于用计算机处理数字信号、更适合多种数字业务，容易实现多网合一等。相应地，数字通信系统也有自己的缺点：它需要更宽的传输带宽和较复杂的同步系统等。2．所设计子系统的原理本小组所设计的子系统是数字通信系统，内容是研究不同信道条件下不同的调制解调方法。数字调制技术的两种方法：(1)利用模拟调制的方法去实现数字式调制，即把数字调制看成是模拟调制的一个特例，把数字基带信号当做模拟信号的特殊情况处理；(2)利用数字信号的离散取值特点通过开关键控载波，从而实现数字调制。这种方法通常称为键控法，比如对载波的相位进行键控，便可获得相移键控（PSK）基本的调制方式。图42PSK通信系统调制解调模型信道产生随机基带信号调制解调抽样判决码元再生高斯白噪声5数字通信的主要特点：a.抗干扰能力强b.易加密通信c.差错可控d.适应飞速发展的计算机通信的要求六、详细设计与编码1．设计方案本次仿真用2PSK对信号进行调制解调。①2PSK的调制原理如果两个频率相同的载波同时开始振荡，这两个频率同时达到正最大值，同时达到零值，同时达到负最大值，它们应处于同相状态；如果其中一个开始得迟了一点，就可能不相同了。如果一个达到正最大值时，另一个达到负最大值，则称为反相。一般把信号振荡一次（一周）作为360度。如果一个波比另一个波相差半个周期，我们说两个波的相位差180度，也就是反相。当传输数字信号时，1码控制发0度相位，0码控制发180度相位。载波的初始相位就有了移动，也就带上了信息。相移键控是利用载波的相位变化来传递数字信息，而振幅和频率保持不变。在2PSK中，通常用初始相位0和π分别表示二进制“1”和“0”。因此，2PSK信号的时域表达式为（1）式中，n表示第n个符号的绝对相位：)cos(A)(2PSKnctte6（2）因此，上式可以改写为（3）由于两种码元的波形相同，极性相反，故BPSK信号可以表述为一个双极性全占空矩形脉冲序列与一个正弦载波的相乘：（4）式中（5）这里s(t)为双极性全占空(非归零)矩形脉冲序列，g(t)是脉宽为Ts的单个矩形脉冲，而an的统计特性为（6）BPSK信号的调制原理框图如图2所示。与2ASK信号的产生方法相比较，只是对是S(t)的要求不同。在2ASK中S(t)是单极性的，而在BPSK中S(t)是双极性的基带信号。图52PSK信号的调制原理框开关电路S(t)e2psk(t)1800移相coswct0”时发送“”时发送“，1,00nPtPttecc1,cosA,cosA)(2PSK概率为概率为ttsteccos)(2PSKnsnnTtgats)()(PPan1,1,1概率为概率为7②2PSK的解调原理2PSK信号的解调方法是相干解调法。由于PSK信号本身就是利用相位传递信息的，所以在接收端必须利用信号的相位信息来解调信号。下图3中给出了一种2PSK信号相干接收设备的原理框图。图中经过带通滤波的信号在相乘器中与本地载波相乘，然后用低通滤波器滤除高频分量，在进行抽样判决。判决器是按极性来判决的。即正抽样值判为1，负抽样值判为0。图62PSK的相干接收机原理框图③2PSK的仿真框图图7基于MATLAB的2PSK调制解调仿真框图2．编程工具的选择本仿真所用的工具软件是MATLAB.2011b。该软件的功能强大，最擅长矩阵处理，并在系统仿真、数字信号处理、图形图像分析、数理统计、通信及自动控制领域得到广泛应用，同时MATLAB内部有许coswctbe2psk(t)e输出dca带通滤波器相乘器低通滤波器抽样判决器信道产生随机基带信号调制解调抽样判决码元再生高斯白噪声8多与通信有关的函数，这样程序编写方便，也便于观察波形特征。3．编码与测试2PSK主函数main.m:%------------2PSK通信系统仿真主函数----------------clearall;closeall;clc;max=64;g=randint(1,max);%长度为max的随机二进制序列比特率a1=[];b1=[];f=1000;%载波频率1kHzt=0:2*pi/1999:2*pi;%采样频率2kHzforn=1:length(g);ifg(n)==0;A=zeros(1,2000);%每个值2000个点elseg(n)==1;A=ones(1,2000);enda1=[a1A];%s(t),码元宽度2000c=cos(2*pi*f*t);%载波信号b1=[b1c];%与s(t)等长的载波信号,变为矩阵形式endfigure(1);subplot(3,2,1);plot(a1);gridon;axis([02000*length(g)-22]);title('二进制信号序列');9a2=[];b2=[];forn=1:length(g);ifg(n)==0;B=ones(1,2000);%每个值2000个点c=cos(2*pi*f*t);%载波信号elseg(n)==1;B=ones(1,2000);c=cos(2*pi*f*t+pi);%载波信号enda2=[a2B];%s(t),码元宽度2000b2=[b2c];%与s(t)等长的载波信号endtiaoz=a2.*b2;%e(t)调制figure(1);subplot(3,2,2);plot(tiaoz);gridon;axis([02000*length(g)-22]);title('2PSK调制信号');figure(2);subplot(3,2,1);plot(abs(fft(a1)));axis([02000*length(g)0400]);title('原始信号频谱');figure(2);subplot(3,2,2);plot(abs(fft(tiaoz)));axis([02000*length(g)0400]);title('2PSK信号频谱');%------------------带有高斯白噪声的信道--------------------tz=awgn(tiaoz,10);%信号tiaoz中加入白噪声，信噪比为10figure(1);subplot(3,2,3);plot(tz);gridonaxis([02000*length(g)-22]);title('通过高斯白噪声信道后的信号');figure(2);subplot(3,2,3);plot(abs(fft(tz)));axis([02000*length(g)0400]);title('加入白噪声的2PSK信号频谱');10jiet=2*b1.*tz;%同步解调figure(1);subplot(3,2,4);plot(jiet);gridonaxis([02000*length(g)-22]);title('相乘后信号波形')figure(2);subplot(3,2,4);plot(abs(fft(jiet)));axis([02000*length(g)0400]);title('相乘后信号频谱');%----------------------低通滤波器----------------------fp=500;fs=700;rp=3;rs=20;fn=11025;ws=fs/(fn/2);wp=fp/(fn/2);%计算归一化角频率[n,wn]=buttord(wp,ws,rp,rs);%计算阶数和截止频率[b,a]=butter(n,wn);%计算H(z）figure(3);freqz(b,a,1000,11025);subplot(2,1,1);axis([040000-1003])title('LPF幅相频图');jt=filter(b,a,jiet);figure(1);subplot(3,2,5);plot(jt);gridonaxis([02000*length(g)-22]);title('经低通滤波器后信号波形')figure(2);subplot(3,2,5);plot(abs(fft