建模与仿真实验报告（最终五篇）

第一篇：建模与仿真实验报告

重 重 庆 庆

大 大 学 学

学 学 生 生

实 实 验 报 告 实验课程名称

物流系统建模与仿真

开课实验室

物流工程实验室

学 学

院

自动化

年级

专业班

程 物流工程 2 班

学 学 生 生

姓 姓 名

段竞男

学

号

201 ２4 ９12

开 开 课 时 间

2014

至１５

学年第 二

学 学期 期 总 总

成 成 绩 绩

教师签名

自动化 学院制 《 物流系统建模与仿真 》实验报告

开课实验室：

****年**月**日 日 学院 自动化 年级、专业、班 12级物流工程2班 姓名 段竞男 成绩

课程 名称 物流系统建模与仿真 实验项目 名

称 产品测试工艺仿真与分析实验 指导教师 张莹莹 教师评语

成绩

一、实验目得 通过建立单存放区域、单处理工作台得简单模型,了解 5 个基本建模步骤。学习使用统计分析工具.二、实验原理 某工厂车间对三类产品进行检验。这三种类型得产品按照一定得时间间隔方式到达。随后，不同类型得产品被分别送往三台不同得检测机进行检测，每台检测机只检测一种特定得产品类型.其中，类型 1 得产品到第一台检测机检测,类型２得产品到第二台检测机检测,类型 3 得产品到第三台检测机检测。产品检测完毕后,由传送带送往货架区,再由叉车送到相应得货架上存放。类型 1 得产品存放在第 2 个货架上，类型 2 得产品存放在第 3 个货架上，类型 3 得产品存放在第 1 个货架上。

三、使用仪器、材料

一台ＰＣ机,ｆlｅｘsiｍ软件 四、实验步骤

1)创建模型布局

使用鼠标将需要得对象从对象库中拖放到正视图窗口中,根据需要使用鼠标改变对象位置、大小与转角。

2)连接端口

按下键盘上得“A“键，用鼠标拖放在对象间建立输出端口-输入端口连接；方向为从流出实体得对象到流入实体得对象;模型中得对象发出与接收实体需要这种连接。

3)编辑外观、设置对象行为

(1)参数窗口（Parameteｒs Ｗindｏw）

双击对象（或在右键菜单选择 Pａramｅters）

；用于对各种对象得自身特性得设置、编辑。

(2)属性窗口(Ｐropertieｓ Windｏｗ）

右键单击对象,在弹出菜单中选择 Ｐrｏpertieｓ;用于编辑与查瞧所有对象都拥有得一般性信息。

(3)模型树视图(Model Tree Vｉew)模型中得所有对象都在层级式树结构中列出;包含对象得底层数据结构;所有得信息都包含在此树结构中。

4)重置运行

(1)重置模型并运行

(2)控制仿真速度（不会影响仿真结果）

检测机器 1 检测机器 2 检测机器 3 传送带 传送带 传送带 货架 1 货架 2 货架 3 产品 1 产品 2 产品 3

(3)设置仿真结束时间

5)观察结果

(1)使用“Statiｓtics”(统计)菜单中得 Reｐｏrts anｄ Statiｓtics(报告与统计）生成所需得各项数据统计报告.(2)其她报告功能包括：对象属性窗口得统计项;记录器对象;可视化工具对象;通过触发器记录数据到全局表。

五、实验过程原始记录((数据、图表、计算等))

1、运行结果得平面视图:

2、运行结果得立体视图

３、运行结果得暂存区数据分析结果图:

第一个暂存区

第二个暂存区 由报表分析可知 5 次实验中，第一个暂存区得平均等待时间为 1１、４6，而第二个暂存区得平均等待时间为１3、0２，略大于第一个暂存区,由此可见,第二个暂存区得工作效率基本上由第一个暂存区决定。

４、运行结果三个检测台得数据分析结果图,三个检测台得 statｅ饼图：

（１)处理器一：

由实验结果分析可得,处理器一只有５3%得时间处于工作状态，有３２、3％得时间就是处于闲置状态,并且该处理器得准备时间较长，占总时间得 14、7％,这些数据表明该处理器得运行速度完全能满足,甚至超过系统得要求，可以适当得选择更处理速度慢一点得处理器来降低系统成本。

（2）处理器二:

由实验结果分析可得,处理器二只有 16、９％得时间处于工作状态,有６6％得时间就是处于闲置状态，并且有 1７、1％得时间处于准备时间,以上数据说明处理器二闲置时间过长，工作效率低，不能很好地配合物料二得到达速度。

(3）处理器三

由实验结果分析可得，处理器三得只有１６、9％得时间处于工作状态，有 66、1%得时间就是处于闲置状态,并且有 16、9%得时间处于准备时间，以上数据说明处理器三得工作效率低,不能很好地配合物料三得到达速度，可以适当得降低处理器三得处理速度。

六、实验结果及分析

１))对得到得数据做简单分析，提出改进措施。

答:通过对实验数据分析,发现现有得机器设备得设置基本能满足系统得要求，但工作效率低，大大浪费了设备得工作能力,可以适当得提高物料得到达速度，或者降低三台处理器得处理速度.2 2）

讨论：本实验根据三个处理器得统计信息,通过状态图分析各处理器得工作状态，通过暂存区材料得平均等待时间，分析这个检测流程得效率如何？就是否存在瓶颈？如果存在,怎样才能改善整个系统得绩效呢？ 答：由报表分析可知５次实验中，第一个暂存区得平均等待时间为 11、46，而第二个暂存区得平

均等待时间为 13、02,略大于第一个暂存区,由此可见，第二个暂存区得工作效率基本上由第一个暂存区决定。处理器一只有５3%得时间处于工作状态，有 32、３％得时间就是处于闲置状态,并且该处理器得准备时间较长,占总时间得１４、7%,这些数据表明该处理器得运行速度完全能满足,甚至超过系统得要求,可以适当得选择更处理速度慢一点得处理器来降低系统成本。由实验结果分析可得，处理器二只有 16、９%得时间处于工作状态,有 66%得时间就是处于闲置状态,并且有 17、1％得时间处于准备时间，以上数据说明处理器二闲置时间过长，工作效率低，不能很好地配合物料二得到达速度。由实验结果分析可得,处理器三得只有 1６、9%得时间处于工作状态，有 66、1%得时间就是处于闲置状态，并且有１6、9%得时间处于准备时间，以上数据说明处理器三得工作效率低,不能很好地配合物料三得到达速度，可以适当得降低处理器三得处理速度。整体来瞧,整个检测流程效率不高，但就是并不存在瓶颈，能够达到系统得要求,但工作效率低，大大浪费了设备得工作能力,可以适当得提高物料得到达速度，或者降低三台处理器得处理速度.《物流系统建模与仿真》实验报告

开 ：

开 课 实 验 室 ：

年 年

月 月

日 日 学院 自动化 年级、专业、班 １２级物流工程２班 姓名 段竞男 成绩

课程 名称 物流系统建模与仿真 实验项目 名

称 多产品多阶段制造系统仿真与分析实验 指导教师 张莹莹 教师评语

成绩

一、实验目得 假定在保持车间逐日连续工作得条件下，对系统进行 365 天得仿真运行（每天按８ 小时计算)，计算每组机器队列中得平均产品数以及平均等待时间。通过仿真运行,找出影响系统得瓶颈因素，并对模型加以改进。

二、实验原理 某制造车间由 5 组机器组成，第 1,2,3，4,5 组机器分别有３，２,4,3，1 台相同得机器。这个车间需要加工三种原料，三种原料分别要求完成４、3 与５ 道工序，而每道工序必须在指定得机器组上处理,按照事先规定好得工艺顺序进行。

概念模型参考如下：

三、使用仪器、材料

一台ＰC 机，flexsｉm 软件 四、实验步骤

1)创建模型布局 使用鼠标将需要得对象从对象库中拖放到正视图窗口中，根据需要使用鼠标改变对象位置、大小与转角。

2)连接端口 按下键盘上得“A”键，用鼠标拖放在对象间建立输出端口—输入端口连接;方向为从流出实体得对象到流入实体得对象；模型中得对象发出与接收实体需要这种连接。本实验此步骤为关键环节，连接线较复杂，故连接端口时应注意保持清晰思路.机器 1 机器 1 机器 1 暂存区 机器 2 机器 2 暂存区 机器 3 机器 3 机器 3 机器 3

暂存区 机器 4 机器 4 机器 4 暂存区 机器 5

暂存区 1 2 3

3)定义对象参数(1)定义Ｓource 在模型中，共有 3 个Ｓource 实体，每个Ｓource 对应一类原料,也就就是说，一个 Source 生成一类原料.我们需要设置每个 Souｒce 实体，使得每类原料得到达间隔时间满足系统得要求。

(2)分别按照要求定义机器组１、机器组 2、机器组 3、机器组 4、机器组 5、暂存区参数。

4)模型运行

(1)设置 Expeｒiｍenter

(2)重置模型并运行 5)观察结果 仿真进行过程中，可以瞧到红、黄、蓝三种不同颜色得原料从系统中流过,经过不同机器组得加工,最后离开系统。

仿真运行到１７52０0 单位时间得时候，自动停止。

五、实验过程原始记录((数据、图表、计算等）

1、运行结果得平面视图: 2、提供运行结果得立体视图：

3、提供运行结果五个暂存区数据分析结果图

由数据图可知第一个暂存区得为 物料输入量为 888 ８,平均数量为 0、６3, 最大等待时间达到了 241、、94，平均等待时间为 12、5，由此可见暂存区一得货物 状态较为空闲，不就是 系统得瓶颈。

为 由数据图可知第二个暂存区得物料输入量为 1117 ３,平均容量为 1 ９67、0 ４, 最大等待时间为 630１ １7、1４ ４为，平均等待时间为 30３ ３ ８８、42, 由此可见暂存区二得货物周转率较高, 处于比较忙碌得状态，非常 容易成为系统得瓶颈。

由数据图可知第三个暂存区得物料输入量为５ ５8 ３５,平均数量为 0、11，最大等待时间达到了 1２ ２4、25,平均等待时间为 3、４0, 由此可见暂存区三得效率较高，较为适应系统得要求。

为 由数据图可知第四个暂存区得物料输入量为 3604，平均数量为 5、6２ ２，最大等待时间为１ １64 ２、37，平均等待时间为 273、55，由此可见暂存区四货物周转率高，并且运行效率高, 能够适应系统得要求，安排合理。

由数据图可知第五 个暂存区得物料为 输入量为 7 ２08,平均数量为２ ２43、75, 最大等待时间为 1 ０９3６、85,平均等待时间为 5956, 由此可见暂存区五 容易造成货物得堆积与等待，不 能够适应系统得要求，安排 不够 合理, 为该系统得瓶颈.六、实验结果及分析

对得到得数据做简单分析，提出改进措施：

暂存区 1 2 3 4 5平均等待时间 12、50 30388、42 3、40 273、55 5956

平均数量 0、63 1967、04 0、11 5、62 243、75 由数据图可知第一个暂存区得物料输入量为 8８88，平均数量为 0、63，最大等待时间达到了241、94，平均等待时间为１2、5,由此可见暂存区一得货物状态较为空闲,不就是系统得瓶颈.由数据图可知第二个暂存区得物料输入量为１１173,平均容量为 1967、04，最大等待时间为 63017、14,平均等待时间为 30388、42，由此可见暂存区二得货物周转率较高，处于比较忙碌得状态，非常容易成为系统得瓶颈.由数据图可知第三个暂存区得物料输入量为 5835,平均数量为０、11，最大等待时间达到了 1２4、２5,平均等待时间为 3、4０，由此可见暂存区三得效率较高,较为适应系统得要求。由数据图可知第四个暂存区得物料输入量为 36０4,平均数量为 5、62,最大等待时间为 1642、３7，平均等待时间为 2７3、５5，由此可见暂存区四货物周转率高,并且运行效率高，能够适应系

统得要求,安排合理.由数据图可知第五个暂存区得物料输入量为 72０8，平均数量为 243、75,最大等待时间为 1０936、85,平均等待时间为 595６，由此可见暂存区五容易造成货物得堆积与等待,不能够适应系统得要求,安排不够合理,为该系统得瓶颈.由数据分析可知暂存区二、五就是系统得瓶颈所在，应该提高暂存区二、暂存区五得处理前一道工序速度,即增加工序一与工序二得机器数量，以此来提高工序 1、2 得处理速度，或者调整物料加工得顺序,从而提高整个系统得运营效率。

《物流系统建模与仿真》实验报告

开课实验室：

****年**月**日 日 学院 自动化 年级、专业、班 １2 级物流工程 2班 姓名段竞男 成绩

课程 名称 物流系统建模与仿真 实验项目 名

称 混合流水线系统仿真与分析实验 指导教师 张莹莹 教师评语

成绩

一、实验目得 主要掌握单台处理器在处理多种产品时对于处理顺序以及时间参数得设置.熟悉先进先出这种存储模式得控制方法.二、实验原理 多对象流水线生产有两种基本形式。一种就是可变流水线，其特点就是:在计划期内，按照一定得间隔期,成批轮番生产多种产品；在间隔期内，只生产一种产品,在完成规定得批量后,转生产另一种产品。另一种就是混合流水线,其特点就是：在同一时间内,流水线上混合生产多种产品。按固定得混合产品组组织生产,即将不同得产品按固定得比例与生产顺序编成产品组。一个组一个组地在流水线上进行生产。

三、使用仪器、材料

ＰC 机一台,flexsiｍ 软件 四、实验步骤

1)创建模型布局

从左边得实体库中依次拖拽出所有实体（一个 Source，５个 Quｅｕe,1２个 Pｒｏceｓｓoｒ，一个Ｃonveｙor，一个Ｓｉnｋ)放在右边模型视图中，调整至适当得位置,如图所示:

2)连接端口

根据流动实体得路径来连接不同实体得端口.按住键盘上得“A”键,与前面章节得操作一样,按上图中得箭头所指向依次连接各个实体.分别(注意方向)从 Sｏurｃe 连到 GeneｒalQueｕe，GeneｒalQuｅue 连到ＧenｅｒaｌMachine1，GｅnｅｒalＱuｅuｅ连到 GenerａlMａchine2, GeneralＱueue 连到 GenerａlMaｃｈine3，ＧｅｎｅralMａｃｈiｎe1 连到 DrillinｇQueue，GeneralMacｈine2 连到 DriｌlingQｕeue，GｅｎｅrａlＭａcｈinｅ3 连到 DrillｉngQueuｅ，DrillｉｎgQueue 连到 DｒillinｇＭａchine1，DriｌlinｇQueｕｅ连到 DrillingMａchｉne２,DrillinｇQueue连到DriｌlingMachine３,DriｌlｉngＭachine1连到ＭillｉｎgQｕeuｅ,ＤriｌlinｇMaｃhｉｎe2 连到 MilliｎgQｕeue,DrｉlｌinｇMacｈine3 连到ＭilｌｉngQｕeue,MiｌｌinｇＱueｕe 连到 MilｌｉnｇＭachine１，MillingＱueｕｅ连到ＭilｌｉngＭaｃhine2,ＭillinｇMaｃhiｎｅ1 连到 GriｎdiｎgQuｅue,ＭiｌlingＭａchine２连到 GrindｉngＱｕeue，ＧrｉndｉｎgＱｕｅｕｅ连到 GrｉnｄingＭａchiｎe1，GrinｄinｇＱueue 连到 GｒindinｇMａｃｈｉne２,GrｉndingＱｕｅue 连到ＧｒindｉnｇMachｉｎｅ3，GｒiｎdiｎｇMachｉｎe1 连到 T

ｅｓｔｉngQueｕe,GrｉｎdingＭachｉne2连到TｅｓtiｎgＱueue,GrindｉｎｇMachｉｎe3连到TesｔingＱueue,TestingQueue 连到 TｅsｔingＭachine，ＴeｓtｉngＭachine 连到 Cｏｎveyor，Ｃoｎveyor 连到Ｓｉnk。完成后，如图所示：

3)定义对象参数

分别按照要求定义 Sourｃe、各机器工位、暂存区参数。

4)模型运行

(1)重置模型并运行(2)加快仿真模型运行速度 如果我们只就是关心仿真结果,而对仿真得过程不感兴趣，则我们可以加快仿真速度,迅速得到结果.图 3、１仿真速度控制比例条 如图 3、1，鼠标左键一直按住比例尺，移动到合适得比例位置,以便迅速得到结果。

5)观察结果

仿真进行过程中，可以瞧到红、绿、蓝三种不同颜色得产品从系统中流过,经过不同机器组得加工，最后离开系统,如下图:

6)结果分析

当仿真运行自动结束后，我们打开 Flexｓｉｍ得工具栏里得 Stａts 目录下得Ｓtandａrd Report 选项，如图所示:

我们通过来增加、以及来减少需要输出得报告内容,使得报告包含以上所列得 5 个部分得数据:iｄlｅ就是空闲时间，pｒoｃessiｎg 就是工作时间，blockeｄ就是产品在设备等待时间，ｓｔatｓ_staytiｍｅavｇ就是平均停留时间，stats_input 就是输入产品数,stats_output 就是输出产品数,设置完成后,如图所示:

点击后生成如图所示表格:

从上表中我们可以很方便得瞧到总运行时间就是 7728min，以及各个设备得输入输出产品数,处理时间等信息.录记始原程过验实、五ﻬ 五、实验过程原始记录((数据、图表、计算等）

1、运行结果得二维平面图

２、运行结果得三维立体图

4、投产按照 1、2、3 得顺序进行得仿真报告 Flexsim Summary Report Time: 7728

Object Class stats_output stats_staytimemin stats_input idle blocked processing

Source1 Source 1700 0 0 0 7720 0 Queue2 Queue 1700 0 1700 0 0 0 Processor3 Processor 500 4 500 594 0 2400 Processor4 Processor 600 4 600 297 0 2700 Processor5 Processor 600 4 600 297 0 2700 ConCon17010 1700 0 0

veyor6 veyor 0 0 Processor7 Processor 600 4 600 203 0 2800 Processor8 Processor 600 4 600 203 0 2800 Processor9 Processor 500 4 500 599 0 2400 Queue10 Queue 1700 0 1700 0 0 0 QueQue1700 1700 0 0

ue11 ue 0 0 Processor12 Processor 850 3 850 10 0 3050 Processor13 Processor 850 3 850 11 0 3050 Processor14 Processor 567 4 567 796 0 2268 ProcessorProcessor 567 4 567 797 0 2268 Processor16 Processor 566 4 566 798 0 2264 Queue17 Queue 1700 0 1700 0 0 0 Queue18 Queue 1700 0 1700 0 0 0 Processor19 Processor 1700 1 1700 18 0 7700 Sink20 Sink 0 0 1700 0 0 0 由报表分析可知，按照 1 1、2 2、３得顺序投产得总得处理时间为７7 7 ２8 8。

投产按照１、3、2 得顺序进行得仿真报告:

Flexsim Summary Report Time: 7728

Object Class stats_input stats_output stats_staytimeavg idle processing blocked Source1 Source 0 1700 0 0 0 7720 Queue2 Queue 1700 1700 4、588235 0 0 0 ProcProc502 502 4、79586 2408 0

essor3 essor 6813 Processor4 Processor 600 600 4、5 298 2700 0 Processor5 Processor 598 598 4、501672 306 2692 0 Queue6 Queue 1700 1700 0、411765 0 0 0 Processor7 Processor 600 600 4、666667 203 2800 0 PrPr60604、20280

ocessor8 ocessor 0 0 666667 3 00 Processor9 Processor 500 500 4、8 599 2400 0 Queue10 Queue 1700 1700 28、32294 0 0 0 Processor11 Processor 850 850 3、588235 10 3050 0 ProcessorProcessor 850 850 3、588235 11 3050 0 Queue13 Queue 1700 1700 0 0 0 0 Processor14 Processor 567 567 4 796 2268 0 Processor15 Processor 567 567 4 797 2268 0 Processor16 Processor 566 566 4 798 2264 0 QueuQueu1700 1700 23290 0 0

e17 e、648 Conveyor19 Conveyor 1700 1700 10 0 0 0 Sink20 Sink 1700 0 0 0 0 0 Processor21 Processor 1700 1700 4、529412 18 7700 0 由报表分析可知，按照 1、3、2 得顺序投产得总得处理时间为 77２８.投产按照 2、1、3 得顺序进行得仿真报告：

Flexsim Summary Report Time:

7725 Object Class stats_input stats_output stats_staytimeavg idle processing blocked Source1 Source 0 1700 0 0 0 7720 Queue2 Queue 1700 1700 4、588235 0 0 0 Processor3 Processor 503 503 4、795229 588 2412 0 Processor4 Processor 600 600 4、5 299 2700 0 Processor5 Processor 597 597 4、502513 311 2688 0

Queue6 Queue 1700 1700 0、408824 0 0 0 Processor7 Processor 600 600 4、666667 205 2800 0 Processor8 Processor 600 600 4、666667 205 2800 0 Processor9 Processor 500 500 4、8 605 2400 0 Queue10 Queue 1700 1700 28、16529 0 0 0 Processor11 Processor 851 851 3、587544 11 3053 0 Processor12 Processor 849 849 3、588928 15 3047 0 Queue13 Queue 1700 1700 0 0 0 0 Processor14 Processor 567 567 4 796 2268 0 Processor15 Processor 567 567 4 798 2268 0 Processor16 Processor 566 566 4 804 2264 0 Queue17 Queue 1700 1700 2317、985 0 0 0 Conveyor19 Conveyor 1700 1700 10 0 0 0 Sink20 Sink 1700 0 0 0 0 0 Processor21 Processor 1700 1700 4、529412 15 7700 0 由报表分析可知，按照 2 2、１、３得顺序投产得总得处理时间为 772 ５.析分及果结验实、六ﻬ 六、实验结果及分析1）

实验结果: :由以上报表分析可知 6 种投产顺序中，按照 2、1、3 得顺序投产得总得处理时间最短为７7２5，故以这种投产顺序作为最佳投产方式。

２）讨论: :系统还存在很多可以改善得地方，请指出有哪些地方就是有待改善得，理由就是什么.答：可以适当得降低处理器得运行速度，除了第三组得处理器空闲时间较短外，其余各组处理器得空闲时间过长。尤其就是，其中暂存区 17 得平均等待时间过长,可以适当得调整运行速度，以提高整个系统得运行效率。

实验报告打印格式说明 1.标题:三号加粗黑体 2.开课实验室:5 号加粗宋体 3.表中内容:(1)标题:５号黑体(2)正文：5 号宋体

4.纸张:１6 开(20cm×2６.5ｃm)5.版芯 上距：2cm 下距：2cm 左距:2.8ｃｍ 右距:２．8cｍ 说明: 1、“年级专业班”可填写为“00 电子 1 班",表示 2000 级电子工程专业第１班.２、实验成绩可按五级记分制（即优、良、中、及格、不及格)，或者百分制记载，若需要将实验成绩加入对应课程总成绩得，则五级记分应转换为百分制。

第二篇：仿真与建模实验报告

实验报告

13工业工程2班 李伟航 13工业工程2班

实验10

一、实验目的: 1.学习库存系统查库与订货处理的结构建模方法

2.学习用Equation模块、Equation（I）模块读写数据库的方法 3.学习用Equation模块、Equation（I）模块进行编程计算的方法

二、实验问题

1.打开上次实验你保存的文件（这个文件要保存好，下次实验还要使用），然后根据以上视频，进行操作实验。

2.简述用Equation模块计算订货量的程序逻辑。

3.简述用Equation（I）模块计算并累加总订货成本到数据库中的程序逻辑。4.Equation模块与Equation（I）模块有何不同？

5.在本案例的假设前提下，在一笔订货的在途货物运输期间（即提前期期间），会不会再次发出订货指令？或者换句话说，每次查库时，会不会有已订但未到的货？为什么？

三.实验过程

1.生成查库员（查库信号）

用Create模块每天生成一个库存检查员实体（实际代表一个查库信号）Create模块具体设置如下图：

2.判断是否需要订货

利用Select Item Out模块、Equation模块和Simulation Variable模块检查库存，并判断是否需要订货。若需要，就将库存检查员实体发送到Select Item Out模块的上端口输出进行后续处理；若不需要订货，就将库存检查员实体发送到下端口输出，简单地离开系统。其中，Equation模块的设置如下图。

3.无需订货的处理

由上一步Equation中设置可知，当s=1时，即无需订货的情况下，直接将库存检查员实体从Select Item Out模块下端口输出，通过Exit模块离开系统。Select Item Out模块设置如下，当s=0时从上端口输出，否则从下端输出。

4.订货处理-建立模型

当s=0时，即需要订货的情况下，库存检查员实体有Select Item Out模块上端口输出，后续订货处理模型如下图：

5.计算订货量

用一个Equation模块获取数据库中的当前库存（kc）和最大库存(ds)，计算订货量dh。Equation模块设置如下：

6.设置订货量属性

将第5步Equation模块计算得到的订货量（dh）输出给Set模块，将订货量赋值给检查实体dh属性。Set属性设置如下图。

7.用Equation（I）模块计算并累加总订货成本到数据库中

因此用Equation（I）模块获取数据库中固定费（gdf）和可变订货费（kbf），再根据订货量算出本次订货成本，把它累加到总订货费用（toc）并写入数据库中。具体Equation（I）模块设置如下图。

8.设置提前期

从下订单到所订货物入库的这段时间为订货提前期，这里用Activity模块表示，设置该提前期服从0.5-1天的均匀分布。订货提前期的设置如下。

9.到货后更新数据库中的库存数据

库存检查实体离开Activity，表示提前期结束，货物到达，到货后用Get模块获取订货量，用Write模块将订货量累加到数据库中的sc表的当前库存（kc）中，根据订货动态增加当前库存。Write模块设置如下。

10.观察运行结果 运行模型，观察数据库sc表中的库存（kc）由变为13，累计订货费用达到了11125。

四.实验问题的分析解答

1.上述过程即为实验步骤。本实验实验模型如下：

2.答：用Equation模块计算订货量的程序逻辑：是当检测到库存小于订货点xs时，用最大库存ds减去目前库存，即为订货量，即dh=ds-kc；过程如实验步骤5所示。

3.答：用Equation（I）模块计算并累加总订货成本到数据库中的程序逻辑是：因为每订货一次的成本=固定费（gdf）+可变订货费（kbf）*订货量（dh），将每次的订货成本累加到总订货成本当中去，即toc=toc+gdf+kbf*dh，即更新了总订货成本。具体步骤如步骤7所示。

4.答：Equation模块是由value库中添加的值模块，是用来读取传递系统产生的值（如产生的随机数）并进行相应的逻辑运算，而对于实体传来的属性则不可以读取；而相反的，Equation（I）模块是由Item库中添加的实物模块，用来读取和传递实体属性并进行相应的逻辑运算。两者都可以读取数据可中的数据。

5.答：一笔订货的在途货物运输期间，不会再次发出订货指令。因为在Create模块中设置了检查实体到来的间隔为1天，而订货提前期服从0.5-1的均匀分布，所以不会出现每次查库有已订但未到的货的情况。

实验11

一、实验目的

1.通过实验理解各种库存性能指标的含义 2.学习用Read模块读取数据库数据的方法

3.学习用Max&Min模块、Mean模块、Equation模块以及各种计算模块计算和采集库存系统性能指标的方法

二.实验问题

1.打开上次实验你保存的文件（这个文件要保存好，下次实验还要使用），然后根据以上视频，进行操作实验。

2.本实验中，用read模块读数据库数据时，在其option页要做何设置？ 3.用Mean Variance模块计算平均每天总成本的均值和置信区间时，模块对话框要做何设置？

4.利用本次实验建立的模型，实验比较以下各组（s，S）下的平均每天总成本：(20, 40)(20, 60)(20, 80)(20, 100)(40, 60)(40, 80)(40, 100)(60,80)(60,100)(80,100)对每种情况重复运行20次，写出各项配置下的平均每天总成本的均值和置信区间。并写出以上最优的（使得平均每天总成本最低）的（s，S）。

5.添加必要的模块，分别计算平均每天缺货成本、平均每天储存成本、平均每天订货成本的均值和置信区间。6.用Plotter，Discrete Event模块绘制当前库存水平的波动曲线，并同时绘制一条高度为20（即订购点）的水平直线，和一条高度为0的水平直线（连接constant模块到plotter DE模块），观察并库存曲线的波动情况以及和两条直线的关系，根据你的观察，缺货情况经常发生吗？

三.实验过程

1.在库存发生变化时读取数据库中当前库存水平

每当库存数据发生变化时，用Read模块读取库存（kc）的值，Read设置如下：

2.计算平均每天储存成本 当库存大于0时，通过max模块与0比较得到当前库存值，并通过Mean&Varience模块（运行10次）计算平均每天的库存，再用Equation模块计算平均每天储存成本。计算平均每天储存成本的模型与设置如下所示：

3.计算平均每天缺货成本 当库存小于0时，通过max模块与0比较得到当前缺货数，并通过Mean&Varience模块（运行10次）计算平均每天的缺货数，用Equation模块（缺货数取反）计算平均每天缺货成本。计算平均每天缺货成本的模型与设置如下所示：

4.计算平均每天订货成本

利用Read模块读取运行完后数据库中的总订货成本，再用Math模块除以系统运行当前时间，得到平均每天订货成本。平均每天订货成本模型和设置如下：

5.计算平均每天总成本

利用Math模块将第2、3、4步所计算出的平均每天储存成本，平均每天缺货成本和平均每天订货成本相加即得到平均每天总成本。

6.计算平均每天总成本的均值和置信区间

在Math模块后用Mean&Varience模块计算运行10次之后的平均每天总成本的均值和置信区间。

四.实验问题的分析解答

1.上述过程即为实验步骤。本实验实验模型如下：

2.答：由于本实验要求在库存发生变化时读取数据库的库存水平，因此在Read模块中的Option选项中，选择Discrete event-read data during when 以及勾中data sources change，意思即为数据发生变化时读取数据库。设置如下：

3.答：用Mean&Varience模块计算平均每天总成本的均值和置信区间时，要勾中Calculate for multiplte simulations，意思是计算出运行多次情况（这里为10）次时平均每天总成本的均值和置信区间。

4.答： 10种情况的运行结果如下，其中最优的（s，S）为（20,60）。

（20，40）（20，60）（20，80）

均值：125.8870元 均值：120.3812元 均值：121.4371元

置信区间（125.8870±2.2602）置信区间（120.3812±1.8750）置信区间（121.4371±1.0904）

（20，100）（40,60）（40，80）

均值：128.1696元 均值：127.3233元 均值：127.5943元

置信区间（128.1696±1.1935）置信区间（127.3233±1.2890）置信区间（127.5943±0.8358）

（40,100）（60,80）（60，100）

均值：134.0826元 均值：146.2447元 均值：147.6999元

置信区间（134.0826±1.2611）置信区间（146.2447±1.3544）置信区间（147.6999±1.1276）

（80,100）

均值：167.5711元 置信区间（167.5711±1.2291）

5.答：用3个Mean&Varience模块计算平均每天存储成本、缺货成本和订货成本的均值和置信区间，Mean&Varience模块勾中Calculate for multiplte simulations。添加模块后的模型以及运行结果如下：

平均每天存储成本：平均每天缺货成本

均值约为9.5140元 均值约为16.2608元 置信区间为（9.5140±0.2234）置信区间为（16.2608±0.9608）

平均每天订货成本：

均值约为98.7903元

置信区间为（98.7903±1.3140）

6.答： Plotter Discrete Event模块设置和运行后的库存水平波动曲线如下图所示，由曲线图可以看出库存水平在0和20之间以及上下变动。由图可以看出，位于0以下的曲线分布比较多，这意味着缺货的情况发生的比较频繁；曲线在0和20之间的分布也较多，这意味着需要订货的情况也发生的比较多。

实验12 一.实验目的: 1.学习库存仿真优化方法

2.学习如何设置模块表格的内容为决策变量

二.实验问题

1.打开上次实验你保存的文件，然后根据以上视频，建立优化模型。

2.在优化模块（Optimizer）中，整数型决策变量和连续型（实数型）的输入方法有何不同？

3.请在你的模型上，使用Extendsim软件提供的优化器来寻找（s，S）的最佳设置（总成本最小）。令s在1和99之间取值（步长为1，即为整数），S在2和100之间取值（步长为1，即为整数）。要注意s和S必须是整数而且满足s < S。a.优化器参数由先选择Quicker Defaults，写出优化结果（s，S）和平均每天总成本。

b.再将优化器参数由选择Better Defaults，写出优化结果（s，S）和平均每天总成本。

4.在上题（第2题）的基础上，通过将库存检查间隔（Evaluation Interval，目前为1天）作为变量加入到优化变量集中，来研究在每天开始工作时查看并补充（当需要时）库存是否是最佳方案，让该值在半天到5天之间连续取值，s和S的取值情况与上题相同。应用优化器求取最优设置。

a.先优化器参数选择Quicker Defaults，写出优化结果（s，S）、Evaluation Interval和平均每天总成本。

b.再将优化器参数选择Better Defaults，写出优化结果（s，S）、和平均每天总成本。

三.实验过程

1.用Data Init模块初始化最小最大库存(s，S)

2.设置Optimizer模块目标函数和决策变量 从value库中将Optimizer模块放进模型中，在查库与订货处理模型中将决策变量订货点(xs，1~99)和最大库存(ds，2~100)和输出变量平均每天总成本（toc）克隆拖放到Optimizer模块上，然后建立目标函数方程和约束条件。Optimizer模块设置如下：

3.设置Optimizer模块约束方程

在Optimizer模块中的Constraints中添加约束方程，约束条件为s和S必须是整数而且满足s < S，约束条件的代码设置如下：

4.设置Optimizer模块运行参数

在Optimizer模块下的Run Parameters中，由于模型为随机模型，所以单机Random mode下的Quicker Defaults按钮快速设置所有优化参数（速度快精度低），然后点击New Run。运行结束后选择Better Default按钮设置规模更大的优化参数（耗时长精度高），并比较两次的结果。

5.运行优化,查看结果

通过Optimizer模块下的Results可以看MinCost数值的变化，运行结束时最顶行会给出最优解。分别用Quicker Defaults参数和Better Defaults参数运行两次比较结果。结果见实验结论3。

四.实验问题的分析解答

1.实验过程如上述实验步骤所示。

2.答：整数型决策变量的输入方法为输入范围的时候不要输入小数点，如1；而实数型的输入方法为输入范围时输入带有小数点的数字，如1.0。

3.答：选择Quicker Defaults参数，（s，S）优化结果为（27,43），平均每天总成本约为122.37627元。

Quicker Defaults

选择Better Defaults参数，（s，S）优化结果为（20,58），平均每天总成本约为130.8615元。

Better Defaults

4.答：将库存检查间隔ei克隆拖入Optimizer模块上，在订货模型create模块后加队列防止其堵塞，ei的范围设置如下。分别运行Quicker Defaults参数Better Defaults进行优化。

 Quicker Defaults做法

第一步：选择Quicker Defaults参数,并开始运行

第二步：分析结果。

（s，S）优化结果为（9,74），Evaluation Interval（ei）为1.0186，平均每天总成本约为123.0515元。

Quicker Defaults

 Better Defaults做法

第一步：选择Better Defaults参数

第二步：分析结果。（s，S）优化结果为（38,54），Evaluation Interval（ei）为1.0714，平均每天总成本约为,120.19126元。

Better Defaults

第三篇：仿真实验报告

仿真软件实验

实验名称：基于电渗流的微通道门进样的数值模拟

实验日期：2013.9.4一、实验目的1、对建模及仿真技术初步了解

2、学习并掌握Comsol Multiphysics的使用方法

3、了解电渗进样原理并进行数值模拟

4、运用Comsol Multiphysics建立多场耦合模型，加深对多耦合场的认识

二、实验设备

实验室计算机，Comsol Multiphysics 3.5a软件。

三、实验步骤

1、建立多物理场操作平台

打开软件，模型导航窗口，“新增”菜单栏，点击“多物理场”，依次新增：“微机电系统模块/微流/斯 托 克 斯 流（mmglf）”

“ACDC模块/静态，电/传导介质DC（emdc）”

“微 机 电 系 统 模 块/微流/电动流（chekf）”

2、建立求解域

工作界面绘制矩形，参数设置：宽度6e-5，高度3e-6，中心（0,0）。复制该矩形，旋转90°。两矩形取联集，消除内部边界。5和9两端点取圆角，半径1e-6。求解域建立完毕。

3、网格划分

菜单栏，网格，自由网格参数，通常网格尺寸，最大单元尺寸：4e-7。

4、设置求解域参数

求解域模式中，斯托克斯流和传导介质物理场下参数无需改动，电动流物理场下，D各向同性，扩散系数1e-8，迁移率2e-11，x速度u，y速度v，势

能V。

5、设置边界条件

mmglf—入口1和7边界“进口/层流流进/0.00005”

出口5和12边界“出口/压力，粘滞应力/0”；

emdc—入口1和7边界“电位能/10V”

出口5和12边界“接地”

其余边界“电绝缘”；

chekf—入口1“浓度/1”，7“浓度/0”

出口5和12“通量/向内通量-nmflux_c_chekf”

其余边界“绝缘/对称”。

6、样品预置

（1）求解器参数默认为稳态求解器，不用修改。

（2）求解器管理器设置求解模式：初始值/初始值表达式，点变量值不可解和线

性化/从初始值使用设定。

（3）首先求解流体，对斯托克斯流求解，观察求解结果，用速度场表示。

（4）再求解电场，改变求解模式，点变量值不可解和线性化/当前解，对传导介

质DC求解，观察求解结果，用电位能表示。

（5）再求解电动流，不改变求解模式，观察求解结果，用电动流浓度表示。

7、样品上样

（1）改变emdc进口，边界7电位能由10改为3。对传导介质DC求解，结果用

电位能表示。

（2）改变chekf进口，7边界改为“通量/向内通量-nmflux_c_chekf”

；求解域

中x速度和y速度改为0去除载流作用；求解器设置改为瞬态求解器，时间改为“0:0.00001:0.00001”。求解模式全部使用当前解，对电动流求解，结果用浓度表示。

再求两次解，完成上样。

8、分离样品

（1）改变chefk进口，7边界“浓度/0”，1边界“浓度/-nmflux_c_chekf”。

（2）改变cmdc进口，7边界“电位能/10”，1边界“电位能/3”。

（3）重新求解电场。求解模式为初始值表达式和当前解，对传到介质DC求解，结果用电位能表示。

（4）样品分离求解。求解模式全部为当前解，对电动流求解，结果用浓度表示。

四、实验结果

五、讨论

在本次试验中，每一步操作都必须严格正确，而且参数的把握也一定要

到位，只有对每一步的设置做到精确无误，才能保证最后的实验结果。我在样品上样时一直未能获得良好的上样结果，发现对瞬态求解器的时间比例进行修改，可以获得良好上样结果，同时，在样品分离改变chefk左进口浓度时发现修改数值导致结果错误，遂未修改浓度，得到了正确结果。因此，一定要在实验时对参数正确设置。

通过对仿真实验课程的学习，及本次试验，我体会到仿真技术对于实验的帮助非常巨大，使得实验室进行的许多实验可以通过计算机模拟直接完成，节省了资源消耗，并极大地提高了实验效率。本课程的学习也让我了解到了仿真及建模技术的要领。我也基本掌握了Comsol Multiphysics

这款软件，我相信在今后我会将我对本课程的学习运用到实际中。

第四篇：动态系统建模仿真实验报告四旋翼仿真

动态系统建模仿真

实验报告（2）

四旋翼飞行器仿真

2012

1实验内容

基于Simulink建立四旋翼飞行器的悬停控制回路，实现飞行器的悬停控制；

建立UI界面，能够输入参数并绘制运动轨迹；

基于VR

Toolbox建立3D动画场景，能够模拟飞行器的运动轨迹。

2实验目的通过在Matlab

环境中对四旋翼飞行器进行系统建模，使掌握以下内容：

四旋翼飞行器的建模和控制方法

在Matlab下快速建立虚拟可视化环境的方法。

3实验器材

硬件：PC机。

工具软件：操作系统：Windows系列；软件工具：MATLAB及simulink。

4实验原理

4.1四旋翼飞行器

四旋翼飞行器通过四个螺旋桨产生的升力实现飞行，原理与直升机类似。

四个旋翼位于一个几何对称的十字支架前，后，左，右四端，如图

所示。旋翼由电机控制；整个飞行器依靠改变每个电机的转速来实现飞行姿态控制。

图1四旋翼飞行器旋转方向示意图

在图

中，前端旋翼

和后端旋翼

逆时针旋转，而左端旋翼

和右端的旋翼

顺时针旋转，以平衡旋翼旋转所产生的反扭转矩。

由此可知，悬停时，四只旋翼的转速应该相等，以相互抵消反扭力矩；同时等量地增大或减小四只旋翼的转速，会引起上升或下降运动；增大某一只旋翼的转速，同时等量地减小同组另一只旋翼的转速，则产生俯仰、横滚运动；增大某一组旋翼的转速，同时等量减小另一组旋翼的转速，将产生偏航运动。

4.2建模分析

四旋翼飞行器受力分析,如图

所示

图2四旋翼飞行器受力分析示意图

旋翼机体所受外力和力矩为：

重力mg,机体受到重力沿方向；

四个旋翼旋转所产生的升力

(i=

1,2,3,4)，旋翼升力沿方向；

旋翼旋转会产生扭转力矩

(i=

1,2,3,4)。垂直于叶片的旋翼平面，与旋转矢量相反。

力模型为：，旋翼通过螺旋桨产生升力。是电机转动力系数，可取，为电机转速。旋翼旋转产生旋转力矩Mi(i=1,2,3,4)，力矩Mi的旋向依据右手定则确定。力矩模型为，其中是电机转动力系数，可取为电机转速。当给定期望转速后，电机的实际转速需要经过一段时间才能达到。实际转速与期望转速之间的关系为一阶延迟：响应延迟时间可取0.05s(即)。期望转速则需要限制在电机的最小转速和最大转速之间，范围可分取[1200rpm，7800rpm]。

飞行器受到外界力和力矩的作用，形成线运动和角运动。线运动由合外力引起，符合牛顿第二定律：

r为飞机的位置矢量。

角运动由合力矩引起。四旋翼飞行器所受力矩来源于两个方面：1）旋翼升力作用于质心产生的力矩；2）旋翼旋转产生的扭转力矩。角运动方程如下式所示。其中，L

为旋翼中心建立飞行器质心的距离，I

为惯量矩阵。

4.3控制回路设计

控制回路包括内外两层。外回路由Position

Control

模块实现。输入为位置误差，输出为期望的滚转、俯仰和偏航角。内回路由Attitude

Control

模块实现，输入为期望姿态角，输出为期望转速。Motor

Dynamics

模块模拟电机特性，输入为期望转速，输出为力和力矩。Rigid

Body

Dynamics

是被控对象，模拟四旋翼飞行器的运动特性。

图3包含内外两个控制回路的控制结构

（1）内回路：姿态控制回路

对四旋翼飞行器，我们唯一可用的控制手段就是四个旋翼的转速。因此，这里首先对转速产生的作用进行分析。假设我们希望旋翼1的转速达到，那么它的效果可分解成以下几个分量：

：使飞行器保持悬停的转速分量；

：除悬停所需之外，产生沿ZB轴的净力；

：使飞行器负向偏转的转速分量；

：使飞行器正向偏航的转速分量；

因此，可以将期望转速写成几个分量的线性组合：

其它几个旋翼也可进行类似分析，最终得到：

在悬浮状态下，四个旋翼共同的升力应抵消重力，因此：

此时，可以把旋翼角速度分成几个部分分别控制，通过“比例-微分”控制律建立如下公式：

综合以上三式可得到期望姿态角-期望转速之间的关系，即内回路。

外回路：位置控制回路

外回路采用以下控制方式：通过位置偏差计算控制信号（加速度）；建立控制信号与姿态角之间的几何关系；得到期望姿态角，作为内回路的输入。期望位置记为。可通过PID

控制器计算控制信号：

是目标悬停位置是我们的目标悬停位置(i=1,2,3)，是期望加速度，即控制信号。注意：悬停状态下线速度和加速度均为0，即。

通过俯仰角和滚转角控制飞行器在XW和YW平面上的运动，通过控制偏航角，通过控制飞行器在ZB轴上的运动。可得：

根据上式可按照以下原则进行线性化：

（1）将俯仰角、滚转角的变化作为小扰动分量，有，，（2）偏航角不变，有，其中初始偏航角，为期望偏航角（3）在悬停的稳态附近，有

根据以上原则线性化后，可得到控制信号（期望加速度）与期望姿态角之间的关系：

则内回路的输入为：

5实验步骤与结果

（1）

根据控制回路的结构建立simulink模型；

（2）

为了便于对控制回路进行参数调整，利用Matlab软件为四旋翼飞行器创建GUI参数界面；

（3）

利用Matlab的VR

Toolbox建立四旋翼飞行器的动画场景

（4）

根据系统的结构框图，搭建Simulink模块以实现模拟飞行器在指定位置的悬停。使用默认数据，此时xdes=3，ydes=4，zdes=5，开始仿真，可以得到运动轨迹x、y、z的响应函数，同时可以得到在xyz坐标中的空间运动轨迹。然后点击GUI中的VR按钮使simulink的工作空间中载入系统仿真所需的参数，把x、y、z的运动轨迹和Roll，Pitch，Yaw输入至VR中的模拟飞行器中，观察飞行器的运动轨迹和运动姿态，然后再使用一组新的参数xdes=-8，ydes=3，zdes=6进行四旋翼飞行器运动进行仿真模拟，可以看出仿真结果和动画场景相吻合。

6实验总结与心得

此次MATLAB实验综合了SIMULINK、GUI和VR场景等多个部分，对四旋翼飞行器运动进行了仿真模拟。由仿真结果可以看出，四旋翼飞行器最终位置达到了期望给定的位置，三个方向的响应曲线最终平稳，对应飞行器悬停在期望位置，达到了控制要求。

本次试验收获很多，学习到了很多知识，首先是熟悉了SIMULINK由简至繁搭建系统的过程，学习了利用VR建立虚拟模型，并在SIMULINK中连接。其次是熟悉了MATLAB

GUI界面的编写和搭建过程。Matlab提供了强大的用户图形界面，以帮助用户不必编写底层程序而直接在软件包基础上进行自行开发，这点在诸多软件中都有所体现。另外通过实验，对四旋翼飞行器的受力分析、模型建立、控制回路设计等有了较为细致的了解。

第五篇：信号与系统仿真实验报告

信号与系统

仿真

实 验 报 告

班级： 学号： 姓名： 学院：

实验一

一、实验者姓名：

二、实验时间：

三、实验地点：

四、实验题目：

5(s25s6)求三阶系统H(s)3的单位阶跃响应，并绘制响应波形图。

s6s210s8

五、解题分析：要知道求单位阶跃响应需知道所用函数，以及产生波形图所需要用到的函数。

六、试验程序：

num=[5 25 30];den=[1 6 10 8];step(num,den,10);title(‘Step response’)

七、实验结果：

实验所得波形图如下：

Step response4.543.53Amplitude2.521.510.50012345Time(sec)678910

八、实验心得体会：通过本次试验了解学会了一些新的函数的应用。了解到了N阶系统的单位阶跃响应的计算方法，和系统的响应波形图的函数应用和绘制方法。为后面的实验打下基础，并对信号仿真和《信号与系统》这门课程之间的联系有所增加，对《信号与系统》这门课里的问题也有了更加深入地了解。

九、实验改进想法：无。

实验二

一、实验者姓名：

二、实验时间：

三、实验地点：

四、实验题目:

一个因果线性移不变系统y(n)0.81y(n2)x(n)x(n2)，求：（1）H(z)；（2）冲激响应h(n)；（3）单位阶跃响应u(n)；（4）H(ej)，并绘出幅频和相频特性。

五、解题分析：离散卷积是数字信号处理中的一个基本运算，MTLAB提供的计算两个离散序列卷积的函数是conv，其调用方式为 y=conv(x,h)。其中调用参数x,h为卷积运算所需的两个序列，返回值y是卷积结果。

MATLAB函数conv的返回值y中只有卷积的结果，没有y的取值范围。由离散序列卷积的性质可知，当序列x和h的起始点都为k=0时，y的取值范围为k=0至length(x)+length(h)-2。

许多离散LTI都可用如下的线性常系数的差分方程描述

ay[kn]bx[kn]

nnn0n0NN其中x[k]、y[k]分别系统的输入和输出。在已知差分方程的N个初始状态y[k]，和输入x[k]，就可由下式迭代计算出系统的输出

y[k](an/a0)y[kn](bn/b0)x[kn]

n1n0NM利用MATLAB提供的filter函数，可方便地计算出上述差分方程的零状态响应。filter函数调用形式为 y=filter(b,a,x)。其中 a[a0,a1,...,aN]，b[b0,b1,...,bM]，分别表示差分方程系数。X表示输入序列，y表示输出序列。输出序列的长度和序列相同。

当序列的DTFT可写成ej的有理多项式时，可用MATLAB信号处理工具箱提供的freqz函数计算DTFT的抽样值。另外，可用MATLAB提供的abs、angle、real、imag等基本函数计算 DTFT的幅度、相位、实部、虚部。若X（ej）可表示为

b0b1ej...bMejMB(ej)X(e)jjjNA(e)a0a1e...aNe则freqz的调用形式为 X=freqz(b,a,w)，其中的b和 a分别是表示前一个

j式子中分子多项式和分母多项式系数的向量，即a[a0,a1,...,aN]，w为抽样的频率点，向量w的长度至少为2。返回值X就是DTFTb[b0,b1,...,bM]。在抽样点w上的值。注意一般情况下，函数freqz的返回值X是复数。

六、实验程序:

clc;clear;close;b=[1 0-1];a=[1 0-0.81];figure(1);subplot(2,1,1);dimpulse(b,a,20)subplot(2,1,2);dstep(b,a,50)w=[0:1:512]*pi/512;figure(2);freqz(b,a,w)

七、实验结果：

冲击响应图及阶跃响应图:

Impulse Response1Amplitude0.50-0.50246810Time(sec)Step Response12141618201Amplitude0.500510152025Time(sec)3035404550 100Magnitude(dB)0-100-200-30000.10.20.30.40.50.60.70.8Normalized Frequency( rad/sample)0.91100Phase(degrees)500-50-10000.10.20.30.40.50.60.70.8Normalized Frequency( rad/sample)0.91

八、实验心得体会：通过实验我们知道了使用Matlab来绘出出一个线性移不变系统的幅频和相频曲线。并知道了在《信号与系统》中得一些差分方程和各种响应，譬如零输入相应、零状态响应、全响应、自由响应、强迫响应、冲击响应、单位阶跃响应等等各种响应在Matlab中的函数表达方式和他们的求法，以及系统的幅频和相频曲线的绘制都有了一定深刻的认识。

九、实验改进想法：无。

实验三

一、实验者姓名：

二、实验时间：

三、实验地点：

四、实验题目：

模拟信号x(t)2sin(4t)5cos(8t)，求N64的DFT的幅值谱和相位谱。

五、解题分析：在MATLAB信号处理工具箱中，MATLAB提供了4个内部函数用于计算DFT和IDFT，它们分别是：fft(x)，fft(x,N)，ifft(X)，ifft(X,N)。

fft(x)计算M点的DFT。M是序列x的长度，即M=length(x)。

fft(x,N)计算N点的DFT。若M>N，则将原序列截短为N点序列，再计算其N点DFT；若M

ifft(X)计算M点的IDFT。M是序列X的长度。

ifft(X,N)计算N点IDFT。若M>N，则将原序列截短为N点序列，再计算其N点IDFT；若M

六、实验程序：

clc;clear;close;N=64;n=0:63;t=d*n;q=n*2*pi/N;x=2*sin(4*pi*t)+5*cos(8*pi*t);y=fft(x,N);subplot(3,1,1);plot(t,x);title(‘source signal’);subplot(3,1,2);plot(q,abs(y));title(‘magnitude’);subplot(3,1,3);plot(q,angle(y));title(‘phase’);

七、实验结果：

***0100806040200|F(k)|05101520Frequency253035

***0100806040200|F(k)|05101520Frequency253035 4321|jW|0-1-2-3-405101520Frequency253035Step Response400020000-2000 Amplitude-4000-6000-8000-10000-12000-1400001234n(samples)5678

八、实验心得体会：通过本次试验我知道了求取模拟信号在N等于一定值时的的DFT的幅值谱和相位谱的求法。通过本次实验，对幅值谱和相位谱有了更深的了解，并与课程《信号与系统》里的一些相关知识连接到了一起，使得学到的只是更加深刻、有意义。

九、实验改进想法：无。

实验四

一、实验者姓名：

二、实验时间：

三、实验地点：

四、实验题目：

将信号x(t)sin(240t)做离散序列，比较原序列与经过FFT和IFFT变换后的序列，并做出说明。

五、解题分析：此题需要对信号做离散序列，还要做FFT和IFFT变换，然后得到图像进行比较。连续时间函数与离散时间函数在编程中的区别主要体现在如下两个方面：第一，自变量的取值范围不同，离散时间函数的自变量是整数，而连续时间函数的自变量为一定范围内的实数；第二，绘图所用的函数不同，连续函数图形的绘制不止一个。本实验中要求绘制离散时间信号图，可以应用MATLAB中的函数来实现。用MATLAB表示一离散序列，可用两个向量来表示。其中一个向量表示自变量的取值范围，另一个向量表示序列的值。之后画出序列波形。当序列是从0开始时，可以只用一个向量x来表示序列。由于计算机内寸的限制，MATLAB无法表示一个无穷长的序列。对于典型的离散时间信号，可用逻辑表达式来实现不同自变量时的取值。

六、实验程序：

t=0:1/255:1;x=sin(2*pi*120*t);y=real(ifft(fft(x)));subplot(2,1,1);plot(t,x);title(‘原波形’);subplot(2,1,2);plot(t,y);

七、实验结果：

原波形10.50-0.5-100.10.20.30.40.50.60.70.80.91恢复的波形10.50-0.5-100.10.20.30.40.50.60.70.80.91

八、实验心得体会：通过对做信号的离散序列以及经FFT和IFFT的变换，了解了相关特性。通过计算机做出的信号波形图，我们能够很直白的看出原波形和经过变换后的波形的差别。

九、实验改进想法：无。

实验五

一、实验者姓名：

二、实验时间：

三、实验地点：

四、实验题目：

2s，激励信号22(s1)100x(t)(1cot)sco1s0(t)0，求（1）带通滤波器的频率响应；（2）输出稳态响应并绘制图形。已知带通滤波器的系统函数为H(s)

五、解题分析：需要知道求频率响应的方法，并绘制图形。

六、实验程序：

clear;t=linspace(0,2*pi,1001);w=[99,100,101];U=[0.5,1,0.5];b=[2,0];a=[1,2,10001];u1=U*cos(w’*t+angle(U’)*ones(1,1001));H=polyval(b,j*w)./polyval(a,j*w);H=freqs(b,a,w);subplot(2,1,1),plot(w,abs(H)),grid;subplot(2,1,2),plot(w,angle(H)),grid;u21=abs(U(1)*H(1))*cos(99*t+angle(U(1)*H(1)));u22=abs(U(2)*H(2))*cos(100*t+angle(U(2)*H(2)));u23=abs(U(3)*H(3))*cos(101*t+angle(U(3)*H(3)));u2=u21+u22+23;figure(2);subplot(2,1,1),plot(t,u1);subplot(2,1,2),plot(t,u2);

七、实验结果：

10.90.80.79910.50-0.5-19999.299.499.699.8100100.2100.4100.6100.810199.299.499.699.8100100.2100.4100.6100.8101

210-1-***222101234567

八、实验心得体会：通过本次试验，了解了频率响应求法，加深了对输出稳态响应的印象。

九、实验改进想法：无。