AM及SSB调制与解调详解.pdf

成绩：1通信原理课程设计通信原理课程设计设计题目：设计题目：AM 及及 SSB 调制与解调及抗噪声性能分析调制与解调及抗噪声性能分析班班 级：级：学生姓名：学生姓名：学生学号：学生学号：指导老师：指导老师：成绩：2目目 录录一、引言.31.1 概述.31.2 课程设计的目的.31.3 课程设计的要求.3二、AM 调制与解调及抗噪声性能分析.42.1 AM 调制与解调.42.1.1 AM 调制与解调原理.42.1.2 调试过程.62.2 相干解调的抗噪声性能分析.92.2.1 抗噪声性能分析原理.92.2.2 调试过程.10三、SSB 调制与解调及抗噪声性能分析.123.1 SSB 调制与解调原理.123.2 SSB 调制解调系统抗噪声性能分析.133.3 调试过程.15四、心得体会.19五、参考文献.19 成绩：3一、引言一、引言1.11.1 概述概述通信原理是通信工程专业的一门极为重要的专业基础课，但内容抽象，基本概念较多，是一门难度较大的课程，通过 MATLAB 仿真能让我们更清晰地理解它的原理，因此信号的调制与解调在通信系统中具有重要的作用。本课程设计是 AM 及 SSB 调制解调系统的设计与仿真，用于实现AM 及 SSB 信号的调制解调过程，并显示仿真结果，根据仿真显示结果分析所设计的系统性能。在课程设计中，幅度调制是用调制信号去控制高频载波的振幅，使其按调制信号的规律变化，其他参数不变。同时也是使高频载波的振幅载有传输信息的调制方式。1.21.2 课程设计的目的课程设计的目的在此次课程设计中，我需要通过多方搜集资料与分析：(1)掌握模拟系统 AM 和 SSB 调制与解调的原理；(2)来理解并掌握 AM 和 SSB 调制解调的具体过程和它在 MATLAB 中的实现方法；(3)掌握应用 MATLAB 分析系统时域、频域特性的方法，进一步锻炼应用 MATLAB 进行编程仿真的能力。通过这个课程设计，我将更清晰地了解 AM 和 SSB 的调制解调原理，同时加深对 MATLAB 这款通信原理辅助教学操作的熟练度。1.31.3 课程设计的要求课程设计的要求(1)熟悉 MATLAB 的使用方法，掌握 AM 信号的调制解调原理，以此为基础用 MATLAB 编程实现信号的调制解调；(2)设计实现 AM 调制与解调的模拟系统，给出系统的原理框图，对系统的主要参数进行设计说明；(3)采用 MATLAB 语言设计相关程序，实现系统的功能，要求采用一种方式进行仿真，即直接采用 MATLAB 语言编程的静态方式。要求采用两种以上调制信号源进行仿真，并记录各个输出点的波形和频谱图；(4)对系统功能进行综合测试，整理数据，撰写课程设计论文。成绩：4二、二、AMAM 调制与解调及抗噪声性能分析调制与解调及抗噪声性能分析2.12.1 AMAM 调制与解调调制与解调2.1.1 AM 调制与解调原理幅度调制是由调制信号去控制高频载波的幅度，使正弦载波的幅度随着调制信号而改变的调制方案，属于线性调制。AM信号的时域表示式：频谱：调制器模型如图所示：m t mstcosct0A图1-1 调制器模型AM的时域波形和频谱如图所示：成绩：5时域 频域 图1-2 调制时、频域波形AM信号的频谱由载频分量、上边带、下边带三部分组成。它的带宽是基带信号带宽的2倍。在波形上，调幅信号的幅度随基带信号的规律而呈正比地变化，在频谱结构上，它的频谱完全是基带信号频谱在频域内的简单搬移。所谓相干解调是为了从接受的已调信号中，不失真地恢复原调制信号，要求本地载波和接收信号的载波保证同频同相。相干载波的一般模型如下：将已调信号乘上一个与调制器同频同相的载波，得 twtmAtmAtwtmASccAM2cos)(21)(21cos)(tcoswt)(0020c 由上式可知，只要用一个低通滤波器，就可以将第1项与第2项分离，无失真的恢复出原始的调制信号 )(21)(00TMATM 相干解调的关键是必须产生一个与调制器同频同相位的载波。如果同频同相位的条件得不到满足，则会破坏原始信号的恢复。LPF mst pst dst coscc tt 成绩：62.1.2 调试过程t=-1:0.00001:1;%定义时长 A1=6;%调制信号振幅A2=10;%外加直流分量f=3000;%载波频率w0=2*f*pi;%角频率 Uc=cos(w0*t);%载波信号subplot(5,2,1);plot(t,Uc);%画载波信号title(载波信号);axis(0,0.01,-1,1);%坐标区间T1=fft(Uc);%傅里叶变换subplot(5,2,2);plot(abs(T1);%画出载波信号频谱title(载波信号频谱);axis(5800,6200,0,200000);%坐标区间mes=A1*cos(0.002*w0*t);%调制信号subplot(5,2,3);plot(t,mes);%画出调制信号 title(调制信号);T2=fft(mes);%傅里叶变换subplot(5,2,4);plot(abs(T2);%画出调制信号频谱title(调制信号频谱);axis(198000,202000,0,1000000);%坐标区间Uam1=A2*(1+mes/A2).*cos(w0).*t);%AM 已调信号subplot(5,2,5);plot(t,Uam1);%画出已调信号title(已调信号);T3=fft(Uam1);%已调信号傅里叶变换subplot(5,2,6);plot(abs(T3);%画出已调信号频谱title(已调信号频谱);成绩：7axis(5950,6050,0,900000);%坐标区间sn1=20;%信噪比db1=A12/(2*(10(sn1/10);%计算对应噪声方差n1=sqrt(db1)*randn(size(t);%生成高斯白噪声Uam=n1+Uam1;%叠加噪声后的已调信号 Dam=Uam.*cos(w0*t);%对AM已调信号进行解调subplot(5,2,7);plot(t,Dam);%滤波前的AM解调信号title(滤波前的AM解调信号波形);T4=fft(Dam);%求AM信号的频谱subplot(5,2,8);plot(abs(T4);%滤波前的AM解调信号频谱title(滤波前的AM解调信号频谱);axis(187960,188040,0,600000);Ft=2000;%采样频率fpts=100 120;%通带边界频率fp=100Hz 阻带截止频率fs=120Hzmag=1 0;dev=0.01 0.05;%通带波动1%，阻带波动5%n21,wn21,beta,ftype=kaiserord(fpts,mag,dev,Ft);%kaiserord估计采用凯塞窗设计的FIR滤波器的参数b21=fir1(n21,wn21,Kaiser(n21+1,beta);%由fir1设计滤波器z21=fftfilt(b21,Dam);%FIR低通滤波subplot(5,2,9);plot(t,z21,r);%滤波后的AM解调信号title(滤波后的AM解调信号波形);axis(0,1,-1,10);T5=fft(z21);%求AM信号的频谱subplot(5,2,10);plot(abs(T5),r);%画出滤波后的AM解调信号频谱title(滤波后的AM解调信号频谱);axis(198000,202000,0,500000);运行结果:成绩：8 成绩：92.22.2 相干解调的抗噪声性能分析相干解调的抗噪声性能分析2.2.1 抗噪声性能分析原理AM线性调制系统的相干解调模型如下图所示。图3.5.1 线性调制系统的相干解调模型图中可以是AM调幅信号，带通滤波器的带宽等于已调信号带宽10。下面讨论AM调制系)(mtS统的抗噪声性能11。AM信号的时域表达式为 tcosw)t(m)t(c0 ASAM 通过分析可得AM信号的平均功率为 2)t(m2)(220iASAM又已知输入功率,其中B表示已调信号的带宽。BN0in由此可得AM信号在解调器的输入信噪比为 HAMAMABANSfn4)t(mn2)t(m)(02200220iiAM信号经相干解调器的输出信号为 )t(m21)t(m0因此解调后输出信号功率为 )t(m41)t(m)(2200AMS在上图中输入噪声通过带通滤波器之后，变成窄带噪声，经乘法器相乘后的输出噪声为)t(nipicccscccccscn(t)n(t)cosw tn(t)cosw t-n(t)sinw tcosw t11n(t)n(t)cos2w t-n(t)sin2w t22 经LPF后，)t(n21)t(nc0因此解调器的输出噪声功率为 i2c20041)t(n41)t(nNN可得AM信号经过解调器后的输出信噪比为 成绩：10 HAMBNSfn2)t(mn)t(m)(020200由上面分析的解调器的输入、输出信噪比可得AM信号的信噪比增益为 )t(m)t(m22202ii00ANSNSGAM2.2.2 调试过程clf;%清除窗口中的图形t=0:0.01:2;%定义变量区间fc=50;%给出相干载波的频率A=10;%定义输入信号幅度fa=5;%定义调制信号频率mt=A*cos(2*pi*fa.*t);%输入调制信号表达式xzb=5;%输入小信躁比(dB)snr=10.(xzb/10);db=A2./(2*snr);%由信躁比求方差nit=sqrt(db).*randn(size(mt);%产生小信噪比高斯白躁声psmt=(A+mt).*cos(2*pi*fc.*t);%输出调制信号表达式psnt=psmt+nit;%输出叠加小信噪比已调信号波形xzb1=30;%输入大信躁比(dB)snr1=10.(xzb1/10);db1=A2./(2*snr1);%由信躁比求方差nit1=sqrt(db1).*randn(size(mt);%产生大信噪比高斯白躁声psnt1=psmt+nit1;%输出已调信号波形subplot(2,2,1);%划分画图区间plot(t,nit,g);%画出输入信号波形title(小信噪比高斯白躁声);xlabel(t);ylabel(nit);subplot(2,2,2);plot(t,psnt,b);title(叠加小信噪比已调信号波形);成绩：11xlabel(时间);ylabel(输出调制信号);subplot(2,2,3);plot(t,nit1,r);%length用于长度匹配title(大信噪比高斯白躁声);%画出输入信号与噪声叠加波形xlabel(t);ylabel(nit);subplot(2,2,4);plot(t,psnt1,k);title(叠加大信噪比已调信号波形);%画出输出信号波形xlabel(时间);ylabel(输出调制信号);运行结果：由上图可见，当输入信号一定时，随着噪声的加强，接收端输入信号被干扰得越严重。而相应的输出波形相对于发送端的波形误差也越大。而当噪声过大时信号几难分辨。这信噪比变小导致的在实际的信号传输过程中当信道噪声过大将会导致幅度相位等各失真当然由于非线性元件如滤波器等的存在。非线性失真也会随噪声加大而变大。成绩：12三、三、SSBSSB 调制与解调及抗噪声性能分析调制与解调及抗噪声性能分析 3.13.1 SSBSSB 调制与解调原理调制与解调原理单边带调制信号是将双边带信号中的一个边带滤掉而形成的。根据方法的不同，产生 SSB 信号的方法有：滤波法和相移法。由于滤波法在技术上比较难实现,所以在此我们将用相移法对 SSB 调制与解调系统进行讨论与设计。相移法和 SSB 信号的时域表示设单频调制信号为载波为则其双边带信号 DSB 信号的时域表示式为若保留上边带，则有若保留下边带，则有将上两式合并得：由希尔伯特变换故单边带信号经过希尔伯特变换后得:把上式推广到一般情况，则得到 式中tAtmmmcos)(ttcccos)(tAtAttAtsmcmmcmcmmDSB)cos(21)cos(21coscos)(1()cos()2USBmCmstAt11coscossinsin22mmcmmcAtAt11coscossinsin22mmcmmcAttAtt1()cos()2LSBmCmstAttAtAmmmmsinso c上上上上上上上上)()(tmtm上上上上上上上)()(Mtm 成绩：13若M()是m(t)的傅里叶变换，则上式中的-jsgn可以看作是希尔伯特滤波器传递函数，即移相法 SSB 调制器方框图相移法是利用相移网络，对载波和调制信号进行适当的相移，以便在合成过程中将其中的一个边带抵消而获得 SSB 信号。相移法不需要滤波器具有陡峭的截止特性，不论载频有多高，均可一次实现 SSB 调制。SSB 信号的解调 SSB 信号的解调不能采用简单的包络检波，因为 SSB 信号是抑制载波的已调信号，它的包络不能直接反映调制信号的变化，所以仍需采用相干解调。SSB 信号的性能SSB 信号的实现比 AM、DSB 要复杂，但 SSB 调制方式在传输信息时，不仅可节省发射功率，而且它所占用的频带宽度比 AM、DSB 减少了一半。它目前已成为短波通信中一种重要的调制方式。3.23.2 SSBSSB 调制解调系统抗噪声性能分析调制解调系统抗噪声性能分析 噪声功率这里，B=fH 为 SSB 信号的带通滤波器的带宽。信号功率SSB 信号sgn)()(jMMsgn)(/)()(jMMHho01144iNNn B 成绩：14与相干载波相乘后，再经低通滤波可得解调器输出信号因此，输出信号平均功率输入信号平均功率为单边带解调器的输入信噪比为单边带解调器的输出信噪比为制度增益因为在 SSB 系统中，信号和噪声有相同表示形式，所以相干解调过程中，信号和噪声中的正交分量均被抑制掉，故信噪比没有改善。o1()()4m tm t22oo1()()16Sm tm t()()m tm t上上上上上上上上上上上上上上上上上上上上上上上)(412tmSiBntmBntmNSii02024)()(4122oo001()()16144m tSm tNn Bn Boo/1/SSBiiSNGSN 成绩：153.33.3 调试过程调试过程先建立 3 个 M 文件1.afd_buttfunction b,a=afd_butt(Wp,Ws,Rp,As);if Wp=0 error(Passband edge must be larger than 0)endif Ws=Wp error(Stopband edge must be larger than Passband edge)endif(Rp=0)|(As t0=0.1;fs=12000;%t0 采样区间，fs 采样频率 fc=1000;Vm0=2.5;ma=0.25;%fc 载波频率，Vm0 输出载波电压振幅，ma 调幅度 n=-t0/2:1/fs:t0/2;%定义变量区间 N=length(n);A=4;%定义调制信号幅度 x1=A*cos(150*pi*n);%调制信号 x2=hilbert(x1,N);%对 x1 做希尔伯特变换 y=(Vm0*x1.*cos(2*pi*fc*n)-Vm0*x2.*sin(2*pi*fc*n)/2;%保留上边带的已调波信号 xzb=2;%输入小信噪比(dB)snr=10.(xzb/10);h,l=size(x1);%求调制信号的维度 fangcha=A*A./(2*snr);%由信噪比求方差 nit=sqrt(fangcha).*randn(h,l);%产生高斯白噪声 成绩：16 yn=y+nit;%叠加小信噪比噪声的已调波信号 xzb=10;%输入小信噪比(dB)snr1=10.(xzb/10);h,l=size(x1);%求调制信号的维度 fangcha=A*A./(2*snr1);%由信噪比求方差 nit1=sqrt(fangcha).*randn(h,l);%产生高斯白噪声 yn1=y+nit1;%叠加小信噪比噪声的已调波信号 figure(1)subplot(2,2,1)%划分画图区间 plot(n,x1)%画出调制信号的波形 title(调制信号);subplot(2,2,2)%划分画图区间 plot(n,y)%画出已调波信号波形 title(已调波信号);subplot(2,2,3)%划分画图区间 plot(n,yn)%画出叠加噪声的已调波信号波形 title(叠加小信噪比噪声的已调波信号);subplot(2,2,4)%划分画图区间 plot(n,yn1)%画出叠加噪声的已调波信号波形 title(叠加大信噪比噪声的已调波信号);X=fft(x1);%调制信号 x1 的傅里叶变换 Y=fft(y);%已调信号 y 的傅里叶变换 Yn=fft(yn);%叠加小信噪比噪声的已调信号 yn 的傅里叶变换 Yn1=fft(yn1);%叠加大信噪比噪声的已调信号 yn 的傅里叶变换 w=0:2*pi/(N-1):2*pi;%定义变量区间 figure(2)subplot(2,2,1)%划分画图区间 plot(w,abs(X)%画出调制信号频谱波形 axis(0,pi/4,0,3000);%给出横纵坐标的范围 title(调制信号频谱);subplot(2,2,2)%划分画图区间 plot(w,abs(Y)%画出已调波信号频谱 axis(pi/6,pi/4,0,2500);%给出横纵坐标的范围 title(已调波信号频谱);subplot(2,2,3)%划分画图区间 plot(w,abs(Yn)%画出叠加小信噪比噪声的已调波信号频谱 axis(pi/6,pi/4,0,2500);%给出横纵坐标的范围 title(叠加小信噪比噪声的已调波信号频谱);subplot(2,2,4)%划分画图区间 plot(w,abs(Yn1)%画出叠加大信噪比噪声的已调波信号频谱 axis(pi/6,pi/4,0,2500);%给出横纵坐标的范围 title(叠加大信噪比噪声的已调波信号频谱);y1=y-2*cos(1500*pi*n);y2=Vm0*y1.*cos(2*pi*fc*n);%将已调幅波信号的频谱搬移到原调制信号的频谱处 wp=40/N*pi;ws=60/N*pi;Rp=1;As=15;T=1;%滤波器参数设计 OmegaP=wp/T;OmegaS=ws/T;cs,ds=afd_butt(OmegaP,OmegaS,Rp,As);b,a=imp_invr(cs,ds,T);y=filter(b,a,y2);yn=y+nit;figure(3)成绩：17 subplot(2,1,1)%划分画图区间 plot(n,y)%画出解调波波形 title(解调波);Y=fft(y);%解调波 y 的傅里叶变换 subplot(2,1,2)%划分画图区间 plot(w,abs(Y)%画出解调信号频谱 axis(0,pi/6,0,2500);%给出横纵坐标的范围 title(解调信号频谱);运行结果:*Butterworth Filter Order=6-0.0500.05-4-2024信 信 信 信-0.0500.05-10-50510信 信 信 信 信-0.0500.05-20-1001020信 信 信 信 信 信 信 信 信 信 信 信 信 信-0.0500.05-10-50510信 信 信 信 信 信 信 信 信 信 信 信 信 信 成绩：1800.20.40.60100020003000信 信 信 信 信 信0.550.60.650.70.7505001000150020002500信 信 信 信 信 信 信0.550.60.650.70.7505001000150020002500信 信 信 信 信 信 信 信 信 信 信 信 信 信 信 信0.550.60.650.70.7505001000150020002500信 信 信 信 信 信 信 信 信 信 信 信 信 信 信 信-0.05-0.04-0.03-0.02-0.0100.010.020.030.040.05-10-50510信 信 信00.050.10.150.20.250.30.350.40.450.505001000150020002500信 信 信 信 信 信可以清晰地看出，加大噪声后，信号的波形杂乱无章，起伏远大于加小噪声时的波形。造成此现象的原因是当信噪比较小时，噪声的功率在解调信号中所占比重较大，所以会造成杂 成绩：19波较多的情况；而信噪比很大时，噪声的功率在解调信号中所占比重就很小了，噪声部分造成的杂乱波形相对就不是很明显，甚至可以忽略。综上所述，叠加噪声会造成解调信号的失真，信噪比越小，失真程度越大。所以当信噪比低于一定大小时，会给解调信号带来严重的失真，导致接收端无法正确地接收有用信号。所以在解调的实际应用中，应该尽量减少噪声的产生。四、心得体会四、心得体会这次的课程设计时间虽短但收获很多。我们用 MATLAB 进行了 AM 及 SSB 调制与解调的研究。不但又加深了课本的知识，而且也对 matlab 的基本知识有了一定掌握。本次课程设计中实现了通信基本知识与 MATLAB 的结合，并在实际中设计并仿真 AM 及 SSB 调制与解调的过程。这次课程设计中我们不得不对 AM 原理以及 SSB 原理进行更深一层次的理解，对书中原来学到的只知其果不懂其因的理论，在设计中也有了更深刻的认识。这次程序需要自己写而我们 MATLAB 的基础不是很好。这次课程设计虽然很简单，由于没有基础，查找了很多相关的资料，而且我们在编译和调试过程中除了很多次错误，这是过程，学习就是在过程中进行的，经过自己几天的劳动，再加上同学们和老师的帮助，不仅对读程序有了很大提高，而且自己的编译水平也上了一个新台阶，更加熟系了 MATLAB 的应用，也对其中的函数有了大概的了解总之这次课程设计收获很大，是一场难得的锻炼机会。五、参考文献五、参考文献1樊昌信，曹丽娜。通信原理（第六版）。国防工业出版社。2仲麟，王峰。MATLAB 仿真技术与应用教程M。北京：国防工业出版社。20043 陶亚雄，刘南平，王坚。现代通信原理M（第三版）。北京：电子工业出版社，2009.44 解相吾，解文博。通信电子电路M。北京：人民邮电出版社 2010.65 陈磊，唐晓辉。现代通信原理实验指导书。桂林:桂林航天工业高等专科学校，电子工程系。2009.7