模拟调制技术性能比较综述

1、课题名称 ： 通信系统仿真设计 课 题： 模拟调制技术性能比较 院 （系）： 专 业： 班 级： 学生姓名： 学 号： 指导教师： 年 12 月 8 日 摘要 调制，在通信系统中的作用至关重要。所谓调制，就是把信号转换成适合在信道中传 输的形式的一种过程。而模拟调制技术是非常重要的一种调制方式。本文主要运用 MATLAB 软件，通过写程序完成模拟调制技术(本文主要侧重线性调制)的仿真，并对仿真结果进 行分析，比较模拟调制技术性能。 关键字：线性调制 AM调制 DSB 调制 SSB调制 MATLAB仿真 、通信系统仿真设计的目的和要求 1、目的：以 MATLAB作为编程和通信仿真的工具， 通过仿

2、真， 加深对通信原理知识的理 解，同时，掌握利用 MATLAB来进行通信仿真设计的基本方法。 2、具体要求如下：模拟调制技术性能比较的仿真设计 给定单音信号 m(t) 2cos2 t ，载波 s(t) A cos 20 t ；请按原理及给定条件和要 求，借助 MATLAB设计一模拟通信系统： 当采用 AM调制时， A取 2； 当采用 DSB和SSB调制时， A取 1； 解调时，加入经接收端带通滤波后，功率为 0.1 的窄带高斯白噪声。 AM解调时要用相干解调与包络检波； 包络检波时：取 R=50，C分别为 2F，20F和0.2 F；电路结构见下图。 按要求仿真并画出各种调制及解调后的信号波形

3、图； 分析比较各种技术的性能； 调用的子程序必须能运行并与主程序能有效连接。 、幅度调制(线性调制)的原理 幅度调制是用调制信号去控制高频正弦载波的幅度，使其按调制信号的规律变化的过 程。幅度调制器的一般模型如图 1所示。 图中， 为调制信号， 为已调信号， 为滤波器的 冲激响应，则已调信号的时域和频域一般表达式分别为 式中， 为调制信号 的频谱， 为载波角频率。 由以上表达式可见，对于幅度调制信号，在波形上，它的幅度随基带信号规律而变化； 在频谱结构上， 它的频谱完全是基带信号频谱在频域内的简单搬移。 由于这种搬移是线性的， 因此幅度调制通常又称为线性调制，相应地，幅度调制系统也称为线性调制

4、系统。 在图 1的一般模型中， 适当选择滤波器的特性，便可得到各种幅度调制信号， 例如： 常规双边带调幅（ AM）、抑制载波双边带调幅（ DSB-SC）、单边带调制（ SSB）和残留边带调 制（ VSB）信号等。 1）常规双边带调幅（ AM） 1. AM 信号的表达式、频谱及带宽 AM 调制器模型如图 2所示。假设滤波器为全通网络 （ 1）， 调制信号 叠加直流 后再与载波相乘，则输出的信号就是常规 双边带调幅（ AM ）信号。 AM 信号的时域和频域表示式分别为 式中， 为外加的直流分量； 可以是确知信号也可以是随机信号，但通常认为其平均值 为0，即。 AM 信号的典型波形和频谱分别如图 3

5、所 示，图中假定调制信号 的上限频率为 。 显然，调制 信号 的带宽为 。 由图可知， AM 信号波形的包络与输 入基带信号 成正比，故用包络检波的 方法很容易恢复原始调制信号。 但为了 保证包络检波时不发生失真，必须满足 ，否则将出现过调幅现象而 带来失真。 由频谱可以看出， AM 信号的频谱 是由载频分量和上、下两个边带组成（通常称 频谱中画斜线的部分为上边带，不画斜线的部分为下边带） 。上边带的频谱与原调制信号的 频谱结构相同，下边带是上边带的镜像。因此， AM 信号是带有载波的双边带信号，它的带 宽为基带信号带宽的两倍，即 式中， 为调制信号 的带宽， 为调制信号的最高频率。 2. A

6、M 信号的解调 调制过程的逆过程叫做解调。 AM 信号的解调是把接收到的已调信号 还原为调制 信号 。 AM 信号的解调方法有两种：相干解调和包络检波解调。 （1）相干解调 由 AM 信号的频谱可知，如果将已调信号的频谱搬回到 原点位置， 即可得到原始的调制信号频谱， 从而恢复出原始信 号。解调中的频谱搬移同样可用调制时的相乘运算来实现。 相 干解调的原理框图如图 4 所示。将已调信号乘上一个与调制器 同频同相的载波，得 由上式可知， 只要用一个低通滤波器， 就可以将第 1项与第 2项分离， 无失真的恢复出原始 的调制信号 想干解调的关键是必须产生一个与调制器同频同相位的的载波。 如果同频同相

7、位条件得不到 满足，则会破坏原始信号的回复。 （2）包络检波法 由 的波形可见， AM 信号波形的包络与输入基带信号 成正比，故可以用包络 检波的方法恢复原始调制信号。 包络检波器一般由半波或全波整流器和低通滤波器组成。 D 、电阻R和电容 C组成。当RC满足条件 个二极管峰值包络检波器如图 5所示，电路由二极管 时，包络检波器的输出与输入信号的包络十分相近， 即 隔去直流后即可得到原信号。 包络检波法属于非相干解调法， 其特点是： 解调效率高， 解调器输出近似为相干解调的 2倍；解调电路简单，特别是接收端不需要与发送端同频同相位的载波信号，大大降低实现 难度。故几乎所有的调幅（ AM ）式接

8、收机都采用这种电路。 综上所述， 可以看出， 采用常规双边带幅度调制传输信息的好处是解调电路简单， 可采 用包络检波法。 缺点是调制效率低， 载波分量不携带信息，但却占据了大部分功率，白白浪 费掉。如果抑制载波分量的传送，则可演变出另一种调制方式，即抑制载波的双边带调幅 （DSB-SC ）。 2）抑制载波的双边带调幅（ DSB） 1. DSB信号的表达式、频谱及带宽 在幅度调制的一般模型中， 若假设滤波器为全通网络 1），调制信号中无直 或称抑制载波双边带 （DSB ） 流分量，则输出的已调信号就是无载波分量的双边带调制信号 调制信号，简称双边带（ DSB ）信号。 DSB 调制器模型如图 6

9、所示。可见 DSB 信号实质上就是基带 信号与载波直接相乘，其时域和频域表示式分别为 DSB 信号的包络不再与 的变化规律一致， 因此不能进行包络检波， 需采用相干解调； 除不再含有载频分量离散谱外， DSB信号的频谱与 AM 信号的完全相同，仍由上下对称的两 个边带组成。故 DSB信号是不带载波的双边带信号，它的带宽与AM 信号相同，也为基带信 号带宽的两倍， 即 式中， 为调制信号带宽， 为调制信号的最高频率。 2. DSB信号的功率分配及调制效率 由于不再含有载波成分， DSB 信号的功率就等于边带功率，是调制信号功率 的一半，即 式中， 为边带功率， 为调制信号功率。 显然，DSB 信

10、号的调制效率为 100%。 3. DSB信号的解调 DSB信号只能采用相干解调，其模型与 AM 信号相干解调时完全相同，如图 4所示。此 时，乘法器输出 经低通滤波器滤除高次项，得 即无失真地恢复出原始电信号。 抑制载波的双边带幅度调制的好处是， 节省了载波发射功率， 调制效率高； 调制电路简 单，仅用一个乘法器就可实现。缺点是占用频带宽度比较宽，为基带信号的2倍。 3）单边带调制（ SSB） 由于 DSB信号的上、下两个边带是完全对称的，皆携带了调制信号的全部信息，因此， 从信息传输的角度来考虑，仅传输其中一个边带就够了。这就又演变出另一种新的调制方 式 单边带调制（ SSB）。 1. SS

11、B信号的产生 法 中 产生 SSB 信号的方法很多，其中最基本的方 有滤波法和相移法。 （1）用滤波法形成 SSB 信号 用滤波法实现单边带调制的原理图如图 7所示， 图 的 为单边带滤波器。产生 SSB 信号最直观方法的是，将 设计成具有理想高 通特性 或理想低通特性 的单边带滤波器， 从而只让所需的一个边带通过， 而滤 除另一个边带。产生上边带信号时 即为 ，产生下边带信号时 即为 。 SSB 信号的频谱可表示为 用滤波法形成 SSB 信号，原理框图简洁、直观，但存在的 一个重要问题即单边带滤波器不易制作。理想特性的滤波器是 不可能做到的，实际滤波器从通带到阻带总有一个过渡带。滤 波器的实

12、现难度与过渡带相对于载频的归一化值有关，过渡带 的归一化值愈小，分割上、下边带就愈难实现。而一般调制信 号都具有丰富的低频成分，经过调制后得到的 DSB 信号的上、 下边带之间的间隔很窄，要想通过一个边带而滤除另一个，要 求单边带滤波器在 附近具有陡峭的截止特性（即很小的过渡 带），这就使得滤波器的设计与制作很困难，有时甚至难以实现。为此，实际中往往采用多 级调制的办法，目的在于降低每一级的过渡带归一化值，减小实现难度。 2）用相移法形成 SSB 信号 SSB 信号的时域表示式为 式中， “ ”对应上边带信号， “+”对应下边带信号；表示把 的所有频率成分均相 移 ，称 是 的希尔伯特变换。

13、根据上式可得到用相移法形成 SSB 信号的一般模型，如图 8所示。图中，为希尔 伯特滤波器，它实质上是一个宽带相移网络，对 中的任意频率分量均相移 。 该网络要对调制信号的所有频率 相移法形成 SSB 信号的困难在于宽带相移网络的制作， 分量严格相移 ，这一点即使近似达到也是困难的。 2. SSB信号的带宽、功率和调制效率 从 SSB 信号调制原理图中可以清楚地看出， SSB 信号的频谱是 DSB 信号频谱的一个边 带，其带宽为 DSB 信号的一半，与基带信号带宽相同，即 式中， 为调制信号带宽， 为调制信号的最高频率。 由于仅包含一个边带，因此 SSB信号的功率为 DSB 信号的一半，即 显

14、然，因 SSB 信号不含有载波成分，单边带幅度调制的效率也为100% 。 3. SSB信号的解调 从 SSB 信号调制原理图中不难看出， SSB 信号的包络 不再与调制信号 成正比，因此 SSB 信号的解调也不能采 用简单的包络检波，需采用相干解调，如图9所示。 此时，乘法器输出 经低通滤波后的解调输出为 从而可得到无失真的调制信号。 综上所述， 单边带幅度调制的好处是，节省了载波发射功率， 调制效率高； 频带宽度只有双 边带的一半，频带利用率提高一倍。缺点是单边带滤波器实现难度大。 三、MATLAB程序仿真与分析 1、AM 调制 给定单音信号 m(t) 2cos2 t ，载波 s(t) A

15、cos20 t ， A取2，解调时，加入经接 收端带通滤波后，功率为 检波时：取 R=50 ，C分别为 2F，20F和0.2 子程序见附录) 0.1 的窄带高斯白噪声， AM解调时要用相干解调与包络检波，包络 F。其 MATLAB程序见附录 包络检波 ，其仿真结果如图 10和图 11： 400 t/s 度 幅 B.单音信号频谱 D.加噪声后的AM已调信号频谱 300 200 100 500 400 300 度 幅 w/(rad*-s1) F.载波信号频谱 E.载波信号波形 200 100 t/s 度 幅 w/(rad*-s1) J.AM已调信号频谱 I.AM已调信号波形 1000 500 t/

16、s 度 幅 0 0 0.2 0.4 0.6 0.8 1 w/(rad*-s1) C.加噪声后的AM已调信号波形 1000 500 t/s 度 幅 0 0.2 0.4 0.6 0.8 1 1 3 w/(rad*-s1) 通过低通滤波器后的波形 0 0 通过低通滤波器后的频谱 0 0 0.2 0.4 0.6 0.8 1 2 t/s 度 幅 200 100 w/(rad*-s1) 度 幅 300 图 10 AM 调制与相干解调过程信号波形频谱图 从图10的A 图、E图和 I图可看出， 在波形上， 已调信号的幅度随基带信号的规律呈正 比地变化；由 B图、F图和 J 图可看出，在频谱结构上，它的频谱完全

17、是基带信号频谱在频 域内的简单搬移，若调制信号频率为 ，载波频率 c ，调制后信号频率搬移至 c 处； 从 D 图和 H图，可看出解调时的频谱亦是在频域内的简单搬移。因此，可验证，AM调制是 线性调制。 A. 包 络 检 波 后 的 AM 信 号 波 形 （ 其 中 R=50， C=2*10 -5） 从 I 图和 C图、J图和 D图可以得出加入噪声是小噪声信号，因为其波形和频谱几乎不 变。从 A图和 K 图可以看出，经调制解调的单音信号波形与原波形比较虽存在一定的失真， 但整体上还是基本满足原信号特征的。 10 度5 幅5 0 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 -1.

18、5 -1 -0.5 0 0.5 1 1.5 t/s C.包络检波 后的AM信号波 形（其 中R=50， C=2*10 -7） t/s B. 包 络 检 波 后 的 AM 信 号 波 形 ( 其 中 R=50， C=2*10 -6) 度 幅 10 5 0 -1.5 -1 -0.5 0.5 1.5 2 度 幅 10 5 0 -2 0 t/s D.包络检波 后的AM信号波 形(其 中R=50， C=6*10 -3) 0 -2 2 -1.5 -1 -0.5 0.5 1 1.5 度 幅 0 t/s 图 11 AM 包络检波 由前面的原理分析可知，包络检波时要合理选择 RC 的参数，使其满足关系式 ，防止

19、拖尾。对于给定的参数，经计算可发现不满足条件，其仿真结果如 图11的 A 图、 B 图、和 C 图所示，出现了严重失真。根据计算，选取一个比较符合条件的 参数，其仿真结果如图 11的 D 图所示。 2、DSB 调制 给定单音信号 m(t) 2cos2 t ，载波 s(t) A cos20 t ， A取1，解调时，加入经接 收端带通滤波后， 功率为 0.1 的窄带高斯白噪声。 其 MATLAB程序见附录， 其仿真结果如下： A.调制信号波形 B.调制信号频谱 C.加噪声后的 DSB已 调信号波 形 200 加噪声后的DSB已调信号 频谱 150 度 100 幅 50 0 1 2 3 t/s 2

20、1 0 -1 -2 0 t/s 度 幅 度 幅 I.DSB已调信号波形 w/(rad*s 1) J.DSB已调信 号频谱 500 400 300 200 100 w/(rad*s 1) 0 0 0.2 0.4 0.6 0.8 t/s G.加噪声后的 DSB解 调信号波 形 t/s 300 200 度 幅 100 0 0 1 23 w/(rad*s -1) H.加噪声 后的DSB解调 信号频谱 K.加噪 DSB解 调信号 度 幅 w/(rad*s 1) L.加噪DSB解调 信号 t/s w/(rad*s 1) 图 12 DSB 调制与解调过程波形和频谱图 由图 12 的频谱图可以看出 DSB的调

21、制与解调同 AM一样是线性的， DSB调制信号与 AM 调制信号的带宽一样。加入噪声后， 已调信号出现了相位和幅度的失真，相干解调后， 得到 了携带高频成分的信号， 之后通过低通滤波器，滤除了高频成分， 保留了低频成分，恢复了 原基带信号。从图 12的 A和 K可以得出 DSB调制信号的调制效率较 AM调制来说较高 2、SSB调制 给定单音信号 m(t) 2cos2 t ，载波 s(t) A cos20 t ， A取1，解调时，加入经接 收端带通滤波后， 功率为 0.1 的窄带高斯白噪声。 其 MATLAB程序见附录， 其仿真结果如下： A.调制信号波形 B.调 制信号 频 谱 w/(rad*

22、s 1) F.载波信号频谱 度 幅 t/s 度 幅 300 200 100 0 0.2 0.4 0.6 0.8 1 C.加 噪声 后的SSB已 调信 号波形 -2 -1 0 1 2 t/s D.加 噪 声后的SSB已 调信 号频谱 200 150 度幅 100 50 0123 E.载波 信号 波 形 t/s w/(rad*s -1) G.加 噪声 后的SSB解 调信 号波形 度 幅 0 w/(rad*s -1) I.SSB下边带调制 J.SSB下 边带调制信号 频 谱 t/s w/(rad*s -1) t/s 度 幅 t/s 加噪 声后的SSB解 调信号 w/(rad*s -1) 加噪 声 后

23、 的 SSB解 调 信 号 通过低通滤波器后的频谱 100 度幅 50 0123 w/(rad*s -1) 0 图13 SSB 调制与解调过程波形和频谱图 由图 13 可知，与 DSB调制信号相比，由于 SSB调制信号保留的是上边带，所以其带宽 只有 AM的一半。加入噪声后，出现了相位和幅度的失真，相干解调后，得到了携带高频成 分的信号，之后通过低通滤波器，滤除了高频成分，保留了低频成分，恢复了原基带信号。 总之，三种线性调制方法各有优缺点：AM调制的优点是接收设备简单，缺点是功率利 用率低，抗干扰能力差； DSB调制的优点是功率利用率高，且带宽与AM相同，但接收要求 同步解调，设备较复杂；

24、SSB调制的优点是功率利用率和频带利用率都较高，抗干扰能力优 于 AM，且带宽只有 AM的一半，缺点是发送和接收设备都较复杂。 四、小结 在这次基于 MATLAB的通信系统仿真设计中，主要是通过分析、查阅、掌握模拟调制技 术的相关原理，在 MATLAB环境下设计了对幅度调制( AM调制、 DSB调制和 SSB调制)的仿 真。我最大的收获是对 MATLAB软件的使用有了更深的了解， 尤其是懂得了新的函数， 同时， 提高了我查阅相关图书的能力，此外，对模拟调制技术原理有了更加深刻的认识。 五、附录 AM调制与解调程序 % 调AM制的参数设置 % clear all;clc; Am=sqrt(2);

25、 % 调制信号振幅 fm=1; % 调制信号频率 fc=10; % 载波信号频率 Ac=2; % 载波信号振幅 ts=0.01; t0=2; t=-t0+0.0001:ts:t0; N=length(t); w=0:2*pi/(N-1):2*pi; % 载波和调制信号的波形 % mt=Am*cos(2*pi*fm.*t); ct=Ac*cos(2*pi*fc.*t); figure(1); subplot(3,4,1);plot(t,mt); xlabel(t/s),ylabel( 幅度 ),title(A. 单音信号波形 ); % axis(-2,2,-10,10); subplot(3,4

26、,2);plot(w,real(fft(mt); xlabel(w/(rad*s-1),ylabel( 幅度 ),title(B. 单音信号频谱 ); axis(0,1,0,400); subplot(3,4,5);plot(t,ct); xlabel(t/s),ylabel( 幅度 ),title(E. 载波信号波形 ); axis(0,2,-3,3); subplot(3,4,6);plot(w,real(fft(ct); xlabel(w/(rad*s-1),ylabel( 幅度 ),title(F. 载波信号频谱 ); axis(0,1,0,500); % AM 信号调制程序 % A0

27、=2; % 直流偏量 amt=(A0+mt).*ct; subplot(3,4,9);plot(t,real(amt); xlabel(t/s),ylabel( 幅度 ),title(I.AM 已调信号波形 ); axis(-2,2,-10,10); subplot(3,4,10);plot(w,real(fft(amt); xlabel(w/(rad*s-1),ylabel( 幅度 ),title(J.AM 已调信号频谱 ); axis(0,1,0,1000); % 在信道中加入窄带高斯白噪声 % fangcha=0.1; h,l=size(mt); nit=sqrt(fangcha).*r

28、andn(h,l); 10 amtn=amt+nit; subplot(3,4,3);plot(t,real(amtn); xlabel(t/s),ylabel( 幅度 ); title(C. 加噪声后的 AM已调信号波形 ); subplot(3,4,4);plot(w,real(fft(amtn); xlabel(w/(rad*s-1),ylabel( 幅度 ); title(D. 加噪声后的 AM已调信号频谱 ); axis(0,pi,0,1000); %对 加噪声后的 AM已调信号进行“相干解调”% amtm=amtn.*ct; subplot(3,4,7);plot(t,real(a

29、mtm); xlabel(t/s),ylabel( 幅度 ), title(G. 加噪声后的 AM相干解调信号波形 ); subplot(3,4,8);plot(w,real(fft(amtm); xlabel(w/(rad*s-1),ylabel( 幅度 ); title(H. 加噪声后的 AM相干解调信号频谱 ); axis(0,pi,0,1000); % 解调信号通过低通滤波器 % Rp=0.1;Rs=50;% 设置通带波纹( dB)和阻带最小衰减 (dB) Wp=5/100;Ws=10/100;%设置通带截止频率和阻带起始角频率 n,Wn=ellipord(Wp,Ws,Rp,Rs);%

30、 椭圆滤波器的最小阶数 n 和截止频率 Wn b,a=ellip(n,Rp,Rs,Wn,low); X1=filtfilt(b,a,amtm); X=X1/2-A0; subplot(3,4,11);plot(t,real(X); axis(-2,2,-2,2); xlabel(t/s),ylabel( 幅度 ), title(K. 加噪声后的 AM解调信号 , 通过低通滤波器后的波形 ); subplot(3,4,12);plot(w,real(fft(X); xlabel(w/(rad*s-1),ylabel( 幅度 ); title(L. 加噪声后的 AM解调信号 , 通过低通滤波器后的

31、频谱 ); % “包络检波” % amtm1=BaoLuoJianBo(50,2*10-5,N,amtn); amtm2=BaoLuoJianBo(50,2*10-6,N,amtn); amtm3=BaoLuoJianBo(50,2*10-7,N,amtn); amtm4=BaoLuoJianBo(50,6*10-3,N,amtn); figure(2); subplot(4,1,1); plot(t,real(amtm1); xlabel(t/s),ylabel( 幅度 ); hold on;plot(t,A0+mt,r-); 11 title(A. 包络检波后的 subplot(4,1,2

32、); plot(t,real(amtm2); xlabel(t/s),ylabel( title(B. 包络检波后的 hold on;plot(t,A0+mt,r-); subplot(4,1,3); plot(t,real(amtm3); xlabel(t/s),ylabel( title(C. 包络检波后的 hold on;plot(t,A0+mt,r-); subplot(4,1,4); plot(t,real(amtm4); xlabel(t/s),ylabel( title(D. 包络检波后的 hold on;plot(t,A0+mt,r-); AM信号波形（其中 幅度 ）; AM信

33、号波形（其中 幅度 ）; AM信号波形（其中 幅度 ）; AM信号波形（其中 R=50，C=2*10-5 ) R=50，C=2*10-6 ) R=50，C=2*10-7 ) R=50，C=6*10-3 ) 包络检波子程序 function out=BaoLuoJB(R,C,N,st) w=-pi:10000*pi/(N-1):pi; for k=1:N if (st(k)0) st(k)=0; end end H=1./(1+1i.*w*R*C); out=conv(st,ifft(H); DSB调制与解调程序 % DS调B 制参数设置 % clear all;clc; Am=sqrt(2);

34、 % 调制信号振幅 fm=1; % 调制信号频率 fc=10; % 载波信号频率 Ac=1; % 载波信号振幅 ts=0.01; t0=2; t=-t0+0.0001:ts:t0; N=length(t); w=0:2*pi/(N-1):2*pi; % 载 波和调制信号波形 % mt=Am*cos(2*pi*fm.*t); ct=Ac*cos(2*pi*fc.*t); subplot(3,4,1);plot(t,mt); 12 xlabel(t/s),ylabel( 幅度 ),title(A. 调制信号波形 ); axis(-2,2,-2,2); subplot(3,4,2);plot(w,r

35、eal(fft(mt); xlabel(w/(rad*s-1),ylabel( 幅度 ),title(B. 调制信号频谱 ); axis(0,1,0,400); subplot(3,4,5);plot(t,ct); xlabel(t/s),ylabel( 幅度 ),title(E. 载波信号波形 ); axis(-2,2,-2,2); subplot(3,4,6);plot(w,real(fft(ct); xlabel(w/(rad*s-1),ylabel( 幅度 ),title(F. 载波信号频谱 ); axis(0,1,0,400); % DSB信号调制程序 % dsbt=mt.*ct;

36、subplot(3,4,9); plot(t,real(dsbt);hold on;plot(t,mt,r-); xlabel(t/s),ylabel( 幅度 ),title(I.DSB 已调信号波形 ); axis(-2,2,-2,2); subplot(3,4,10);plot(w,real(fft(dsbt); xlabel(w/(rad*s-1),ylabel( 幅度 ),title(J.DSB 已调信号频谱 ); axis(0,1,0,500); % 在 信道中加入窄带高斯白噪声% fangcha=0.1; h,l=size(mt); nit=sqrt(fangcha).*randn

37、(h,l); dsbtn=dsbt+nit; subplot(3,4,3);plot(t,real(dsbtn); xlabel(t/s),ylabel( 幅度 ); title(C. 加噪声后的 DSB已调信号波形 ); subplot(3,4,4);plot(w,real(fft(dsbtn); xlabel(w/(rad*s-1),ylabel( 幅度 ); title( 加噪声后的 DSB已调信号频谱 ); axis(0,pi,0,200); %对 加窄带高斯白噪声后的 DSB已调信号解调 % dsbtm=(dsbtn.*ct)*2; subplot(3,4,7);plot(t,rea

38、l(dsbtm); xlabel(t/s),ylabel( 幅度 ), title(G. 加噪声后的 DSB解调信号波形 ); subplot(3,4,8);plot(w,real(fft(dsbtm); xlabel(w/(rad*s-1),ylabel( 幅度 ); title(H. 加噪声后的 DSB解调信号频谱 ); axis(0,pi,0,300); %解 调信号通过低通滤波器 % 13 Rp=0.1;Rs=50;Wp=5/100;Ws=10/100; n,Wn=ellipord(Wp,Ws,Rp,Rs); b,a=ellip(n,Rp,Rs,Wn,low); X=filtfilt(

39、b,a,dsbtm); subplot(3,4,11);plot(t,real(X);axis(-2,2,-2,2); xlabel(t/s),ylabel( 幅度 ), title(K.加噪 DSB解调信号 , 通过低通滤波器的波形 ); subplot(3,4,12);plot(w,real(fft(X); xlabel(w/(rad*s-1),ylabel( 幅度 ); title(L.加噪 DSB解调信号 , 通过低通滤波器的频谱 ); axis(0,pi,0,100); SSB调制与解调程序 % SSB调制参数设置 % clear all;clc; Am=sqrt(2); % 调制信

40、号振幅 fm=1; % 调制信号频率 fc=10; % 载波信号频率 Ac=1; % 载波信号振幅 ts=0.01; t0=2; t=-t0+0.000001:ts:t0; N=length(t); w=0:2*pi/(N-1):2*pi; % 载 波和调制信号波形 % mt=Am*cos(2*pi*fm.*t); ct=Ac*cos(2*pi*fc.*t); subplot(3,4,1);plot(t,mt); xlabel(t/s),ylabel( 幅度 ),title(A. 调制信号波形 ); % axis(-2,2,-2,2); subplot(3,4,2);plot(w,real(fft(mt); xlabel(w/(rad*s-1),ylabel( 幅度 ),title(B. 调制信号频谱 ); axis(0,1,0,400); subplot(3,4,5);plot(t,ct); xlabel(t/s),ylabel( 幅度 ),title(E. 载波信号波形 ); % axis(-2,2,-2,2); subplot(3,4,6);plot(w,real(fft(ct); xlabel(w/(rad*s-1),ylabel( 幅度 ),title(F. 载