文档线性幅度调制系统的仿真设计课程设计报告-毕设论文.doc

通信系统课程设计报告通信系统课程设计报告课题名称 线性幅度调制系统的仿真设计 学生姓名 班 级 学 号 指导教师 设计地点 13年 12月 1日目 录序言 1第 1 章MATLAB软件介绍 2第2章 AM调制解调 32.1 AM调制指标要求 32.2 整体设计 32.2.1 AM调制设计 32.2.2 AM解调设计 4 2.3 AM调制流程图 52.4 AM调制仿真结果及分析 6第3章 DSB调制解调 93.1 DSB调制指标要求 93.2 整体设计 93.2.1 DSB调制设计 93.2.2 DSB解调设计 103.3 DSB调制流程图 113.4 DSB调制仿真结果及分析12第4章 SSB调制解调 144.1 SSB 调制指标要求 144.2 整体设计 144.2.1 SSB调制设计 144.2.2 SSB解调设计 154.3 SSB调制流程图 164.4 SSB调制仿真结果及分析 17参考文献 19体会与建议 20附录 21序 言这次课程设计的重点就是模拟通信系统中的调制解调系统的基本原理以及仿真，并在MATLAB软件平台上的仿真实现几种常见的模拟调制方式。幅度调制是用调制信号去控制高频载波的振幅，使其按调制信号的规律变化的过程 ，常分为标准调幅（AM）、抑制载波双边带调制（DSB）、单边带调制（SSB）和残留边带调制（VSB）等。在通信系统中，信道的频段往往是有限的，而原始的通道信号的频段与信道要求的频段是不匹配的，这就要求将原始信号进行调制再进行发送。相应的在接收端对调制的信号进行调解，恢复原始的信号，而且调制解调还可以在一定的程度上抑制噪声对通信型号的干扰。调制解调技术按照通信信号是模拟的还是数字的可分为模拟调制调制和数字调制解调。但不论是模拟调制解调还是数字调制解调，都是通过将通信信号的信息加载到载波的幅度、频率或者相位上，因此调制解调可有分为幅度调制、频率调制和相位调制。通信系统有不同的分类的方法，根据是否采用调制，将通信系统分为基带传输和频带传输。基带传输是将未经频带调制的信号直接传送，调制的方式有很多。调制方式按照传输特性，调制方式可分为线性调制和非线性调制。广义的线性调制，是指已调波中被调参数随调 制信号成线性变化的调制过程。狭义的线性调制,是指把调制信号的频谱搬移到载波频率两侧而成为上、下边带的调制过程。此时只改变频谱中各分量的频率，但不改变各分量振幅的相对比例，使上边带的频谱结构与调制信号的频谱相同，下边带的频谱结构则是调制信号频谱的镜像。狭义的线性调制有调幅(AM)、抑制载波的双边带调制(DSB-SC)和单边带调制(SSB)、残留边带调制（VSB）。调制和解调在通信系统中是一个极为重要的组成部分，采用什么样的调制与解调方式将直接影响通信系统的性能。而模拟调制技术的原理还可以推广到数字调制中，因此我们有必要对模拟调制技术进行研究。 第 1 章MATLAB 软件介绍MATLAB是matrix & laboratory两个词的组合，意为矩阵工厂（矩阵实验室）。是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB和Mathematca、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C+，JAVA的支持。可以直接调用，用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。第 2 章 AM调制信号21AM调制指标要求（1）信源为fm=7Hz,Am=7V的余弦信号，载波fc=70Hz（2）根据线性幅度调制原理，确定调制系统设计方案（3）画出AM调制信号时域波形和频谱图（4）对数据结果进行分析2.2 整体设计 2.2.1 AM调制设计标准调幅就是常规双边带调制，简称调幅（AF）。假设调制信号m(t)的平均值为0，将其叠加一个直流偏量后与载波相乘(图2-1),即可形成调幅信号。 载波乘法器加法器m(t) S(t)图2-1 AM调制模型AM信号的时域和频域表达式分别为 SAM=(t)A0+m(t)coswc(t)=A0coswc(t)+m(t)coswc(t) (1) SAM()= A0(+C)+ (-)+1/2M（+）+M(-C) （2）式中，A0为外加的直流分量；m(t)可以是确知信号也可以是随机信号，但通常认为其平均值为0，即 m(t) =0 。AM信号的频谱SAM(t)是由载频分量和上、下两个边带组成（通常称频谱中画斜线的部分为上边带，不画斜线的部分为下边带）。上边带的频谱与原调制信号的频谱结构相同，下边带是上边带的镜像。显然，无论是上边带还是下边带，都含有原调制信号的完整信息。故AM信号是带有载波的双边带信号，它的带宽为基带信号带宽的两倍，即BAM=2Bm=2fH (3)式中，Bm=fH为调制信号m(t)的带宽，fH为调制信号的最高频率。2.2.2 AM解调设计（1） AM信号的解调原理及方式 解调是将位于载波的信号频谱再搬回来，并且不失真的恢复出原始基带信号。解调的方式有两种：相干解调和非相干解调。相干解调适用于各种线性调制系统，非线性解调一般适用于幅度调制（AM）信号。（2）相干解调相干解调也叫同步检波。解调与调制的实质一样，均是频谱搬移。解调是把在载频位置的以调信号的谱搬回到原始基带位置。由AM信号的频谱可知，如果将已调信号的频谱搬回到原点位置，即可得到原始的调制信号频谱，从而恢复出原始信号。解调中的频谱搬移同样可用调制时的相乘运算来实现。 载波低通滤波器乘法器S（t） m(t)图2-2 相干解调原理图SAM(t)coswct=A0+m(t)cos2wct=1/2A0+m(t)+1/2A0+m(t)cos2wct(1)2.3AM调制流程图 开始 为各个变量赋初值生成源信号m(t)对m(t)进行AM调制计算带宽画出AM信号波形画出AM相干解调波形画出AM频谱图 结束图2-3 AM调制流程图2.4 AM调制仿真及分析 仿真程序close all;clear all;dt=0.0001; %时间采样间隔fm=7; %信源最高频率fc=70; %载波中心频率T=1 %信号时长N=T/dt;t=0:N-1*dt;mt=7*cos(2*pi*fm*t); %信源，调制信号m(t)表达式%AM modulationA=7; %不失真条件，A要大于m(t)的最大值s_am=(A+mt).*cos(2*pi*fc*t); %AM时域表达式B=2*fm; %带通滤波器带宽figure(1) %表示显示在一张图片上subplot(311) %3表示上下显示三个，1表示一行显示一个，1表示显示的位置plot(t,s_am);hold on; %画出AM信号波形 ，hold on 继续 plot(t,A+mt,r-); %标出AM的包络title(AM调制信号)xlabel(t); %横坐标显示tgrid on%AM demodulationrt=s_am.*cos(2*pi*fc*t); %相干解调 表达式rt=rt-mean(rt);f,rf=FFT_SHIFT(t,rt);t,rt=RECT_LPF(f,rf,2*fm); %低通滤波subplot(312)plot(t,rt);hold on;title(AM相干解调后的信号波形)xlabel(t)grid onsubplot(313)f,sf=FFT_SHIFT(t,s_am); %调制信号频谱plot(f,sf);hold on; %画出频谱图形axis(-2*fc 2*fc 0 1.5*max(sf); title(AM信号频谱)xlabel(f); %横坐标显示fgrid on 图2-4AM仿真图形（不过载）A=7。由图2-4可以看出AMD 的调制信号幅值的为-20到+20，周期为0.13图2-5AM仿真图形（过载）A=1由图2-5可以看出AMD 的调制信号的幅值为-10到+10周期为0.13（1）当未满足|m(t)|maxA0的条件时，就会出现“过调幅”现象，这时包络检波将会发生失真。本题过调幅时的参数（A0=10，|m(t)|=7）表达式为f(t)=(8+7)coswct。（2）可以看出AM调制占用频带较宽，已调信号的频带宽度是调制信号的频带的两倍。（3）AM信号是带有载波的双边带信号，它的带宽为基带信号带宽的两倍，所以信道利用率不高，传输效率很低。（4）由图2-4看出带宽为BAM=2Bm=2fH=14hz，与理论值一致。第3章 DSB调制信号3.1 DSB调制指标要求（1）信源为fm=7Hz,Am=7V的余弦信号，载波fc=70Hz（2）根据线性幅度调制原理，确定调制系统设计方案（3）画出DSB调制解调信号时域波形和频谱图（4）对数据结果进行分析3.2 整体设计3.2.1 DSB调制设计在AM信号中，载波分量并不携带信息，信息完全由边带传送。如果将载波抑制，只需在将直流去掉，即可输出抑制载波双边带信号，简称双边带信号（DSB）。DSB调制器模型如图3-1所示。m(t)SDSB(t)cosc(t)图3-1 DSB调制模型 其中，设正弦载波为c(t)=Acos(ct+0) (1)式中，为载波幅度；为载波角频率；为初始相位（假定0为0）。调制过程是一个频谱搬移的过程，它是将低频信号的频谱搬移到载频位置。而解调是将位于载频的信号频谱再搬回来，并且不失真地恢复出原始基带信号。3.2.2 DSB解调设计双边带解调通常采用相干解调的方式，它使用一个同步解调器，即由相乘器和低通滤波器组成。在解调过程中，输入信号和噪声可以分别单独解调。相干解调的原理框图如图2-2所示 c(t)=Acos(ct+0) (2)乘法器低通滤波器 (t)Sm(t)m0(t)图3-2 相干解调器原理图信号传输信道为高斯白噪声信道，其功率为2。DSB信号只能运用相干解调，SDSB(t).coswct=m(t)cos2wct=1/2m(t)+1/2m(t)cos2wct (3)经低通滤波器滤除高次项，得M0(t)=1/2m(t) (4) 无失真的恢复原始电信号。3.3 DSB调制流程图 开始 为各个变量赋初值生成源信号m(t)对m(t)进行DSB调制计算带宽画出DSB信号波形画出DSB相干解调波形画出DSB频谱图 结束图3-3 DSB调制流程图3.4 DSB调制仿真及分析仿真程序close all;clear all;dt=0.0001; %时间采样间隔fm=7; %信源最高频率fc=70; %载波中心频率T=1 %信号时长N=T/dt;t=0:N-1*dt;mt=7*cos(2*pi*fm*t); %信源，调制信号m(t)表达式%DSB modulations_dsb=mt.*cos(2*pi*fc*t); %DSB时域表达式B=2*fm; %带通滤波器带宽figure(1) %表示显示在一张图片上subplot(312) %3表示上下显示三个，1表示一行显示一个，2表示显示的位置plot(t,s_dsb);hold on; %画出DSB信号波形 ，hold on 继续 plot(t,mt,r-); %标出mt的波形 title(DSB调制信号)xlabel(t); %横坐标显示tgrid on%DSB demodulationrt=s_dsb.*cos(2*pi*fc*t); %相干解调 表达式rt=rt-mean(rt);f,rf=FFT_SHIFT(t,rt);t,rt=RECT_LPF(f,rf,2*fm); %低通滤波subplot(312)plot(t,rt); hold on;title(DSB相干解调后的信号波形)xlabel(t);grid onsubplot(313)f,sf=FFT_SHIFT(t,s_dsb); %调制信号频谱plot(f,sf);hold on; %画出频谱图形axis(-2*fc 2*fc 0 1.5*max(sf); title(DSB信号频谱)xlabel(f); %横坐标显示fgrid on图3-4 DSB仿真图形由图3-4可以看出，DSB的调制信号幅值为-10到+10周期为0.13，有上边带和下边带（1）DSB的频谱不含载频，只有两个边带，仍是线性频谱搬移（频谱线性搬移即将调制信号频谱线性搬移至载频附近）。（2）DSB的带宽和AM的带宽相同，所以频带利用率都比较低。（3）与AM信号频谱相比，因为不存在载波分量，DSB信号的频谱少掉了载频部分的频谱；（4）DSB信号所需的传输带宽仍是调制信号的2倍。（5）由图3-4看出带宽为BAM=2Bm=2fH=14hz，与理论值一致。第4章 SSB调制信号4.1 SSB调制指标要求（1）信源为fm=7Hz,Am=7V的余弦信号，载波fc=70Hz（2）根据线性幅度调制原理，确定调制系统设计方案（3）画出SSB调制信号时域波形和频谱图（4）对数据结果进行分析4.2 整体设计4.2.1 SSB调制设计用滤波法实现单边带调制的原理图如图4-1所示，图中的HSSB()为单边带滤波器。产生SSB信号最直观方法的是，将HSSB()设计成具有理想高通特性HH()或理想低通特性HL()的单边带滤波器，从而只让所需的一个边带通过，而滤除另一个边带。产生上边带信号时HSSB()即为HH()，产生下边带信号时HSSB()即为HL()。 载波 Hssb(w)乘法器 Sdsb(t)m(t)Sssb(t)图4-1 SSB信号调制模型显然，SSB信号的频谱可表示为SSSB(w)=SDSB(w)HSSB(w)=1/2M(w+wc)+M(w-wc)HSSB(w) (1)用滤波法实现SSB信号，原理框图简洁、直观，但存在的一个重要问题是单边带滤波器不易制作。这是因为，理想特性的滤波器是不可能做到的，实际滤波器从通带到阻带总有一个过渡带。滤波器的实现难度与过渡带相对于载频的归一化值有关，过渡带的归一化值愈小，分割上、下边带就愈难实现。而一般调制信号都具有丰富的低频成分，经过调制后得到的DSB信号的上、下边带之间的间隔很窄，要想通过一个边带而滤除另一个，要求单边带滤波器在附近具有陡峭的截止特性即很小的过渡带，这就使得滤波器的设计与制作很困难，有时甚至难以实现。为此，实际中往往采用多级调制的办法，目的在于降低每一级的过渡带归一化值，减小实现难度。4.2.2 SSB解调设计从SSB信号调制原理图中不难看出，SSB信号的包络不再与调制信号成正比，因此SSB信号的解调也不能采用简单的包络检波，需采用相干解调低通滤波器乘法器 载波 Sp(t)Sssb(t)m(t)图4-2 SSB信号相干解调原理图乘法器输出 Sp(t)=Sssb(t)*cosct=1/2m(t)cosctm(t)sinct cosct=1/2m(t) cos2ct1/2m(t)cosct sinct =1/4m(t)+1/4m(t) cos2ct1/4m(t)sin2ct(2)经低通滤波后的解调输出为m0(t)=1/4m(t)(3)因而可得到无失真的调制信号。4.3 SSB调制流程图 开始 为各个变量赋初值生成源信号m(t)对m(t)进行DSB调制计算带宽画出DSB信号波形画出DSB相干解调波形画出DSB频谱图 结束图4-3 SSB调制流程图4.4 SSB调制仿真及分析仿真程序close all;clear all;dt=0.0001; %时间采样间隔fm=7; %信源最高频率fc=70; %载波中心频率T=1 %信号时长N=T/dt;t=0:N-1*dt;mt=7*cos(2*pi*fm*t); %信源，调制信号m(t)表达式%SSB modulation 调幅s_ssb=real(hilbert(mt).*exp(j*2*pi*fc*t); %SSB时域表达式B=fm; %带通滤波器带宽figure(1) %表示显示在一张图片上subplot(313) %3表示上下显示三个，1表示一行显示一个，3表示显示的位置plot(t,s_ssb);hold on; %画出SSB信号波形 ，hold on 继续 plot(t,mt,r-); %标出mt的包络title(SSB调制信号)xlabel(t);grid on%SSB demodulation 解调，检波rt=s_ssb.*cos(2*pi*fc*t); %相干解调f,rf=FFT_SHIFT(t,rt);t,nr=RECT_LPF(f,rf,2*fm); %低通滤波subplot(312)plot(t,nr); hold on;title(SSB相干解调后的信号波形)xlabel(t);grid onsubplot(313)f,sf=FFT_SHIFT(t,s_ssb); %单边带信号频谱plot(f,sf);hold on; %画出频谱图形axis(-2*fc 2*fc 0 1.5*max(sf); title(SSB信号频谱)xlabel(f);grid on图4-4 SSB仿真图形由图4-4可以看出，SSB的调制信号幅值为-10到+10，周期为0.13，只有一个下边带（1）由频谱仿真结果可以看出，SSB信号的实现比AM，DSB要复杂，SSB调制方式在传输时，传输的是一个边带，节省发射功率；（2）SSB的传输带宽比AM、DSB压缩一半，频带利用率高。（3）和DSB频谱图相比，SSB的频谱图出现的只有一个边带，但仍是频谱线性搬移。 （4）由图4-4看出带宽为BAM=2Bm=2fH=14hz，与理论值一致。参考文献1 樊昌信.曹丽娜等.通信原理（第六版）.北京.国防工业出版社。2 王兴亮.通信系统原理教程【M】.西安电子科技大学出版。3 刘颖. 数字通信原理与技术 北京.北京邮电大学出版社，20024 徐金明.MATLAB应用基础.北京.北京交通大学出版社。5 赵鸿图、茅艳等.通信原理MATLAB仿真教程【M】.北京.人民邮电出版社。6 汪浩 软件无线电调制解调系统的仿真与实现 航空电子技术 2005（36）57-59体会与建议通过这次的课程设计，我更进一步地了解了AM、DSB和SSB的调制与相干解调原理，特别是对其调制解调电路用Matlab来实现并进行分析。这次的课程设计让我们在原来的理论基础上加以实践，让我们更深刻地了解了模拟调制的过程，这对我们以后写毕业论文和找工作方面都有巨大的益处。依照以前所学MTLAB的书上介绍的MATLAB软件，实验准备之时我对MATLAB的相关知识及其应用也一定的了解，回顾了一下MATLAB编程语言。然后根据线性调制幅度原理确定设计了方案。通过上网查找与相关线性调制幅度相关的实验程序，结合本次试验的任务一点一滴的开始编写最终的程序。程序写完后发现图形总是出现失真，通过询问同学，知道了问题所在并加以改正，从而达到了各类调制解调系统的仿真实现。经过这次设计中我接触了很多新的知识，扩展了我的知识面，更加锻炼了我的动手能力，使我受益匪浅。然而，验收的时候老师提出的问题我却没能很好的回答出来，这主要归结于自己平时的疏忽与不努力，在验收前不能把所有的问题都想的更全面，尤其对于原理部分更是模糊不清。所以这次的验收让我非常深刻地反省了我自己，虽然没有回答出问题，但是认识到了自身的不足，以后我会努力改正这些缺点！这次的课程设计对我来说是一次深刻的教训，以后我会更加认真地去完成这珍贵的课程设计，更好的把理论运用于实践，争取把问题想得更深层次，尽力地发现缺点并完善自我。最后，感谢指导老师对我此次课程设计提出的建议与批评，我相信经过这次设计可以学到很多知识并且认识到自己的不足，对于以后的毕业论文以及工作会有很大的帮助。附录调用函数function f,sf=FFT_SHIFT(t,st)T=t(end);df=1/T; %频率间隔N=length(t); f=-N/2:N/2-1*df;sf=fft(st); %傅里叶变化sf=T/N*fftshift(sf); %将零频分量移至频谱的中心functiont,st=IFFT_SHIFT(f,Sf)df=f(2)-f(1);fmax=(f(end)-f(1)+df);dt=1/fmax;N=length(f);t=0:N-1*dt;Sf=fftshift(Sf);st=fmax*ifft(Sf);st=real(st); %求实部function t,st=RECT_LPF(f,Sf,B)df=f(2)-f(1); %频率间隔fN=length(f);RectH=zeros(1,fN); %矩阵BN=floor(B/df);%将矩阵中的内容从小到大排列BN_SHIFT=-BN:BN-1+floor(fN/2);RectH(BN_SHIFT)=1;Yf=RectH.*Sf;t,st=IFFT_SHIFT(f,Yf);第 21 页 共 23 页