基于MATLAB 的模拟调制实验报告
一、实验目的
1.进一步学习调制的知识，掌握调频与调角两种模拟调制技术。
2.进一步学习MATLAB 的编程，熟练使用MATLAB 进行作图。
二、实验原理
1.调制的概念
调制（modulation ）就是对信号源的信息进行处理加到载波上，使其变为适 合 于信道传输的形式的过程，是使载波随信号而改变的技术。
一般，用来传送消息的信号()t u c 叫作载波或受调信号，代表所欲传送消息的信
号叫作调制信号，调制后的信号()t u 叫作已调信号。用调制信号()t u Ω控制载波
的某些参数，使之随()t u Ω而变化，就可实现调制。
2.调制的目的
频谱变换
当所要传送的信号的频率或者太低，或者频带很宽，对直接采用电磁波的形 式进行发送很不利，需要的天线尺寸很大，而且发射和接受短的天线与谐振回路的参数变化范围很大。为了信息有效与可靠传输，往往需要将低频信号的基带频谱搬移到适当的或指定的频段。这样可以提高传输性能，以较小的发送功率与较短的天线来辐射电磁波。
实现信道复用
为了使多个用户的信号共同利用同一个有较大带宽的信道，可以采用各种复用技术。如模拟电话长途传输是通过利用不同频率的载波进行调制。将各用户话音每隔4 kHz 搬移到高频段进行传输。
提高抗干扰能力
不同的调制方式，在提高传输的有效性和可靠性方面各有优势。如调频广播
系统，它采用的频率调制技术，付出多倍带宽的代价，由于抗干扰性能强，其音质比只占10 kHz 带宽的调幅广播要好得多。扩频通信就是以大大扩展信号传输带宽，以达到有效抗拒外部干扰和短波信道多径衰落的特殊调制方式。
3.调制的种类
根据()t u Ω和()t u c 的不同类型和完成调制功能的调制器传递函数不同，调制
分为以下多种方式：
（1）．按调制信号()t u Ω的类型分为：
● 模拟调制：调制信号()t u Ω是连续变化的模拟量，如话音与图像信号。
● 数字调制：调制信号是数字化编码符号或脉冲编码波形。 （2）．按载波信号()t u c 的类型分：
● 连续波调制：载波信号为连续波形，通常以正弦波作为载波。
● 脉冲调制：载波信号是脉冲波形序列。
（3）．按调制器的不同功能分：
● 幅度调制：以调制信号去控制载波的幅度变化，如模拟调幅，脉冲幅度调制
（PAM ），幅移键控（ASK ）。
● 频率调制：以调制信号去控制载波信号的频率变化，如模拟调频（FM ），频
移键控（FSK ），脉宽调制（PDM ）。
● 相位调制：以调制信号去控制载波信号的相位变化，如模拟调相（PM ），相
移键控（PSK ），脉位调制（PPM ）。
（4）．按调制器的传输函数分：
● 线性调制：已调信号的频谱与调制信号频谱是线性的频谱位移关系。如各种
幅度调制，幅移键控（ASK ）。
● 非线性调制：已调信号的频谱与调制信号频谱没有线性关系，即调制后派生
出大量不同于调制信号的新的频率成份。如调频（FM ），调相（PM ），频移键
控（FSK ）。
三、实验过程
这次实验主要利用MATLAB 实现幅度调制与角度调制。
1.幅度调制原理
AM 调制
AM 是用调制信号去控制高频正弦载波的幅度，使其按调制信号的规律变化
的过程。如图1所示。
图1
调制信号()t m 叠加直流0A 后再与载波相乘，则输出的信号就是常规双边带
（AM ）调幅波。其时域表达式为：
0()[()]*cos AM c S t A m t t ω=+
AM 信号波形的包络与输入基带信号()t m 成正比，故用包络检波的方法很容
易恢复原始调制信号。但为了保证包络检波时不发生失真，必须满足
()max 0t m A ≥，否则将出现过调幅现象而带来失真。AM 信号的频谱是由载频分量
和上、下两个边带组成。上边带的频谱与原调制信号的频谱结构相同，下边带是
上边带的镜像。AM 信号是带有载波的双边带信号，它的带宽信号带宽的两倍2AM H B f =。
下面利用matlab 实现AM 的调制。为方便观察，令输入信号为单一正弦波信号。
-----------------AM调制输出信号和频谱---------------------- dt=0.0001; %时间采样频率
fH=10; %调制信号最高频率
fc=1000; %载波中心频率
T=0.2; %信号时长
N=T/dt;
t=[0:N-1]*dt;
mt=sqrt(2)*cos(2*pi*fH*t); %调制信号
A=2;
s_am=(A+mt).*cos(2*pi*fc*t);
[f,Xf]=FFT_SHIFT(t,s_am); %已调信号频谱
subplot(211);
plot(t,s_am);hold on; %画出AM信号波形
plot(t,A+mt,'r--'); %表示AM上包络
plot(t,-(A+mt),'r--'); %表示AM下包络
title('AM调制信号及其包络');
xlabel('时间t');
subplot(212); %画出AM波频谱
plot(f,abs(Xf));
axis([700 1300 0 max(Xf)]);
title('AM信号功率谱');
xlabel('频率f');
图2
结果分析：
从图2中频谱图可以看出，调幅过程实际就是频谱搬移，是一个线性过程。
经调制后，调制信号的频谱被搬移到载频附近，成为下边频（带）和下边频（带）。
信号功率正比于频谱的平方，所以可以看出载波信号的功率占整个调幅波功率的
绝大部分。
同样由公式计算可以得到 2(1)2a O OT m P P =+
其中，OT P 为载波功率。
所以在发送中，我们可以只发送上、下边带，此时为DSB 调制；只发送上边带或
者下边带，此时为SSB 。
DSB 调制
在幅度调制的一般模型中，若假设滤波器为全通网络（＝1），，调制信
号m(t)中无直流分量，则输出的已调信号就是无载波分量的双边带调制信号，
或称抑制载波双边带（DSB-SC ）调制信号，简称双边带（DSB ）信号。DSB 调制
器模型如下图3所示。
t c ωcos )
(t m )
(t s DSB ×图3-7 D SB 调制器模型
图3
下面利用matlab实现DSB的调制。为方便观察，令输入信号为单一正弦波信号。
------------------DSB调制信号输出和频谱-------------------- dt=0.0001; %时间采样频率
fH=10; %调制信号最高频率
fc=1000; %载波中心频率
T=0.1; %信号时长
N=T/dt;
t=[0:N-1]*dt;
mt=sqrt(2)*cos(2*pi*fH*t); %调制信号
s_dsb=mt.*cos(2*pi*fc*t);
[f,sf]=FFT_SHIFT(t,s_dsb); %已调信号频谱
subplot(211)
plot(t,s_dsb);hold on; %画出DSB信号波形
plot(t,mt,'r--'); %标示mt波形
plot(t,-mt,'r--');
title('DSB调制信号及其包络');
xlabel('时间t');
subplot(212)
plot(f,PSD);
axis([750, 1250, 0 ,0.015]);
title('DSB信号频谱');
xlabel('频率f');
图4
结果分析：
从图4可以看出，DSB 信号时域图与频谱中不含有载波分量，仅有上下边带，发
射功率大大降低。
SSB
由于DSB 信号的上、下两个边带是完全对称的，皆携带了调制信号的全部信息，
因此，从信息传输的角度来考虑，仅传输其中一个边带就够了。用滤波法实现单
边带调制的原理图如下图5所示
t c ωcos )
(t m )
(t s SSB ×图3-9 S SB信号的滤波法产生)(ωSSB H
图5
下面利用matlab实现SSB的调制。为方便观察，令输入信号为单一正弦波信号。
------------------SSB调制信号输出和频谱-------------------- dt=0.0001; %时间采样频率
fH=10; %调制信号最高频率
fc=1000; %载波中心频率
T=0.1;
N=T/dt;
t=[0:N-1]*dt;
mt=sqrt(2)*cos(2*pi*fH*t); %调制信号
s_ssb=real(hilbert(mt).*exp(j*2*pi*fc*t));
[f,sf]=FFT_SHIFT(t,s_ssb); %单边带信号频谱
subplot(211)
plot(t,s_ssb);hold on; %画出SSB信号波形
plot(t,mt,'r--'); %标示mt 的包络
plot(t,-mt,'r--'); %标示mt 的包络
title('SSB调制信号');
xlabel('时间t');
subplot(212)
plot(f,abs(sf));
axis([750 1250 0 max(abs(sf))]);
title('SSB信号频谱');
xlabel('频率f');
图6
结果分析：
从图6的调制信号频谱可以看出，输出信号只有上（或者下）边带，大大减
小了发射功率。而且SSB 信号频带可节约一半，这对于日益拥挤的短波波段（3
—30MHz ）来说有重大意义，因为这样就能在同一波段中，使所容纳的频道数目
增加一倍，大大提高了短波波段的利用率。
总结：
在以上模拟线性调制中，AM 波发送信号的上下边带和载波，带宽为调信号
的最大频率的两倍，发射功率较大；DSB 波只发送信号的上下边带，发射功率大
大减小，但带宽仍然为调信号的最大频率的两倍；SSB 波只发送信号上边带或者
下边带，带宽是AM 与DSB 波的一半。但是，在解调过程中，AM 波采用非相干解
调（包络检波），电路十分简单，采用二极管检波电路就可以恢复出调制信号。
DSB 、SSB 需采用相干解调，电路较复杂。所以在广播通信中广泛采用AM 调制。
2.角度调制
在调制时，若载波的频率随调制信号变化，称为频率调制，简称调频；若载
波的相位随调制信号变化，称为相位调制，简称调相。在这两种调制中，载波的
幅度保持恒定不变，而频率和相位的变化都表现为载波的瞬时相位变化，故把调
频和调相统称为角度调制或调角。
调角的一般表达式为 ()cos[()]m c s t A t t ω?=+
式中：A 为载波的恒定振幅，[()]c t t ω?+为信号的瞬时相位，极为()t θ；()t ?为
相对于载波的相位c t ω的瞬时相位偏移；[()]/c d t t t ω?+是信号的瞬时角频率，记
为()t ω；而()/d t t ?成为相对于载频c ω的瞬时频偏。
摘要 本次设计主要是以Matlab为基础平台，对FM信号进行仿真。介绍了FM信号，及其调制和解调的基本原理，并设计M文件，分析在混入噪声环境下的波形失真，以及分析FM的抗噪声性能。本设计的主要目的是对Matlab的熟悉和对模拟通信理论的更深化理解。 关键词：Matlab；FM；噪声
前言 (2) 1 设计基础 (3) 1.1 Matlab及M文件的简介 (3) 1.2模拟调制概述 (4) 1.2.1模拟调制系统各个环节分析 (5) 1.2.2 模拟调制的意义 (6) 2 FM基本原理与实现 (7) 2.1 FM的基本原理 (7) 2.1.1调制 (7) 2.1.2解调 (8) 2.2 FM的实现 (8) 2.2.1 FM调制的实现 (8) 2.2.2 FM解调的实现 (9) 2.3 调频系统的抗噪声性能 (10) 2.3.1 高斯白噪声信道特性 (10) 3 FM的仿真实现与分析 (14) 3.1 未加噪声的FM解调实现 (14) 3.2 叠加噪声时的 FM解调 (16) 总结 (20) 致谢 (21) 参考文献 (22) 附录 (23)
通信按照传统的理解就是信息的传输。在当今高度信息化的社会，信息和通信已成为现代社会的命脉。信息作为一种资源，只有通过广泛传播与交流，才能产生利用价值，促进社会成员之间的合作，推动社会生产力的发展，创造出巨大的经济效益。而通信作为传输信息的手段或方式，与传感技术、计算机技术相融合，已成为21世纪国际社会和世界经济发展的强大动力。可以预见，未来的通信对人们的生活方式和社会的发展将会产生更加重大和意义深远的影响。 在通信系统中，从消息变换过来的原始信号所占的有效频带往往具有频率较低的频谱分量(例如语音信号)，如果将这种信号直接在信道中进行传输，则会严重影响信息传送的有效性和可靠性，因此这种信号在许多信道中均是不适宜直接进行传输的。在通信系统的发射端通常需要有调制过程，将调制信号的频谱搬移到所希望的位置上，使之转换成适于信道传输或便于信道多路复用的已调信号；而在接收端则需要有解调过程，以恢复原来有用的信号。调制解调方式常常决定了一个通信系统的性能。随着数字化波形测量技术和计算机技术的发展，可以使用数字化方法实现调制与解调过程。 调制在通信系统中具有重要的作用。通过调制，不仅可以进行频谱搬移，把调制信号的频谱搬移到所希望的位置上，从而将调制信号转换成适合于信道传输或便于信道多路复用的已调信号，而且它对系统的传输有效性和传输可靠性有着很大的影响。调制方式往往决定了一个通信系统的性能。调制技术是指把基带信号变换成传输信号的技术。基带信号是原始的电信号，一般是指基本的信号波形，在数字通信中则指相应的电脉冲。在无线遥测遥控系统和无线电技术中调制就是用基带信号控制高频载波的参数（振幅、频率和相位），使这些参数随基带信号变化。用来控制高频载波参数的基带信号称为调制信号。未调制的高频电振荡称为载波（可以是正弦波，也可以是非正弦波，如方波、脉冲序列等）。被调制信号调制过的高频电振荡称为已调波或已调信号。已调信号通过信道传送到接收端，在接收端经解调后恢复成原始基带信号。
数学实验报告 班级: 学号： 姓名： 实验序号：1 日期：年 月 日 实验名称：特殊函数与图形 ◆ 问题背景描述：绘图是数学中的一种重要手段，借助图形，可以使抽象的对象得到 明白直观的体现，如函数的性质等。同时，借助直观的图形，使初学者更容易接受新知识，激发学习兴趣。 ◆ 实验目的：本实验通过绘制一些特殊函数的图形，一方面展示这些函数的特点属性， 另一方面，就 Matlab 强大的作图功能作一个简单介绍。 实验原理与数学模型： 1、 球2222x y z R ++= ，x=Rsin φcos θ, y= Rsin φsin θ, z= cos φ, 0≤θ≤2π , 0≤φ≤π 环面 222222222()4(),(cos )cos ,x y z a r a x y x a r φθ+++-=+=- (cos )sin ,sin ,02,02y a r z r φθφφπθπ=-=≤≤≤≤ 2、 平面摆线：2 22 31150,(sin ),(1cos ),0233 x y x a t t y a t t π+-==-=-≤≤ 3、 空间螺线：（圆柱螺线）x=acost , y=asint , z=bt ;（圆锥螺线）22 cos ,sin ,x t t y t t z t === 4、 椭球面sin cos ,sin sin ,cos ,02,0x a y b z c φθφθφθπφπ===≤<≤≤ 双叶双曲面3 tan cos ,tan sin ,sec ,02,22 x a y b z c π φθφθφθπφπ===≤<- << 双曲抛物面2 sec ,tan 2 u x au y bu z θθ=== 实验所用软件及版本：mathematica(3.0) 主要内容（要点）： 1、 作出下列三维图形（球、环面） 2、 作出下列的墨西哥帽子 3、 作出球面、椭球面、双叶双曲面，单叶双曲面的图形 4、 试画出田螺上的一根螺线 5、 作出如图的马鞍面
Matlab通信原理仿真 学号： 2142402 姓名：圣斌
实验一Matlab 基本语法与信号系统分析 一、实验目的： 1、掌握MATLAB的基本绘图方法； 2、实现绘制复指数信号的时域波形。 二、实验设备与软件环境： 1、实验设备：计算机 2、软件环境：MATLAB R2009a 三、实验内容： 1、MATLAB为用户提供了结果可视化功能，只要在命令行窗口输入相应的命令，结果就会用图形直接表示出来。 MATLAB程序如下： x = -pi::pi; y1 = sin(x); y2 = cos(x); %准备绘图数据 figure(1); %打开图形窗口 subplot(2,1,1); %确定第一幅图绘图窗口 plot(x,y1); %以x，y1绘图 title('plot(x,y1)'); %为第一幅图取名为’plot(x,y1)’ grid on; %为第一幅图绘制网格线 subplot(2,1,2) %确定第二幅图绘图窗口 plot(x,y2); %以x，y2绘图 xlabel('time'),ylabel('y') %第二幅图横坐标为’time’,纵坐标为’y’运行结果如下图： 2、上例中的图形使用的是默认的颜色和线型，MATLAB中提供了多种颜色和线型，并且可以绘制出脉冲图、误差条形图等多种形式图： MATLAB程序如下： x=-pi:.1:pi; y1=sin (x); y2=cos (x); figure (1); %subplot (2,1,1); plot (x,y1); title ('plot (x,y1)'); grid on %subplot (2,1,2); plot (x,y2);
南昌航空大学 数学与信息科学学院 实验报告 课程名称：数学实验 实验名称： MATLAB基本操作 实验类型：验证性■综合性□ 设计性□ 实验室名称：数学实验室 班级学号： 10 学生姓名：钟 X 任课教师（教师签名）： 成绩： 实验日期： 2011-１０- 1０
一、实验目的 1、熟悉MATLAB基本命令与操作 2、熟悉MATLAB作图的基本原理与步骤 3、学会用matlab软件做图 二、实验用仪器设备、器材或软件环境 计算机MATLAB软件 三、实验原理、方案设计、程序框图、预编程序等 问题1：在区间【0,2π】画sinx 实验程序: >> x=linspace(0,2*pi,30); >> y=sin(x); >> plot(x,y) 问题2：在【0,2π】用红线画sinx,用绿圈画cosx,实验程序:
>> x=linspace(0,2*pi,30); >> y=sin(x); >> z=cos(x); >> plot(x,y,'r',x,z,'co') >> 问题3：在【0，π】上画y=sinx的图形。 实验程序: >> ezplot('sin(x)',[0,pi]) >> 问题4：在【0，π】上画x=cos3t,y=sin3t星形图形。
实验程序: >> ezplot('cos(t).^3','sin(t).^3',[0,pi]) >> 问题5：[-2,0.5]，[0,2]上画隐函数 实验程序: >> ezplot('exp(x)+sin(x*y)',[-2,0.5,0,2]) >> 问题6：在[－２，２]范围内绘制tanh的图形。实验程序: >> fplot('tanh',[-2,2])
闽江学院 《通信原理设计报告》 题目：基于MATLAB的模拟调制系统仿真与测试学院：计算机科学系 专业：12通信工程 组长：曾锴（3121102220） 组员：薛兰兰（3121102236） 项施旭（3121102222） 施敏（3121102121） 杨帆（3121102106） 冯铭坚（3121102230） 叶少群（3121102203） 张浩（3121102226） 指导教师：余根坚 日期：2014年12月29日——2015年1月4日
摘要在通信技术的发展中，通信系统的仿真是一个重点技术，通过调制能够将信号转化成适用于无线信道传输的信号。 在模拟调制系统中最常用最重要的调制方式是用正弦波作为载波的幅度调制和角度调制。在幅度调制中，文中以调幅、双边带和单边带调制为研究对象，从原理等方面阐述并进行仿真分析；在角度调制中，以常用的调频和调相为研究对象，说明其调制原理，并进行仿真分析。利用MATLAB下的Simulink工具箱对模拟调制系统进行仿真，并对仿真结果进行时域及频域分析，比较各个调制方式的优缺点，从而更深入地掌握模拟调制系统的相关知识，通过研究发现调制方式的选取通常决定了一个通信系统的性能。 关键词模拟调制；仿真；Simulink 目录 第一章绪论 (1) 1.1 引言 (1) 1.2 关键技术 (1) 1.3 研究目的及意义 (2) 1.4 本文工作及内容安排 (2) 第二章模拟调制原理 (3) 2.1 幅度调制原理 (3) 2.1.1 AM调制 (4) 第三章基于Simulink的模拟调制系统仿真与分析 (6) 3.1 Simulink工具箱简介 (6) 3.2 幅度调制解调仿真与分析 (8) 3.2.1 AM调制解调仿真及分析 (8) 第四章总结 (12) 4.1 代码 (13) 4.2 总结 (14)
实验一MATLAB操作基础 实验目的和要求： 1、熟悉MATLAB的操作环境及基本操作方法。 2、掌握MATLAB的搜索路径及设置方法。 3、熟悉MATLAB帮助信息的查阅方法 实验内容： 1、建立自己的工作目录，再设置自己的工作目录设置到MA TLAB搜索路径下，再试 验用help命令能否查询到自己的工作目录。 2、在MA TLAB的操作环境下验证课本；例1-1至例1-4，总结MATLAB的特点。 例1-1
例1-2 例1-3 例1-4
3、利用帮助功能查询inv、plot、max、round等函数的功能。 4、完成下列操作： （1）在matlab命令窗口输入以下命令： x=0:pi/10:2*pi; y=sin(x); （2）在工作空间窗口选择变量y，再在工作空间窗口选择回绘图菜单命令或在工具栏中单击绘图命令按钮，绘制变量y的图形，并分析图形的含义。
5、访问mathworks公司的主页，查询有关MATLAB的产品信息。 主要教学环节的组织： 教师讲授实验目的、开发环境界面、演示实验过程，然后同学上机练习。 思考题： 1、如何启动与退出MA TLAB集成环境？ 启动： （1）在windows桌面，单击任务栏上的开始按钮，选择‘所有程序’菜单项，然后选择MA TLAB程序组中的MA TLABR2008b程序选项，即可启动 MATLAB系统。 （2）在MA TLAB的安装路径中找到MA TLAB系统启动程序matlab.exe，然后运行它。 （3）在桌面上建立快捷方式后。双击快捷方式图标，启动MA TLAB。 退出： （1）在MA TLAB主窗口file菜单中选择exitMATLAB命令。 （2）在MA TLAB命令窗口中输入exit或quit命令。 （3）单击MATLAB主窗口的关闭按钮。 2、简述MATLAB的主要功能。 MATLAB是一种应用于科学计算领域的数学软件，它主要包括数值计算和符 号计算功能、绘图功能、编程语言功能以及应用工具箱的扩展功能。 3、如果一个MATLAB命令包含的字符很多，需要分成多行输入，该如何处理？
MATLAB 仿真实验报告 课题名称：MATLAB 仿真——图像处理 学院：机电与信息工程学院 专业：电子信息科学与技术 年级班级：2012级电子二班 一、实验目的 1、掌握MATLAB处理图像的相关操作，熟悉相关的函数以及基本的MATLAB语句。 2、掌握对多维图像处理的相关技能，理解多维图像的相关性质 3、熟悉Help 命令的使用，掌握对相关函数的查找，了解Demos下的MATLAB自带的原函数文件。 4、熟练掌握部分绘图函数的应用，能够处理多维图像。 二、实验条件
MATLAB调试环境以及相关图像处理的基本MATLAB语句，会使用Help命令进行相关函数查找 三、实验内容 1、nddemo.m函数文件的相关介绍 Manipulating Multidimensional Arrays MATLAB supports arrays with more than two dimensions. Multidimensional arrays can be numeric, character, cell, or structure arrays. Multidimensional arrays can be used to represent multivariate data. MATLAB provides a number of functions that directly support multidimensional arrays. Contents ： ●Creating multi-dimensional arrays 创建多维数组 ●Finding the dimensions寻找尺寸 ●Accessing elements 访问元素 ●Manipulating multi-dimensional arrays操纵多维数组 ●Selecting 2D matrices from multi-dimensional arrays从多维数组中选择二维矩 阵 (1)、Creating multi-dimensional arrays Multidimensional arrays in MATLAB are created the same way as two-dimensional arrays. For example, first define the 3 by 3 matrix, and then add a third dimension. The CAT function is a useful tool for building multidimensional arrays. B = cat(DIM,A1,A2,...) builds a multidimensional array by concatenating（联系起来）A1, A2 ... along the dimension DIM. Calls to CAT can be nested（嵌套）. (2)、Finding the dimensions SIZE and NDIMS return the size and number of dimensions of matrices. (3)、Accessing elements To access a single element of a multidimensional array, use integer subscripts（整数下标）. (4)、Manipulating multi-dimensional arrays
实验一 MATLAB 基本操作及运算 一、 实验目的 1、 理解Matlab 数据对象的特点； 2、 掌握基本Matlab 运算规则； 3、 掌握Matlab 帮助的使用方法； 二、 实验的设备及条件 计算机一台（带有MATLAB7.0以上的软件环境）。 三、 实验内容 要求建立一个名为experiment01.m 的，把与实验内容1-7相关的实验命令都放入该文件中，题与题之间用相应注释分割。注意对实验中出现的相关函数或变量，请使用help 或doc 查询相关帮助文档，学习函数的用法。 1、 建立以下标量： 1） a=10 2） b=2.5×1023 3） c=2+3i ，（i 为虚数单位） 4） d=3/2πj e ,（j 为虚数单位，这里要用到exp ，pi ） 2、 建立以下向量： 1） aVec=[3.14 15 9 26] 2） bVec=????? ???????18228871.2 3） cVec=[5 4.8 … -4.8 -5 ] （向量中的数值从5到-5，步长为-0.2） 4） dVec=[100 100.01 … 100.99 101] （产生1到10之间的等对数间隔向量，参考logspace ，注意向量的长度） 3、 建立以下矩阵： 1）???? ??????=2222 aMat aMat 一个9×9的矩阵，其元素全为2；（参考ones 或zeros ）
2）??????? ?????????=1000005000001 bMat bMat 是一个9×9的矩阵，除主对角上的元素为[1 2 3 4 5 4 3 2 1]外，其余元素均为0。（参考diag ）。 3）100 20109212291111 =cMat cMat 为一个10×10的矩阵，可有1:100的向量来产生（参考reshape ） 4）???? ??????=NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN dMat dMat 为3×4的NaN 矩阵，（参考nan ） 5）?? ????---=8710225113eMat 6）产生一个5×3随机整数矩阵fMat ，其值的范围在-3到3之间。（参考rand 和floor 或ceil ） 4、 使用题1中的变量计算下列等式的x,y,z 的值： 1） ) 6/)15((11--+=a e x 2） g g h h b a y /121,)(=+=提示π,参考sqrt 。 3） c c a d c d c R z ))3/sin()]))([(log(π-+= ，其中R 表示取括号内复数的实数部分，c 表示c 的共轭复数，log 是自然对数。（参考real ，conj ，log ） 5、 使用题2中的向量求解一下等式： 1）)25.2/(22 25.221 cVec e xVec -=π， 其中cVec 指的是题2 中定义的向量cVec ，一下雷同。 2）22)(bVec aVec yVec T +=,T aVec 表示aVec 的转置 3） )/1(log 10dVec zVec =，10log 表示已10为底的对数，参考log10 6、 使用题2和题3中所产生的向量和矩阵计算以下等式，注意本题的操作
课程设计任务书 学生姓名：杨刚专业班级：电信1302 指导教师：工作单位：武汉理工大学 题目：信号分析处理课程设计 －基于MATLAB的模拟信号频率调制（FM）与解调分析 初始条件： 1.Matlab6.5以上版本软件； 2.先修课程：通信原理等； 要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求） 1、利用MATLAB中的simulink工具箱中的模块进行模拟频率（FM）调制与解调，观 察波形变化 2、画出程序设计框图，编写程序代码，上机运行调试程序，记录实验结果（含计算结 果和图表等），并对实验结果进行分析和总结； 3、课程设计说明书按学校统一规范来撰写，具体包括： ⑴目录；⑵理论分析； ⑶程序设计；⑷程序运行结果及图表分析和总结； ⑸课程设计的心得体会（至少800字，必须手写。）； ⑹参考文献（不少于5篇）。 时间安排： 周一、周二查阅资料，了解设计内容； 周三、周四程序设计，上机调试程序； 周五、整理实验结果，撰写课程设计说明书。 指导教师签名： 2013 年 7月 2 日 系主任（或责任教师）签名： 2013年 7月 2日
目录 1 Simulink简介 (1) 1.1 Matlab简介······················································错误！未定义书签。 1.2 Simulink介绍 ···················································错误！未定义书签。 2 原理分析 ·····························································错误！未定义书签。 2.1通信系统 ·························································错误！未定义书签。 2.1.1通信系统的一般模型 ···································错误！未定义书签。 2.1.2 模拟通信系统 (3) 2.2 FM调制与解调原理···········································错误！未定义书签。 3 基于Matlab方案设计 (6) 3.1 Matlab代码 (6) 3.2 Matlab仿真 (8) 4 基于Simulink方案设计 (12) 4.1 使用Simulink建模和仿真的过程 (12) 4.1.1 Simulink模块库简介 (12) 4.1.2 调制解调模块库简介 (13) 4.2 FM调制与解调电路及仿真 (14) 4.3 仿真结果分析 (17) 5 心得体会 ·····························································错误！未定义书签。 6 参考文献 (20) 本科生课程设计评定表
实验一 MATLAB 及仿真实验（控制系统的时域分析） 一、实验目的 学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性； 二、预习要点 1、 系统的典型响应有哪些？ 2、 如何判断系统稳定性？ 3、 系统的动态性能指标有哪些？ 三、实验方法 （一） 四种典型响应 1、 阶跃响应： 阶跃响应常用格式： 1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。 2、),(Tn sys step ；表示时间范围0---Tn 。 3、),(T sys step ；表示时间范围向量T 指定。 4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。 2、 脉冲响应： 脉冲函数在数学上的精确定义：0 ,0)(1)(0 ?==?∞ t x f dx x f 其拉氏变换为：) ()()()(1)(s G s f s G s Y s f === 所以脉冲响应即为传函的反拉氏变换。 脉冲响应函数常用格式： ① )(sys impulse ； ② ); ,();,(T sys impulse Tn sys impulse ③ ),(T sys impulse Y = （二） 分析系统稳定性 有以下三种方法： 1、 利用pzmap 绘制连续系统的零极点图； 2、 利用tf2zp 求出系统零极点； 3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析 Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.
实验一matlab基本操作 一、实验目的 熟悉matlab的安装与启动；熟悉matlab用户界面；熟悉matlab功能、建模元素；熟悉matlab优化建模过程。 二、实验设备与工具 1.计算机 2.matlab软件 三、实验步骤 1. 了解matlab的硬件和软件必备环境； 2. 启动matlab； 3. 学习优化建模过程。 四、实验报告要求 1. 写出matlab系统界面的各个构成；以及系统布局区的组成；以及每一部 分的功能； 2. 优化建模过程应用举例 五、实验内容 （一）、Matlab操作界面 1.命令窗口（command window） 2.命令历史窗口（command history） 3.工作空间管理窗口（workspace） 4.当前路径窗口（current directory） （二）、优化建模过程应用举例 1、简单矩阵 123 456 789 A ?? ?? =?? ?? ?? 的输入步骤。 （1）在键盘上输入下列内容 A = [1,1,3; 4,5,6; 7,8,9] （2）按【Enter】键，指令被执行。 （3）在指令执行后，MATLAB指令窗中将显示以下结果： A = 1 2 3 4 5 6 7 8 9
2、矩阵的分行输入。 A=[1,2,3 4,5,6 7,8,9] A = 1 2 3 4 5 6 7 8 9 3、指令的续行输入 S=1-1/2+1/3-1/4+ ... 1/5-1/6+1/7-1/8 S = 0.6345 4、画出衰减振荡曲线t e y t 3sin 3-=及其它的包络线3 0t e y -=。t 的取值范围是]4,0[π。 t=0:pi/50:4*pi; y0=exp(-t/3); y=exp(-t/3).*sin(3*t); plot(t,y,'-r',t,y0,':b',t,-y0,':b') grid 5、画出2222) sin(y x y x z ++=所表示的三维曲面。y x ,的取值范围是]8,8[-。 clear;x=-8:0.5:8; y=x';
1 线性模拟调制 1.1模拟调制原理 模拟调制是指用来自信源的基带模拟信号去调制某个载波，而载波是一个确知的周期性波形。模拟调制可分为线性调制和非线性调制，本文主要研究线性调制。 线性调制的原理模型如图1.1所示。设c(t)=Acos2t f o π,调制信号为m(t),已调信号为s(t)。 图1.1 线性调制的远离模型 调制信号m(t)和载波在乘法器中相乘的结果为：t A t m t s w o cos )()('=，然后通过一个传输函数为H(f)的带通滤波器，得出已调信号为。 从图1.1中可得已调信号的时域和频域表达式为： （1-1） 式（1-1）中，M(f)为调制信号m(t)的频谱。 由于调制信号m(t)和乘法器输出信号之间是线性关系，所以成为线性调制。带通滤波器H(f)可以有不同的设计，从而得到不同的调制种类。 1.2双边带调制DSB 的基本原理 在幅度调制的一般模型中，若假设滤波器为全通网络，调制信号m(t)中无直流分量，则输出的已调信号就是无载波分量的双边带调制信号，或称抑制载波双边带（DSB ）调制信号，简称双边带（DSB ）信号。 设正弦型载波c(t)=Acos( t) ，式中：A 为载波幅度， 为载波角频率。 根据调制定义，幅度调制信号（已调信号）一般可表示为： (t)=Am(t)cos(t) （1-2） ?? ???-++==) ()]()([21)()(*]cos )([)(f H f f M f f M f s t h t t m t s o o o w m(t) H(t) A os t w o c s(t) )(' t s
其中，m(t)为基带调制信号。 设调制信号m(t)的频谱为M()，则由公式2-2不难得到已调信号 (t)的频谱： )]()([2 )(c c m M M A s ωωωωω-++= （1-3） 由以上表示式可见，在波形上，幅度已调信号随基带信号的规律呈正比地变化；在频谱结构上，它的频谱完全是基带信号频谱在频域的简单搬移。 标准振幅就是常规双边带调制，简称调幅（AM ）。假设调制信号m(t)的平均值为0，将其叠加一个直流偏量 后与载波相乘，即可形成调幅信号。其时域表达式为: )cos())(()(0t t m t c AM A s ω+= （1-4） 式中： 为外加的直流分量；m(t)可以是确知信号，也可以是随机信号。 若为确知信号，则AM 信号的频谱为: （1-5） AM 信号的频谱由载频分量、上边带、下边带三部分组成。AM 信号的总功率包括载波功率和边带功率两部分。只有边带功率才与调制信号有关，也就是说，载波分量并不携带信息。因此，AM 信号的功率利用率比较低。 AM 调制器模型如下图所示。 图1.2 AM 调制器模型 AM 信号的时域和频域表达式分别为 （1-6） （1-7） 式中，A o 为外加的直流分量；m(t)可以是确知信号也可以是随机信号，但通常认为其平均值为0，即0)(=t m — 。 由频谱可以看出，AM 信号的频谱由载波分量、上边带、下边带三部分组成。上边带的频谱结构与原调制信号的频谱结构相同，下边带是上边带的镜像。因此，AM 信号是带有载波 分量的双边带信号，他的带宽是基带信号带宽 的2倍，即 ) (cos )()(cos ) (cos )]([)(t w c t m t w c A t w c t m A o t s o AM +=+=)]()([2 1)]()([)(w c w M w c w M w c w w c w A o t s AM -+++-++=δδπ)] ()([2 1)]()([)(0 ω ω ω ω ωωωδωδπωc c c c m M M A s -+++-++=f H
实验一小球做自由落体运动内容：一小球竖直方向做自由落体，并无损做往返运动。程序： theta=0:0.01:2*pi x=cos(theta) y=sin(theta) l=1 v=1 while l<10 for t=1:10 y=y+(-1)^l*v*t plot(x,y,[-1,1],[-56,2],'.') axis equal pause(0.1) end l=l+1 end 结果：
-50 -40 -30 -20 -10 收获：通过运用小球自由落体规律，及（-1）^n 来实现无损往 返运动！ 实验二 旋转五角星 内容：一个五角星在圆内匀速旋转 程序：x=[2 2 2 2 2 2] y=[0 4/5*pi 8/5*pi 2/5*pi 6/5*pi 0] y1=2*sin(y) x1=2*cos(y) theta=0:4/5*pi:4*pi
x2=2*cos(theta) y2=2*sin(theta) plot(x,y,x1,y1,x2,y2) axis equal theta1=theta+pi/10 x2=2*cos(theta1) y2=2*sin(theta1) plot(x2,y2) axis equal theta=0:4/5*pi:4*pi for rot=pi/10:pi/10:2*pi x=2*cos(theta+rot) y=2*sin(theta+rot) plot(x,y) pause(0.1) end 结果：
-2 -1.5-1-0.500.51 1.52 -2-1.5-1-0.500.511.5 2 收获：通过theta1=theta+pi/10，我们可以实现五角星在圆内匀速 旋转！ 实验三 转动的自行车 内容：一辆自行车在圆内匀速转动 程序：x=-4:0.08:4; y=sqrt(16-x.^2); theta1=-pi/2:0.01*pi:3*pi/2; x3=0.5*cos(theta1); y3=0.5*sin(theta1); theta=-pi/2+0.02*pi for k=1:100
MATLAB/Simulink 电力系统建模与仿真 实验报告 姓名：****** 专业：电气工程及其自动化 班级：******************* 学号：*******************
实验一无穷大功率电源供电系统三相短路仿真 1.1 无穷大功率电源供电系统仿真模型构建 运行MATLAB软件，点击Simulink模型构建，根据电路原理图，添加下列模块： （1）无穷大功率电源模块（Three-phase source） （2）三相并联RLC负荷模块（Three-Phase Parallel RLC Load） （3）三相串联RLC支路模块（Three-Phase Series RLC Branch） （4）三相双绕组变压器模块（Three-Phase Transformer (Two Windings)） （5）三相电压电流测量模块（Three-Phase V-I Measurement） （6）三相故障设置模块（Three-Phase Fault） （7）示波器模块（Scope） （8）电力系统图形用户界面（Powergui） 按电路原理图连接线路得到仿真图如下： 1.2 无穷大功率电源供电系统仿真参数设置 1.2.1 电源模块 设置三相电压110kV，相角0°，频率50Hz，接线方式为中性点接地的Y形接法，电源电阻0.00529Ω，电源电感0.000140H，参数设置如下图：
1.2.2 变压器模块 变压器模块参数采用标幺值设置，功率20MVA，频率50Hz，一次测采用Y型连接，一次测电压110kV，二次侧采用Y型连接，二次侧电压11kV，经过标幺值折算后的绕组电阻为0.0033，绕组漏感为0.052，励磁电阻为909.09，励磁电感为106.3，参数设置如下图： 1.2.3 输电线路模块 根据给定参数计算输电线路参数为：电阻8.5Ω，电感0.064L，参数设置如下图： 1.2.4 三相电压电流测量模块 此模块将在变压器低压侧测量得到的电压、电流信号转变成Simulink信号，相当于电压、电流互感器的作用，勾选“使用标签（Use a label）”以便于示波器观察波形，设置电压标签“Vabc”，电流标签“Iabc”，参数设置如下图：
MATLAB基本操作 实验报告 课程名称： 院系： 专业班级： 学号： 学生姓名： 指导教师： 开课时间：至学年第学期
一、学生撰写要求 按照实验课程培养方案的要求，每门实验课程中的每一个实验项目完成后，每位参加实验的学生均须在实验教师规定的时间内独立完成一份实验报告，不得抄袭，不得缺交。 学生撰写实验报告时应严格按照本实验报告规定的内容和要求填写。字迹工整，文字简练，数据齐全，图表规范，计算正确，分析充分、具体、定量。 二、教师评阅与装订要求 1.实验报告批改要深入细致，批改过程中要发现和纠正学生实验报告中的问题，给出评语和实验报告成绩，签名并注明批改日期。实验报告批改完成后，应采用适当的形式将学生实验报告中存在的问题及时反馈给学生。 2.实验报告成绩用百分制评定，并给出成绩评定的依据或评分标准（附于实验报告成绩登记表后）。对迟交实验报告的学生要酌情扣分，对缺交和抄袭实验报告的学生应及时批评教育，并对该次实验报告的分数以零分处理。对单独设课的实验课程，如学生抄袭或缺交实验报告达该课程全学期实验报告总次数三分之一以上，不得同意其参加本课程的考核。 3.各实验项目的实验报告成绩登记在实验报告成绩登记表中。本学期实验项目全部完成后，给定实验报告综合成绩。 4.实验报告综合成绩应按课程教学大纲规定比例（一般为10-15%）计入实验课总评成绩；实验总评成绩原则上应包括考勤、实验报告、考核（操作、理论）等多方面成绩； 5.实验教师每学期负责对拟存档的学生实验报告按课程、学生收齐并装订，按如下顺序装订成册：实验报告封面、实验报告成绩登记表、实验报告成绩评定依据、实验报告（按教学进度表规定的实验项目顺序排序）。装订时统一靠左侧按“两钉三等分”原则装订。
实验一 Matlab基础知识 一、实验目的： 1.熟悉启动和退出Matlab的方法。 2.熟悉Matlab命令窗口的组成。 3.掌握建立矩阵的方法。 4.掌握Matlab各种表达式的书写规则以及常用函数的使 用。 二、实验内容： 1.求[100，999]之间能被21整除的数的个数。（rem） 2.建立一个字符串向量，删除其中的大写字母。（find） 3.输入矩阵，并找出其中大于或等于5的元素。（find） 4.不采用循环的形式求出和式 63 1 2i i= ∑ 的数值解。(sum) 三、实验步骤： ●求[100，199]之间能被21整除的数的个数。（rem） 1.开始→程序→Matlab 2.输入命令： ?m=100:999; ?p=rem(m,21); ?q=sum(p==0) ans=43 ●建立一个字符串向量，删除其中的大写字母。（find） 1.输入命令：
?k=input('’,’s’)； Eie48458DHUEI4778 ?f=find(k>=’A’&k<=’Z’)； f=9 10 11 12 13 ?k(f)=[ ] K=eie484584778 ●输入矩阵，并找出其中大于或等于5的元素。（find） 1.输入命令： ?h=[4 8 10;3 6 9; 5 7 3]； ?[i,j]=find(h>=5) i=3 j=1 1 2 2 2 3 2 1 3 2 3 ●不采用循环的形式求出和式的数值解。（sum） 1.输入命令： ?w=1:63； ?q=sum(2.^w) q=1.8447e+019
实验二 Matlab 基本程序 一、 实验目的： 1. 熟悉Matlab 的环境与工作空间。 2. 熟悉M 文件与M 函数的编写与应用。 3. 熟悉Matlab 的控制语句。 4. 掌握if,switch,for 等语句的使用。 二、 实验内容： 1. 根据y=1+1/3+1/5+……+1/（2n-1）,编程求：y<5时最大n 值以及对应的y 值。 2. 编程完成，对输入的函数的百分制成绩进行等绩转换，90~100为优，80~89为良，70~79为中，60~69为及格。 3. 编写M 函数文件表示函数 ，并分别求x=12和56时的函数值。 4. 编程求分段函数 2226;03 56;0532 1;x x x x y x x x x x x x +-<≠=-+≤<≠≠-+且且及其它，并求输入x=[-5.0,-3.0,1.0,2.0,2.5,3.0,3.5]时的输出y 。 三、 实验步骤： 根据y=1+1/3+1/5+……+1/（2n-1）,编程求：y<5时最大n 值以及对应的y 值。 1. 打开Matlab ，新建M 文件 2. 输入命令： 51022-+x
实验报告 课程名称：控制理论（乙） 指导老师：林峰 成绩：__________________ 实验名称：MATLAB 仿真实验 实验类型：________________同组学生姓名：__________ 一、实验目的和要求（必填） 二、实验内容和原理（必填） 三、主要仪器设备（必填） 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析（必填） 七、讨论、心得 实验九 控制系统的时域分析 一、 实验目的： 1．用计算机辅助分析的办法，掌握系统的时域分析方法。 2．熟悉Simulink 仿真环境。 二、实验原理及方法： 系统仿真实质上就是对系统模型的求解，对控制系统来说，一般模型可转化成某个微分方程或差分方程表示，因此在仿真过程中，一般以某种数值算法从初态出发，逐步计算系统的响应，最后绘制出系统的响应曲线，进而可分析系统的性能。控制系统最常用的时域分析方法是，当输入信号为单位阶跃和单位冲激函数时，求出系统的输出响应，分别称为单位阶跃响应和单位冲激响应。在MATLAB 中，提供了求取连续系统的单位阶跃响应函数step ，单位冲激响应函数impulse ，零输入响应函数initial 等等。 二、实验内容： 二阶系统，其状态方程模型为 ? 1x -0.5572 -0.7814 1x 1 = + u ? 2x 0.7814 0 2x 0 1x y = [1.9691 6.4493] +[0] u 2x 四、实验要求： 1．编制MATLAB 程序，画出单位阶跃响应曲线、冲击响应曲线、系统的零输入响应、斜坡输入响应； （1）画出系统的单位阶跃响应曲线； A=[-0.5572 -0.7814;0.7814 0 ]; B=[1;0];
Tutorial 5 实验报告 实验名称：Matlab 绘图操作 实验目的： 1、 掌握绘制二维图形的常用函数； 2、 掌握绘制三维图形的常用函数； 3、 掌握绘制图形的辅助操作。 实验内容： 1. 设sin .cos x y x x ?? =+ ??+?? 23051，在x=0~2π区间取101点，绘制函数的曲线。 2. 已知： y x =21，cos()y x =22，y y y =?312，完成下列操作： （1） 在同一坐标系下用不同的颜色和线性绘制三条曲线； （2） 以子图形式绘制三条曲线； （3） 分别用条形图、阶梯图、杆图和填充图绘制三条曲线。 3. 已知：ln(x y x x ≤=??+>??0102 ，在x -≤≤55区间绘制函数曲线。 4. 绘制极坐标曲线sin()a b n ρθ=+，并分析参数a 、b 、n 对曲线形状的影响。 5．在xy 平面内选择区域[][],,-?-8888， 绘制函数z =的三种三维曲面图。 6. 用plot 函数绘制下面分段函数的曲线。 ,(),,x x f x x x x ?+>? ==??+
8. 在同一坐标轴中绘制下列两条曲线。 （1）.y x =-205 （2）sin()cos ,sin()sin x t t t y t t π=?≤≤? =?303 实验结果： 1. 2. （1）
（2）
（3）
实验一 模拟通信的MATLAB 仿真 姓名：左立刚 学号：031040522 简要说明： 实验报告注意包括AM ，DSB ，SSB ，VSB ，FM 五种调制与解调方式的实验原理，程序流程图，程序运行波形图，simulink 仿真模型及波形，心得体会，最后在附录中给出了m 语言的源程序代码。 一．实验原理 1．幅度调制（AM ） 幅度调制（AM ）是指用调制信号去控制高频载波的幅度，使其随调制信号呈线性变化的过程。AM 信号的数学模型如图3-1所示。 图2-1 AM 信号的数学模型 为了分析问题的方便，令 δ =0， 1.1 AM 信号的时域和频域表达式 ()t S AM =[A 0 +m ()t ]cos t c ω （2-1） ()t S AM =A 0 π[()()ωωωωδC C ++-]+()()[]ωωωωc c M M ++-2 1 （2-2）
AM 信号的带宽 2 =B AM f H （2-3） 式中， f H 为调制信号的最高频率。 2.1.3 AM 信号的功率P AM 与调制效率 η AM P AM =()222 2 t m A +=P P m c + （2-4） 式中，P C =2 A 为不携带信息的载波功率；()2 2 t m P m =为携带信息的边带 功率。 ()() t t m A m P P AM C AM 2 2 2+= = η （2-5） AM 调制的优点是可用包络检波法解调，不需要本地同步载波信号，设备简单。AM 调制的最大缺点是调制效率低。 2.2、双边带调制（DSB ） 如果将在AM 信号中载波抑制，只需在图3-1中将直流 A 0 去掉，即可输出 抑制载波双边带信号。 2.2.1 DSB 信号的时域和频域表达式 ()()t t m t c DSB S ωcos = （2-6） ()()()[]ωωωωωC C DSB M M S ++-=2 1 （2-7） DSB 信号的带宽 f B B H AM DSB 2 == （2-8）