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=U m cos(w0t),载波为u c=U c cos(w c t),那么调幅信号（已调波）可表示为：u am=U am(t)cos(w c t),式中，Uam(t)为已调波的瞬时振幅值（也称为调幅波的包络函数）。
由于调幅信号的瞬时振幅与调制信号成线性关系，既有：U am(t)=U cm+K a U am cos(w0t)=U cm(1+(K a U0m/U cm)cos(w0t))=U cm(1+m a cos(w0t))式子中，Ka为比例常数，一般由调制电路的参数决定；m a=k a U0m/U cm为调制系数（无单位），ma反映了调幅波振幅的改变量，常用百分数表示。
把上述两式可以得出单频信号调幅波的表达式为：u am=U cm(1+m a cos(w0t))cos(w c t);以上分析是在单一正弦信号作为调制信号的情况下进行的。
实际传送的调制信号往往并非单一频率的信号，而是一个具有连续频谱的限带信号。
如果将某一连续信号的限带信号U0(t)=f(t)作为调制信号，那么调幅波可表示为：Uam=[U cm+k a f(t)]cos(w c t)将其f(t)利用傅立叶级数展开为：F(t)=∑ U0n cos(0n t)将上面两式联合，则调幅波的表达式为：U am=U cm[1+∑m n cos(w0n t)]cos(w c t)式子中：m n=k a U0n/U cm则我们可以根据上述Uam式子，来进行调制的数学表达，让Matlab来实现计算和绘图。
另外，上式中的m a为调制系数，它反映了信号调制的强弱程度，一般ma的值越大调幅度越深。
在调幅波信号的分析中常用贫与分析法（即采用频谱图）来表达振幅调制的特。
在上述单频调幅信号的频谱中，对最后的u am表达式利用三角函数的公式展开为：u am=U cm(1+m a cos(w0t))cos(w c t)=U cm[cos(w c t)+1/2macos(w c+w0)t+1/2macos(w c-w0)t] 可见，单频调幅波并不是一个简单的正弦波，其中包含有三个频率分量，即载波分量w c,上边频（USF）分量w u=w c+w0和下边频（LSF）分量w1=w c-w0,上下边频分量相当于载波是对称的，每个边频分量的振幅是调幅波包络振幅的一半。
对于限带调幅信号的频谱，将U am的表达式展开：U am=U cm[1+∑mn cosw0t]cos(w c t)=Ucm{cos(wct)+∑[1/2mncos(wc-w0)t+1/2mncos(wc+w0)]}可见，经调制后限带信号的各个频率都会产生各自的上变频和下边频，叠加后就形成了所谓的上边频带和下边频带。
因为上、下边频幅度相等且成对出现，所以上下边频带的频谱分布相对于载波是镜像对称的。
第 4 章课程设计项目进度表第 5 章课程设计任务分配表第 6 章达到的效果6.1 程序设计思想首先验证了标准调幅调制的设计原理，另外可以得出AM调幅波的特点：（1）调幅波大振幅（包络）随调制信号变化，而且包络的变化规律与调制信号波形一致，表明调制信号（信息）记载在调幅波的包络中。
（2）由上面的U AM(t)的表达式: u am=U cm(1+m a cos(w0t))cos(w c t);得出调幅波的包络函数为：U am(t)=U cm(1+m a cos(w0t)),则得出了调幅波包络波峰值为：U am|max=U am(1+m a)波谷值为：U am|min=U am(1-m a)（3）分析时域图和频谱图可见：载波分量并不包含调制信息，调制信息只包含在上，下边频分量，边频的振幅反映了调制信号幅度的大小。
并且单频调幅波的频谱实质上是把低频调制信号的频谱线性搬移到载波的上下边频，调幅过程实质上就是一个频谱的线性搬移过程。
（4）由五种不同的调制系数m a得到的不同的调至结果图可以看出，调制系数m a反映了调幅的强弱程度，一般ma的值越大调幅越深（图a,b,c）。
当m a=0时，表示未调幅，即无调幅作用；当m a=1时，调制系数的百分比达到100%，U m=U cm,此时的包络振幅的最小值U am|min=0;当m a>1时（图e），已调波的包络形状与调制信号不一样，产生了严重的包络失真。
这种情况称之为过量调幅[3]。
实际应用中必须尽力避免。
因此，在振幅调制过程中为了避免产生过量调幅失真，保证已调波的包络真实的反映出调制信号的变化规律，要求调制系数m a必须满足：0<m a<1.（5）对于多频率的调制信号进行调幅时，由于各频率信号的幅度不同，因而调制系数mn也不相同。
长引用合成调制系数m’=√(错误!未找到引用源。
)。
一般在调制过程中，要保证所有的频率信号的调制都不会引起过量调幅失真。
6.2 程序最终实现结果a)当调制信号是单一频率时：1.载波信号的分析：图6.1 2.调制信号的分析：图6.23.已调波信号的分析图6.3放大x轴图6.44.改变调制深度m（上例中m=0.5）当m=0.3时：图6.5 当m=0.8时：图6.6当m=1时：6.7 当m=1.5时（即m>1）:图6.8b)当调制信号不是单一频率时：1.调制信号：图6.9 2.已调波信号：图6.103.当K=5时，此时的m都大于1（上面的k=1.5时，m<1）:图6.11第 7 章源程序7.1 当调制信号是单一频率时1. 载波信号：t=-1:0.000001:1; %定义时域t的围和步进A0=10; %载波信号的振幅A0f1=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; %定义时域t的围和步进A1=5; %调制信号的振幅A1f2=30; %调制信号的频率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的值m1=k*A1/A0; %分别计算m值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 章设计心得本设计完成了标准振幅调制的模拟，以及对调制系数影响的设计实现。