AM模拟调制系统的设计与仿真

1、摘要调幅，英文是Amplitude?Modulation(AM)。调幅也就是通常说的中波，范围在503-1060KHZ。调幅是用声音的高低变为幅度的变化的电信号。？本课程设计主要研究了AM模拟调制系统的设计和仿真。在本次通信系统仿真训练中，我主要通过了解模拟幅度调制和解调的原理和其实现方法， 然后根据其模拟幅度调制系统的原理给出了调制和解调的框图。其次弄懂了AM模拟调制的基本原理。最后利用Matlab软件仿真模拟幅度调制系统，实现AM调制和相干解调， 给出了调制信号、 载波信号及已调信号及解调信号的波形图和频谱图， 并计算了该系统的信噪比。关键词：调制解调AM模拟调制信噪比目录前言调制在通信系

2、统中的作用是至关重要的。所谓调制，就是把信号转换成适合在信道中传输的形式的一种过程。广义的调制分为基带调制和带通调制（也称载波调制）。在大多数场合，调制一般指载波调制。载波调制，就是用调制信号去控制载波的参数的过程，使载波的某一个或某几个参数按照调制信号的规律而变化。调制信号是指来自信源的信息信号（基带信号），这些新号可以是模拟的，也可以是数字的。未接受调制的周期性振荡信号称为载波，它可以是正弦波，也可以是非正弦波。载波调制后称为已调信号，它包含有调制信号的全部特征。解调则是调制的逆过程，其作用是将已调信号中的调制信号恢复出来。此次设计主要进行模拟调至系统的模拟和仿真，最常用和最重要的模拟调制

3、方式是用正弦波作为载波的幅度调制和角度调制。常见的调幅AM双边带DSB单边带VSB等调制就是幅度调制的几个典型实例；而频率调制FM是角度调制中被广泛采用的一种。本文主要分析了AM在高斯白噪声影响下的波形变化，通过对有无噪声解调信号波形的对比分析，估计AM调制解调系统的性能。一、调制及解调原理调制：将各种数字基带信号转换成适于信道传输的数字调制信号（已调信号或频带信号）；时域定义：调制就是用基带信号去控制载波信号的某个或几个参量的变化，将信息荷载在其上形成已调信号传输，而解调是调制的反过程，通过具体的方法从已调信号的参量变化中将恢复原始的基带信号。频域定义：调制就是将基带信号的频谱搬移到信道通带

4、中或者其中的某个频段上的过程，而解调是将信道中来的频带信号恢复为基带信号的反过程.根据所控制的信号参量的不同，调制可分为：调幅，使载波的幅度随着调制信号的大小变化而变化的调制方式。调频，使载波的瞬时频率随着调制信号的大小而变，而幅度保持不变的调制方式。调相，利用原始信号控制载波信号的相位。调制的目的是把要传输的模拟信号或数字信号变换成适合信道传输的信号，这就意味着把基带信号（信源）转变为一个相对基带频率而言频率非常高的代通信号。该信号称为已调信号，而基带信号称为调制信号。调制可以通过使高频载波随信号幅度的变化而改变载波的幅度、相位或者频率来实现。调制过程用于通信系统的发端。在接收端需将已调信号

5、还原成要传输的原始信号，也就是将基带信号从载波中提取出来以便预定的接受者（信宿）处理和理解的过程。该过程称为解调。调制的方式有很多。根据调制信号是模拟信号还是数字信号，载波是连续波(通常是正弦波)还是脉冲波，相应的调制方式有模拟连续波调制(简称模拟调制)、数字连续波调制(简称数字调制)、模拟脉冲波调制和数字脉冲波调制等。计算机内的信息是由“0”和“1”组成数字信号，而在电话线上传递的却只能是模拟电信号(模拟信号为连续的，数字信号为间断的)。于是，当两台计算机要通过电话线进行数据传输时，就需要一个设备负责数模的转换。这个数模转换器就是我们这里要讨论的Modem计算机在发送数据时，先由Modem巴

6、数字信号转换为相应的模拟信号，这个过程称为“调制”，也成D/A转换。经过调制的信号通过电话载波传送到另一台计算机之前，也要经由接收方的Modena责把模拟信号还原为计算机能识别的数字信号，这个过程我们称“解调”，也称A/D转换。正是通过这样一个“调制”与“解调”的数模转换过程，从而实现了两台计算机之间的远程通讯。解调原理调制过程的逆过程叫做解调。AM信号的解调是把接收到的已调信号还原为调制信号国。)。AM号的解调方法有两种：相干解调和包络检波解调。(1)相干解调(2)由AM信号的频谱可知，如果将已调信号的频谱搬回到原点位置，即可得到原始的调制信号频谱，从而恢复出原始信号。解调中的频谱搬移同样可

7、用调制时的相乘运算来实现。相干解调的原理框图如图1所示。SmBPF一ZMLPFm0*图1相干解调器一般模型将已调信号乘上一个与调制器同频同相的载波，得1.1.一SAM(t)cosctAom(t)cosctAom(t)A。m(t)cos2J(1)22由上式可知，只要用一个低通滤波器，就可以将第1项与第2项分离，无失真的恢复出原始的调制信号1.m0(t)-Aom(t)(2)2相干解调的关键是必须产生一个与调制器同频同相位的载波。如果同频同相位的条件得不到满足，则会破坏原始信号的恢复。(2)包络检波法由狈的波形可见，AM信号波形的包络与输入基带信号册Q)成正比，故可以用包络检波的方法恢复原始调制信号

8、。包络检波器一般由半波或全波整流器和低通滤波器组成。包络检波法属于非相干解调法，其特点是：解调效率高，解调器输出近似为相干解调的2倍；解调电路简单，特别是接收端不需要与发送端同频同相位的载波信号，大大降低实现难度。故几乎所有的调幅(AM式接收机都采用这种电路。二、模拟调制模拟调制是指用来自信源的基带模拟信号去调制某个载波，而载波是一个确知的周期性波形。模拟调制可分为线性调制和非线性调制，本文主要研究线性调制。线性调制的原理模型如图所示。设c(t)Acos2fot,调制信号为m，已调信号为s。图2线性调制的远离模型调制信号m(t)和载波在乘法器中相乘的结果为：s(t)m(t)Acosot,然后通

9、过一个传输函数为H的带通滤波器，得出已调信号为。s(t)m(t)cosot*h(t)1(3)s(f)-M(ffo)M(ffo)H(f)式(2-1)中，M(f)为调制信号m(t)的频谱。由于调制信号m和乘法器输出信号之间是线性关系，所以成为线性调制。带通滤波器H(f)可以有不同的设计，从而得到不同的调制种类。在幅度调制的一般模型中，若假设滤波器为全通网络，调制信号m(t)中无直流分量，则输出的已调信号就是无载波分量的双边带调制信号，或称抑制载波双边带(DSB调制信号，简称双边带(DSB信号。设正弦型载波c(t)Acos(用,式中：A为载波幅度，稣为载波角频率。根据调制定义，幅度调制信号(已调信号

10、)一般可表示为：Sm(t)Am(t)cos(ct)(4)其中，m为基带调制信号。设调制信号m(t)的频谱为M(w),则由公式2-2不难得到已调信号Sm(t)的频谱：A.Sm()M(c)M(c)(5)2由以上表示式可见，在波形上，幅度已调信号随基带信号的规律呈正比地变化；在频谱结构上，它的频谱完全是基带信号频谱在频域内的简单搬移。标准振幅就是常规双边带调制，简称调幅(AM)。假设调制信号m(t)的平均值为0,将其叠加一个直流偏量A。后与载波相乘，即可形成调幅信号。其时域表达式为：SAM(t)(A0m(t)cos(ot)(6)式中：Ao为外加的直流分量；m(t)可以是确知信号，也可以是随机信号。若

11、为确知信号，则AMJ号的频谱为：Sm()A)(c)(c)1M(c)M(JAM信号的频谱由载频分量、上边带、下边帝人部分组成。AM信号的总功率包括载波功率和边带功率两部分。只有边带功率才与调制信号有关，也就是说，载波分量并不携带信息。因此，AM信号的功率利用率比较低。AM调制器模型如下图所示。图3AM调制器模型AM1号的时域和频域表达式分别为SAM(t)Am(t)COSc(t)e(8)Aocosc(t)m(t)cosc(t)1SAM(t)Ao(c)(c)M(c)M(c)(9)式中，Ao为外加的直流分量；m(t)可以是确知信号也可以是随机信号， 但通常认为其平均值为0,即m(t)0o由频谱可以看出

12、，AM信号的频谱由载波分量、上边带、下边带三部分组成。上边带的频谱结构与原调制信号的频谱结构相同，下边带是上边带的镜像。因此，AM信号是带有载波分量的双边带信号，他的带宽是基带信号带宽fH的2倍，即AMM制典型波形和频谱如图1-1所示：图4AM调制典型波形和频谱如果在AM调制模型中将直流A去掉，即可得到一种高调制效率的调制方式一抑制载波双边带信号(DSB-SC),简称双边带信号。其时域表达式为SDSB(t)m(t)cos(ct)(10)式中，假设的平均值为0。DSB的频谱与AM的谱相近，只是没有了在3 处的函数6,即1 .Sm-M(c)M(c)(11)2AMB 调原理与抗噪性能解调是调制的逆过

13、程，其作用是从接收的已调信号中恢复原基带信号(即调制信号)。解调的方法可分为两类：相干解调和非相干解调(包络检波)相干解调，也称同步检波，为了无失真地恢复原基带信号，接收端必须提供一个与接收的已调载波严格同步(同频同相)的本地载波(称为相干载波)，它与接受的已调信号相乘后,经低通滤波器取出低频分量，即可得到原始的基带调制信号。包络检波器就是直接从已调波的幅度中提取原调制信号，通常由半波或全波整流器和低通滤波器组成。AM因此，解调器输出端的有用信号功率为：212TSOm0(t)m(t)(15)4解调AM信号时，接收机中的带通滤波器的中心频率?与调制频率”相同，此解调器输入端的窄带噪声：n(t)n

14、c(t)cos(ct)ns(t)sin(ct)(16)它与相干载波cosset)相乘后，得：,、，、1,、1_,、一.ni(t)cos(ct)-nc(t)一nc(t)cos(2ct)ns(t)sin(2&(17)22经低通滤波器后，解调器最终输出噪声为：1n0(t)n/t)(18)2故输出噪声功率为：2111Non2(t)n；(t)-NiNB(19)444式中，B23，为AM的带通滤波器的带宽，no为噪声单边功率谱密度。解调器输入信号平均功率为：12,Si-m2(t)(20)2可得解调器的输入信噪比：n(t)图5AM相干解调性能分析模型设解调器输入信号为:与相干载mco，、2，、m(t

15、)cos(ct)m(t)cos(1ct)低通滤波器-m(t)cos(ct)2经低通滤波器后，输出信号为：m(t)cos(221cosctct)mo(t)m。俨no(t)(13)(14)NiSinB法三种，其中窗函数设计法是最常用的，其次是频率采样法，但这两种方法在设计中还会存在一些不足之处，所以需要优化的设计方法，而等波纹逼近法很好的弥补了窗函数法和频率采样法的不足。对于数字高通、带通滤波器的设计，通用方法为双线性变换法。可以借助于模拟滤波器的频率转换设计一个所需类型的过渡模拟滤波器，再经过双线性变换将其转换策划那个所需的数字滤波器。具体设计步骤如下：(1)确定所需类型数字滤波器的技术指标。(

16、2)将所需类型数字滤波器的边界频率转换成相应的模拟滤波器的边界频率，转换公式为QCD)(3)将相应类型的模拟滤波器技术指标转换成模拟低通滤波器技术指标。(4)设计模拟低通滤波器。(5)通过频率变换将模拟低通转换成相应类型的过渡模拟滤波器。(6)采用双线性变换法将相应类型的过渡模拟滤波器转换成所需类型的数字滤波器。我们知道，脉冲响应不变法的主要缺点是会产生频谱混叠现象，使数字滤波器的频响偏离模拟滤波器的频响特性。为了克服之一缺点，可以采用双线性变换法。下面我们介绍用窗函数法设计FIR滤波器的步骤。如下：(1)根据对阻带衰减及过渡带的指标要求，选择用窗数类型(矩形窗、三角窗、汉宁窗、哈明窗、凯塞窗

17、等)，并估计窗口长度No先按照阻带衰减选择窗函数类型。原则是在保证阻带衰减满足要求的情况下，尽量选择主瓣的窗函数。(2)构造希望逼近的频率响应函数。(3)计算h(n)。(4)加窗得到设计结果解调器的输出信噪比：So点因此制度增益为：GAM也就是说，AM信号的解调器使信噪比改善一倍。FIR数字滤波器的设计方法主要有窗函数设计充m2(t)，BSo/NoSi/Ni(22)(23)频率采样设计法以及等波纹逼近设计接下来，我们根据语音信号的特点给出有关滤波器的技术指标：低通滤波器的性能指标：通带边界频率fp=300Hz,阻带截止频率fc=320Hz,阻带最小衰减As=100db,通带最大衰减Ap=1dB

18、在Matlab中，可以利用函数firl设计FIR滤波器，利用函数butter,cheby1和ellip设计IIR滤波器，利用Matlab中的函数freqz画出各步步器的频率响应。hn=fir1(M,wc,window),可以指定窗函数向量window。如果缺省window参数，则fir1默认为哈明窗。其中可选的窗函数有RectangularBarlrttHammingHannBlackman窗，其相应的都有实现函数。MATLABJ号处理工具箱函数buttpbuttorbutter是巴特沃斯滤波器设计函数，具有5种调用格式，本课程设计中用到的是：N,wc=butter(N,wc,Rp,As,s)

19、,该格式用于计算巴特沃斯模拟滤波器的阶数N和3dB截止频率wc。MATLA第号处理工具箱函数cheblap,cheblord和cheeby1是切比雪夫I型滤波器设计函数。我们用至U的是cheeby1函数，其调用格式如下：B,A=cheby1(N,Rp,wpo,ftypr)B,A=cheby1(N,Rp,wpo,ftypr,s)函数butter,cheby1和ellip设计IIR滤波器时都是默认的双线性变换法， 所以在设计滤波器时只需要代入相应的实现函数即可。三、AM调制解调系统的MATLA昉真及其分析MATLA 篁现信号DSB调制采用MATLAB函数modulate实现，其函数格式为：Y=MO

20、DULATE(X,Fc,Fs,METHOD,OPT为基带调制信号，Fc为载波频率，Fs为抽样频率，METHOD为调制方式选择，DSBWI制时为am,OPMDSBiM制时可不选，Fs需满足Fs2*Fc+BWBW的调制信号带宽。DSB信号解调采用MATLAB函数demod实现， 其函数使用格式为：X=DEMOD(Y,Fc,Fs,METHOD,OPT讷DSB已调信号，Fc为载波频率，Fs为抽样频率，METHOD解调方式选择，DSB解调时为ani,OPTEDSBiM制时可不选。观察信号频谱需对信号进行即可得到原始的基带调制信号。图9为经过相干解调后的滤波前的解调信号波形图及其频谱滤波后的解调信号波形1

21、050度幅-1000.10.20.30.40.50.60.70.80.91tx105滤波后的解调信号频谱54度幅10050100150200250300350400450500频率1050-5-10C1050-5图。图10滤波后0.10.20.30.40.50.60.70.80.9时间加噪声2的解调信号波形的解调信号的-1000.10.20.30.40.50.60.70.80.9傅里叶变换，采用MATLAB数fft实现，其函数常使用格式为：Y=FFT(X,N),X为时域函数，N为傅里叶变换点数选择，一般取值20频域变换后，对频域函数取模，格式：Y1=ABS(Y)再进行频率转换，转换方法：f=(

22、0:length(Y)-1)*Fs/length(Y)分析解调器的抗噪性能时，在输入端加入高斯白噪声，采用MATLAB1数awgn实现，其函数使用格式为：Y=AWGN(X,SNR)加高斯白噪声于X中，SNM信噪比，单位为dB,其值在假设X的功率为0dBM的情况下确定。信号的信噪比为信号中有用的信号功率与噪声功率的比值，根据信号功率定义，采用MATLAB函数var实现，其函数常使用格式为：Y=VAR(X),返回向量的方差，则信噪比为:SNR=VAR(X1)/VAR(X2)绘制曲线采用MATLABS数plot实现，其函数常使用格式：PLOT(X,Y),X为横轴变量，Y为纵轴变量，坐标范围限定AXI

23、S(x1x2y1y2),轴线说明XLABEL和YLABEL程序设计流程图见附录。MATLAB 仿真及其分析图6载波信号的波形及其频谱图7调制信号的波形及其频谱载波是指被调制以传输信号的波形，一般为正弦波。载波信号，就是把普通信号加载到一定频率的高频信号上，在没有加载普通信号的高频信号时，高频信号的波幅是固定的，加载之后波幅就随着普通信号的变化而变化。本设计采用频率为8000Hz振幅为10的余弦载波信号，其波形及频谱图如图6所示。调制信号是指来自信源的消息信号(基带信号)，这些信号可以是模拟的，也可以是数字的。本设计中产生的基带信号频率为10H乙振幅为5的余弦信号。该基带信号的波形及其频谱如图7

24、所示。图8AM调制信号的波形及其频谱图9滤波前的解调信号的波形及其频谱调制就是由调制信号去控制高频载波的幅度，使之随调制信号作线性变化的过程，本设计的AM已调信号的波形及其频谱如图8所示，从图中可以看出已调信号幅度随基带信号的幅度变化而变化，且把频率为10HZ的基带信号搬移到频率为8000Hz的载频上，实现了信道的多路复用，提高了信道利用率。由式(4)可知，其解调后的信号中含有高频及噪声成分，经低通滤波器取出低频分量,加噪声1的解调信号波形时间波形及其频谱图11加入不同噪声的解调信号波形图10为经过低通滤波器，滤除了噪声及高频成分之后的解调信号的波形图及其频谱图，其波形与基带信号波形基本一致。

25、将信道中存在的不需要的电信号统称为噪声， 通信系统的噪声是叠加在信号上的， 没有传输信号时通信系统中也有噪声，噪声永远存在于通信系统中。图11为在信道中加入不同程度的高斯白噪声经过低通滤波器后输出的解调信号波形，可以看出其解调出的信号波形图发生明显的失真。图12输入/输出信噪比图13低通滤波器的幅度函数图12为该系统采用相干解调时输入、输出信噪比的对比图。总结在课程设计过程中，着重研究了DSB号调制与解调原理和MATLA瞅拟实现，熟悉了信号波形、频谱的和系统性能的分析方法，了解了数字滤波器的设计与使用方法，综合提高了自己的专业技能。通过对程序的设计， 我进一步熟悉了MATLAB1发环境， 对M

26、ATLAB勺一些具体操作和应用有了更深入的了解。如：有要求的正弦信号的产生，基本图形的绘制和各种的函数的使用等。同时，这次设计我对数字信号处理和通信原理课本上学到的知识点有了更深入的理解和掌握。比如对信号的调制和解调过程有了更深层的理解，学会了如何使用MATLAB寸信号进行SS琳制和解调，了解了低通滤波器的MATLABJ计方法。还有很重要的一点是，我学会了如何安排设计所需的时间及合理利用网络资源等普遍实用的学习方法，通过和同学探讨，拓宽了我的眼界。参考文献1,李建新.现代通信系统分析与仿真一MATLAB!信工具箱.西安：西安电子科技大学出版社，2000.2,樊昌信,通信原理,北京：国防工业出版

27、社，2002.3,刘敏,MATLAB,通信仿真与应用,北京：国防工业出版社.5,吴伟陵等着,移动通信原理,北京：电子工业出版社，2005附录输DSB周制信号加入白噪声对已调信号解调显示解调后，瑚耒制牖M第i|岫鸠尚波形及频谱结束总体设计结构图程序代码：载波心A0=10;%载波信号振幅A1=5;%调制信号振幅fs=8000;%载波彳言号频率fc=1000;ws=fs*pi;wc=fc*pi;y=A0*cos(ws*t);%载波心figure(1);subplot(2,2,1);plot(t,y);xlabel(t);title（载频信号波形）;gridon;subplot(2,2,2);Y=ff

28、t(y);%对y的傅里叶变换plot(abs(Y);%频谱取模xlabel(所)；title（载波信号频谱）;gridon;调制信号A0=10;%载波信号振幅A1=5;%调制信号振幅fs=8000;%载波彳言号频率fc=1000;ws=fs*pi;wc=fc*pi;%调制信号figure(2);subplot(2,1,1);plot(t,x);xlabel(t);title(调制信号);axis(0,1,-10,10);%gridon;subplot(2,1,2);X=fft(x);%plot(abs(X);%xlabel(频率)；title(调制信号频谱);gridon;AM已调信号A0=1

29、0;A1=5fs=8000;fc=1000;ws=fs*pi;wc=fc*pi;%调制信号%调制度%AM已调信号figure(3);subplot(2,1,1);plot(t,y1);gridon;xlabel(t);title(AM调制信号波形);axis(0,1,-10,10);subplot(2,1,2);Y1=fft(y1);%plot(abs(Y1);%gridon;xlabel(频率)；定义坐标轴的范围对x的傅里叶变换取模值傅里叶变换取模title(AM调制信号频谱);axis(0,15000,0,500000);gridon;解调信号A0=10;A1=5fs=8000;fc=10

30、00;ws=fs*pi;wc=fc*pi;%AM调制信号已调信号调制度y11=y1.*cos(ws*t);%对AM已调信号进行解调figure(4);subplot(2,1,1);plot(t,y11);gridon;xlabel(t);title(未滤波的解调倍);axis(0,1,0,10);subplot(2,1,2);Y11=fft(y11);%傅里叶变换plot(abs(Y11);%取模gridon;xlabel(所)；title(未滤波的解调倍);axis(0,500,0,500000);gridon;AMB调彳S号FIR滤波A0=10;A1=5fs=8000;%AMFt=2000

31、;%fpts=100,120;%阻带截止频率fs=120Hz;mag=10;n21,wn21,beta,ftype=kaiserord(fpts,mag,dev,Ft);%计的FIR滤波器的参数b21=fir1(n21,wn21,Kaiser(n21+1,beta);%h,w=freqz(b21,1);%figure(3);plot(w/pi,abs(h);gridontitle(FIR低通滤波器);figure(4);subplot(2,1,1);y111=fftfilt(b21,y11);%plot(t,y111);gridon;xlabel(t);title(,滤波后的解调信号波形,);axis(0,1,-10,10);subplot(2,1,2);Y111=fft(y111);%plot(abs(Y111);%gridon;xlabel(频率)；title(,滤波后的解调信号频谱,);axis(0,500,0,500000);gridon;加噪声A0=10;A1=5fs=8000;fc=1000;ws=fs*pi;wc=fc*pi;%基带信号调制度fc=1000;ws=fs*pi;wc=fc*pi;调制信号y11=y1.*cos(ws*t);%对AME调信号进行解调已调信号