xx 学 院
课程设计报告
课程名称: 设计题目: 学生班级: 学生: 指导教师: 完成日期:
通信原理
基于matlab的AM调制仿真 13通信(1) xxxxxxx 王小x 2015-12-24
数学与计算机学院
.下载可编辑.
. .. .. .
.下载可编辑.
. .. .. .
课程设计项目研究报告
目 录
第 1 章 项目简介 ........................................................ 1
1.1 项目名称 .......................................................... 1 1.2 开发人员 .......................................................... 1 1.3 指导教师 .......................................................... 1 第 2 章 项目研究意义 .................................................... 1
2.1 课程设计概述 ...................................................... 1 2.2 研究意义 .......................................................... 2 第 3 章 采用的技术 ...................................................... 2
3.1 课程设计的方案设计原理 ............................................ 2 第 4 章 课程设计项目进度表 .............................................. 4 第 5 章 课程设计任务分配表 .............................................. 4 第 6 章 达到的效果 ...................................................... 4
6.1 程序设计思想 ...................................................... 4 6.2 程序最终实现结果 .................................................. 5 第 7 章 源程序 ......................................................... 11
7.1 当调制信号是单一频率时 ........................................... 11 7.2 当调制信号不是单一频率时 ......................................... 13 第 9 章 设计心得 ....................................................... 15 第 10 章 参考文献 ...................................................... 16
.下载可编辑.
. .. .. .
第 1 章 项目简介
1.1 项目名称
基于matlab的AM调制仿真
1.2 开发人员
阮文添、清海
1.3 指导教师
王小文
第 2 章 项目研究意义
2.1 课程设计概述
正弦载波幅度随调制信号而变化的调制,简称调幅(AM)。调幅的技术和设备比较简单,频谱较窄,但抗干扰性能差,广泛应用于长中短波广播、小型无线、电报等电子设备中。早期的无线电报机采用火花式放电器产生高频振荡。传号时火花式发报机发射高频振荡波,空号时发报机没有输出。这种电报信号的载波不是纯正弦波,它含有很多谐波分量,会对其他信号产生严重干扰。理想的模拟正弦波调幅是:载波幅度与调制信号瞬时值成线性关系。
MATLAB是美国MathWorks公司出品的商业数学软件,用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境.MATLAB是matrix&laboratory两个词的组合,意为矩阵工厂(矩阵实验室)。是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境[1]。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式,代表了当今国际科学计算软件的先进水平。
.下载可编辑.
. .. .. .
2.2 研究意义
现在的社会越来越发达,科学技术不断地更新,在信号和模拟电路里面经常要用到调
制与解调,而AM的调制与解调是最基础的,也是经常用到的。用AM调制与解调可以在电路里面实现很多功能,制造出很多有用又实惠的电子产品,为我们的生活带来便利。 造我们日常生活中用的收音机也是采用了AM调制方式,而且在军事和民用领域都有十分重要的研究课题。
第 3 章 采用的技术
3.1 课程设计的方案设计原理
信号通过一定的传输介质在发射机和接收机之间进行传送时,信号的原始形式一般不适合传输。因此,必须转化他们的形式。将低频信号加到高频载波的过程,或者说把信息加载到信息载体上以便阐述的处理过程,伟调制。所谓“加载”,其实质是使高频载波信号(信息载体)的某个特性参数随信息信号的大小呈现性变化的过程。通常称代表信息的信号为调制信号,称信息载体信号为载波信号,称调制后的频带信号为已调波信号。在多种调制中,最先应用的一种就是标准振幅调制(AM)。标准振幅调制是一种相对便宜的,质量不高的调制形式。在频域中已调波频谱是基带调制信号频谱的线性位移;在时域中,已调波包络与调制信号波形呈线性关系[2]。
对于单频信号的调制情况,如果设单品调制信号为u0=Umcos(w0t),载波为uc=Uccos(wct),那么调幅信号(已调波)可表示为:uam=Uam(t)cos(wct),式中,Uam(t)为已调波的瞬时振幅值(也称为调幅波的包络函数)。由于调幅信号的瞬时振幅与调制信号成线性关系,既有:
Uam(t)=Ucm+KaUamcos(w0t)
=Ucm(1+(KaU0m/Ucm)cos(w0t)) =Ucm(1+macos(w0t))
式子中,Ka为比例常数,一般由调制电路的参数决定;ma=kaU0m/Ucm为调制系数(无单位),ma反映了调幅波振幅的改变量,常用百分数表示。把上述两式可以得出单频信号调幅波的表达式为:
.下载可编辑.
uam=Ucm(1+macos(w0t))cos(wct);
以上分析是在单一正弦信号作为调制信号的情况下进行的。实际传送的调制信号往
. .. .. .
往并非单一频率的信号,而是一个具有连续频谱的限带信号。如果将某一连续信号的限带信号U0(t)=f(t)作为调制信号,那么调幅波可表示为:
Uam=[Ucm+kaf(t)]cos(wct)
将其f(t)利用傅立叶级数展开为:
F(t)=∑ U0ncos(0nt)
将上面两式联合,则调幅波的表达式为:
Uam=Ucm[1+∑mncos(w0nt)]cos(wct)
式子中:mn=kaU0n/Ucm
则我们可以根据上述Uam式子,来进行调制的数学表达,让Matlab来实现计算和绘图。
另外,上式中的ma为调制系数,它反映了信号调制的强弱程度,一般ma的值越大调幅度越深。
在调幅波信号的分析中常用贫与分析法(即采用频谱图)来表达振幅调制的特。在上述单频调幅信号的频谱中,对最后的uam表达式利用三角函数的公式展开为:
uam=Ucm(1+macos(w0t))cos(wct)
=Ucm[cos(wct)+1/2macos(wc+w0)t+1/2macos(wc-w0)t]
可见,单频调幅波并不是一个简单的正弦波,其中包含有三个频率分量,即载波分量wc,上边频(USF)分量wu=wc+w0和下边频(LSF)分量w1=wc-w0,上下边频分量相当于载波是对称的,每个边频分量的振幅是调幅波包络振幅的一半。对于限带调幅信号的频谱,将Uam的表达式展开:
Uam=Ucm[1+∑mncosw0t]cos(wct)
=Ucm{cos(wct)+∑[1/2mncos(wc-w0)t+1/2mncos(wc+w0)]}
可见,经调制后限带信号的各个频率都会产生各自的上变频和下边频,叠加后就形成了所谓的上边频带和下边频带。因为上、下边频幅度相等且成对出现,所以上下边频带的频谱分布相对于载波是镜像对称的。
.下载可编辑.
. .. .. .
第 4 章 课程设计项目进度表
日期 2015-12.21~12.22 2015-12.22~12.23 2015-12.23~12.24 完成的工作 项目可行性研究,研究报告 写AM仿真程序 AM仿真程序的调试 2015-12.24~12.25 撰写课程设计总结报告 2015-12.25 交课程设计纸质和电子版材料 第 5 章 课程设计任务分配表
成员 座号 项目容 1、单一频率AM调制信号程序的编写 2、实验心得的编写 3、AM项目意义的编写 4、报告的复查 1、非单一频率AM调制信号程序的编写 2、AM设计思想的编写 3、AM设计心得的编写 4、报告的整理 序号 清海 1号 01 阮文添 24号 02 第 6 章 达到的效果
6.1 程序设计思想
首先验证了标准调幅调制的设计原理,另外可以得出AM调幅波的特点: (1)调幅波大振幅(包络)随调制信号变化,而且包络的变化规律与调制信号波形一致,表明调制信号(信息)记载在调幅波的包络中。
.下载可编辑.
. .. .. .
(2)由上面的UAM(t)的表达式: uam=Ucm(1+macos(w0t))cos(wct);得出调幅波的包络函数为:Uam(t)=Ucm(1+macos(w0t)),则得出了调幅波包络 波峰值为:Uam|max=Uam(1+ma)
波谷值为:Uam|min=Uam(1-ma)
(3)分析时域图和频谱图可见:载波分量并不包含调制信息,调制信息只包含在上,下边频分量,边频的振幅反映了调制信号幅度的大小。并且单频调幅波的频谱实质上是把低频调制信号的频谱线性搬移到载波的上下边频,调幅过程实质上就是一个频谱的线性搬移过程。
(4)由五种不同的调制系数ma得到的不同的调至结果图可以看出,调制系数ma反映了调幅的强弱程度,一般ma的值越大调幅越深(图a,b,c)。当ma=0时,表示未调幅,即无调幅作用;当ma=1时,调制系数的百分比达到100%,Um=Ucm,此时的包络振幅的最小值Uam|min=0;当ma>1时(图e),已调波的包络形状与调制信号不一样,产生了严重的包络失真。这种情况称之为过量调幅[3]。实际应用中必须尽力避免。因此,在振幅调制过程中为了避免产生过量调幅失真,保证已调波的包络真实的反映出调制信号的变化规律,要求调制系数ma必须满足:0 m’=√(错误!未找到引用源。)。一般在调制过程中,要保证所有的频率信号的调制都不会引起过量调幅失真。 6.2 程序最终实现结果 a) 当调制信号是单一频率时: 1. 载波信号的分析: .下载可编辑. . .. .. . 图6.1 2. 调制信号的分析: 图6.2 .下载可编辑. . .. .. . 3. 已调波信号的分析 图6.3 放大x轴 图6.4 4. 改变调制深度m(上例中m=0.5) 当m=0.3时: .下载可编辑. . .. .. . 图6.5 当m=0.8时: 图6.6 当m=1时: .下载可编辑. . .. .. . 6.7 当m=1.5时(即m>1): 图6.8 .下载可编辑. . .. .. . b) 当调制信号不是单一频率时: 1. 调制信号: 图6.9 2. 已调波信号: 图6.10 .下载可编辑. . .. .. . 3. 当K=5时,此时的m都大于1(上面的k=1.5时,m<1): 图6.11 第 7 章 源程序 7.1 当调制信号是单一频率时 1. 载波信号: t=-1:0.000001:1; A0=10; %定义时域t的围和步进 %载波信号的振幅A0 %载波信号的频率 f1=3000; w0=f1*pi*2; .下载可编辑. . .. .. . Uc=A0*cos(w0*t); figure(1); %载波信号 %绘图 subplot(2,1,1); plot(t,Uc);title('高频载波'); axis([0,0.01,-15,15]); subplot(2,1,2); Y1=fft(Uc); %对载波信号进行快速傅立叶变换 %定义图像显示的横纵坐标围 plot(abs(Y1));title('高频载波频谱'); axis([5000,7000,0,20000000]); 2. 调制信号: t=-1:0.000001:1; A1=5; f2=30; %调制信号 %绘图 %定义时域t的围和步进 %调制信号的振幅A1 %调制信号的频率 w1=f2*pi*2; U0=A1*cos(w1*t); figure(1); subplot(2,1,1); plot(t,U0);title('调制信号'); axis([0,1,-15,15]); subplot(2,1,2); Y2=fft(U0); %对调制信号进行快速傅立叶变换 %定义图像显示的横纵坐标围 plot(abs(Y2));title('调制信号频谱'); axis([0,1000,0,10000000]); 3. 已调波信号: t=-1:0.00001:1; A0=10;A1=5; %给定调制信号和载波的振幅 %给定调制信号和载波的频率 %频率向角频率的转换 f0=3000;f1=30; w0=2*f0*pi; .下载可编辑. . .. .. . w1=2*f1*pi; m=0.5; %调制系数 Uam=A1*(1+m*cos(w1*t)).*cos((w0).*t); %已调波信号 subplot(2,1,1); plot(t,Uam); axis([0,0.25,-10,10]); grid on; title('AM调制信号波形'); subplot(2,1,2); Y3=fft(Uam); %对已调波进行快速傅里叶变换 %绘图 plot(abs(Y3)),grid; title('AM调制信号频谱'); axis([5000,7000,0,1000000]); 7.2 当调制信号不是单一频率时 1. 调制信号: t=-1:0.00001:1; A1=5; A2=4; A3=3; f1=100; f2=200; f3=300; w1=2*f1*pi; w2=2*f2*pi; w3=2*f3*pi; U1=A1*cos(w1*t)+A2*cos(w2*t)+A3*cos(w3*t) %合成的调幅波形式 %绘图 %频率向角速度的转换 %三种不同频率调幅波的频率 %三种不同频率调幅波的幅值 subplot(2,1,1); .下载可编辑. . .. .. . plot(t,U1); axis([0,0.03,-30,30]); grid on; title('调制信号波形'); subplot(2,1,2); Y3=fft(U1); %对调幅波进行快速傅立叶变换 plot(abs(Y3)),grid; title('调制信号频谱'); axis([0,1000,0,1000000]); 2. 已调波信号: t=-1:0.00001:1; A0=10; A1=5; A2=4; A3=3; f0=3000; f1=100; f2=200; f3=300; w0=2*f0*pi; w1=2*f1*pi; w2=2*f2*pi; w3=2*f3*pi; k=1.5; %设定比例常数K的值 %分别计算m值 %频率向角速度的转换 %载波频率 %三种不同频率调幅波的频率 %载波的幅值 %三种不同频率调幅波的幅值 m1=k*A1/A0; m2=k*A2/A0; m3=k*A3/A0; U1=A1*cos(w1*t)+A2*cos(w2*t)+A3*cos(w3*t); .下载可编辑. . .. .. . %合成的调幅波形式 Uam=A0*(1+m1*cos(w1*t)+m2*cos(w2*t)+m3*cos(w3*t)).*cos((w0).*t); %已调波 %绘图 subplot(2,1,1); plot(t,Uam); axis([0,0.03,-30,30]); grid on; title('AM已调波信号波形'); subplot(2,1,2); Y3=fft(Uam); %对已调波作快速傅里叶变换 plot(abs(Y3)),grid; title('AM已调波信号频谱'); 第 9 章 设计心得 本设计完成了标准振幅调制的模拟,以及对调制系数影响的设计实现。Matlab方面相对于以前对它的了解,学会了使用快速傅里叶变换,并对其实现绘图分析,实验中暂时没有用到自己编写的M函数文件[1],直接用到了一个置的fft变换,得出了信号的频谱图;采用控制变量法,进行对调制系数m的影响的设计探究。在试验中,绘图的坐标围影响了图形显示的细节与否,将围调小,可以看到其实际波形变化,大围时可以轻易分析其轮廓,即外部特征。 调制技术是通信电子线路中很基础,很重要的一个环节,是实现通信必不可少的一项技术。这里选取调制中最简单最基础的标准振幅调制进行设计,因为其他的振幅调制方式如:双边带调幅(DSB),单边带(SSB),AM残边带等形式,这些调制的数学基础,即函数表达式基本类似。实验中运用Matlab这一强大的数学处理软件,很容易就能得到要想的傅里叶变换和精准的绘图,然后根据输出结果进行理论分析,对于学习Matlab软件非常有利。 .下载可编辑. . .. .. . 第 10 章 参考文献 [1] 建平,戴永夏,潘玲玲,王睿韬.,数字信号处理实验教程——基于MATLB,DSP 和SOPC实现[M], 清华大学,2010年 [2] 王卫东、高频电子线路、电子工业,2009年 [3] 王华 有军,matlab电子仿真与应用教程,国防工业,2010年 .下载可编辑. . .. .. . 学校地址:省武夷山市武夷大道16号 设计单位:数学与计算机学院 版本号:WyuKcsj Ver2007 .下载可编辑. 因篇幅问题不能全部显示,请点此查看更多更全内容