ASK调制解调通信系统

1、信号与通信系统课程设计说明书题目：设计ASK调制解调通信系统系部：信息与控制工程学院专 业： 电子信息工程班级：XXXX级X班学生姓名：XX学号：XXXXXXXXXX指导教师：XXX2018年 6月12日信号与通信系统课程设计说明书目 录 TOC o 1-5 h z HYPERLINK l bookmark4 o Current Document 1设计任务与要求2 HYPERLINK l bookmark6 o Current Document 设计任务2 HYPERLINK l bookmark8 o Current Document 设计要求2 HYPERLINK l bookmark1

2、0 o Current Document 2设计方法与内容3 HYPERLINK l bookmark12 o Current Document MATLAB 简介3 HYPERLINK l bookmark14 o Current Document ASK信号调制原理3 HYPERLINK l bookmark16 o Current Document ASK解调原理4 HYPERLINK l bookmark18 o Current Document 3仿真实现过程5 HYPERLINK l bookmark20 o Current Document ASK信号的产生5 HYPERLINK

3、 l bookmark22 o Current Document 载波信号波形5 HYPERLINK l bookmark24 o Current Document ASK调制解调实现6 HYPERLINK l bookmark26 o Current Document 叠加噪声的ASK调制解调7 HYPERLINK l bookmark28 o Current Document 4结论10 HYPERLINK l bookmark30 o Current Document 5附录11 HYPERLINK l bookmark32 o Current Document 参考文献18信号与通信系

4、统课程设计说明书1设计任务与要求设计任务.根据题目查阅有关资料，掌握数字带通调制技术。.学习MATLAB软件，掌握MATLAB各种函数的使用。.据数字带通调制原理，运行MATLAB进行编辑，仿真调制过程，记录并分析仿真 结果。设计要求.掌握ASK调制解调原理2.绘制出ASK信号解调前后在时域和频域中的波形，观察解调前后频谱的变化理 解ASK信号解调原理。信号与通信系统课程设计说明书2设计方法与内容MATLAB 简介Matlab是MathWorks公司于1982年推出的一套高性能的数值计算和可视化 软件。它集数值分析、矩阵运算、信号处理和图形显示于一体，构成了一个方便、 界面良好的用户环境。它还

5、包括了 Toolbox(工具箱)的各类问题的求解工具，可 用来求解特定学科的问题。其特点是：(1)可扩展性：Matlab最重要的特点是易于扩展，它允许用户自行建立指定 功能的M文件。对于一个从事特定领域的工程师来说，不仅可利用Matlab所提供 的函数及基本工具箱函数，还可方便地构造出专用的函数。从而大大扩展了其应 用范围。当前支持Matlab的商用Toolbox(工具箱)有数百种之多。而由个人开发 的Toolbox则不可计数。(2)易学易用性：Matlab不需要用户有高深的数学知识和程序设计能力，不 需要用户深刻了解算法及编程技巧。(3)高效性：Matlab语句功能十分强大，一条语句可完成十

6、分复杂的任务。 如fft语句可完成对指定数据的快速傅里叶变换，这相当于上百条C语言语句的 功能。它大大加快了工程技术人员从事软件开发的效率。ASK信号调制原理数字信号对载波信号的振幅调制称为振幅键控，即 ASK(Amplitude Shift Keying)。2ASK就是调制信号为二进制数字基带信号时的振幅键控。简单的说，振幅键控就是利用载波的幅度变化来传递数字信息，而其频率和初始 相位保持不变。在2ASK中，载波的幅度只有两种变化状态，分别对应二进制信息的 “0” 或 “1”。2ASK已调信号可表示为e 0 = s(t) cos 3 ct式中，3c为载波角频率，s(t)为单极性NRZ矩形脉冲

7、序列s(t) = San g(t-n Ts)其中，g(t)是持续时间为Ts、高度为1的矩形脉冲，口为二进制数字(1,出现概率为p=i = (0，出现概率为1-p信号与通信系统课程设计说明书ASK解调原理2ASK/OOK信号有两种基本的解调方法：非相干解调（包络检波法）和相干解调 （同步检测法）。本课程设计要求的是相干解调，如图2-1:e2 ASK图2-1相干解调信号与通信系统课程设计说明书3仿真实现过程ASK信号的产生图3-1二进制基带信号时域谱和频域谱先将源程序创建M文件，自定义路径;编写语句x=ceil(rand(1,100000)-0.5)生 成一段随机的二进制基带信号，其中rand产生

8、随机矩阵，ceil为取整函数；再编写 语句FFT1=fft(x,128); FFT1=abs(FFT1)对随机序列进行傅里叶变换并取绝对值。其 时域谱和频域谱如图3-1。载波信号波形图3-2载波信号时域谱和频域谱信号与通信系统课程设计说明书在调制解调系统中，载波信号的频率一般要大于信号源的频率。信号源频率为 12 Hz，所以将载波频率设置为36 Hz，编写正弦函数carry=cos(2*pi*Fc*t)并进行 频域转换FFT2=fft(carry,256); FFT2=abs(FFT2)。其时域谱和频域谱如图3-2。ASK调制解调实现编写数字带通调制函数y=dmod(x,Fc,Fd,Fs,as

9、k,2)进行2ASK调制，其中x为 输入信号，Fc为载波频率，Fd为码速率，Fs为采样频率，Fs/Fd必须是一个正整数。 再进行频域转换FFT3=fft(y,256); FFT3=abs(FFT3)，其时域谱和频域谱如图3-3：图3-3已调信号时域谱和频域谱图3-4解调信号解调时域谱和频域谱编写数字带通调制函数z=ddemod(y,Fc,Fd,Fs,ask,2)进行2ASK解调，其中x信号与通信系统课程设计说明书为输入信号，Fc为载波频率，Fd为码速率，Fs为采样频率。再进行频域转换 FFT4=fft(z,64);FFT4=abs(FFT4),由于调制到解调的信道中没有加入噪声，所以调 制前信

10、号(即原始信号)与调制后信号相同，其时域谱和频域谱如图3-4。叠加噪声的ASK调制解调图3-5叠加小噪声调制信号时域谱和频域谱图3-6叠加大噪声调制信号时域谱和频域谱信号与通信系统课程设计说明书由于信道中的噪声是叠加在信号上的，噪声始终是存在的，通常称它为加性噪声 或加性干扰。对已调信号y分别叠加高斯小噪声(信噪比为6)和大噪声(信噪比为 -2)则编写Ynt1=awgn(y,6);Ynt2=awgn(y,-2)并进行频域转换，其时域谱和频域谱 分别如图3-5和图3-6。图3-7叠加小噪声解调信号时域谱和频域谱图3-8叠加大噪声解调信号时域谱和频域谱对小噪声Ynt1和大噪声信号Ynt2分别解调，

11、编写数字带通调制函数 z1=ddemod(Ynt1,Fc,Fd,Fs,ask,2);z2=ddemod(Ynt2,Fc,Fd,Fs,ask,2);并进行信号与通信系统课程设计说明书频域转换，其时域谱和频域谱如分别如图3-7和图3-8。T 下 wtr日 P H | S | X |LElI e ZE.di t eft In= ert Tool 3 Luxk七口尸 ndciw HLp 口8 UI 与 医应 C 1?* S H O信噪1=匕与谓码率的美军-10 -e -S -4-20246 B 10f吉噪 bL SNR(r/dB)图3-9误差数和误差率误码率是衡量一个数字通信系统性能的重要指标。在信道

12、高斯白噪声的干扰下， 二进制2ASK数字调制误码率取决于解调器输入信噪比，编写误码率函数 br,Pe1=symerr(x,z1);br,Pe2=symerr(x,z2)，其中 br 为符号误差数，Pe1 为符 号误差率。绘制如图3-9。调用函数semilogy(SNR,Pe);绘制信噪比与误码率的关系曲线如图3-10，由此可 得出：与无噪声时(误码率为零)相比较，当信噪比较大时，噪声小误码率低；反之， 信噪比较小时噪声大误码率高。信号与通信系统课程设计说明书4结论通过此次的课程设计，使我收获了很多。不仅让我对专业方面的知识认识更加深 刻，同时也掌握了学习方法，懂得了无论做什么事，都要先弄懂原理

13、，从根本出发， 将一切问题分块分部解决的方法。让我今后都受益颇深。在本次课程设计的过程中，也遇到了很多问题和难点。比如说当我做到解调部分 的时候，程序都写完了，但就是波形出不来。不能够正确的仿真出波形和原信号。这 让我很头疼，调了好久也没有输出。最后，通过上网搜索资料，并且与同学们一起讨 论研究。终于一点点发现了问题，最后调出了波形。虽然过程很坎坷，但结果却令人欣慰。我觉得在这次的课程设计中，收获的不仅 是我们自己的方面，更多的也有团队的合作和与别人的交流和沟通。这让我们能够更 好的与别人分享和交流，为我们在今后的学习和工作道路上都有所帮助。10信号与通信系统课程设计说明书5附录%程序名称：A

14、SK.m%程序功能：ASK解调解调程序代码clc;%清除命令工作窗里的内容clear;%清除内存空间变量 %产生二进制随机序列x=ceil(rand(1,100000)-0.5) %产生二进制随机序列并取大于x的最小整数figure(1)%窗口 1，包含时域谱和频域谱subplot(2,1,1)stairs(x);xlabelC时间 t);figure(1)%窗口 1，包含时域谱和频域谱subplot(2,1,1)stairs(x);xlabelC时间 t);ylabel(序列值);%分块图函数subplot，图形窗口分成2块子窗口的第1个图像%第2个图像%x轴标注%y轴标注title(二进制

15、随机序列);添加图像标题axis(1 21 -1 2)%控制坐标轴的范围grid on%图像中添加栅格grid on%对随机序列进行频谱分析FFT1=fft(x,128);%对随机序列进行傅里叶变换FFT1=abs(FFT1);FFT1=fft(x,128);%对随机序列进行傅里叶变换FFT1=abs(FFT1);%对傅里叶变换取绝对值figure(1)subplot(2,1,2)%subplot(2,1,2)%第2个图像plot(FFT1);xlabelC频率 f);ylabel(幅度 FFT1);title(随机序列频谱);axis(0 128 0 50) grid on%载波信号t=1/

16、360:1/360:20;%t=1/360:1/360:20;%载波时间范围Fc=36;%Fc=36;%载波频率carry=cos(2*pi*Fc*t);%正弦载波信号carry=cos(2*pi*Fc*t);%正弦载波信号11信号与通信系统课程设计说明书figure(2)%窗口 2,包含时域谱和频域谱subplot(2,1,1)plot(carry);xlabel(时间 t);ylabel(幅度 carry);title(载波信号)；axis(1 600 -2 2) grid on%对载波信号进行频谱分析FFT2=fft(carry,256);%对载波信号进行傅里叶变换FFT2=abs(FF

17、T2);%对傅里叶变换取绝对值figure(2)subplot(2,1,2)plot(FFT2);xlabelC频率 f);ylabelC幅度 FFT2);title(载波信号频谱)；axis(0 256 0 100)grid on%ASK的调制Fd=12;%Fd为码速率，Fs为采样频率Fs=360;y=dmod(x,Fc,Fd,Fs,ask,2);%调用数字带通调制函数dmod进行2ASK调制 for i=1:20if x(i)=0yy(30*(i-1)+1:30*i)=0;elseyy(30*(i-1)+1:30*i)=y(30*(i-1)+1:30*i);endend%对20个随机码元进

18、行判别，若码元为0则该码元周期内调制信号为零figure(3)subplot(2,1,1)plot(yy);12信号与通信系统课程设计说明书xlabelC时间 t);丫匕51(幅度 y);title(已调信号)；axis(1 600 -2 2) grid on%对已调信号进行频谱分析FFT3=fft(y,256);%对已调信号进行傅里叶变换FFT3=abs(FFT3);%对傅里叶变换取绝对值figure(3)subplot(2,1,2)plot(FFT3);xlabel(频率 f);ylabel(幅度 FFT3);title(已调信号频谱)；axis(0 256 0 50)grid on%AS

19、K的解调z=ddemod(y,Fc,Fd,Fs,ask,2); %调用数字带通调制函数dmod进行2ASK解调%对傅里叶变换取绝对值subplot(2,1,1)stairs(z);xlabel(时间 t);y1abe1(幅度 z);title(解调信号)；axis(1 21 -1 2) grid on%对解调信号进行频谱分析FFT4=fft(z,64);%对解调信号进行傅里叶变换FFT4=abs(FFT4);%对傅里叶变换取绝对值figure(4)subplot(2,1,2)plot(FFT4);13信号与通信系统课程设计说明书xlabelC频率 f);ylabel(幅度 FFT4);titl

20、e(解调信号频谱)；axis(0 64 0 50)grid on%加入高斯小噪声,SNR为6Ynt1=awgn(y,6);%加入高斯小噪声，信噪比为6subplot(2,1,1)plot(Ynt1);xlabelC时间 t);ylabelC幅度 Ynt1);title(加小噪声信号)；axis(1 600 -2 2)grid on%对加小噪声信号进行频谱分析FFT5=fft(Ynt1,256);%对加入小噪声的调制信号进行傅里叶变换FFT5=abs(FFT5);%对傅里叶变换取绝对值figure(5)subplot(2,1,2)plot(FFT5);xlabelC频率 f);ylabelC幅度

21、 FFT5);title(加小噪声信号频谱)axis(0 256 0 50)grid on%ASK加小噪声信号的解调及误码率z1=ddemod(Ynt1,Fc,Fd,Fs,ask,2);%调用数字带通调制函数dmod对加小噪声信号进行解调br,Pe1=symerr(x,z1)%对解调后加小噪声信号误码分析，br为符号误差数，Pe1为符号误差率figure(6)subplot(2,1,1)14信号与通信系统课程设计说明书stairs(zl);xlabelC时间 t);ylabel(幅度 zl);title(加小噪声解调信号)；axis(1 21 -1 2)grid on%对加小噪声解调信号进行频

22、谱分析FFT6=fft(z1,64);%对加入小噪声的解调信号进行傅里叶变换FFT6=abs(FFT6);%对傅里叶变换取绝对值subplot(2,1,2)plot(FFT6);xlabel(频率 f);ylabel(幅度 FFT6);title(加小噪声解调信号频谱)；axis(0 64 0 50)grid on%加入高斯大噪声,SNR为-2Ynt2=awgn(y,3);%加入高斯大噪声，信噪比为-2subplot(2,1,1)plot(Ynt2);xlabel(时间 t);ylabel(幅度 Ynt2);title(加大噪声信号)；axis(1 600 -2 2)grid on%对加大噪声

23、信号进行频谱分析FFT7=fft(Ynt2,256);%对加入大噪声的调制信号进行傅里叶变换FFT7=abs(FFT7);%对傅里叶变换取绝对值figure(7)subplot(2,1,2)plot(FFT7);15信号与通信系统课程设计说明书xlabelC频率 f);ylabel(幅度 FFT5);title(加大噪声信号频谱)axis(0 256 0 50) grid on%ASK加大噪声信号的解调及误码率z2=ddemod(Ynt2,Fc,Fd,Fs,ask,2);%调用数字带通调制函数dmod对加大噪声信号进行解调br,Pe2=symerr(x,z2)%对解调后加大噪声信号误码分析，br为符号误差数，Pe1为符号误差率subplot(2,1,1)stairs(z2);xlabelC时间 t);ylabelC幅度 z2);title(加大噪声解调信号)；axis(1 21 -1 2) grid on%对加大噪声解调信号进行频谱分析