物流系统建模与仿真实验
实验报告
实验学期至学年第学期年级专业班级
学生姓名学号
一、实验名称
产品测试工艺仿真与分析实验
二、实验内容
（1）实验背景
某制造车间由5 组机器组成，第1，2，3，4，5 组机器分别有3，2，4，3，1 台相同的机器。这个车间需要加工三种原料，三种原料分别要求完成4、3 和5 道工序，而每道工序必须在指定的机器组上处理，按照事先规定好的工艺顺序进行。 （2）概念模型如下：
（2）模型实体设计需求元素列表
模型元素 系统元素 备注
Flowitem 原料 不同实体类型代表不同类型的原料，分别标为1、2、3
Processor 机器 进行不同的参数定义以表征不同机器组中的机器 Conveyor 传送带 Rack 货架
Operator 操作员 可以进行搬运或加工等操作的人 Dispatcher 调度器 给操作员进行任务分配的控制器 Transporter 叉车 进行搬运操作的小车 Queue 机器组暂存区
Source 原材料库 原材料的始发处 Sink
成品库
原料加工后的最终去处
三、实验结果分析（模型运行图、数据统计图表）
1.模型平面图
2.模型立体图
检测机器1
检测机器2 检测机器3 传送带
传送带 传送带 货架1
货架2
货架3
产品1
产品2
产品3
3. 检测台1数据分析结果与state饼图
4. 检测台2数据分析结果图与state饼图
5. 检测台数据
分析结果图与
state饼图
6. 暂存
区1数据分析结果图
7. 暂存区2数据分析结果图
四、实验过程中遇到的问题及解决方
1. 数据结果
从仿真结果中可以看出三个检测机的调整时间分别为
21.9%、21.2%、24.6%，工作时间分别为58%、67.2%、69%，
空闲时间分别为20%、12.6%、8.7%。检测机几乎有1/5的
时间时再调整状态，检测机1有1/5的时间时在空闲的状
态。产品在三个检测机的平均逗留时间分别为36.49S、43.22S、38.67S。
暂存区1平均有9个产品，平均逗留时间为134S；暂存区2平均有20个产品，平均逗留时间为314S。
2. 模型瓶颈分析
从上面的分析得出三个检测机的工作效率分别为58%、67.2%、69%，其工作效率不是很高，而且暂存区2中产品数量多，平均逗留时间也长。
检测机工作效率不高，会造成利润的降低；暂存区产品多需要更大面积的暂存区去存放产品，增加产品的存储费用，产品在暂存区逗留时间越长，不确定因素造成的牛鞭效应就越大。
所以检测机工作效率不高，暂存区2库存量大都是该系统的瓶颈问题。应该着力解决系统的瓶颈，提高检测机的工作效率，降低暂存区的库存，减少产品在暂存区的逗留时间。
3. 模型改进
将检测机的准备时间缩短，可以抬高检测机的工作效率；增加叉车的工作效率，降低暂存区2的库存，减少产品在暂存区的逗留时间。
五、实验收获与体会
在本次实验中通过建立单存放区域、单处理工作台的简单模型，了解5个基本建模步骤。初步熟悉了Flexsim仿真软件如何应用。了解了每个实体的含义，实体属性的设置方法，实体之间的连接方法，如何通过观察仿
真结果找到系统的瓶颈等。
一、实验名称
多产品多阶段制造系统仿真与分析实验
二、实验内容
（1）实验背景
某制造车间由5 组机器组成，第1，2，3，4，5 组机器分别有3，2，4，3，1 台相同的
机器。这个车间需要加工三种原料，三种原料分别要求完成4、3 和5 道工序，而每道工序必须
在指定的机器组上处理，按照事先规定好的工艺顺序进行。 （2）概念模型参考如下：
、（3
）模型实体设计需求元素列表
模型元素 系统元素 备注
Flowitem 原料 不同实体类型代表不同类型的原料，分别标为1、2、3 Processor 机器 进行不同的参数定义以表征不同机器组中的机器 Queue 机器组暂存区
Source 原材料库 原材料的始发处，每天连续8小时提供原料，不同的原料有不同的间隔到达时间，因此用不同的Source 生成不同类型的原料 Sink
成品库
原料加工后的最终去处
三、实验结果分析（模型运行图、数据统计图表）
1. 仿真模型平面图
机器1
机器1
机器1
暂存区
机器2
机器2
暂存区
机器3
机器3
机器3
机器3
暂存区
机器4
机器4
机器4
暂存区
机器5
暂存区
1
2
3
2. 仿真模型立体图
暂存区1数据分析结果图暂存区2数据分析结果图
暂存区3数据分析结果图暂存区4数据分析图
暂存区5数
据分析结果
1. 仿真结果分析
通过仿真结果数据得出5个暂存区的等待加工的产品的平均数量分别为644、3、0.06、12.37、87.61，平均等待时间分别为9239、52、1.56、293、1241。从以上数据可以看出暂存区1产品数量过大，而且等待时间也过长；暂存区5的产品数量偏大，等待时间偏长。
所以此系统中暂存区1和暂存区5是该系统的瓶颈。
2.模型改进
解决瓶颈问题我们首先要查看与发生瓶颈处处相连接环节。与暂存区1相连接的是机器2、3、4分别对应着产品1、2、3的加工。要想解决暂存区1的瓶颈问题，必须要加快机器2、3、4加工产品1、2、3的加工效率。
与暂存区5相连接的是机器是机器1，对应的产品是3，要想解决暂存区5的瓶颈问题，就要加快机器1加工产品3的加工效率。
五、实验收获与体会
通过此次实验了解了多产品多阶段产制造系统运行过程，要想保证系统中不存在瓶颈必须协调好各种产品的加工时间与加工顺序。
一、实验名称
混合流水线系统仿真与分析实验
二、实验内容
1. 实验背景
多对象流水线生产有两种基本形式。一种是可变流水线，其特点是：在计划期内，按照一定的间隔期，成批轮番生产多种产品；在间隔期内，只生产一种产品，在完成规定的批量后，转生产另一种产品。另一种是混合流水线，其特点是：在同一时间内，流水线上混合生产多种产品。按固定的混合产品组组织生产，即将不同的产品按固定的比例和生产顺序编成产品组。一个组一个组地在流水线上进行生产。
2. 概念模型
3. 实验说明
（1）一个工厂有5个不同车间（普通车间，钻床车间，铣床车间，磨床车间，检测车间），加工3种类型产品。每种产品都要按工艺顺序在5个不同的车间完成5道工序。
（2）假定在保持车间逐日连续工作条件下，仿真在多对象标准化生产采用不同投产顺序
来生产给定数量的3种产品。通过改变投产顺序使产量、品种、工时和负荷趋于均衡，来减少时间损失。
（3）如果一项作业在特定时间到达车间，发现该组机器全都忙着，该作业就在该组机器处排（4）入一个FIFO规则的队列的暂存区，如果有前一天没有完成的任务，第二天继续加工。
表3.1: 车间配备（单位：台）
普通车间钻床车间铣床车间磨床车间检测车间
机器数量 3 3 2 3 1
表3.2: 加工时间（单位：min）
普通机床钻床铣床磨床检测产品1 5 5 4 4 6
产品2 4 4 3 4 3
产品3 4 5 3 4 1
表3.3：产品数量
总数（个）每批量（个）时间间隔（min）产品1 1000 10 30
产品2 500 5 30
产品3 200 2 30
三、实验结果分析（模型运行图、数据统计图表）
仿真模型平面图
仿真模型立体图
投产顺序为1、2、3时的仿真报告
投产顺序为1、3、2时的仿真报告
仿真顺序为2、1、3时的仿真报告
投产顺序为2、3、1时的仿真报告
投产顺序为3、1、2时的仿真报告
投产顺序为3、2、1时的仿真报告
四、实验过程中遇到的问题及解决方
1.模型仿真结果分析
从上面的仿真报告中看出，不同的投产顺序，生产该批产品所需时间不同。其中投产顺序为3、1、2时为最佳方案，因为此种方案的生产时间最短为6833S。
从仿真报告中得出，虽然投产顺序为3、1、2时是最佳方案，但是此系统产品在暂存区3和暂存区5中逗留时间过长，分别为142.91S和1909.21S。产品在暂存区中逗留时间过长会造成产品堆积，增加库存，库存费用也会随着增加。同时产品逗留时间过长会增加不确定因素造成牛鞭效应。所以暂存区3和暂存区5是该系统的瓶颈。
2.模型改进
要想减少暂存区3和暂存区5库存，我们需要增加暂存区后面环节的工作效率，即增加铣床
和检测机的工作效率，减少每件产品的加工时间和检测时间。或者是增加铣床和检测机的数量用以满足生产以满足生产需求。
五、实验收获与体会
通过此次仿真实验熟悉了多产品多阶段制造系统的两种生产方式，可变流水线生产和混合流水线生产。了解到在可变流水线生产方式中不同投产顺序会有不同的生产效率，在工作使用此种方式生产产品时我们要充分考虑产品的投产顺序，以确保为最佳生产方式，使获得最大利润。
一、实验名称
垃圾回收场仿真与分析实验
二、实验内容
1. 实验背景
近几十年来，由于人类的滥砍、滥伐，无情的破坏我们的大自然，地球上能用的资产和能源逐渐地减少，环保团体发现如果我们不再注重保护环境，终有一天我们会失去地球这个美好的家园。所以近年来，环保团体大力的提倡垃圾回收，位于某地的一家垃圾回收站，把回收来的资源分成铁铝罐、保持瓶和塑胶三大类后存储起来。下面这个模型就是对该资源回收站的仿真。
2. 概念模型
3.模型实体设计需求元素列表
模型元素系统元素备注
Flowitem 回收到的垃圾
Source 产生垃圾
Queue 垃圾暂存区
Separator 垃圾分拣装置将回收到的垃圾分拣后输出
Conveyor 传送设备
FlowNode 传送设备
Reservoir 垃圾存储容器
Rack 垃圾存储容器
Crane 垃圾搬运设备
Transporter 垃圾搬运设备
三、实验结果分析（模型运行图、数据统计图表）
仿真模型平面图
仿真模型立体图
分拣器数据分析结果图分拣器饼图
分拣器数据报告
四、实验过程中遇到的问题及解决方
1. 仿真结果分析
从分拣器饼图得出该系统的分拣器有50.3%的时间处于空闲状态，5.6%的时间处于阻塞状态，工作时间还不到一半；从仿真数据报告得出该系统中垃圾在暂存区2和暂存区13中的逗留时间较长，但是库存量不是很多只有8个。
产品在暂存区逗留时间长主要是由于垃圾搬运设备速度慢，工作效率低造成的。
2.模型改进
解决该系统的瓶颈问题，首先要加快搬运设备和传送带的速度，协调好传送带与其他搬运设别之间的速度。
五、实验收获与体会
通过此次实验我们了解到分拣器、流节点、储液罐等仿真工具在仿真模型中的应用，同时也了解到了垃圾场垃圾处理过程，在垃圾处理过程中要提高垃圾的搬运速度，提高垃圾处理效率。
在日常生活中我们也要注意保护好环境，不乱丢垃圾，提高垃圾分类的意识，方便工作人员的垃圾回收。
一、实验名称
配送中心仿真与分析实验
二、实验内容
1. 实验背景
配送中心时从事货物配送并组织对用户的送货，以实现销售和供应服务的现代流通设施。它不同于传统的仓储设施，在现代商业社会中，配送中心已经成为连锁企业的商流中心、物流中心、信息流中心，是连锁经营得以正常运转的关键设施。该实验是一个典型的配送中心建模过程，该配送中心从三个供应商进货，向三个生长商发货。
2. 概念模型
3. 模型实体设计需求元素列表
模型元素系统元素备注
Flowitem 产品
Source 发生产品3个Source发生产品的速度相同且快于供货商供应速度
模型前面的三个
Processor（按模型流程）供货商
3个Processor加工速率不同，按照模型的系统数
据进行设定
Rack 配送中心3个Rack分别对应3个供货商
Queue 生产商仓库3个Queue订货条件不同，根据模型的系统数据进行设定
模型后面的三个Processor（按模型流程）生产商3个Processor加工速率不同，按照模型的系统数据进行设定
Sink 产品收集装置产品的最终去处
三、实验结果分析（模型运行图、数据统计图表）
仿真模型平面图
仿真模型立体图
配送中心1数据配送中心2数据
配送中心3数据
四、实验过程中遇到的问题及解决方
1.计算
配送中心的进货总成本
进货成本=产品进货数量*产品进货单价
从仿真结果中得到：产品进货数量为1149 + 2189 + 1933 = 5271
进货单价为3元每件
进货成本=5279*3=15837元
供货总收入
供货总收入=供货数量*供货价格
供货数量=1142+2185+1930=5257
供货价格5元/件
供货总收入=5257*5=26285元
存货总成本
存货总成本=平均库存数量*平均库存时间（小时）/100
平均存货量分别为5.54、2.39、2.76
平均存货时间为42.32小时、9.58小时、12.52小时
每件产品在配送中心存货100小时费用1元。
库存成本=5.54*42.32/100+2.39*9.58/100+2.76*12.52/100=2.92元
利润
利润=供货总收入-进货总成本-存货成本=10445.08元
2. 模型仿真结果分析
从仿真结果中得出产品1在配送中心的平均逗留时间为42.32小时，逗留时间明显比其他两个产品要长，这样会造成库存成本的增加，而且产品长时间的逗留在配送中心，对市场的应对能力会变差，影响企业的应变性，不确定性造成的牛鞭效应加剧。
应该减少产品在配送中心的逗留时间，增产品的流转速度。多批次少批量的进货策略可以减少库存与库存时间。
五、实验收获与体会
通过此次实验我们学到如何按照订单来安排生产、存储与发货。研究该配送中心的即时库存成本和利润。配送中心应该尽量减少库存与库存时间，逐渐逼近零库存的趋势，这样可以大大降低成本。也应该规划好合适进行发货，何时进行进货。
了解到配送中心的作用，配送中心是连接生产商与供应商的桥梁，它可以通过联合运输降低运输成本，降低供应链成本，配送中心是企业发展的需要。
一、实验名称
流体建模与仿真分析实验
二、实验内容
1. 实验背景
一个操作员把两种不同类型的实体运送到Item To Fluid，通过Item To Fluid转化成两种流体，并通过管道运输至储液罐。然后两个储液罐中的流体物质通过管道运至Mixer混合，生成一种新产物，新产物经过流体处理器处理，并通过Fluid To Item转化为实体（flow items），最后经过输送机运至吸收器离开模型。
2. 概念模型
3.模型实体设计需求元素列表
模型元素系统元素备注
Flowitem 回收到的垃圾
Source 产生垃圾
Queue 垃圾暂存区
Separator 垃圾分拣装置将回收到的垃圾分拣后输出
Conveyor 传送设备
FlowNode 传送设备
Reservoir 垃圾存储容器
Rack 垃圾存储容器
Crane 垃圾搬运设备
Transporter 垃圾搬运设备
三、实验结果分析（模型运行图、数据统计图表）
《计算方法》实验报告 姓名： 班级： 学号： 实验日期: 2011年10月26日
一、实验题目： 数值积分 二、实验目的： 1．熟悉matlab 编写及运行数值计算程序的方法。 2．进一步理解数值积分的基础理论。 3．进一步掌握应用不同的数值积分方法求解给定的积分并给出数据结果及误差分析。 三、实验内容： 1．分别用复合梯形求积公式及复合辛普森求积公式计算积分xdx x ln 10 ? , 要求计算精度达到410-,给出计算结果并比较两种方法的计算节点数. 2．用龙贝格求积方法计算积分dx x x ?+3 021,使误差不超过510-. 3．用3=n 的高斯-勒让德公式计算积分?3 1 sin x e x ,给出计算结果. 4．用辛普森公式(取2==M N ) 计算二重积分.5 .00 5 .00 dydx e x y ? ? - 四、实验结果： 1.（1）复合梯形法: 将区间[a,b]划分为n 等份，分点n k n a b h kh a x k ,2,1,0,,=-=+=在每个区间[1,+k k x x ](k=0,1,2,···n-1)上采用梯形公式，则得 )()]()([2)()(1 11 1 f R x f x f h dx x f dx x f I n n k k k b a n k x x k k ++===∑?∑? -=+-=+ 故)]()(2)([21 1 b f x f a f h T n k k n ++=∑-=称为复合梯形公式 计算步长和划分的区间 Eps=1E-4 h1=sqrt(Eps/abs(-(1-0)/12*1/(2+1))) h1 =0.0600 N1=ceil(1/h1) N1 =17 用复合梯形需要计算17个结点。 复合梯形： function T=trap(f,a,b,n) h=(b-a)/n;
信号与系统仿真实验报告1．实验目的 了解MATLAB的基本使用方法和编程技术，以及Simulink平台的建模与动态仿真方法，进一步加深对课程内容的理解。 2．实验项目 信号的分解与合成，观察Gibbs现象。 信号与系统的时域分析，即卷积分、卷积和的运算与仿真。 信号的频谱分析，观察信号的频谱波形。 系统函数的形式转换。 用Simulink平台对系统进行建模和动态仿真。 3．实验内容及结果 3.1以周期为T，脉冲宽度为2T1的周期性矩形脉冲为例研究Gibbs现象。 已知周期方波信号的相关参数为：x(t)=∑ak*exp(jkω),ω=2*π/T,a0=2*T1/T,ak=sin(kωT1)/kπ。画出x(t)的波形图（分别取m=1,3,7,19,79,T=4T1）,观察Gibbs现象。 m=1; T1=4; T=4*T1;k=-m:m; w0=2*pi/T; a0=2*T1/T; ak=sin(k*w0*T1)./(k*pi); ak(m+1)=a0; t=0:0.1:40; x=ak*exp(j*k'*w0*t); plot(t,real(x)); 3.2求卷积并画图 (1)已知：x1(t)=u(t-1)-u(t-2), x2(t)=u(t-2)-u(t-3)求：y(t)=x1(t)*x2(t)并画出其波形。 t1=1:0.01:2; f1=ones(size(t1)); f1(1)=0; f1(101)=0; t2=2:0.01:3; f2=ones(size(t2)); f2(1)=0; f2(101)=0; c=conv(f1,f2)/100;
t3=3:0.01:5; subplot(311); plot(t1,f1);axis([0 6 0 2]); subplot(312); plot(t2,f2);axis([0 6 0 2]); subplot(313); plot(t3,c);axis([0 6 0 2]); (2）已知某离散系统的输入和冲击响应分别为：x[n]=[1,4,3,5,1,2,3,5], h[n]=[4,2,4,0,4,2].求系 统的零状态响应，并绘制系统的响应图。 x=[1 4 3 5 1 2 3 5]; nx=-4:3; h=[4 2 4 0 4 2]; nh=-3:2; y=conv(x,h); ny1=nx(1)+nh(1); ny2=nx(length(nx))+nh(length(nh)); ny=[ny1:ny2]; subplot(311); stem(nx,x); axis([-5 4 0 6]); ylabel('输入') subplot(312); stem(nh,h); axis([-4 3 0 5]); ylabel('冲击效应') subplot(313); stem(ny,y); axis([-9 7 0 70]); ylabel('输出'); xlabel('n'); 3.3 求频谱并画图 (1) 门函数脉冲信号x1(t)=u(t+0.5)-u(t-0.5) N=128;T=1; t=linspace(-T,T,N); x=(t>=-0.5)-(t>=0.5); dt=t(2)-t(1); f=1/dt; X=fft(x); F=X(1:N/2+1); f=f*(0:N/2)/N; plot(f,F)
模拟电路仿真实验报告 张斌杰生物医学工程141班 MUltiSim软件使用 一、实验目的 1、掌握MUltiSim软件的基本操作和分析方法。 二、实验内容 1、场效应管放大电路设计与仿真 2、仪器放大器设计与仿真 3、逻辑电平信号检测电路设计与仿真 4、三极管Beta值分选电路设计与仿真 5、宽带放大电路设计与仿真 三、MUItiSim软件介绍 MUItiSim是美国国家仪器(NI)有限公司推出的以WindOWS为基础的仿真工具，适用于板级的模拟/数字电路板的设计工作。它包含了电路原理图的图形输入、电路硬件描述语言输入方式，具有丰富的仿真分析能力。工程师们可以使用MUItiSinl交互式地搭建电路原理图，并对电路进行仿真。MUltiSiIn提炼了SPICE 仿真的复杂内容,这样工程师无需懂得深入的SPlCE技术就可以很快地进行捕获、仿真和分析新的设计，这也使其更适合电子学教育。通过MUItiSiIn和，PCB设计工程师和电子学教育工作者可以完成从理论到原理图捕获与仿真再到和测试这样一个完整的综合设计流程。 实验名称:
仪器放大器设计与仿真 二、实验目的 1、 掌握仪器放大器的设计方法 2、 理解仪器放大器对共模信号的抑制能力 3、 熟悉仪器放大器的调试功能 4、 掌握虚拟仪器库中关于测试模拟电路仪器的使用方法，如示波器，毫伏 表信 号发生器等虚拟仪器的使用 三、设计实验电路图： 四、测量实验结果: 出为差模放大为399mvo 五、实验心得： 应用MUIti S im 首先要准备好器件的PSPiCe 模型，这是最重要的，没有这个 东西免谈，当然SPiCe 高手除外。下面就可以利用MUItiSinl 的元件向导功 能制作 差模分别输入信号InW 第二条线与第三条线: 共模输入2mv 的的电压，输出为2mv 的电压。 第一条线输
本科实验报告 课程名称：计算机数值方法 实验项目：方程求根、线性方程组的直接解 法、线性方程组的迭代解法、代数插值和最 小二乘拟合多项式 实验地点：行勉楼 专业班级： ******** 学号： ********* 学生姓名： ******** 指导教师：李誌，崔冬华 2016年 4 月 8 日
y = x*x*x + 4 * x*x - 10; return y; } float Calculate(float a,float b) { c = (a + b) / 2; n++; if (GetY(c) == 0 || ((b - a) / 2) < 0.000005) { cout << c <<"为方程的解"<< endl; return 0; } if (GetY(a)*GetY(c) < 0) { return Calculate(a,c); } if (GetY(c)*GetY(b)< 0) { return Calculate(c,b); } } }; int main() { cout << "方程组为：f(x)=x^3+4x^2-10=0" << endl; float a, b; Text text; text.Getab(); a = text.a; b = text.b; text.Calculate(a, b); return 0; } 2.割线法： // 方程求根（割线法）.cpp : 定义控制台应用程序的入口点。// #include "stdafx.h" #include"iostream"
心得体会 使用不同的方法，可以不同程度的求得方程的解，通过二分法计算的程序实现更加了解二分法的特点，二分法过程简单，程序容易实现，但该方法收敛比较慢一般用于求根的初始近似值，不同的方法速度不同。面对一个复杂的问题，要学会简化处理步骤，分步骤一点一点的循序处理，只有这样，才能高效的解决一个复杂问题。
测试装置动态特性仿真实验 班级：7391 学号：2009301828 姓名：张志鹏 一、实验目的 1、加深对一阶测量装置和二阶测量装置的幅频特性与相频特性的理解； 2、加深理解时间常数变化对一阶系统动态特性影响； 3、加深理解频率比和阻尼比变化对二阶系统动态特性影响； 4、使学生了解允许的测量误差与最优阻尼比的关系。 二、实验原理 1、 一阶测量装置动态特性 一阶测量装置是它的输入和输出关系可用一阶微分方程描述。一阶测量装置的频率响应函数为： 式中：S S 为测量装置的静态灵敏度；τ为测量装置的时间常数。 一阶测量装置的幅频特性和相频特性分别为： 可知，在规定S S =1的条件下，A （ω）就是测量装置的动态灵敏度。 当给定一个一阶测量装置，若时间常数τ确定，如果规定一个允许的幅值误差ε，则允许测量的信号最高频率ωH 也相应地确定。 为了恰当的选择一阶测量装置，必须首先对被测信号的幅值变化范围和频率成分有个初步了解。有根据地选择测量装置的时间常数τ，以保证A （ω）≥1-ε 能够满足。 2、二阶测量装置动态特性 二阶测量装置的幅频特性与相频特性如下： 幅频特性202220)/(4))/(1(/1)(ωωξωωω--=A 相频特性2200))/(1/()/(2()(ωωωωξφ--=arctg w Α（ω）是ξ和ω/0ω的函数，即具有不同的阻尼比ξ的测试装置当输入信??????ωτ+ωτ-ωτ+=ωτ+=ω22s s )(1j ) (11S j 11S )j (H ()()2 11 A ωτ+=ω()ωτ -=ωφarctan
号频率相同时，应具有不同的幅值响应，反之，当不同的频率的简谐信号送入同一测试装置时它们的幅值响应也不相同，同理具有不同的阻尼比ξ的测试装置当输入信号频率相同时，应有不同的相位差。 (1).当ω=0时，Α（ω）=1;(2).当ω→∞，A （ω）=0;(3).当ξ≥0.707时随着输入信号频率的加大，Α（ω）单调的下降, ξ＜0.707时Α（ω）的特性曲线上出现峰值点;(4)如果ξ=0，))/(1/(1))/(1(/1)(202 20ωωωωω-=-=A ，显然，其峰值点出现在ω=0ω处。其值为“∞”，当ξ从0向0.707变化过程中随着的加大其峰值点逐渐左移，并不断减小。 对以上二阶环节的幅频特性的结论论证如下： (1).当ω=0时A(ω)=1 (2).当ω→∞时,A(ω)=0 (3).要想得到A(ω)的峰值就要使202220)/(4))/(1(/1)(A ωωξ-ωω-=ω 中的202220)/(4))/(1(ωωξωω--取最小值。 令:t=20)/(ωω t t t f 224)1()(ξ+-= 对其求导可得t=1-22ξ时,f(t)取最小值.由于t=20)/(ωω≥0,所以1-22ξ≥0, 2ξ必须小于1/2时,f(t)才有最小值,即ξ>2/2时，A(ω)不出现峰值点;当ξ＜2/2时4244)(ξξ-=t f ,f(t)对ξ求导得)21(82ξξ-，可以看出f(t): ξ属于[0， 2/2]时单调递增，于是得A(ω)的峰值点A 为4244/1)(/1ξξ-=t f ; 在ξ属于[0，2/2]递减。 (4).当ξ=0时 A=∞，t=20)/(ωω，ω/0ω=1,即ξ=0时A(ω)的峰值为∞，且必出现在ω/0ω=1时，当ξ=2/2时,t=0→ω=0，A(ω)=1. 还可以看出，在ξ属于[0，2/2]增大时t=1-22ξ就减小，即f(t)的峰值左平移。 （二）阻尼比的优化 在测量系统中，无论是一阶还是二阶系统的幅频特性都不能满足将信号中的所有频率都成比例的放大。于是希望测量装置的幅频特性在一段尽可能宽的范围内最接近于1。根据给定的测量误差，来选择最优的阻尼比。
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)体系设计： 程序是一个简单的计算器，能正确输入数据，能实现加、减、乘、除等算术运算，运算结果能正确显示，可以清楚数据等。 (2)设计思路： 1)先在Visual C++ 6.0中建立一个MFC工程文件，名为 calculator. 2)在对话框中添加适当的编辑框、按钮、静态文件、复选框和单 选框 3)设计按钮，并修改其相应的ID与Caption. 4)选择和设置各控件的单击鼠标事件。 5)为编辑框添加double类型的关联变量m_edit1. 6)在calculatorDlg.h中添加math.h头文件，然后添加public成 员。 7)打开calculatorDlg.cpp文件，在构造函数中，进行成员初始 化和完善各控件的响应函数代码。 (3)程序清单：
●添加的public成员： double tempvalue; //存储中间变量 double result; //存储显示结果的值 int sort; //判断后面是何种运算：1.加法2.减法3. 乘法 4.除法 int append; //判断后面是否添加数字 ●成员初始化： CCalculatorDlg::CCalculatorDlg(CWnd* pParent /*=NULL*/) : CDialog(CCalculatorDlg::IDD, pParent) { //{{AFX_DATA_INIT(CCalculatorDlg) m_edit1 = 0.0; //}}AFX_DATA_INIT // Note that LoadIcon does not require a subsequent DestroyIcon in Win32 m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME); tempvalue=0; result=0; sort=0; append=0; }
大学物理仿真实验报告一一塞曼效应 一、实验简介 塞曼效应是物理学史上一个著名的实验。荷兰物理学家塞曼（Zeeman）在1896年发现把产生光谱的光源置于足够强的磁场中，磁场作用于发光体，使光谱发生变化，一条谱线即会分裂成几条偏振化的谱线，这种现象称为塞曼效应。 塞曼效应是法拉第磁致旋光效应之后发现的又一个磁光效应。这个现象的发现是对光的 电磁理论的有力支持，证实了原子具有磁矩和空间取向量子化，使人们对物质光谱、原子、分子有更多了解。 塞曼效应另一引人注目的发现是由谱线的变化来确定离子的荷质比的大小、符号。根据 洛仑兹（H.A?Lorentz）的电子论，测得光谱的波长，谱线的增宽及外加磁场强度，即可称得离子的荷质比。由塞曼效应和洛仑兹的电子论计算得到的这个结果极为重要，因为它发表在J、 J汤姆逊（J、J ThomSOn）宣布电子发现之前几个月，J、J汤姆逊正是借助于塞曼效应由洛仑 兹的理论算得的荷质比，与他自己所测得的阴极射线的荷质比进行比较具有相同的数量级，从而得到确实的证据，证明电子的存在。 塞曼效应被誉为继X射线之后物理学最重要的发现之一。 1902年，塞曼与洛仑兹因这一发现共同获得了诺贝尔物理学奖（以表彰他们研究磁场对光的效应所作的特殊贡献）。至今，塞曼效应依然是研究原子内部能级结构的重要方法。 本实验通过观察并拍摄Hg（546.1 nm）谱线在磁场中的分裂情况，研究塞曼分裂谱的特征，学习应用塞曼效应测量电子的荷质比和研究原子能级结构的方法。 二、实验目的 1?学习观察塞曼效应的方法观察汞灯发出谱线的塞曼分裂； 2?观察分裂谱线的偏振情况以及裂距与磁场强度的关系； 3?利用塞曼分裂的裂距，计算电子的荷质比 e m e数值。 三、实验原理 1、谱线在磁场中的能级分裂 设原子在无外磁场时的某个能级的能量为E0,相应的总角动量量子数、轨道量子数、 自旋量子数分别为J、L、S。当原子处于磁感应强度为B的外磁场中时，这一原子能级将 分裂为2J 1层。各层能量为 E = E o MgJ B B（1） 其中M为磁量子数，它的取值为J , J -1 ,…，-J共2J 1个；g为朗德因子；J B为 hc 玻尔磁矩（A B= ）；B为磁感应强度。 4兀m 对于L-S耦合
中南大学系统仿真实验报告 指导老师胡杨 实验者 学号 专业班级 实验日期 2014.6.4 学院信息科学与工程学院
目录 实验一MATLAB中矩阵与多项式的基本运算 (3) 实验二MATLAB绘图命令 (7) 实验三MATLAB程序设计 (9) 实验四MATLAB的符号计算与SIMULINK的使用 (13) 实验五MATLAB在控制系统分析中的应用 (17) 实验六连续系统数字仿真的基本算法 (30)
实验一MATLAB中矩阵与多项式的基本运算 一、实验任务 1．了解MATLAB命令窗口和程序文件的调用。 2．熟悉如下MATLAB的基本运算： ①矩阵的产生、数据的输入、相关元素的显示； ②矩阵的加法、乘法、左除、右除； ③特殊矩阵：单位矩阵、“1”矩阵、“0”矩阵、对角阵、随机矩阵的产生和运算； ④多项式的运算：多项式求根、多项式之间的乘除。 二、基本命令训练 1．eye(m) m=3; eye(m) ans = 1 0 0 0 1 0 0 0 1 2．ones(n)、ones(m,n) n=1;m=2; ones(n) ones(m,n) ans = 1 ans = 1 1
3．zeros(m,n) m=1,n=2; zeros(m,n) m = 1 ans = 0 0 4．rand(m,n) m=1;n=2; rand(m,n) ans = 0.8147 0.9058 5．diag(v) v=[1 2 3]; diag(v) ans = 1 0 0 0 2 0 0 0 3 6．A\B 、A/B、inv(A)*B 、B*inv(A) A=[1 2;3 4];B=[5 6;7 8]; a=A\B b=A/B c=inv(A)*B d=B*inv(A) a = -3 -4 4 5 b = 3.0000 -2.0000 2.0000 -1.0000
西安交通大学实验报告 课程：数据结构实验实验名称：利用单摆测量重力加速度 系别：实验日期： 专业班级：实验报告日期： 姓名：学号： 第 1页 / 共3页 一、实验简介 单摆实验是个经典实验，许多著名的物理学家都对单摆实验进行过细致的研究。本实验的目的是学习进行简单设计性实验的基本方法，根据已知条件和测量精度的要求，学会应用误差均分原则选用适当的仪器和测量方法，学习累积放大法的原理和应用，分析基本误差的来源及进行修正的方法。 二、实验原理 单摆的结构参考图1单摆仪,一级近似的周期公式为 由此通过测量周期摆长求重力加速度。 三、实验内容 1、设计要求: (1) 根据误差均分原理,自行设计试验方案,合理选择测量仪器和方法. (2) 写出详细的推导过程,试验步骤. (3) 用自制的单摆装置测量重力加速度g,测量精度要求△g/g < 1%. 2、可提供的器材及参数: 游标卡尺、米尺、千分尺、电子秒表、支架、细线(尼龙线)、钢球、摆幅测量标尺(提供硬白纸板自制)、天平(公用).
假设摆长l≈70.00cm;摆球直径D≈2.00cm;摆动周期T≈1.700s; 米尺精度△米≈ 0.05cm;卡尺精度△卡≈0.002cm;千分尺精度△千≈0.001cm;秒表精度△秒≈0.01s;根据统计分析,实验人员开或停秒表反应时间为0.1s左右,所以实验人员开,停秒表总的反应时间近似为△人≈0.2s. 3、对重力加速度g的测量结果进行误差分析和数据处理,检验实验结果是否达到设计要求. 4、自拟实验步骤研究单摆周期与摆长,摆角,悬线的质量和弹性系数,空气阻力等因素的关系,试分析各项误差的大小. 5、自拟试验步骤用单摆实验验证机械能守恒定律. 四、实验仪器 单摆仪，摆幅测量标尺，钢球，游标卡尺 五、实验操作 1. 用米尺测量摆线长度； 2. 用游标卡尺测量小球直径； 3. 把摆线偏移中心不超过5度，释放单摆，开始计时，单摆摆过50个周期后停止计时，记录所用时间； 六、实验结果
计算方法实验报告格式 小组名称: 组长姓名(班号)： 小组成员姓名(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一个完整的实验，应包括数据准备、理论基础、实验内容及方法，最终对实验结果进行分析，以达到对理论知识的感性认识，进一步加深对相关算法的理解，数值实验以实验报告形式完成，实验报告格式如下： 一、实验名称 实验者可根据报告形式需要适当写出． 二、实验目的及要求 首先要求做实验者明确，为什么要做某个实验，实验目的是什么，做完该实验应达到什么结果，在实验过程中的注意事项，实验方法对结果的影响也可以以实验目的的形式列出． 三、算法描述（实验原理与基础理论） 数值实验本身就是为了加深对基础理论及方法的理解而设置的，所以要求将实验涉及到的理论基础，算法原理详尽列出． 四、实验内容 实验内容主要包括实验的实施方案、步骤、实验数据准备、实验的算法以及可能用到的仪器设备． 五、程序流程图 画出程序实现过程的流程图，以便更好的对程序执行的过程有清楚的认识，在程序调试过程中更容易发现问题． 六、实验结果 实验结果应包括实验的原始数据、中间结果及实验的最终结果，复杂的结果可以用表格
形式列出，较为简单的结果可以与实验结果分析合并出现． 七、实验结果分析 实验结果分析包括对对算法的理解与分析、改进与建议． 数值实验报告范例 为了更好地做好数值实验并写出规范的数值实验报告，下面给出一简单范例供读者参考． 数值实验报告 小组名称: 小组成员(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一、实验名称 误差传播与算法稳定性． 二、实验目的 1．理解数值计算稳定性的概念． 2．了解数值计算方法的必要性． 3．体会数值计算的收敛性与收敛速度． 三、实验内容 计算dx x x I n n ? += 1 10 ，1,2,,10n = ． 四、算法描述 由 dx x x I n n ? += 1 10 ，知 dx x x I n n ?+=--101110，则
《物流仿真实验》 实验报告书 实验报告题目: 物流仿真实验学院名称: 管理学院 专业: 物流管理 班级: 物流1303 姓名: 孟颖颖 学号: 2 成绩: 2016年7月 实验报告 一、实验名称 物流仿真实验 二、实验要求 ⑴根据模型描述与模型数据对配送中心进行建模;
⑵分析仿真实验结果,进行利润分析,找出利润最大化的策略。 三、实验目的 1、掌握仿真软件Flexsim的操作与应用,熟悉通过软件进行物流仿真建模。 2、记录Flexsim软件仿真模拟的过程,得出仿真的结果。 3、总结Flexsim仿真软件学习过程中的感受与收获。 三、实验设备 (1)硬件及其网络环境 服务器一台:PII400/10、3G/128M以上配置、客户机100台、局域网或广域网。 (2)软件及其运行环境 Flexsim,Windows 2000 Server、SQL Server 7、0以上版本、IIS 5、0、SQL Server 数据库自动配置、IIS 虚拟目录自动配置 四、实验步骤 1 概念模型 1个Sink到操作区,如图:
第二步:连接端口 根据配送流程,对模型进行适宜的连接,所有端口连接均用A连接,如图: 第三步:Source的参数设置 为使Source产生实体不影响后面Processor的生产,尽可能的将时间间隔设置尽可能的小,并对三个Source做出同样的设定。 打开Source参数设置窗口,将时间到达间隔设置为常数1,同时为对三个实体进行区别,进行设置产品颜色,点击触发器,打开离开触发的下拉菜单,点击设置临时实体类型,设置不同实体类型,颜色自然发生变化。并对另外两个Source 进行同样的设置,如图:
一、系统描述 1.1.系统背景 本系统将基于下面的卫星屏幕快照创建一个模型。当前道路网区域的两条道路均为双向，每个运动方向包含一条车道。Tapiolavagen路边有一个巴士站，Menninkaisentie路边有一个带五个停车位的小型停车场。 1.2.系统描述 （1）仿真十字路口以及三个方向的道路，巴士站，停车点；添加小汽车、公交车的三维动画，添加红绿灯以及道路网络描述符； （2）创建仿真模型的汽车流程图，三个方向产生小汽车，仿真十字路口交通运行情况。添加滑条对仿真系统中的红绿灯时间进行实时调节。添加分析函数，统计系统内汽车滞留时间，用直方图进行实时展示。 二、仿真目标 1、timeInSystem值：在流程图的结尾模块用函数统计每辆汽车从产生到丢弃的，在系统中留存的时间。 2、p_SN为十字路口SN方向道路的绿灯时间，p_EW为十字路口EW方向道路的绿灯时间。 3、Arrival rate：各方向道路出现车辆的速率（peer hour）。
三、系统仿真概念分析 此交通仿真系统为低抽象层级的物理层模型，采用离散事件建模方法进行建模，利用过程流图构建离散事件模型。 此十字路口交通仿真系统中，实体为小汽车和公交车，可以源源不断地产生；资源为道路网络、红绿灯时间、停车点停车位和巴士站，需要实施分配。系统中小汽车（car）与公共汽车（bus）均为智能体，可设置其产生频率参数，行驶速度，停车点停留时间等。 四、建立系统流程 4.1.绘制道路 使用Road Traffic Library中的Road模块在卫星云图上勾画出所有的道路，绘制交叉口，并在交叉口处确保道路连通。 4.2.建立智能体对象 使用Road Traffic Library中的Car type模快建立小汽车（car）以及公共汽车（bus）的智能体对象。 4.3.建立逻辑 使用Road Traffic Library中的Car source、Car Move To、Car Dispose、
XXXXX 实验报告 学院（部）XX学院 课程名称生产系统仿真实验 学生姓名 学号 专业 2012年9月10日
《生产系统仿真》实验报告 年月日 学院年级、专业、班实验时间9月10日成绩 课程名称生产系统仿真 实训项目 名称 系统仿真软件的基础应 用 指导 教师 一、实验目的 通过对Flesim软件进一步的学习，建立模型，运用Flesim软件仿真该系统，观察并分析运行结果，找出所建模型的问题并进行改进，再次运行循环往复，直到找出构建该系统更为合理的模型。 二、实验内容 １、建立生产模型。 该模型生产三种产品，产品到达速率服从均值为20、方差为2的正态分布；暂存器的最大容量为25个；检测器的检测时间服从均值为30的指数分布，预制时间为10s；传送带的传送速率为1m/s，带上可容纳的最大货件数为10个。 ２、运行生产模型。 ３、对运行结果进行分析，提出改进方案在运行，直到找到更为合理的模型。 三、实验报告主要内容 １、根据已有数据建立生产模型。 将生产系统中所需实体按组装流程进行有序的排列，并进行连接如图１所示
图1 ２、分别对发生器、暂存器、检验台和传送带进行参数设置。 （１）发生器的产品到达速率服从均值为20、方差为2的正态分布。如图2所示。 （２）暂存器的最大容量设置为25件。如图3所示。 （3）设置检验台的检测时间服从均值为30s的指数分布，预制时间为10s.如图4所示。 （4）传送带的传送速率为1m/s，最大容量为10件。如图5所示 图2 图3 图4 图5 3、对发生器及暂存器进一步设置。 （1）发生器在生成产品时设置三种不同类型的产品，通过颜色区分。如图6所示。 （2）暂存器在输出端口通过设置特定函数以使不同颜色的产品在不同的检验台检验。如图7所示。
哈尔滨理工大学实验报告 控制系统仿真 专业：自动化12-1 学号：1230130101 姓名：
一．分析系统性能 课程名称控制系统仿真实验名称分析系统性能时间8.29 地点3# 姓名蔡庆刚学号1230130101 班级自动化12-1 一．实验目的及内容： 1. 熟悉MATLAB软件的操作过程； 2. 熟悉闭环系统稳定性的判断方法； 3. 熟悉闭环系统阶跃响应性能指标的求取。 二．实验用设备仪器及材料： PC, Matlab 软件平台 三、实验步骤 1. 编写MATLAB程序代码； 2. 在MATLAT中输入程序代码，运行程序； 3.分析结果。 四．实验结果分析： 1.程序截图
得到阶跃响应曲线 得到响应指标截图如下
2.求取零极点程序截图 得到零极点分布图 3.分析系统稳定性 根据稳定的充分必要条件判别线性系统的稳定性最简单的方法是求出系统所有极点,并观察是否含有实部大于0的极点，如果有系统不稳定。有零极点分布图可知系统稳定。
二．单容过程的阶跃响应 一、实验目的 1. 熟悉MATLAB软件的操作过程 2. 了解自衡单容过程的阶跃响应过程 3. 得出自衡单容过程的单位阶跃响应曲线 二、实验内容 已知两个单容过程的模型分别为 1 () 0.5 G s s =和5 1 () 51 s G s e s - = + ，试在 Simulink中建立模型，并求单位阶跃响应曲线。 三、实验步骤 1. 在Simulink中建立模型，得出实验原理图。 2. 运行模型后，双击Scope，得到的单位阶跃响应曲线。 四、实验结果 1．建立系统Simulink仿真模型图，其仿真模型为
（1）机能学基础性虚拟实验教学软件包含四个相对独立的操作实验：家兔的基本实验虚拟操作、蟾蜍的基本实验虚拟操作、大鼠的基本实验虚拟操作、小鼠的基本实验虚拟操作。所有内容全部采用人机互动的虚拟仿真操作来完成，同时配合动画演示，相关仪器设备的使用和操作知识。 我们以大小鼠和蟾蜍的基本实验虚拟操作举例说明： 《大、小鼠基本操作综合实验》介绍了大、小鼠在实验中经常用到的几种基本操作，通过虚拟操作的演示和互动，把实验中的重点、难点表示出来，使学生通过该虚拟实验，熟悉大小鼠实验的各项基本操作，掌握实验的重点。 虚拟实验操作流程及技术点描述： 大小鼠的捉持主要采用动画演示的形式，生动体现了捉持的要点。 大小鼠的固定，又分为徒手固定，固定板固定，头部固定以及固定器固定。学生可以自行选择固定方式，对大小鼠进行固定。 大小鼠的分组与编号；分组演示了如何使用Excel软件取得随机数字后分组。编号着重介绍了背毛单色标记法。 常用给药方法的虚拟操作：灌胃法，皮下注射法，皮内注射法，肌肉注射法，腹腔注射法，静脉注射法.部分采用透视或同步放大局部让学生更直观更系统的学习以上的给药方式及注意事项。 常用麻醉方法的虚拟操作：通过虚拟实验——吸入麻醉和腹腔注射麻醉，让学生熟悉并掌握常用麻药的使用及配制方法。 大小鼠取血的虚拟操作：分为摘眼球取血法，眼眶后静脉丛穿刺取血法，心脏取血，腹主动脉采血法。 大鼠处死方法的演示，脊椎脱臼法，急性失血法，麻醉致死法，气体窒息致死法，击打法。 大鼠主要脏器摘取：学生可动手摘取虚拟大鼠的主要脏器，可掌握各主要脏器的位置和摘取后的性状。 家兔的基本实验虚拟操作内容包括： 家兔麻醉方法，颈部手术包含颈部皮肤切开、分离皮下筋膜、气管插管、颈动脉插管、颈外静脉插管、颈部迷走神经、交感神经、降压神经分离等内容，家兔腹部手术包含回盲部肠系膜分离术、输尿管插管术、膀胱插管术等内容，家兔
南京信息工程大学实验（实习）报告 一、实验目的: 用最小二乘法将给定的十个点拟合成三次多项式。 二、实验步骤： 用matlab编制以函数为基的多项式最小二乘拟合程序，并用于对下列数据作三次多项式最小二乘拟合(取权函数wi=1) x -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 y -2.30 -1 -0.14 -0.25 0.61 1.03 1.75 2.75 4.42 6.94 给定直线方程为：y=1/4*x3+1/2*x2+x+1 三、实验结论： 最小二乘法：通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。最小二乘法还可用于曲线拟合。 一般地。当测量数据的散布图无明显的规律时，习惯上取n次代数多项式。 程序运行结果为： a = 0.9731 1.1023 0.4862 0.2238 即拟合的三次方程为：y=0.9731+1.1023x+0.4862*x2+0.2238*x3
-2.5 -2-1.5-1-0.5 00.51 1.52 2.5 -4-20246 81012 x 轴 y 轴 拟合图 离散点 y=a(1)+a(2)*x+a(3)*x.2+a(4)*x.3 结论： 一般情况下，拟合函数使得所有的残差为零是不可能的。由图形可以看出最小二乘解决了残差的正负相互抵消的问题，使得拟合函数更加密合实验数据。 优点：曲线拟合是使拟合函数和一系列的离散点与观测值的偏差平方和达到最小。 缺点：由于计算方法简单，若要保证数据的精确度，需要大量的数据代入计算。
《计算机仿真》上机实验报告 姓名： 学号： 2012104021 专业：测控 班级： 12级
实验一常微分方程的求解及系统数学模型的转换一．实验目的 通过实验熟悉计算机仿真中常用到的Matlab指令的使用方法，掌握常微分方程求解指令和模型表示及转换指令，为进一步从事有关仿真设计和研究工作打下基础。 二. 实验设备 个人计算机，Matlab软件。 三. 实验准备 预习本实验有关内容（如教材第2、3、5章中的相应指令说明和例题），编写本次仿真练习题的相应程序。 四. 实验内容 1. Matlab中常微分方程求解指令的使用 题目一：请用MATLAB的ODE45算法分别求解下列二个方程。要求:1.编写出Matlab 仿真程序；2.画出方程解的图形并对图形进行简要分析；3.分析下列二个方程的关系。 1．2． 1.function fun=funl(t,x) fun=-x^2;
[t,x]=ode45('fun1',[0,20],[1]); figure(1);plot(t,x); grid 2.function fun=fun2(t,x) fun=x^2; [t,x]=ode45('fun2',[0,20],[-1]); figure(2);plot(t,x); grid
题目二：下面方程组用在人口动力学中，可以表达为单一化的捕食者-被捕食者模式（例如，狐狸和兔子）。其中1x 表示被捕食者， 2x 表示捕食者。如果被捕食者有无限的食物，并且不会出现捕食者。于是有1'1x x ，则这个式子是以指数形式增长的。大量的被捕食者将会使捕食者的数量增长；同样，越来越少的捕食者会使被捕食者的数量增长。而且，人口数量也会增长。请分别调用ODE45、ODE23算法求解下面方程组。要求编写出Matlab 仿真程序、画出方程组解的图形并对图形进行分析和比较。 1.ODE45
实验五双极性不归零码 一、实验目的 1.掌握双极性不归零码的基本特征 2.掌握双极性不归零码的波形及功率谱的测量方法 3.学会用示波器和功率谱分析仪对信号进行分析 二、实验仪器 1.序列码产生器 2.单极性不归零码编码器 3.双极性不归零码编码器 4.示波器 5.功率谱分析仪 三、实验原理 双极性不归零码是用正电平和负电平分别表示二进制码１和０的 码型，它与双极性归零码类似，但双极性非归零码的波形在整个码元持续期间电平保持不变．双极性非归零码的特点是：从统计平均来看，该码型信号在１和０的数目各占一半时无直流分量，并且接收时判决电平为０，容易设置并且稳定，因此抗干扰能力强．此外，可以在电缆等无接地的传输线上传输，因此双极性非归零码应用极广．双极性非归零码常用于低速数字通信．双极性码的主要缺点是：与单极性非归零码一样，不能直接从双极性非归零码中提取同步信号，并且１码和０码不等概时，仍有直流成分。 四、实验步骤
1.按照图3.5-1 所示实验框图搭建实验环境。 2.设置参数：设置序列码产生器序列数N=128；观察其波形及功率谱。 3.调节序列数N 分别等于6 4.256，重复步骤2. 图3.5-1 双极性不归零码实验框图 实验五步骤2图 N=128
实验五步骤3图N=64 N=256
六、实验报告 （1）分析双极性不归零码波形及功率谱。 （2）总结双极性不归零码的波形及功率谱的测量方法。 实验六 一、实验目的 1.掌握双极性归零码的基本特征 2.掌握双极性归零码的波形及功率谱的测量方法 3.学会用示波器和功率谱分析仪对信号进行分析 二、实验仪器 1.序列码产生器 2.单极性不归零码编码器 3.双极性归零码编码器
动态系统建模仿真实验报告（2） 四旋翼飞行器仿真 姓名： 学号： 指导教师： 院系： 2014.12.28
1实验容 基于Simulink建立四旋翼飞行器的悬停控制回路，实现飞行器的悬停控制； 建立GUI界面，能够输入参数并绘制运动轨迹； 基于VR Toolbox建立3D动画场景，能够模拟飞行器的运动轨迹。 2实验目的 通过在 Matlab 环境中对四旋翼飞行器进行系统建模，使掌握以下容： 四旋翼飞行器的建模和控制方法 在Matlab下快速建立虚拟可视化环境的方法。 3实验器材 硬件：PC机。 工具软件：操作系统：Windows系列；软件工具：MATLAB及simulink。 4实验原理 4.1四旋翼飞行器 四旋翼飞行器通过四个螺旋桨产生的升力实现飞行，原理与直升机类似。四个旋翼位于一个几何对称的十字支架前，后，左，右四端，如图 1 所示。旋翼由电机控制；整个飞行器依靠改变每个电机的转速来实现飞行姿态控制。 图1四旋翼飞行器旋转方向示意图
在图 1 中， 前端旋翼 1 和后端旋翼 3 逆时针旋转， 而左端旋翼 2 和右端的旋翼 4 顺时针旋转， 以平衡旋翼旋转所产生的反扭转矩。 由此可知， 悬停时， 四只旋翼的转速应该相等，以相互抵消反扭力矩；同时等量地增大或减小四只旋翼的转速，会引起上升或下降运动；增大某一只旋翼的转速，同时等量地减小同组另一只旋翼的转速，则产生俯仰、横滚运动；增大某一组旋翼的转速，同时等量减小另一组旋翼的转速，将产生偏航运动。 4.2建模分析 四旋翼飞行器受力分析,如图 2 所示 图2四旋翼飞行器受力分析示意图 旋翼机体所受外力和力矩为： 重力mg , 机体受到重力沿w z -方向； 四个旋翼旋转所产生的升力i F (i= 1 , 2 , 3 , 4)，旋翼升力沿b z 方向； 旋翼旋转会产生扭转力矩i M (i= 1 , 2 , 3 , 4)。i M 垂直于叶片的旋翼平面，与旋转矢量相反。 力模型为：2i F i F k ω= ，旋翼通过螺旋桨产生升力。F k 是电机转动力系数， 可取826.1110/N rpm -?，i ω为电机转速。旋翼旋转产生旋转力矩Mi(i=1,2,3,4)，
通信系统实验报告——基于SystemView的仿真实验 班级： 学号： 姓名： 时间：
目录 实验一、模拟调制系统设计分析 -------------------------3 一、实验内容-------------------------------------------3 二、实验要求-------------------------------------------3 三、实验原理-------------------------------------------3 四、实验步骤与结果-------------------------------------4 五、实验心得------------------------------------------10 实验二、模拟信号的数字传输系统设计分析------------11 一、实验内容------------------------------------------11 二、实验要求------------------------------------------11 三、实验原理------------------------------------------11 四、实验步骤与结果------------------------------------12 五、实验心得------------------------------------------16 实验三、数字载波通信系统设计分析------------------17 一、实验内容------------------------------------------17 二、实验要求------------------------------------------17 三、实验原理------------------------------------------17 四、实验步骤与结果------------------------------------18 五、实验心得------------------------------------------27