完整word版,ASK调制解调通信系统

1、信口号与通信系统课程设计说明书题目 :设计ASK调制解调通信系统系部：信息与控制工程学院专业：电子信息工程班级：XXXX级X班学生姓名：XXX 学号:XXXXXXXXXX指-导教师：XXX2018年 6月12日信号与通信系统课程设计说明书参考文献222333455567101118设计任务与要求1.1 设计任务1.2 设计要求设计方法与内容 2.1 MATLAB 简介2.2 ASK信号调制原理 2.3 ASK解调原理仿真实现过程3.1 ASK信号的产生3.2 载波信号波形3.3 ASK调制解调实现 3.4 叠加噪声的ASK调制解调结论附录11设计任务与要求1.1设计任务1. 根据题目查阅有关资

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

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

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

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

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

7、y=dmod(x,Fc,Fd,Fs,'ask',2)进行2ASK调制，输入信号，Fc为载波频率，Fd为码速率，Fs为采样频率，Fs/Fd必须是一个正整数。再进行频域转换FFT3=fft(y,256); FFT3=abs(FFT3)，其时域谱和频域谱如图3-3 :口g佶幵规比Xl ± w.J止仁乂 1丄1.4 上_丄. Xi V EJh. Lop RCJtlwAlJoqqq电 m V/ y 口SO30 公U1QO-so100150闻1r200图3-3已调信号时域谱和频域谱盒3 1.呼X. L 宣_4 g 1 13 一 Il 业 gp SftJD m u 耳 4 I 込回

8、<回 u I U 叵3 口IhK晦01111 11t1112-45BID121 斗IS1 aWOMJ O t娇调f吉予Hfiis!1 1111L111111c11"' "L 1r- -i111._ - -r-T-|-11111111111i111111七二二SS11i_. _i_ J _ _L.宀广七A' =厂-八.亠-=1070fiO50 -40 3E>30W%3(11C1阿率r图3-4解调信号解调时域谱和频域谱编写数字带通调制函数z=ddemod(y,Fc,Fd,Fs,'ask',2) 进行2ASK解调，其中x为输入信号，Fc

9、为载波频率，Fd为码速率，Fs为米样频率。再进行频域转换3-4。FFT4=fft(z,64);FFT4=abs(FFT4),由于调制到解调的信道中没有加入噪声，所以调制前信号(即原始信号)与调制后信号相同，其时域谱和频域谱如图22AlDtH100150200350 4QQ 450500600时|'可t3.4 叠加噪声的ASK调制解调-J- F-Lgur e 5口回因Zil * ERi t aCi frwtInnl r 口a Tkt <ip lindrtwA mu号 h 0包心運）迈 S 1 n Bl加叮、噪声信号图3-5叠加小噪声调制信号时域谱和频域谱时同I图3-6叠加大噪声调制

10、信号时域谱和频域谱通常称它为加性噪声 和大噪声(信噪比为 其时域谱和频域谱Ul旦凶21101彗Fix«£3.1 -LIJi twXikser t.Tqo1=D亡玉*t幻pWirtiowH±l-p 03 h食包巒i ®4< 目I厂口加小噪声擀调洁w20-、 垦4a101214时可图3-7叠加小噪声解调信号时域谱和频域谱匚II旦|区|124G0101214 1G 1Q 20H'J par t加大嘩卫孵训佶弓顚谱£.1 Xa Eii t. Ifi匚jiwcirP工4口.五D.«=kt.op比.indovU-cZp ww曰Ife

11、 致包r?同箱I沪 a ra力口夭堞芦斛洞信号5014Ci30LJ_2010LJIXtJHO'£030 频率f4C5050由于信道中的噪声是叠加在信号上的，噪声始终是存在的，或加性干扰。对已调信号y分别叠加高斯小噪声(信噪比为 6) -2)则编写 Ynt1=awgn(y,6);丫nt2=awgn(y,-2) 并进行频域转换, 分别如图3-5和图3-6。图3-8叠加大噪声解调信号时域谱和频域谱对小噪声Ynt1和大噪声信号Ynt2分别解调，编写数字带通调制函数Z仁ddemod(Y nt1,Fc,Fd,Fs,'ask',2);z2=ddemod(Y nt2,Fc,F

12、d,Fs,'ask',2);并进行频域转换，其时域谱和频域谱如分别如图3-7和图3-8。J F-i CUT"n PMile £di t ZKLZns rrt Zcol sli-ndrw tLcLpI 禺曰 <*7 怎 ini I 巴口1010'<=laYi±iriliii1 0,三TT - -卜厂；口_= E! Iin1J.L=!= = = = = = = ZiT 7 Z T Z - -C 二-二 J._i_1fc_-1_d卜 -1卜H-=iri_ .：hriihi- TThiri- I » = =-«- V

13、 = <£=£71 "=2 = *» = » = = ：=T ；i 9 E ±=5訂二*三=三汁 - ” =t E：=£ = ；学=1 =r三甲h = r = = V = I = J z= I: z_ _ _ _ _ _ ± _ _一_ 一'_ _ t Y H -1 iJ三* 十 FryPs_ _ C _ _ _ _ _ _ - : _T_ : ' _ _ _ -1II11r弋-( 1115 = T z = ; 5-= T三 r 7 z WF h/Ih -10 =104”202-a倍口探 比

14、SNRCr/dD)图3-9误差数和误差率误码率是衡量一个数字通信系统性能的重要指标。在信道高斯白噪声的干扰下， 二进制2ASK数字调制误码率取决于解调器输入信噪比，编写误码率函数br,Pe1=symerr(x,z1);br,Pe2=symerr(x,z2)，其中 br 为符号误差数，Pel 为符号误差率。绘制如图3-9。调用函数semilogy(SNR,Pe);绘制信噪比与误码率的关系曲线如图 3-10，由此可 得出：与无噪声时(误码率为零)相比较，当信噪比较大时，噪声小误码率低；反之， 信噪比较小时噪声大误码率高。4结论通过此次的课程设计，使我收获了很多。不仅让我对专业方面的知识认识更加深

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

16、有所帮助。信号与通信系统课程设计说明书155附录%程序名称:ASK.m%程序功能:ASK解调解调程序代码cic;%清除命令工作窗里的内容clear;%清除内存空间变量x的最小整数%产生二进制随机序列x=ceil(rand(1,100000) -0.5)%产生二进制随机序列并取大于figure(1)%窗口 1,包含时域谱和频域谱sub plot(2,1,1)%第2个图像xlabel('时间 t');%x轴标注ylabel('序列值');%y轴标注title('二进制随机序列');%添加图像标题axis(1 21 -1 2)%控制坐标轴的范围grid

17、 on%图像中添加栅格%对随机序列进行频谱分析FFT1=fft(x,128);%对随机序列进行傅里叶变换FFT仁abs(FFT1);%对傅里叶变换取绝对值figure(1)sub plot(2,1,2)%第2个图像plot(FFT1);xlabel('频率 f);ylabel('幅度 FFT1');title('随机序列频谱');axis(0 128 0 50)grid on%载波信号t=1/360:1/360:20;%载波时间范围Fc=36;%载波频率%分块图函数sub plot，图形窗口分成2块子窗口的第1个图像carry=cos(2* pi*Fc*

18、t);%正弦载波信号stairs(x);figure(2)%窗口 2,包含时域谱和频域谱sub plot(2,1,1)pl ot(carry);xIabelC 时间 t');ylabel('幅度 carry');title（'载波信号'）;axis(1 600 -2 2)grid on%对载波信号进行频谱分析FFT2=fft(carry,256);%对载波信号进行傅里叶变换FFT2=abs(FFT2);%对傅里叶变换取绝对值figure(2)sub plot(2,1,2)plot(FFT2);xlabel('频率 f);ylabel('幅

19、度 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);end end%对20个随机码元进行判别，若码元为0则该码元周期内调制信号为零 plot(yy);xlabel(&#

20、39;时间 t');ylabelC 幅度 y');title（'已调信号'）;axis(1 600 -2 2) grid on%对已调信号进行频谱分析FFT3=fft(y,256);%对已调信号进行傅里叶变换FFT3=abs(FFT3);%对傅里叶变换取绝对值figure(3) sub plot(2,1,2) plot(FFT3);xlabel('频率 f);ylabel('幅度 FFT3');title（'已调信号频谱'）;axis(0 256 0 50) grid on%ASK的解调z=ddemod（y,Fc,Fd,F

21、s,'ask',2）; %调用数字带通调制函数dmod进行2ASK解调figure(4)%对傅里叶变换取绝对值sub plot(2,1,1) stairs(z);xlabel('时间 t');ylabel('幅度 z');title（'解调信号'）;axis(1 21 -1 2) grid on%对解调信号进行频谱分析 figure(4) sub plot(2,1,2) plot(FFT4);FFT4=fft(z,64);%对解调信号进行傅里叶变换FFT4=abs(FFT4);%对傅里叶变换取绝对值信号与通信系统课程设计说明书xI

22、abelC 频率 f);ylabelC 幅度 FFT4');titleC解调信号频谱');axis(O 64 0 50) grid on%加入高斯小噪声，SNR为6Yn t1=awg n(y,6);%加入高斯小噪声，信噪比为6figure(6)sub plot(2,1,1)17figure(5) sub plot(2,1,1) plot(Y nt1);xlabel('时间 t');ylabel('幅度 Ynt1');title('加小噪声信号');axis(1 600 -2 2) grid on%对加小噪声信号进行频谱分析FFT5

23、=fft(Y nt1,256);%对加入小噪声的调制信号进行傅里叶变换FFT5=abs(FFT5);%对傅里叶变换取绝对值figure(5) sub plot(2,1,2) plot(FFT5);xlabel('频率 f);ylabel('幅度 FFT5');title('加小噪声信号频谱')axis(0 256 0 50) grid on%ASK加小噪声信号的解调及误码率 z1=ddemod(Y nt1,Fc,Fd,Fs,'ask',2);%调用数字带通调制函数 dmod对加小噪声信号进行解调 br,P e1=symerr(x,z1)%

24、对解调后加小噪声信号误码分析，br为符号误差数，Pe1为符号误差率信号与通信系统课程设计说明书stairs(zl);xIabelC 时间 t');ylabelC 幅度 z1');title（'加小噪声解调信号'）;axis(1 21 -1 2) grid on%对加小噪声解调信号进行频谱分析FFT6=fft(z1,64);%对加入小噪声的解调信号进行傅里叶变换FFT6=abs(FFT6);%对傅里叶变换取绝对值figure(6)sub plot(2,1,2)plot(FFT6);xlabel('频率 f);ylabel('幅度 FFT6'

25、);title（'加小噪声解调信号频谱'）;axis(0 64 0 50)grid on%加入高斯大噪声，SNR为-2Yn t2=awg n(y,3);%加入高斯大噪声，信噪比为-2figure(7)sub plot(2,1,1)plot(Y nt2);xlabel('时间 t');ylabel('幅度 Ynt2');title（'加大噪声信号'）;axis(1 600 -2 2)grid on%对加大噪声信号进行频谱分析FFT7=fft(Y nt2,256);%对加入大噪声的调制信号进行傅里叶变换FFT7=abs(FFT7);%

26、对傅里叶变换取绝对值figure(7)sub plot(2,1,2)#plot(FFT7);信号与通信系统课程设计说明书xIabelC 频率 f);ylabelC 幅度 FFT5');title（'加大噪声信号频谱'） axis(O 256 0 50) grid on%ASK加大噪声信号的解调及误码率 z2=ddemod(Y nt2,Fc,Fd,Fs,'ask',2);%调用数字带通调制函数 dmod对加大噪声信号进行解调 br,P e2=symerr(x,z2)%对解调后加大噪声信号误码分析，br为符号误差数，Pe1为符号误差率 figure(8) sub plot(2,1,1) stairs(z2);xlabel('时间 t');ylabel('幅度 z2');ti