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

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

2、比目录前言 1一、调制及解调原理 .21.1 调制原理 .21.2 解调原理 3二、模拟调制 .42.1 模拟调制原理 42.2 AM 调制的基本原理 .42.3 AM 解调原理与抗噪性能 .62.4 FIR 数字滤波器设计方法 8三、 AM 调制解调系统的 MATLAB 仿真及其分析 .103.1 AM 调制解调分析的 MATLAB 实现 .103.2 MATLAB 仿真及其分析 .10总结 13参考文献 .14附录 151前言调制在通信系统中的作用是至关重要的。所谓调制，就是把信号转换成适合在信道中传输的形式的一种过程。广义的调制分为基带调制和带通调制（也称载波调制） 。在大多数场合，调制

3、一般指载波调制。载波调制，就是用调制信号去控制载波的参数的过程，使载波的某一个或某几个参数按照调制信号的规律而变化。调制信号是指来自信源的信息信号（基带信号） ，这些新号可以是模拟的，也可以是数字的。未接受调制的周期性振荡信号称为载波，它可以是正弦波，也可以是非正弦波。载波调制后称为已调信号，它包含有调制信号的全部特征。解调则是调制的逆过程，其作用是将已调信号中的调制信号恢复出来。此次设计主要进行模拟调至系统的模拟和仿真，最常用和最重要的模拟调制方式是用正弦波作为载波的幅度调制和角度调制。常见的调幅 AM、双边带 DSB、单边带 VSB 等调制就是幅度调制的几个典型实例；而频率调制 FM 是角

4、度调制中被广泛采用的一种。本文主要分析了 AM 在高斯白噪声影响下的波形变化，通过对有无噪声解调信号波形的对比分析， ，估计 AM 调制解调系统的性能。2一、调制及解调原理1.1 调制原理调制: 将各种数字基带信号转换成适于信道传输的数字调制信号(已调信号或频带信号)；时域定义：调制就是用基带信号去控制载波信号的某个或几个参量的变化，将信息荷载在其上形成已调信号传输，而解调是调制的反过程，通过具体的方法从已调信号的参量变化中将恢复原始的基带信号。 频域定义：调制就是将基带信号的频谱搬移到信道通带中或者其中的某个频段上的过程，而解调是将信道中来的频带信号恢复为基带信号的反过程. 根据所控制的信号

5、参量的不同，调制可分为： 调幅，使载波的幅度随着调制信号的大小变化而变化的调制方式。 调频，使载波的瞬时频率随着调制信号的大小而变，而幅度保持不变的调制方式。 调相，利用原始信号控制载波信号的相位。 调制的目的是把要传输的模拟信号或数字信号变换成适合信道传输的信号，这就意味着把基带信号（信源）转变为一个相对基带频率而言频率非常高的代通信号。该信号称为已调信号，而基带信号称为调制信号。调制可以通过使高频载波随信号幅度的变化而改变载波的幅度、相位或者频率来实现。调制过程用于通信系统的发端。在接收端需将已调信号还原成要传输的原始信号，也就是将基带信号从载波中提取出来以便预定的接受者（信宿）处理和理解

6、的过程。该过程称为解调。 调制的方式有很多。根据调制信号是模拟信号还是数字信号，载波是连续波（通常是正弦波）还是脉冲波，相应的调制方式有模拟连续波调制（简称模拟调制） 、数字连续波调制（简称数字调制） 、模拟脉冲波调制和数字脉冲波调制等。计算机内的信息是由“0”和“1”组成数字信号，而在电话线上传递的却只能是模拟电信号（模拟信号为连续的，数字信号为间断的） 。于是，当两台计算机要通过电话线进行数据传输时，就需要一个设备负责数模的转换。这个数模转换器就是我们这里要讨论的Modem。计算机在发送数据时，先由 Modem 把数字信号转换为相应的模拟信号，这个过程称为“调制” ，也成 D/A 转换。经

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

8、 1 所示。 )(tSm )(tZ)(0tmt0cos图 1 相干解调器一般模型将已调信号乘上一个与调制器同频同相的载波，得（1）ttmAtttmAttS cccAM 2os)(21)(21os)(os)( 000 由上式可知，只要用一个低通滤波器，就可以将第 1 项与第 2 项分离，无失真的恢复出原始的调制信号（2）)(21)(00ttm相干解调的关键是必须产生一个与调制器同频同相位的载波。如果同频同相位的条件得不到满足，则会破坏原始信号的恢复。（2）包络检波法由 的波形可见，AM 信号波形的包络与输入基带信号 成正比，故可以用包络检波的方法恢复原始调制信号。包络检波器一般由半波或全波整流器

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

10、果为： ，然后通过一个传tAtms0 cos)(输函数为 H(f)的带通滤波器，得出已调信号为。从图 2.1 中可得已调信号的时域和频域表达式为：（3）式（2-1）中，M(f)为调制信号 m(t)的频谱。由于调制信号 m(t)和乘法器输出信号之间是线性关系，所以成为线性调制。带通滤波器 H(f)可以有不同的设计，从而得到不同的调制种类。2.2 AM 调制的基本原理在幅度调制的一般模型中，若假设滤波器为全通网络，调制信号 m(t)中无直流分量，则输出的已调信号就是无载波分量的双边带调制信号，或称抑制载波双边带（DSB）调制信号，简称双边带（DSB）信号。设正弦型载波 ，式中：A 为载波幅度， 为

11、载波角频率。()cos()tt根据调制定义，幅度调制信号（已调信号）一般可表示为：（4）()cos()mSttt其中，m(t)为基带调制信号。0()()cos*()1)(2oostmtthfMffHf5设调制信号 m(t)的频谱为 M( )，则由公式 2-2 不难得到已调信号 的频谱： )(tSm（5）()(2( ccmMAS由以上表示式可见，在波形上，幅度已调信号随基带信号的规律呈正比地变化；在频谱结构上，它的频谱完全是基带信号频谱在频域内的简单搬移。标准振幅就是常规双边带调制，简称调幅（AM） 。假设调制信号 m(t)的平均值为 0，将其叠加一个直流偏量 后与载波相乘，即可形成调幅信号。其

12、时域表达式为:0A（6）)cos()()00ttmAtSM式中： 为外加的直流分量；m(t)可以是确知信号，也可以是随机信号。0若为确知信号，则 AM 信号的频谱为:（7） AM 信号的频谱由载频分量、上边带、下边带三部分组成。AM 信号的总功率包括载波功率和边带功率两部分。只有边带功率才与调制信号有关，也就是说，载波分量并不携带信息。因此，AM 信号的功率利用率比较低。AM 调制器模型如下图所示。图 3 AM 调制器模型AM 信号的时域和频域表达式分别为（8）（9）式中， 为外加的直流分量；m(t)可以是确知信号也可以是随机信号，但通常认为其Ao平均值为 0，即 。0)(tm由频谱可以看出，

13、AM 信号的频谱由载波分量、上边带、下边带三部分组成。上边带的频谱结构与原调制信号的频谱结构相同，下边带是上边带的镜像。因此，AM 信号是带有载波分量的双边带信号，他的带宽是基带信号带宽 的 2 倍，即00()()cos()cosAMStmtt0 1()()()()()2MccccSt M0 1()()()()()2mccccSAM fHBAM6AM 调制典型波形和频谱如图 1-1 所示：图 4 AM 调制典型波形和频谱如果在 AM 调制模型中将直流 A 去掉，即可得到一种高调制效率的调制方式抑制载波双边带信号(DSBSC)，简称双边带信号。 其时域表达式为（10）)cos()(ttmSDB式

14、中，假设的平均值为 0。DSB 的频谱与 AM 的谱相近，只是没有了在 处的函数，即（11）)()(21ccmMS2.3 AM 解调原理与抗噪性能解调是调制的逆过程，其作用是从接收的已调信号中恢复原基带信号（即调制信号） 。解调的方法可分为两类：相干解调和非相干解调（包络检波） 。相干解调，也称同步检波，为了无失真地恢复原基带信号，接收端必须提供一个与接收的已调载波严格同步（同频同相）的本地载波（称为相干载波） ，它与接受的已调信号相乘7后，经低通滤波器取出低频分量，即可得到原始的基带调制信号。包络检波器就是直接从已调波的幅度中提取原调制信号，通常由半波或全波整流器和低通滤波器组成。AM 相干

15、解调性能分析模型如图 2.5 所示：图 5 AM 相干解调性能分析模型设解调器输入信号为： (12)()cos()mtSt与相干载波 cos( t)相乘后，得(13)2()1cs21)(cos)2 ttttttmcc经低通滤波器后，输出信号为:(14)()(0tt因此，解调器输出端的有用信号功率为:(15)(41)(220tmts解调 AM 信号时，接收机中的带通滤波器的中心频率 与调制频率 相同，此解调器输入端的窄带噪声:（16） )sin()cos()( tttnt ci 它与相干载波 cos( t)相乘后，得：（17）)2i()2()1)(2)cos( ttttttn csccci经低通

16、滤波器后，解调器最终输出噪声为：（18） )()(0tntc故输出噪声功率为：(19) BNttnNic0220 41)(41)(式中， ,为 AM 的带通滤波器的带宽， 为噪声单边功率谱密度。HfB20n解调器输入信号平均功率为：（20）)(21tmSi带 通滤 波 器sm(t) sm(t)n(t)ni(t)mo(t)no(t)低 通滤 波 器cosct8可得解调器的输入信噪比：（21）BntmNSi02)(1解调器的输出信噪比：（22）ntS02)(因此制度增益为:（23）iAMNSG/0也就是说，AM 信号的解调器使信噪比改善一倍。2.4 FIR 数字滤波器设计方法FIR 数字滤波器的设

17、计方法主要有窗函数设计法、频率采样设计法以及等波纹逼近设计法三种，其中窗函数设计法是最常用的，其次是频率采样法，但这两种方法在设计中还会存在一些不足之处，所以需要优化的设计方法，而等波纹逼近法很好的弥补了窗函数法和频率采样法的不足。对于数字高通、带通滤波器的设计，通用方法为双线性变换法。可以借助于模拟滤波器的频率转换设计一个所需类型的过渡模拟滤波器，再经过双线性变换将其转换策划那个所需的数字滤波器。具体设计步骤如下：（1） 确定所需类型数字滤波器的技术指标。（2） 将所需类型数字滤波器的边界频率转换成相应的模拟滤波器的边界频率，转换公式为 =2/T tan(0.5)（3） 将相应类型的模拟滤波

18、器技术指标转换成模拟低通滤波器技术指标。（4） 设计模拟低通滤波器。（5） 通过频率变换将模拟低通转换成相应类型的过渡模拟滤波器。（6） 采用双线性变换法将相应类型的过渡模拟滤波器转换成所需类型的数字滤波器。我们知道，脉冲响应不变法的主要缺点是会产生频谱混叠现象，使数字滤波器的频响偏离模拟滤波器的频响特性。为了克服之一缺点，可以采用双线性变换法。下面我们介绍用窗函数法设计 FIR 滤波器的步骤。如下：9（1）根据对阻带衰减及过渡带的指标要求，选择串窗数类型（矩形窗、三角窗、汉宁窗、哈明窗、凯塞窗等） ，并估计窗口长度 N。先按照阻带衰减选择窗函数类型。原则是在保证阻带衰减满足要求的情况下，尽量

19、选择主瓣的窗函数。（2）构造希望逼近的频率响应函数。（3）计算 h(n)。（4）加窗得到设计结果。接下来，我们根据语音信号的特点给出有关滤波器的技术指标：低通滤波器的性能指标：通带边界频率 fp=300Hz，阻带截止频率 fc=320Hz，阻带最小衰减 As=100db , 通带最大衰减 Ap=1dB。在 Matlab 中，可以利用函数 fir1 设计 FIR 滤波器，利用函数 butter,cheby1 和 ellip设计 IIR 滤波器，利用 Matlab 中的函数 freqz 画出各步步器的频率响应。hn=fir1(M，wc，window)，可以指定窗函数向量 window。如果缺省 w

20、indow 参数，则 fir1默认为哈明窗。其中可选的窗函数有 Rectangular Barlrtt Hamming Hann Blackman 窗，其相应的都有实现函数。MATLAB 信号处理工具箱函数 buttp buttor butter 是巴特沃斯滤波器设计函数，其有5 种调用格式，本课程设计中用到的是：N,wc=butter(N,wc,Rp,As,s),该格式用于计算巴特沃斯模拟滤波器的阶数 N 和 3dB 截止频率 wc。MATLAB 信号处理工具箱函数 cheblap,cheblord 和 cheeby1 是切比雪夫 I 型滤波器设计函数。我们用到的是 cheeby1 函数，其

21、调用格式如下：B,A=cheby1(N,Rp,wpo,ftypr)B,A=cheby1(N,Rp,wpo,ftypr,s)函数 butter,cheby1 和 ellip 设计 IIR 滤波器时都是默认的双线性变换法，所以在设计滤波器时只需要代入相应的实现函数即可。100 0.0010.0020.0030.0040.0050.0060.0070.0080.009 0.01-10010低低低低低低低低低低0 5000 10000 150000510x 105低低低低低低低低低低0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1-10-50510低低低低低低低低0 50

22、 100 150 200 250 300 350 400 450 5000510x 105低低低低低低低低低低三、 AM 调制解调系统的 MATLAB 仿真及其分析3.1 AM 调制解调分析的 MATLAB 实现信号 DSB 调制采用 MATLAB 函数 modulate 实现，其函数格式为： Y = MODULATE(X,Fc,Fs,METHOD,OPT) X 为基带调制信号，Fc 为载波频率，Fs 为抽样频率，METHOD 为调制方式选择，DSB 调制时为am ，OPT 在 DSB 调制时可不选，Fs 需满足 Fs 2*Fc + BW，BW 为调制信号带宽。 DSB 信号解调采用 MATL

23、AB 函数 demod 实现，其函数使用格式为： X = DEMOD(Y,Fc,Fs,METHOD,OPT) Y 为 DSB 已调信号，Fc 为载波频率，Fs 为抽样频率，METHOD为解调方式选择，DSB 解调时为am ，OPT 在 DSB 调制时可不选。 观察信号频谱需对信号进行傅里叶变换，采用 MATLAB 函数 fft 实现，其函数常使用格式为：Y=FFT(X,N)，X 为时域函数，N 为傅里叶变换点数选择，一般取值 2。频域变换后，对频域函数取模，格式：Y1=ABS(Y)，再进行频率转换，转换方法：f=(0:length(Y)-1)*Fs/length(Y) 分析解调器的抗噪性能时，

24、在输入端加入高斯白噪声，采用 MATLAB 函数 awgn 实现，其函数使用格式为：Y =AWGN(X,SNR)，加高斯白噪声于 X 中，SNR 为信噪比，单位为 dB，其值在假设 X 的功率为0dBM 的情况下确定。 信号的信噪比为信号中有用的信号功率与噪声功率的比值，根据信号功率定义，采用MATLAB 函数 var 实现，其函数常使用格式为：Y =VAR(X)，返回向量的方差，则信噪比为:SNR=VAR(X1)/VAR(X2)。绘制曲线采用 MATLAB 函数 plot 实现，其函数常使用格式：PLOT（X，Y） ，X 为横轴变量，Y 为纵轴变量，坐标范围限定 AXIS(x1 x2 y1

25、y2)，轴线说明 XLABEL 和 YLABEL。程序设计流程图见附录。3.2 MATLAB 仿真及其分析图 6 载波信号的波形及其频谱 图 7 调制信号的波形及其频谱110 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1-10-50510低低低低AM低低低低低低0 5000 10000 15000012345x 105低低低低AM低低低低低低0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 10510低低低低低低低低低低低低0 50 100 150 200 250 300 350 400 450 500012345x 105低低低低低低低

26、低低低低低低低0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1-10-50510t低低低低低低低低低低0 50 100 150 200 250 300 350 400 450 500012345x 105低低低低低低低低低低低低0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1-10-50510低低低低低低低1低低低低低低低0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1-10-50510低低低低低低低2低低低低低低低载波是指被调制以传输信号的波形，一般为正弦波。载波信号，就是把普通信号加载到一定频率的高频

27、信号上，在没有加载普通信号的高频信号时，高频信号的波幅是固定的，加载之后波幅就随着普通信号的变化而变化。本设计采用频率为 8000HZ，振幅为 10 的余弦载波信号，其波形及频谱图如图 6 所示。调制信号是指来自信源的消息信号（基带信号） ，这些信号可以是模拟的，也可以是数字的。本设计中产生的基带信号频率为 10HZ，振幅为 5 的余弦信号。该基带信号的波形及其频谱如图 7 所示。图 8 AM 调制信号的波形及其频谱 图 9 滤波前的解调信号的波形及其频谱调制就是由调制信号去控制高频载波的幅度，使之随调制信号作线性变化的过程，本设计的 AM 已调信号的波形及其频谱如图 8 所示，从图中可以看出

28、已调信号幅度随基带信号的幅度变化而变化，且把频率为 10HZ 的基带信号搬移到频率为 8000HZ 的载频上，实现了信道的多路复用，提高了信道利用率。由式（4）可知，其解调后的信号中含有高频及噪声成分，经低通滤波器取出低频分量，即可得到原始的基带调制信号。图 9 为经过相干解调后的滤波前的解调信号波形图及其频谱图。图 10 滤波后的解调信号的波形及其频谱 图 11 加入不同噪声的解调信号波形121 1.5 2 2.5 3 3.5 4 4.5 501002003004005006007008009001000低低低低/低低低低0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0

29、.9 100.20.40.60.811.21.4低低低低FIR低低低低低图 10 为经过低通滤波器，滤除了噪声及高频成分之后的解调信号的波形图及其频谱图，其波形与基带信号波形基本一致。将信道中存在的不需要的电信号统称为噪声，通信系统的噪声是叠加在信号上的，没有传输信号时通信系统中也有噪声，噪声永远存在于通信系统中。 图 11 为在信道中加入不同程度的高斯白噪声经过低通滤波器后输出的解调信号波形，可以看出其解调出的信号波形图发生明显的失真。图 12 输入/输出信噪比 图 13 低通滤波器的幅度函数图 12 为该系统采用相干解调时输入、输出信噪比的对比图。13总结在课程设计过程中，着重研究了 DS

30、B 信号调制与解调原理和 MATLAB 模拟实现，熟悉了信号波形、频谱的和系统性能的分析方法，了解了数字滤波器的设计与使用方法，综合提高了自己的专业技能。通过对程序的设计，我进一步熟悉了 MATLAB 开发环境，对 MATLAB 的一些具体操作和应用有了更深入的了解。如：有要求的正弦信号的产生，基本图形的绘制和各种的函数的使用等。同时，这次设计我对数字信号处理和通信原理课本上学到的知识点有了更深入的理解和掌握。比如对信号的调制和解调过程有了更深层的理解，学会了如何使用 MATLAB 对信号进行 SSB 调制和解调，了解了低通滤波器的 MATLAB 设计方法。还有很重要的一点是，我学会了如何安排

31、设计所需的时间及合理利用网络资源等普遍实用的学习方法，通过和同学探讨，拓宽了我的眼界。14参考文献1.李建新. 现代通信系统分析与仿真MATLAB 通信工具箱.西安：西安电子科技大学出版社，2000.2.樊昌信. 通信原理. 北京：国防工业出版社，2002.3.刘敏. MATLAB. 通信仿真与应用. 北京：国防工业出版社.4.曹志刚等著. 现代通信原理. 北京：清华大学出版社.2001.55.吴伟陵等著. 移动通信原理. 北京：电子工业出版社，200515附录总体设计结构图开始输入 DSB 调制信号加入白噪声显示加入白噪声前后的已调信号波形及频谱图对已调信号解调显示解调后，加入白噪声前后的解

32、调信号波形及频谱图显示输入输出信噪比关系图结束1617程序代码：载波信号t=-1:0.00001:1;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(载频信号波形);axis(0,0.01,-15,15); grid on;subplot(2,2,2);Y=fft(y); %对 y 的傅里叶变换plot(abs(Y); %频谱取模xlabel(频率);title

33、(载波信号频谱);axis(0,15000,0,1000000grid on;调制信号t=-1:0.00001:1;A0=10; %载波信号振幅A1=5； %调制信号振幅fs=8000; %载波信号频率fc=1000;ws=fs*pi;wc=fc*pi;18x=A1*cos(0.01*wc*t); %调制信号figure(2);subplot(2,1,1);plot(t,x);xlabel(t);title(调制信号 ); axis(0,1,-10,10); %定义坐标轴的范围grid on;subplot(2,1,2);X=fft(x); %对x 的傅里叶变换 plot(abs(X); %取

34、模值xlabel(频率);title(调制信号 频谱);axis(0,500,0,1000000);grid on; AM 已调信号t=-1:0.00001:1;A0=10;A1=5fs=8000; fc=1000;ws=fs*pi;wc=fc*pi;x=A1*cos(0.01*wc*t); %调制信号m=0.15; %调制度y1=3*(1+0.15*x).*cos(ws*t); %AM已调信号figure(3);subplot(2,1,1);plot(t,y1);grid on;xlabel(t);title(AM调制信号波形 );19axis(0,1,-10,10);subplot(2,1

35、,2);Y1=fft(y1); %傅里叶变换 plot(abs(Y1); %取模grid on;xlabel(频率);title(AM调制信号频谱 );axis(0,15000,0,500000); grid on;解调信号t=-1:0.00001:1;A0=10;A1=5fs=8000; fc=1000;ws=fs*pi;wc=fc*pi;x=A1*cos(0.01*wc*t); %调制信号m=0.15; %调制度y1=3*(1+0.15*x).*cos(ws*t); %AM 已调信号y11=y1.*cos(ws*t); %对AM已调信号进行解调figure(4);subplot(2,1,1

36、);plot(t,y11);grid on;xlabel(t);title(未滤 波的解调信号 );axis(0,1,0,10);subplot(2,1,2);Y11=fft(y11); %傅里叶变换 plot(abs(Y11); %取模grid on;20xlabel(频率);title(未滤 波的解调信号 频谱 );axis(0,500,0,500000); grid on;AM 解调信号 FIR 滤波t=-1:0.00001:1;A0=10;A1=5fs=8000; fc=1000;ws=fs*pi;wc=fc*pi;x=A1*cos(0.01*wc*t); %调制信号m=0.15; y

37、1=3*(1+0.15*x).*cos(ws*t); %AM 已调信号y11=y1.*cos(ws*t); %对AM已调信号进行解调Ft=2000; %采样频率fpts=100,120; %通带边界频率 fp=100Hz，阻带截止频率 fs=120Hz； mag=1 0;dev=0.01 0.05; %通带波动 1%，阻带波动 5% n21,wn21,beta,ftype=kaiserord(fpts,mag,dev,Ft); % 估计采用凯塞窗设计的 FIR滤波器的参数b21=fir1(n21,wn21,Kaiser(n21+1,beta); %由 fir1 设计滤波器 h,w=freqz(

38、b21,1); %得到频率响应 figure(3);plot(w/pi,abs(h);grid ontitle(FIR 低通滤波器 );figure(4);subplot(2,1,1);y111=fftfilt(b21,y11); %对 y11 低通滤波plot(t,y111);grid on;xlabel(t);title(滤波后的解调信号波形);axis(0,1,-10,10);subplot(2,1,2);Y111=fft(y111); %求 y111 的频谱 plot(abs(Y111); %对 Y111 取模grid on;xlabel(频率);title(滤波后的解调信号频谱);2

39、1axis(0,500,0,500000); grid on;加噪声t=-1:0.00001:1;A0=10;A1=5fs=8000; fc=1000;ws=fs*pi;wc=fc*pi;x=A1*cos(0.01*wc*t); %基带信号m=0.15; %调制度y1=3*(1+0.15*x).*cos(ws*t); %AM 调制信号y2=awgn(y1,5); %加入高斯白噪声y3=awgn(y1,10);y4=awgn(y1,15);y5=awgn(y1,20);y21=y2.*cos(ws*t); %解调y31=y3.*cos(ws*t);y41=y4.*cos(ws*t);y51=y5

40、.*cos(ws*t);y211=fftfilt(b21,y21); %低通滤波y311=fftfilt(b21,y31);y411=fftfilt(b21,y41);y511=fftfilt(b21,y51);figure(5);subplot(2,1,1);plot(t,y211);grid on;xlabel(t);title(加噪声 1 的解调 信号波形);axis(0,1,-10,10);subplot(2,1,2);plot(t,y311);grid on;xlabel(t);title(加噪声 2 的解调 信号波形);axis(0,1,-10,10);figure(6);subp

41、lot(2,1,1);plot(t,y411);grid on;xlabel(t);22title(加噪声 3 的解调 信号波形);axis(0,1,-10,10);subplot(2,1,2);plot(t,y511);grid on;xlabel(t);title(加噪声 4 的解调 信号波形);axis(0,1,-10,10);X1=y2-y1; %输入噪声 x1=var(y1)/var(X1); %输入信噪比Y1=y211-y111; %输出噪声y1=var(y111)/var(Y1); %输出信噪比X2=y3-y1; x2=var(y1)/var(X2); Y2=y311-y111; y2=var(y111)/var(Y2); X3=y4-y1; x3=var(y1)/var(X3); Y3=y411-y111; y3=var(y111)/var(Y3); X4=y5-y1; x4=var(y1)/var(X4);Y4=y511-y111; y4=var(y111)/var(Y4); figure(7);subplot(1,1,1);in=x1,x2,x3,x4; out=y1,y2,y3,y4;plot(in) hold on;plot(out);xlabel(t);ylabel(输入/出信噪比);axis(1,5,-10,1000);grid on;