数字通信系统课程设计

1、三峡大学数字通信系统课程设计基于matlab的2psk调制与解调系统学 院： 计算机与信息学院 学 号： 2010156203 班级序号： 72 姓 名： 蔡红亮 指导老师： 夏平 日 期： 2013年1月2日 数字通信系统课程设计任务书题 目基于matlab的2psk调制与解调系统设 计目 的1、综合应用信号与系统、随机信号分析、通信原理等课程知识，使学生建立通信系统的整体概念；2、培养学生系统设计与系统开发的思想；3、培养学生利用matlab软件平台进行通信系统仿真的能力。设 计要 求1、每人一组，完成指定设计问题的系统详细功能；2、对通信系统有整体的较深入的理解，深入理解设计系统的基本原

2、理，画出组成的每个通信子系统的原理框图；3、提出设计方案；4、完成仿真软件的编制；5、在matlab仿真平台上，实现仿真系统正常运行；6、提交详细的设计报告。主 要设 备计算机1台，matlab仿真软件一套。主要参考文献1 张辉，曹丽娜等. 现代通信原理与技术（第三版）.西安电子科技大学出版社，2009.2 周炯槃等. 通信原理(第三版). 北京邮电大学出版社，2005. 3 樊昌信. 通信原理（第五版）.国防工业出版社，2006.4 曹志刚等.现代通信原理. 清华大学出版社，1992.5 leon w. couch. digital and analog communication syst

3、ems(影印第5版). 清华大学出版社，1999 .6 stephen g. wilson. digital modulation and coding（影印版）.电子工业出版社，1998.课程设计进度安排（起止时间、工作内容）. .目 录一、课程设计主要内容二、系统设计的原理三、课程设计思路及过程四、课程设计实验结果五、课程设计结果分析六、参考文献 七、附录1、 课程设计的主要内容了解并掌握数字频带传输系统2psk的调制与解调系统，通过matlab来实现整个通信系统的设计，进一步加深对2psk系统的理解，并熟练掌握matlab语言。2、 课程设计原理数字信号的传输方式分为基带传输和带通传输，

4、在实际应用中，大多数信道具有带通特性而不能直接传输基带信号。为了使数字信号在带通信道中传输，必须使用数字基带信号对载波进行调制，以使信号与信道的特性相匹配。这种用数字基带信号控制载波，把数字基带信号变换为数字带通信号的过程称为数字调制。数字调制技术的两种方法：利用模拟调制的方法去实现数字式调制，即把数字调制看成是模拟调制的一个特例，把数字基带信号当作模拟信号的特殊情况处理；利用数字信号的离散取值特点通过开关键控载波，从而实现数字调制。这种方法通常称为键控法，比如对载波的相位进行键控，便可获得相移键控（psk）基本的调制方式。 乘法器)(2tspsk双极性不归零t0wcos)(tf码型变换单/双

5、 图-1模拟调制的方法tcwcos)(tf)(2tspsk开关电路移相01800p 图-2相移键控方法调制过程：数字调相：如果两个频率相同的载波同时开始振荡，这两个频率同时达到正最大值，同时达到零值，同时达到负最大值，它们应处于同相状态；如果其中一个开始得迟了一点，就可能不相同了。如果一个达到正最大值时，另一个达到负最大值，则称为反相。一般把信号振荡一次（一周）作为360度。如果一个波比另一个波相差半个周期，我们说两个波的相位差180度，也就是反相。当传输数字信号时，1码控制发0度相位，0码控制发180度相位。载波的初始相位就有了移动，也就带上了信息。相移键控：利用载波的相位变化来传递数字信息

6、，而振幅和频率保持不变。在2psk中，通常用初始相位0和分别表示二进制“1”和“0”。因此，2psk信号的时域表达式为 (t)=acos其中，表示第n个符号的绝对相位：= 因此，上式可以改写为 当码元宽度和载波周期相等时，2psk调制波形如下： 图-3 2psk调制波形解调过程2psk信号的解调方法是相干解调法。由于psk信号本身就是利用相位传递信息的，所以在接收端必须利用信号的相位信息来解调信号。下图2-3中给出了一种2psk信号相干接收设备的原理框图。图中经过带通滤波的信号在相乘器中与本地载波相乘，然后用低通滤波器滤除高频分量，在进行抽样判决。判决器是按极性来判决的。即正抽样值判为1，负抽

7、样值判为0.图-4 2psk相干解调系统框图0101sttabcd0tttte00110 图-5 2psk信号相干解调各点时间波形可见，2psk信号相干解调的过程实际上是输入已调信号与本地载波信号进行极性比较的过程，故常称为极性比较法解调。 由于2psk信号实际上是以一个固定初相的末调载波为参考的，因此，解调时必须有与此同频同相的同步载波。如果同步载波的相位发生变化，如0相位变为相位或相位变为0相位，则恢复的数字信息就会发生“0”变“1”或“1”变“0”，从而造成错误的恢复。这种因为本地参考载波倒相，而在接收端发生错误恢复的现象称为“倒”现象或“反向工作”现象。绝对移相的主要缺点是容易产生相位

8、模糊，造成反向工作。这也是它实际应用较少的主要原因。 三、课程设计思路及过程本课题主要是通过matlab来实现整个2psk系统，在用matlab语言编程的时候主要分为以下几个模块。模块一、通过码型变换产生双极性数字基带信号，首先调用matlab函数产生一组随机序列，然后将离散的随机序列变换成单极性的二进制数字基带信号，将单极性信号取反并与上面信号相减即可得到我们所需要的二进制双极性数字基带信号，生成的二进制信号带宽为1。模块二、载波信号的产生，采用的载频为10hz，所以基带信号周期是载波周期的10倍，载波振幅为1。模块三、通过模拟调制的方法得到2psk信号：（ 注：在这设置载波信号的初相为0）

9、模块四、已调信号通过信道传输，通过matlab函数产生一个随机信号来模拟高斯白噪声，并与已调信号相加，在这里我们只考虑加性噪声，不考虑乘型噪声。模块五、相干解调过程（注：这里没有载波提取过程，而是直接用本地载波去与已调信号相乘）模块六、用matlab来设计理想的巴特沃斯低通滤波器模块七、根据位同步信号进行抽样判决模块八、2psk信号的频谱分析模块九、采用窗函数方法对信号进行功率谱密度分析4、 课程设计的试验结果1、2psk调制、解调系统信号波形图图-6二进制数字基带信号与载波信号 图-7 2psk信号与接收机接收到的信号波形图-8通过解调、低通输出波形图-9判决输出波形将抽样输出的基带信号波形

10、与原信号相比较可以发现两个信号波形完全相同，说明整个2psk的调制、解调系统是正确的。五、课程设计结果分析1、2psk调制解调系统频谱分析数字基带信号为二进制非归零码，所以它的傅里叶变换应该是门函数，信号带宽为1/ts，其频谱图-9。载波信号为正弦信号，载频为10hz，其频谱如下。图-9 数字基带信号与载波信号频谱2psk信号的实质为基带信号与载波相乘的结果，所以2psk的信号频谱是基带信号频谱的线性搬移，中心频率为fc.，带宽为基带信号带宽的两倍，如图-10，2psk信号通过信道后，加上了高斯白噪声，其频谱如下。图-10 2psk信号与接收机接收的信号频谱相干解调实质是将已调信号与提取的载波

11、信号相乘，所以解调输出的信号的频率成分包括基带信号频率和2倍的载频信号和噪声如图-11，通过低通滤波器后，高频成分被虑掉，剩下基带信号和部分噪声。图-11 解调输出与低通输出的信号频谱图-12 判决输出信号频谱与基带信号频谱通过比较判决输出信号的频谱和基带信号频谱，得出结论，2psk的调制解调系统设计是正确的。2、 各阶段信号功率谱密度分析 图-12 2psk功率谱密度3、“倒”现象分析 在相干解调之前，需要根据已调信号来提取载波信号，如果同步载波的相位发生变化，则恢复的数字信息就会发生“0”变“1”或“1”变“0”，从而造成错误的恢复。这种因为本地参考载波倒相。所以在实际应用中多采用2dps

12、k调制解调系统，可以防止“倒”现象的发生。在matlab程序中，只用将提取的载波加上90度的相位偏移即可看到“倒”现象。 六、参考文献1 张辉，曹丽娜等. 现代通信原理与技术（第三版）.西安电子科技大学出版社，2009.2 周炯槃等. 通信原理(第三版). 北京邮电大学出版社，2005. 3 樊昌信. 通信原理（第五版）.国防工业出版社，2006.4 曹志刚等.现代通信原理. 清华大学出版社，1992.5 leon w. couch. digital and analog communication systems(影印第5版). 清华大学出版社，1999 .6 stephen g. wils

13、on. digital modulation and coding（影印版）.电子工业出版社，1998. 七、附录 程序代码实现：clear allclose alli=20;j=5000;t=linspace(0,20,5000);%在0到20之间等分成5000个点% % 功能实现模块一：产生随机二进制数字基带信号 %rem=round(rand(1,i);%调用函数产生随机序列signal1=t;%将离散的的随机序列转变成连续的单极性二进制码元for n=1:20 if rem(n)=1; signal2(k)=0; else signal2(k)=1; endend;signal3=si

14、gnal1-signal2;%将两个单极性数字基带信号相减即可获得双极性二进制数字基带信号subplot(312);plot(250*t,signal3);title(双极性二进制数字基带信号);axis(0,5000,-2,2);grid on;% 功能实现模块二：产生载波信号 % 载波周期与码元宽度相同且均为1 %fc=10;%载波频率为10hzs=sin(2*pi*fc*t); subplot(313);plot(s);title(载波信号波形);axis(0 5000 -1 1);grid on;% 功能实现模块三：2psk调制部分:基带信号与载波相乘 %psk=signal3.*s;

15、figure(2);subplot(211);plot(t*250,psk);title(已调信号波形);axis(0,5000,-2,2);% 功能实现模块四：调制信号通过信道传输产生加性噪声 %noise=rand(1,j);y1=psk+noise;%模拟加性噪声subplot(212);plot(y1);axis(0,5000,-1,2);title(加噪后的波形);% 功能实现模块五：相干解调实现过程 %y2=y1.*s;%接收机收到信号输入乘法器与载波信号相乘figure(3)subplot(211);plot(y2);title(解调输出波形);grid;% 功能实现模块六：通过

16、巴特沃斯低通滤波器 %fp=100;fs=200;rp=3;rs=20;fn=10000;%设计巴特沃斯低通滤波器的参数 ws=fs/(fn/2); wp=fp/(fn/2);%计算归一化角频率n,wn=buttord(wp,ws,rp,rs);%计算阶数和截止频率b,a=butter(n,wn);y3=2*filter(b,a,y2);subplot(212);plot(y3);title(经低通滤波器后的输出信号波形);grid on;% 功能实现模块七：根据位同步信号进行抽样判决 %y4=linspace(0,19);y4=0;for n=1:20; if y3(n*250)0; y4(

17、n)=1; else y4(n)=-1; endendfigure(4)subplot(211);stem(y4);axis(0 20 -2 2);title(抽样后得到的离散采样值);grid on;xlabel(时间/s);y5=t;for n=1:20 if y4(n)1; for m=j/i*(n-1)+1:j/i*n y5(m)=-1; end else for m=j/i*(n-1)+1:j/i*n y5(m)=1; end endendsubplot(212);plot(y5);axis(0 5000 -2 2);grid on;title(抽样判决输出基带信号);xlabel(

18、采样点数);% 功能实现模块八：各阶段信号频谱分析 %fs=250;n=1024; %采样频率和数据点数n=0:n-1;t=n/fs; %时间序列y=fft(signal3,n);%对信号进行快速fourier变换mag=abs(y);%求取fourier变换的振幅f=n*fs/n;figure(5);subplot(211)plot(f(1:n/2),mag(1:n/2);axis(0 50 0 500);xlabel(频率/hz);ylabel(振幅);title(基带信号频谱（n=1024）);grid on; fs=250;n=1024;n=0:n-1;t=n/fs;y=fft(s,n

19、);%对信号进行快速fourier变换mag=abs(y);%求取fourier变换的振幅f=n*fs/n;subplot(212)plot(f(1:n/2),mag(1:n/2);axis(0 50 0 500);xlabel(频率/hz);ylabel(振幅);title(载波信号频谱（n=1024);grid on; fs=250;n=1024;n=0:n-1;t=n/fs;y=fft(psk,n);%对信号进行快速fourier变换mag=abs(y);%求取fourier变换的振幅f=n*fs/n;figure(6)subplot(211)plot(f(1:n/2),mag(1:n/

20、2);axis(0 50 0 500);xlabel(频率/hz);ylabel(振幅);title(2psk信号频谱（n=1024);grid on; fs=250;n=1024;n=0:n-1;t=n/fs;y=fft(y1,n);%对信号进行快速fourier变换mag=abs(y);%求取fourier变换的振幅f=n*fs/n;subplot(212)plot(f(1:n/2),mag(1:n/2);axis(0 50 0 500);xlabel(频率/hz);ylabel(振幅);title(接收机接收到的信号频谱（n=1024);grid on; fs=250;n=1024;n=

21、0:n-1;t=n/fs;y=fft(y2,n);%对信号进行快速fourier变换mag=abs(y);%求取fourier变换的振幅f=n*fs/n;figure(7)subplot(211)plot(f(1:n/2),mag(1:n/2);axis(0 50 0 500);xlabel(频率/hz);ylabel(振幅);title(解调输出信号频谱（n=1024);grid on; fs=250;n=1024;n=0:n-1;t=n/fs;y=fft(y3,n);%对信号进行快速fourier变换mag=abs(y);%求取fourier变换的振幅f=n*fs/n;subplot(21

22、2)plot(f(1:n/2),mag(1:n/2);axis(0 50 0 500);xlabel(频率/hz);ylabel(振幅);title(低通输出信号频谱（n=1024);grid on; fs=250;n=1024;n=0:n-1;t=n/fs;y=fft(y5,n);%对信号进行快速fourier变换mag=abs(y);%求取fourier变换的振幅f=n*fs/n;figure(8)subplot(211)plot(f(1:n/2),mag(1:n/2);axis(0 50 0 500);xlabel(频率/hz);ylabel(振幅);title(抽样判决输出信号频谱（n

23、=1024);grid on; fs=250;n=1024; %采样频率和数据点数n=0:n-1;t=n/fs; %时间序列y=fft(signal3,n);%对信号进行快速fourier变换mag=abs(y);%求取fourier变换的振幅f=n*fs/n;subplot(212)plot(f(1:n/2),mag(1:n/2);axis(0 50 0 500);xlabel(频率/hz);ylabel(振幅);title(基带信号频谱（n=1024）);grid on;% 功能实现模块九：各阶段信号功率谱密度分析 %fs=1000; n=0:1/fs:1; nfft=1024; window=hamming(100); %海明窗noverlap=20; %数据无重叠 range=half;