MATLAB课程设计

MATLAB 课 程 设 计 报 告     目：利用MATLAB实现信号的幅度调制与解调 班  级： 姓  名： 学  号： 姓    名： 学    号： 指导教师： 一、摘要： 现代通信系统要求通信距离远、信道容量大、传输质量好。在信号处理里面经常要用到调制与解调，而信号幅度调制与解调是最基本，也是经常用到的。用AM调制与解调可以实现很多功能，制造出很多的电子产品。本设计主要研究内容是利用MATLAB实现对正弦信 进行双边带幅度调制，载波信号频率为100Hz，在MATLAB中显示调制信号的波形和频谱，已调信号的波形和频谱，比较信号调制前后的变化。并对已调信号解调，比较了解调后的信号与原信号的区别。信号幅度调制与解调及MATLAB中信号示的基本及绘图函数的调用，实现了对连续时间信号的可视化表示。本文采用MATLAB对信号的幅度进行调制和解调。 二、关键词：幅度、 调制、 解调、 MATLAB。 三、设计内容 1.调制信号 调制信号是原始信息变换而来的低频信号。调制本身是一个电信号变换的过程。调制信号去改变载波信号的某些特征值（如振幅、频率、相位等），导致载波信号的这个特征值发生有规律的变化，这个规律是调制信号本身的规律所决定的。 1.1 MATLAB实现调制信号的波形 本设计的调制信号为正弦波信号 ，通过matlab仿真显示出其波形图如图1-1所示 图1-1调制信号 1.2 MATLAB实现调制信号的频谱 连续傅里叶变换是一个特殊的把一组函数映射为另一组函数的线性算子。傅里叶变换就是把一个函数分解为组成该函数的连续频率谱。在数学中，信号f(t)的傅里叶变换被认为是处在频域中的信号。 离散傅里叶变换的一种快速算法,简称FFT。为了节省电脑的计算时间，实现数字信号的实时处理，减少离散傅里叶变换（DFT）的计算量。 快速傅里叶变换（FFT），并不是一种新的傅立叶分析理论，而是减少DFT计算量的算法设计思想和DFT各种快速算，为计算机处理数据节省了时间。对已知信号进行快速傅里叶变换其抽样频率为100hz，matlab仿真频谱如图1-2 图1-2 已知信号的频谱 1.3 MATLAB实现载波的仿真 载波是被调制以传输信号的波形，一般为正弦波。一般要求正弦载波的频率远高于调制信号的带宽，否则会发生混叠，使传输信号失真。 我们一般需要发送的数据的频率是低频的，如果按照本身的数据的频率来传输，不利于接收和同步。使用载波传输，我们可以将数据的信号加载到载波的信号上，接收方按照载波的频率来接收数据信号，有意义的信号波的波幅与无意义的信号的波幅是不同的，将这些信号提取出来就是我们需要的数据信号。 本实验的载波频率为100hz，起波形图如图1-3所示 图1-3 载波波形图 2．信号的幅度调制 本章对信号的幅度调制原理作详细的解释，并通过matlab对已调信号进行仿真，并显示出已调信号的时域与频域的波形。 2.1 信号的调制 调制，就是把信号转化成适合在信道中传输的形式的一种。广义的调制分为基带调制和带通调制（也称载波调制）。载波调制，就是用调制信号区控制载波参数的过程，即使载波的某一个或几个参数按照调制信号的亏率而变化。调制信号时指来自信源的消息信号。未受调制的周期性振荡信号称为载波，载波调制后称为已调信号，它含有调制信号的全部特征。 2.2  幅度调制原理 1． 信号的幅度调制原理 幅度调制是用调制信号去控制高频载波的幅度，使其按调制信号的规律变化的过程。幅度调制器的一般模型如图2-1所示，载波信号 ，其中A0为外加信号。 图2-1 (2.1) (2.2) 由频谱可以看出，AM信号的频谱由载波分量，上边带，下边带三部分组成。上边带的频谱结构与原调制信号的频谱结构相同，下边带是上边带的镜像。因此AM信号时带有载波分量的双边带调制信号，它的带宽是基带信号带宽fH的2倍。 2.  双边带幅度调制 在AM调制模型中将直流A0去掉，即可得到一种高调制效率的调制方式----抑制载波双边带信号（DSB—SC），简称双边带信号DSB（如图2-2）。其时域表达式为： 图2-2  DSB信号的调制 其时域表达式为：               (2.3) 其频域相应为：           (2.4) 与先前幅度调节相比,只是少了直流分量，相应的少了 ，也就是在 的 处少了冲击响应，与AM信号比较，因为不存在载波分量，DSB信号的调制效率为100%。DSB信号解调时需要采用相干解调，也称同步检波。 2.3 MATLAB实现双边带幅度调制 本节主要通过matlab软件把调制信号 调制到载波信号 上，使载波的幅度按照调制信号的规律进行变化。其波形如图2-3所示 图2-3  已调信号时域波形 2.4 MATLAB实现已调信号的频谱图 与调制信号实现频谱图的原理一样，用matlab对已调信号进行快速傅里叶变换抽样频率为200hz，可以得到已调信号的频谱图，如图2-4所示 图3-4  已调信号频谱图 2.5 幅度调制前后的比较 双边带幅度调制是用调制信号去控制高频载波的幅度，使其按调制信号的规律变化的过程。调制后的信号是把调制信号加载到载波上，使载波上呈现调制信号的包络线。对于已调信号的频谱，实质就是把调制信号由低频段搬迁到高频段进行通信。其波形比较如图2-5所示 图2-5  调制前后时域与频域波形 3. 已调信号的解调 3.1 AM信号的解调原理及方式 解调是调制的逆过程，解调是将位于载波的信号频谱再搬回来，并且不失真的恢复出原始基带信号。对于幅度调制来说，解调是从它的幅度变化提取调制信号的过程。 解调的方式有两种：相干解调与非相干解调。相干解调适用于各种线性调制系统，非相干解调一般适用幅度调制（AM）信号。 1.AM信号的相干解调 所谓相干解调是为了从接受的已调信号中，不失真地恢复原调制信号，要求本地载波和接收信号的载波保证同频同相。相干载波的一般模型如下图3-1所示： 图 3-1  相干解调模型 将已调信号乘上一个与调制器同频同相的载波，得 (3.1) 由上式可知，只要用一个低通滤波器，就可以将第1项与第2项分离，无失真的恢复出原始的调制信号 (3.2) 相干解调的关键是必须产生一个与调制器同频同相位的载波。如果同频同相位的条件得不到满足，则会破坏原始信号的恢复。 2.AM信号的非相干解调 所谓非相干解调是在接收端解调信号时不需要本地载波，而是利用已调信号中的包络信号来恢复原基带信号[7]。因此，非相干解调一般只适用幅度调制（AM）系统。忧郁包络解调器电路简单，效率高，所以几乎所有的幅度调制（AM）接收机都采用这种电路。如下为串联型包络检波器的具体电路。其利用二极管的单向导电特性，将调幅高频信号去掉一半，再利用电容器的充放电特性和低通滤波器滤去高频分量，就可以得到与包络线形状相同的音频信号。 图3-2  AM信号的非相干解调原理 图3-2 当RC满足条件 时，包络检波器的输出基本与输入信号的包络变化呈线性关系，即 (3.3) 其中， 。隔去直流后就得到原信号 3.2 MATLAB实现已调信号的解调 图3-3  解调后的调制信号 解调过程是运用相干解调，对已调信号乘上载波，再通过低通滤波器滤波就可以无失真的恢复出原始的调制信号。此滤波器的设置参数为：通带损耗Rp为0.5，阻带衰减不小于60分贝，通带截止频率为Wp为0.03，阻带起始频率为0.1。其解调后的仿真波形如图3-3所示。 3.3 信号解调前后的比较 本设计是采用双边带调幅来实现的信号调制，用相干波解调后，没有直流分量，幅度减为原来的一半，频率不变化。由频谱图可知解调是把已调信号的频谱搬回到原来的位置，从而恢复出原始信号。恢复信号不失真的关键是低通滤波器的设置，尽可能的减小失真度。 图3-4  调制信号与解调信号对比 4.设计总程序如下： t =-2:0.001:2 y1=sin(t*40*pi) figure subplot(3,2,1) plot(t,y1) title('调制信号') xlabel('时间/s') ylabel('幅度') grid on xlim([-0.1,0.1])          %显示调制信号 fm=20                  %调制信号的频频率 fs=100 N=512 yk=fft(y1,N)      %傅里叶变换 %yw=abs(yk) yw=abs(fftshift(yk)) fw=[-255:256]/N*fs subplot(3,2,2) plot(fw,yw) title('调制信号的频谱') xlabel('频率/hz') ylabel('F(w)') grid on xlim([-8,8]) y3=cos(pi*200*t)    %载波频率100hz  w=2f*pi %subplot(3,2,3) subplot(3,2,3) plot(t,y3) title('载波信号') xlabel('时间/s')    %载波信号显示 ylabel('幅度') grid on xlim([-0.1,0.1]) y4=sin(t*40*pi).*cos(pi*200*t) subplot(3,2,4) plot(t,y4) title('已调信号') xlabel('时间/s') ylabel('幅度') grid on xlim([-0.1,0.1])    %已调信显示 ylim([-2,2]) fs1=200 ykdsb=fft(y4,N) ywdsb=abs(yk) fwdsb=[-255:256]/512*fs1 subplot(3,2,5) plot(fwdsb,ywdsb) title('已调信号的频谱') xlabel('频率/hz') ylabel('Fdsb(w)') grid on xlim([-150,150]) y6=y4.*y3 grid on Rp=0.5        %通带损耗不大于Rp分贝 Rs=60          %阻带衰减不小于RS分贝 Wp=0.03        %通带截止频率 Ws=0.1        %阻带起始频率 [n,Wn]=buttord(Wp,Ws,Rp,Rs,'s') [b,a]=butter(n,Wn,'low') X1=filter(b,a,y6) subplot(3,2,6) plot(t,X1) title('滤波后的f（t）信号') xlabel('时间/s') ylabel('幅度') grid on xlim([-2,-1]) 其对应仿真的图形如图4-1所示： 图4-1 四、设计总结 通过本次Matlab的课程设计，我对Matlab软件有了更深入的了解。为了做好这次的课程设计，查阅了大量资料，并在各个网站上搜索与此有关的知识，这个过程使我受益匪浅，了解到了自主探究学习的很多方法。我觉得这个是最重要的，对于今后任何一个领域或者某一方面的学习研究中都是大有益处的。 对于这次课程设计，其中难度最大的莫过于对调制与解调这部分原理的理解，在MATLAB软件的使用上面，自己的确存在很大的缺陷，对于一些语言还只是停留在理解上面，要自己熟练运用的确还有一段距离，在编程过程中，也努力过，努力想去继承前人的智慧，但无奈自己理解参透能力有限。所以最后在程序这方面还是有待加强。 Matlab软件是一个非常实用的软件，它不但可以进行复杂的数值运算，还能够对图像进行各种分析处理，因此，在今后的学习工作中，我还会继续对Matlab进行更加深入的学习和应用，用它来解决今后将要面临的各种繁琐的问题。 五、参考文献 1、陈后金. 信号与系统[M].北京：高等教育出版社，2007.07. 2、张洁.双边带幅度调制及其 MATLAB 仿真[J].科技经济市场，2006.9 3、陈洁，焦振宇. 基于MATLAB7.0 的信号调制与解调分析[J]. 山西电子技术. 2006(5): 4、傅祖芸. 信息论基础理论与应用[M]. 北京：电子工业出版社，2001.