目录
一、背景 (4)
二、基本要求 (4)
三、设计概述 (4)
四、Matlab设计流程图 (5)
五、Matlab程序及仿真结果图 (6)
1、生成m序列及m序列性质 (6)
2、生成50位随机待发送二进制比特序列，并进行扩频编码 (7)
3、对扩频前后信号进行BPSK调制，观察其时域波形 (9)
4、计算并观察扩频前后BPSK调制信号的频谱 (10)
5、仿真经awgn信道传输后，扩频前后信号时域及频域的变化 (11)
6、对比经信道前后两种信号的频谱变化 (12)
7、接收机与本地恢复载波相乘，观察仿真时域波形 (14)
8、与恢复载波相乘后，观察其频谱变化 (15)
9、仿真观察信号经凯萨尔窗低通滤波后的频谱 (16)
10、观察经过低通滤波器后无扩频与扩频系统的时域波形 (17)
11、对扩频系统进行解扩，观察其时域频域 (18)
12、比较扩频系统解扩前后信号带宽 (19)
13、比较解扩前后信号功率谱密度 (20)
14、对解扩信号进行采样、判决 (21)
15、在信道中加入2040~2050Hz窄带强干扰并乘以恢复载波 (24)
16、对加窄带干扰的信号进行低通滤波并解扩 (25)
17、比较解扩后信号与窄带强干扰的功率谱 (27)
六、误码率simulink仿真 (28)
1、直接扩频系统信道模型 (28)
2、加窄带干扰的直扩系统建模 (29)
3、用示波器观察发送码字及解扩后码字 (30)
4、直接扩频系统与无扩频系统的误码率比较 (31)
5、不同扩频序列长度下的误码率比较 (32)
6、扩频序列长度N=7时，不同强度窄带干扰下的误码率比较 (33)
七、利用Walsh码实现码分多址技术 (34)
1、产生改善的walsh码 (35)
2、产生两路不同的信息序列 (36)
3、用两个沃尔什码分别调制两路信号 (38)
4、两路信号相加，并进行BPSK调制 (39)
5、观察调制信号频谱，并经awgn信道加高斯白噪和窄带强干扰 (40)
6、接收机信号乘以恢复载波，观察时域和频域 (42)
7、信号经凯萨尔窗低通滤波器 (43)
8、对滤波后信号分别用m1和m2进行解扩 (44)
9、对两路信号分别采样，判决 (45)
八、产生随机序列Gold码和正交Gold码 (47)
1、产生Gold码并仿真其自相关函数 (48)
2、产生正交Gold码并仿真其互相关函数 (50)
九、实验心得体会 (51)
直接序列扩频系统仿真
一、背景
直接序列扩频通信系统(DSSS)是目前应用最为广泛的系统。在发送端，直扩系统将发送序列用伪随机序列扩展到一个很宽的频带上去，在接受端又用相同的扩频序列进行解扩，回复出原有信息。由于干扰信息与伪随机序列不相关，扩频后能够使窄带干扰得到有效的抑制，提高输出信噪比。系统框图如下图所示：
二、基本要求：
1.通过matlab建模，对直扩系统进行仿真，数据调制方式可以自由选择，可以
使用基带信号，但最好能使用频带信号，信道为高斯白噪信道。要仿真出扩频前的信号的频偏，扩频后的信号频谱，过信道之后的频谱以及解扩之后的频谱。
2.研究并仿真产生m序列，写出生成m序列的算法。
3.验证直扩系统对窄带干扰的抑制能力，在信道中加入一个窄带强干扰，仿真
出加了干扰后的频谱图和解扩后的频谱图，给出误码率等仿真图。
4.在以上基础上仿真实现码分多址技术，使用Walsh码进行复用，实现多个信
号同时传输。(选做)
可选项：
1.在信道中加入多径，使用rake接收来抗多径效应。
2.产生除m序列之外的其他随机序列，如Gold码，正交Gold码等等。
3.对比无扩频的系统的误码率。
三、设计概述
本次课设完成基本要求，并选作了可选项码分多址，Gold码及误码率对比。通过matlab建模仿真了直扩系统BPSK调制的各点频偏及时域信号，并仿真了窄带强干扰对直扩系统的影响以及利用改善的WALSH码实现码分多址技术。另外，通过matlab的simulink工具盒bertool工具仿真对比了直扩系统和无扩频系统的误码率。
四、matlab设计流程图
基本扩频系统仿真流程图
五、matlab程序及仿真结果图
1、生成m序列及m序列性质
实验产生7位m序列，频率100Hz，模拟线性反馈移位寄存器序列，原理图如下：
clear all;
clc;
X1=0;X2=0;X3=1;
m=350; %重复50遍的7位单极性m序列
for i=1:m
Y3=X3; Y2=X2; Y1=X1;
X3=Y2; X2=Y1;
X1=xor(Y3,Y1);
L(i)=Y1;
end
for i=1:m
M(i)=1-2*L(i); %将单极性m序列变为双极性m序列
end
k=1:1:m;
figure(1)
subplot(3,1,1) %做m序列图
stem(k-1,M);
axis([0,7,-1,1]);
xlabel('k');
ylabel('M序列');
title('移位寄存器产生的双极性7位M序列') ;
subplot(3,1,2)
ym=fft(M,4096);
magm=abs(ym); %求双极性m序列频谱
fm=(1:2048)*200/2048;
plot(fm,magm(1:2048)*2/4096);
title('双极性7位M序列的频谱')
axis([90,140,0,0.1]);
[a,b]=xcorr(M,'unbiased');
subplot(3,1,3) %求双极性m序列自相关函数
plot(b,a);
axis([-20,20,-0.5,1.2]);
title('双极性7位M序列的自相关函数');
由上图可以看出，7位m序列为1，-1，-1，-1，1，-1，1。另外，自相关函数的图形比较尖锐，最大值为1，最小值为-1/7，符合理论结果。
2、生成50位随机待发送二进制比特序列，并进行扩频编码
生成的信息码频率为100/7Hz，利用m序列编码后，频率变为100Hz。
N=50;a=0;
x_rand=rand(1,N); %产生50个0与1之间随机数
for i=1:N
if x_rand(i)>=0.5 %大于等于0.5的取1，小于0.5的取0
x(i)=1;a=a+1;
else x(i)=0;
end
end
t=0:N-1;
figure(2) %做信息码图
subplot(2,1,1)
stem(t,x);
title('扩频前待发送二进制信息序列');
tt=0:349;
subplot(2,1,2)
l=1:7*N;
y(l)=0;
for i=1:N
k=7*i-6;
y(k)=x(i);
k=k+1;y(k)=x(i);k=k+1;y(k)=x(i);k=k+1;y(k)=x(i);k=k+1;y(k)=x(i);k=k+1;y(k)=x(i); k=k+1;y(k)=x(i);
end
s(l)=0;
for i=1:350 %扩频后，码率变为100/7*7=100Hz
s(i)=xor(L(i),y(i));
end
tt=0:7*N-1;
stem(tt,s);
axis([0,350,0,1]);
title('扩频后的待发送序列码');
3、对扩频前后信号进行BPSK调制，观察其时域波形
BPSK调制采用2kHz信号cos(2*2000*t)作为载波
figure(3)
subplot(2,1,2)
fs=2000;
ts=0:0.00001:3.5-0.00001;%为了使信号看起来更光滑，作图时采样频率为100kHz % ps=cos(2*pi*fs*ts);
s_b=rectpulse(s,1000); %将冲激信号补成矩形信号
s_bpsk=(1-2.*s_b).*cos(2*pi*fs*ts);%扩频后信号BPSK调制时域波形
plot(ts,s_bpsk);
xlabel('s');
axis([0.055,0.085,-1.2,1.2])
title('扩频后bpsk信号时域波形');
subplot(2,1,1)
s_bb=rectpulse(x,7000);
s_bpskb=(1-2.*s_bb).*cos(2*pi*fs*ts);%无扩频信号BPSK调制时域波形
plot(ts,s_bpskb);
xlabel('s');
axis([0.055,0.085,-1.2,1.2]);
title('扩频前bpsk信号时域波形')
可以看出，100/7Hz的无扩频信号每0.07s时由于序列极性变换产生相位变换，100Hz的扩频后调制信号每0.01s由于序列极性变换产生相位变换。
课题一： 连续时间信号和系统时域分析及MATLAB实现 课题要求： 深入研究连续时间信号和系统时域分析的理论知识。利用MATLAB强大的图形处理功能、符号运算功能以及数值计算功能，实现连续时间信号和系统时域分析的仿真波形。 课题内容： 一、用MATLAB实现常用连续时间信号的时域波形（通过改变参数，分析其时域特性）。 1、单位阶跃信号， 2、单位冲激信号， 3、正弦信号， 4、实指数信号， 5、虚指数信号， 6、复指数信号。 二、用MATLAB实现信号的时域运算 1、相加， 2、相乘， 3、数乘， 4、微分， 5、积分 三、用MATLAB实现信号的时域变换（参数变化，分析波形变化） 1、反转， 2、使移（超时，延时）， 3、展缩， 4、倒相， 5、综合变化 四、用MATLAB实现信号简单的时域分解 1、信号的交直流分解， 2、信号的奇偶分解 五、用MATLAB实现连续时间系统的卷积积分的仿真波形 给出几个典型例子，对每个例子，要求画出对应波形。 六、用MATLAB实现连续时间系统的冲激响应、阶跃响应的仿真波形。 给出几个典型例子，四种调用格式。 七、利用MATLAB实现连续时间系统对正弦信号、实指数信号的零状态响应的仿真波形。 给出几个典型例子，要求可以改变激励的参数，分析波形的变化。 课题二： 离散时间信号和系统时域分析及MATLAB实现。 课题要求： 深入研究离散时间信号和系统时域分析的理论知识。利用MATLAB强大的图
形处理功能、符号运算功能以及数值计算功能，实现离散时间信号和系统时域分析的仿真波形。 课题内容： 一、用MATLAB绘制常用信号的时域波形（通过改变参数分析其时域特性） 1、单位序列， 2、单位阶跃序列， 3、正弦序列， 4、离散时间实指数序列， 5、离散时间虚指数序列， 6、离散时间复指数序列。 二、用MATLAB实现信号的时域运算 1、相加， 2、相乘， 3、数乘。 三、用MATLAB实现信号的时域变换（参数变化，分析波形的变化） 1、反转， 2、时移（超时，延时）， 3、展缩， 4、倒相。 四、用MATLAB实现离散时间系统卷积和仿真波形 给出几个典型例子，对每个例子要求画出e(k)，h(k)，e(i)，h(i)，h(-i)，Rzs(k)波形。 五、用MATLAB实现离散时间系统的单位响应，阶跃响应的仿真波形 给出几个典型例子，四中调用格式。 六、用MATLAB实现离散时间系统对实指数序列信号的零状态响应的仿真波形 给出几个典型例子，要求可以改变激励的参数，分析波形的变化。 课题三： 连续时间信号傅里叶级数分析及MATLAB实现。 课题要求： 深入研究连续时间信号傅里叶级数分析的理论知识，利用MATLAB强大的图形处理功能，符号运算功能以及数值计算功能，实现连续时间周期信号频域分析的仿真波形。 课题内容： 一、用MATLAB实现周期信号的傅里叶级数分解与综合 以周期矩形波信号为例，绘出包含不同谐波次数的合成波形，观察合成波形与原矩形 波形之间的关系及吉布斯现象。
数学实验报告 班级: 学号： 姓名： 实验序号：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);
成都理工大学工程技术学院 《通信仿真课程设计》报告 班级：信息工程1班 姓名：寇路军 学号： 201620101133 指导教师：周玲 成绩： 2019 年 3月 23 日
目录 通信仿真课程设计报告 (2) 一.绪论 (2) 二．课程设计的目的 (2) 三．模拟调制系统的设计 (3) 3.1 二进制相移键控调制基本原理 (3) 3.2 2PSK信号的调制 (3) 3.2.1模拟调制的方法 (3) 3.3 2PSK信号的解调 (4) 3.4 2PSK的“倒∏现象”或“反向工作” (5) 3.5功率谱密度 (5) 四．数字调制技术设计 (7) 4.1 2PSK的仿真 (7) 4.1.1仿真原理图 (7) 4.1.2 仿真数据 (7) 4.1.3 输出结果 (9) 总结 (10) 参考文献 (11)
通信仿真课程设计报告 一.绪论 随着社会的快速发展，通信系统在社会上表现出越来越重要的作用。目前，我们生活中使用的手机，电话，Internet，ATM机等通信设备都离不开通信系统。随着通信系统与我们生活越来越密切，使用越来越广泛，对社会对通信系统的性能也越高。另外，随着人们对通信设备更新换代速度越来越快。不得不缩短通信系统的开发周期以及提高系统性能。针对这两方面的要求，必需要通过强大的计算机辅助分析设计技术和工具才能实现。自从现代以来，计算机科技走上了快速发展道路，实现了可视化的仿真软件。 通信系统仿真，在目前的通信系统工程设计当中。已成为了不可替代的一部分。它表现出很强的灵活性和适应性。为我们更好地研究通信系统性能带来了很大的帮助。本论文主要针对模拟调制系统中的二进制相移键控调制技术进行设计和基于Simulink进行仿真。通过系统仿真验证理论中的结论。本论文设计的目的之一是进一步加强理论知识，熟悉Matlab软件。 Simulink是MATLAB最重要的组件之一，它提供一个动态系统建模、仿真和综合分析的集成环境。在该环境中，无需大量书写程序，而只需要通过简单直观的鼠标操作，就可构造出复杂的系统。Simulink具有适应面广、结构和流程清晰及仿真精细、贴近实际、效率高、灵活等优点，并基于以上优点Simulink 已被广泛应用于控制理论和数字信号处理的复杂仿真和设计。同时有大量的第三方软件和硬件可应用于或被要求应用于Simulink。 二．课程设计的目的 1.掌握模拟系统2PSK调制和解调原理及设计方法。 2.熟悉基于Simulink的通信系统仿真。
课程设计任务书 学生姓名：董航专业班级：电信1006班 指导教师：阙大顺，李景松工作单位：信息工程学院 课程设计名称：Matlab应用课程设计 课程设计题目：Matlab运算与应用设计5 初始条件： 1.Matlab6.5以上版本软件； 2.课程设计辅导资料：“Matlab语言基础及使用入门”、“Matlab及在电子信息课程中的应 用”、线性代数及相关书籍等； 3.先修课程：高等数学、线性代数、电路、Matlab应用实践及信号处理类相关课程等。 要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求） 1.课程设计内容：根据指导老师给定的7套题目，按规定选择其中1套完成； 2.本课程设计统一技术要求：研读辅导资料对应章节，对选定的设计题目进行理论分析， 针对具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表），并对实验结果进行分析和总结。具体设计要求包括： ①初步了解Matlab、熟悉Matlab界面、进行简单操作； ②MATLAB的数值计算：创建矩阵矩阵运算、多项式运算、线性方程组、数值统计； ③基本绘图函数：plot, plot3, mesh, surf等，要求掌握以上绘图函数的用法、简单图形 标注、简单颜色设定等； ④使用文本编辑器编辑m文件，函数调用； ⑤能进行简单的信号处理Matlab编程； ⑥按要求参加课程设计实验演示和答辩等。 3.课程设计说明书按学校“课程设计工作规范”中的“统一书写格式”撰写，具体包括： ①目录； ②与设计题目相关的理论分析、归纳和总结； ③与设计内容相关的原理分析、建模、推导、可行性分析； ④程序设计框图、程序代码（含注释）、程序运行结果和图表、实验结果分析和总结； ⑤课程设计的心得体会（至少500字）； ⑥参考文献（不少于5篇）； ⑦其它必要内容等。 时间安排：1.5周（分散进行） 参考文献： [1](美)穆尔,高会生,刘童娜,李聪聪．MA TLAB实用教程(第二版) . 电子工业出版社,2010. [2]王正林,刘明．精通MATLAB(升级版) .电子工业出版社,2011. [3]陈杰. MA TLAB宝典(第3版) . 电子工业出版社,2011. [4]刘保柱,苏彦华,张宏林. MATLAB 7.0从入门到精通(修订版) . 人民邮电出版社,2010. 指导教师签名：年月日 系主任（或责任教师）签名：年月日
实验一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
目录 第一章课程设计的任务说明 (1) 1.1 课程设计的目的 (1) 1.2 课程设计的要求 (1) 第二章MA TLAB/SIMULINK简介 (2) 第三章通信技术的历史和发展 (3) 3.1通信的概念 (3) 3.2 通信的发展史简介 (4) 3.3通信技术的发展现状和趋势 (4) 第四章2ASK的基本原理和实现 (5) 4.1 2ASK的产生 (5) 4.2 2ASK的功率谱和带宽 (6) 4.3 2ASK信号的解调及抗噪声性能分析 (6) 第五章Smulink的模型建立和仿真 (10) 5.1 建立模型方框图 (10) 5.2参数设置 (11) 5.3仿真波形图 (15) 5.4 不同信噪比下的误码率 (17) 总结 (18) 参考文献 (19)
第一章课程设计的任务说明 1.1 课程设计的目的 （1）通过利用matlab simulink，熟悉matlab simulink仿真工具。 （2）通过课程设计来更好的掌握课本相关知识，熟悉2ASK的调制与解调。 （3）更好的了解通信原理的相关知识，磨练自己分析问题、查阅资料、巩固知识、创新等各方面能力。 1.2 课程设计的要求 （1）掌握课程设计涉汲到的相关知识，相关概念、原理清晰，明了。 （2）仿真图设计合理、能够正确运行。 （3）按照要求撰写课程设计报告。
第二章MATLAB/SIMULINK简介 美国Mathworks公司于1967年推出了矩阵实验室“Matrix Laboratory”（缩写为Matlab）这就是Matlab最早的雏形。开发的最早的目的是帮助学校的老师和学生更好的授课和学习。从Matlab诞生开始，由于其高度的集成性及应用的方便性，在高校中受到了极大的欢迎。由于它使用方便，能非常快的实现科研人员的设想，极大的节约了科研人员的时间，受到了大多数科研人员的支持，经过一代代人的努力，目前已发展到了7.X版本。Matlab是一种解释性执行语言，具有强大的计算、仿真、绘图等功能。由于它使用简单，扩充方便，尤其是世界上有成千上万的不同领域的科研工作者不停的在自己的科研过程中扩充Matlab的功能，使其成为了巨大的知识宝库。可以毫不夸张的说，哪怕是你真正理解了一个工具箱，那么就是理解了一门非常重要的科学知识。科研工作者通常可以通过Matlab来学习某个领域的科学知识，这就是Matlab真正在全世界推广开来的原因。目前的Matlab版本已经可以方便的设计漂亮的界面，它可以像VB等语言一样设计漂亮的用户接口，同时因为有最丰富的函数库（工具箱），所以计算的功能实现也很简单，进一步受到了科研工作者的欢迎。另外，,Matlab和其他高级语言也具有良好的接口，可以方便的实现与其他语言的混合编程，进一步拓宽了Matlab的应用潜力。可以说，Matlab已经也很有必要成为大学生的必修课之一，掌握这门工具对学习各门学科有非常重要的推进作用。 Simulink是MA TLAB中的一种可视化仿真工具，也是目前在动态系统的建模和仿真等方面应用最广泛的工具之一。确切的说，Simulink是一个用来对动态系统进行建模、仿真和分析的软件包，它支持线性和非线性系统，连续、离散时间模型，或者是两者的混合。系统还可以使多种采样频率的系统，而且系统可以是多进程的。Simulink工作环境进过几年的发展，已经成为学术和工业界用来建模和仿真的主流工具包。在Simulink环境中，它为用户提供了方框图进行建模的图形接口，采用这种结构画模型图就如同用手在纸上画模型一样自如、方便，故用户只需进行简单的点击和拖动就能完成建模，并可直接进行系统的仿真，快速的得到仿真结果。它的主要特点在于：1、建模方便、快捷；2、易于进行模型分析；3、优越的仿真性能。它与传统的仿真软件包微分方程和差分方程建模相比，具有更直观、方便、灵活的优点。Simulink模块库（或函数库）包含有Sinks（输出方式）、Sources（输入源）、Linear（线性环节）、Nonlinear（非线性环节）、Connection（连接与接口）和Extra（其他环节）等具有不同功能或函数运算的Simulink库模块（或库函数），而且每个子模型库中包含有相应的功能模块，用户还可以根据需要定制和创建自己的模块。用Simulink创建的模型可以具有递阶结构，因此用户可以采用从上到下或从下到上的结构创建模型。用户可以从最高级开始观看模型，然后用鼠标双击其中的子系统模块，来查看其下一级的内容，以此类推，从而可以看到整个模型的细节，帮助用户理解模型的结构和各模块之间的相互关系。在定义完一个模型后，用户可以通过Simulink的菜单或MATLAB的命令窗口键入命令来对它进行仿真。菜单方式对于交互工作非常方便，而命令行方式对于运行仿真的批处理非常有用。采用Scope模块和其他的显示模块，可以在仿真进行的同时就可立即观看到仿真结果，若改变模块的参数并再次运行即可观察到相应的结果，这适用于因果关系的问题研究。仿真的结果还可以存放到MATLAB的工作空间里做事后处理。模型分析工具包括线性化和整理工具，MATLAB的所有工具及Simulink本身的应用工具箱都包含这些工具。由于MATLAB和SIMULINK的集成在一起的，因此用户可以在这两种环境下对自己的模型进行仿真、分析和修改模型。但是Simulink不能脱离MA TLAB而独立工作。
自控系统仿真软件课程设计报告 MATLAB 设计题目：牛顿摆球 姓名： 学号： 院系： 班级：1203 指导教师： 2014年12月20日
一．课程设计目的 1、熟悉课程设计的基本流程； 2、掌握MATLAB语法结构及调试方法； 3、熟悉MATLAB函数调用，熟练二维画图; 4、掌握MATLAB语言在控制方面的运用； 5、学会用MATLAB进行基本仿真； 6、掌握MATLAB编程技巧，提高编程水平。 二．系统分析 1.题目的描述： （1）牛顿摆球原理描述 五个质量相同的球体由吊绳固定，彼此紧密排列。当摆动最右侧的球并在回摆时碰撞紧密排列的另外四个球，最左边的球将被弹出，并仅有最左边的球被弹出。当然此过程也是可逆的，当摆动最左侧的球撞击其它球时，最右侧的球会被弹出。当最右侧的两个球同时摆动并撞击其他球时，最左侧的两个球会被弹出。同理相反方向同样可行，并适用于更多的球。 为了更接近现实，在这里我将考虑重力及空气阻力的影响，摆球将不会永无止境的运动下去，由于外界因素的影响，摆球运动一段时间后将回归静止状态。（2）通过MATLAB动画程序制作软件，实现下述过程 当运行程序时，把最右边的小球拉到一定的高度放下，让其碰撞其余四个小球，仅让最左边的小球被弹出，当最左边小球回摆碰撞其它球时，最右边小球又被弹出，如此循环。由于是非理想条件下，摆球的摆动幅度会随摆动次数的增加越来越小，直到静止。 时间停顿两秒，把右边两小球一起拉到一定高度放下，让其碰撞其余三个球，同样仅让左边两球被弹出，当球回摆再次碰撞时，最右边两球又被同时弹出，如此循环，因为外界因素的影响，最终五个球都会静止下来。 （3）整个实验看似简单，但要在MATLAB上完成这样一个动画过程，还是需要下点功夫，克服困难的。经过自己的努力，终于实现了整个过程，这也是一种不小的收获。 2.设计要求： （1）能够实现有阻尼摆动，即摆幅随摆动次数增加越来越小，直到静止。（2）能够让摆球弧线摆动。 三．系统设计 1.系统设计过程 （1）通过函数axis建立坐标系 （2）在坐标系范围内通过函数line画各个支架 （3）通过函数title添加标题“动量守恒实验”、函数text添加标注“牛顿摆球” （4）通过函数line画出五个球，并设定其初始位置，颜色,大小,线条的擦拭方式
实验一 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.
课程设计 基于MATLABsimulink的2FSK系统的仿真 电子与信息工程学院 信息与通信工程系
课程设计评分标准
基于MATLAB/simulink的2FSK系统的仿真 一、摘要 本文是基于matlab和simulink环境下对信号的调制与解调过程的仿真，通过仿真，对系统的误码率的分析，以及理论与仿真结果的比较， 二、关键字：
目录 1 背景知识 (1) 1.1通信简介 (1) 1.2仿真系统的简介： (2) 1.32FSK的调制与解调的原理： (3) 1.3.1 2FSK的产生 (4) 1.3.2 2FSK滤波器的解调及抗噪声性能 (6) 1.3.3 由相关调制解调的原理图 (9) 2 仿真系统模型的设计： (9) 2.1仿真框图 (9) 2.2仿真目的和意义： (9) 2.3仿真思路 (10) 2.4M文件和仿真结果 (10) 2.5 SIMULINK仿真模型图： (16) 2.6结果分析： (21) 2.6.1 Matlab仿真结果分析 (21) 2.6.2 (22) 3 心得体会： (22) 4 参考文献 (22)
1 背景知识 1.1 通信简介 通信就是克服距离上的障碍，从一地向另一地传递和交换消息。消息是信息源所产生的，是信息的物理表现，例如，语音、文字、数据、图形和图像等都是消息。消息有模拟消息（如语音、图像等）以及数字消息（如数据、文字等）之分。所有消息必须在转换成电信号（通常简称为信号）后才能在通信系统中传输。所以，信号是传输消息的手段，信号是消息的物质载体。 相应的信号可分为模拟信号和数字信号，模拟信号的自变量可以是连续的或离散的，但幅度是连续的，如电话机、电视摄像机输出的信号就是模拟信号。数字信号的自变量可以是连续的或离散的，但幅度是离散的，如电船传机、计算机等各种数字终端设备输出的信号就是数字信号。 通信的目的是传递消息，但对受信者有用的是消息中包含的有效内容，也即信息。消息是具体的、表面的，而信息是抽象的、本质的，且消息中包含的信息的多少可以用信息量来度量。 通信技术，特别是数字通信技术近年来发展非常迅速，它的应用越来越广泛。通信从本质上来讲就是实现信息传递功能的一门科学技术，它要将大量有用的信息无失真，高效率地进行传输，同时还要在传输过程中将无用信息和有害信息抑制掉。当今的通信不仅要有效地传递信息，而且还有储存、处理、采集及显示等功能，通信已成为信息科学技术的一个重要组成部分。 通信系统就是传递信息所需要的一切技术设备和传输媒质的总和，包括信息源、发送设备、信道、接收设备和信宿(受信者) ，它的一般模型如图1所示。 →→→→ 信息源发送设备信道接收设备受信者 ↑ 噪声源 图1 通信系统一般模型 通信系统可分为数字通信系统和模拟通信系统。数字通信系统是利用数字信号来传递消息的通信系统，其模型如图2所示，
第一类：单位转换 1.长度单位换算的设计与实现 2.面积单位换算的设计与实现 3.体积单位换算的设计与实现 4.容积单位换算的设计与实现 5.质量单位换算的设计与实现 6.时间单位换算的设计与实现 7.温度单位换算的设计与实现 7.压强单位换算的设计与实现 8.角度单位换算的设计与实现 8.功率单位换算的设计与实现 第二类：曲线绘制 1.直线的自动绘制和相关计算 2.椭圆的自动绘制和相关计算 3.双曲线的自动绘制和相关计算 4.抛物线的自动绘制和相关计算 5.心脏线的自动绘制和相关计算 6.渐开线的自动绘制和相关计算 7.滚圆线的自动绘制和相关计算 8.三叶玫瑰线的自动绘制和相关计算9.四叶玫瑰线的自动绘制和相关计 10.阿基米德螺线的自动绘制和相关计算第三类：曲面绘制 1.球面的自动绘制和相关计算 2.椭球面的自动绘制和相关计算 3.单叶双曲面的自动绘制和相关计算 4.双叶双曲面的自动绘制和相关计算 5.抛物面的自动绘制和相关计算 6.双曲抛物面的自动绘制和相关计算 7.双曲柱面的自动绘制和相关计算 8.椭圆柱面的自动绘制和相关计算 9.抛物柱面的自动绘制和相关计算 10.圆锥面的自动绘制和相关计算 第四类：线性回归 1.男士身高体重相关计算经验公式 2.女士身高体重相关计算经验公式 3.男士胖瘦等级的确定 4.女士胖瘦等级的确定 5.男士身高脚长相关计算经验公式 6.女士身高脚长相关计算经验公式 7.父子身高相关性研究 8.母子身高相关性研究 9.父女身高相关性研究 10.母女身高相关性研究 第五类：学习成绩 1.期末总评自动计算的设计与实现 2.成绩等级自动评定的设计与实现 3.成绩分段自动统计的设计与实现 4.成绩分布折线自动绘制的设计与实现 5.成绩自动统计分析的设计与实现 6.试卷分布自动分析的设计与实现 7.试卷难度自动分析的设计与实现 8.考试成绩名次自动生成的设计与实现
实验一小球做自由落体运动内容：一小球竖直方向做自由落体，并无损做往返运动。程序： 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”，参数设置如下图：
电子科技大学 UNIVERSITY OF ELECTRONIC SCIENCE AND TECHNOLOGY OF CHINA 激光原理与技术 课程设计 课程教师： 作者姓名： 学号：
题目一： 编程计算图示谐振腔的稳定性与光焦度1/F的关系。可取R1=∞, R2=∞, l1=250mm, l2=200mm。，用matlab程序计算光线在腔内的轨迹，演示腔稳定和不稳定时光线在腔内往返次数增加时光线轨迹。初始光线参数可以任意选择。 利用matlab编程如下： clear,clc L1=250;L2=200; R1=inf;R2=inf; syms d; T=[1,L1;0,1]*[1,0;-d,1]*[1,L2;0,1]*[1,0;-2/R2,1]*[1,L2;0,1]*[1 ,0;-d,1]*[1,L1;0,1]*[1,0;-2/R1,1]; A=T(1,1); B=T(1,2); C=T(2,1); D=T(2,2); h=(A+D)/2; ezplot(h,[0,0.012]) title('谐振腔的稳定性');xlabel('透镜光焦度D（/mm）');ylabel('等效 g1g2') 运行结果：
题目二： 和透镜上的模式半径与光焦度1/F的关系。 计算输出镜M 2 利用matlab编程如下： clear,clc L1=250;L2=200;R1=inf;R2=inf;w1=0.5*10^-3; syms d T1=[1,L2;0,1]*[1,0;-2/R2,1]*[1,L2;0,1]*[1,0;-d,1]*[1,L1;0,1]*[1,0;-2/ R1,1]*[1,L1;0,1]*[1,0;-d,1]; A1=T1(1,1);B1=abs(T1(1,2));C1=T1(2,1);D1=T1(2,2);h1=(A1+D1)/2; W1=((w1*B1/pi)^(1/2))/((1-h1^2)^(1/4)); T2=[1,0;-2/R2,1]*[1,L2;0,1]*[1,0;-d,1]*[1,L1;0,1]*[1,0;-2/R1,1]*[1,L1 ;0,1]*[1,0;-d,1]*[1,L2;0,1]; A2=T2(1,1);B2=abs(T2(1,2));C2=T2(2,1);D2=T2(2,2);h2=(A2+D2)/2; W2=((w1*B2/pi)^(1/2))/((1-h2^2)^(1/4)); figure (1) ezplot(W1,[0,0.012]); title('透镜上的光斑半径');xlabel('透镜光焦度D（/mm）');ylabel('光束半径') figure (2); ezplot(W2,[0,0.012]) title('输出镜上的光斑半径');xlabel('透镜光焦度D（/mm）');ylabel('光束半径') figure (3);
Matalab课后作业 学院：电气信息工程及其自动化 班级： 学号： 姓名： 完成日期： 2012年12月23日
1、 matlab 软件主要功能是什么？电气工程及其自动化专业本科生主要用到哪 些工具箱，各有什么功能？ 答：（1）主要功能：工业研究与开发； 数学教学，特别是线性代数；数值分析和科学计算方面的教学与研究；电子学、控制理论和物理学等工程和科学学科方面的教学与研究； 经济学、化学和生物学等计算问题的所有其他领域中的教学与研究；符号计算功能；优化工具；数据分析和可视化功能；“活”笔记本功能；工具箱；非线性动态系统建模和仿真功能。 （2）常用工具箱： （a ） MATLAB 主工具箱：扩充matlab 的数值计算、符号运算功能、图形建模仿真功能、文字处理功能以及与硬件实时交互功能。 （b ）符号数学工具箱：符号表达式、符号矩阵的创建；符号可变精度求解；因式分解、展开和简化；符号代数方程求解；符号微积分；符号微分方程。 （c ） SIMULINK 仿真工具箱： Simulink 是用于动态系统和嵌入式系统的多领域仿真和基于模型的设计工具。对各种时变系统，包括通讯、控制、信号处理、视频处理和图像处理系统，Simulink 提供了交互式图形化环境和可定制模块库来对其进行设计、仿真、执行和测试。 （d ）信号处理工具箱：数字和模拟滤波器设计、应用及仿真；谱分析和估计；FFT 、DCT 等 变换；参数化模型。 （e ）控制系统工具箱：连续系统设计和离散系统设计；状态空间和传递函数以及模型转换；时域响应（脉冲响应、阶跃响应、斜坡响应）；频域响应（Bode 图、Nyquist 图）；根轨迹、极点配置。 2、设y=23e t 4-sin(43t+3 ),要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。 程序如下： clc clear x=0:0.01:1.5; y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3); y1=diff(y); subplot(2,1,1) plot(x,y) subplot(2,1,2) plot(x(1:150),y1) 曲线如下图所示：
实验一 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