系统辨识与建模实验报告电加热炉动态特性辨识实验
姓名学号：
张春燕312102332
同组同学：沈剑312102331
序号：81
指导老师：郭毓
实验时间：2013年5月
系统辨识与建模实验电加热炉动态特性辨识实验
目录
一．实验目的 (3)
2.2 数据获取 (3)
2.3 离线辨识 (3)
3.1 数据预处理 (3)
3.2 结构辨识 (4)
3.3.1 RLS辨识参数 (6)
3.3.2 RELS辨识参数 (7)
3.3.3 RIV辨识参数 (8)
3.3.4 RML辨识参数 (9)
3.4 模型验证 (10)
3.4.1 输入阶跃响应比较 (10)
3.4.2 比较残差 (11)
四．实验结果分析 (12)
五．实验心得 (12)
附录1： (13)
附录2： (13)
附录3： (14)
附录4： (16)
一．实验目的
通过实验了解辨识方法在工程应用中的一些实际问题；了解数据获取和数据处理的各种方法和手段，掌握各种辨识方法的应用特点。
二．实验内容及其步骤
2.1 编写M 序列的产生程序
在实验参数设定时选择加热电压60V ，采样周期为3S ，所以加入的M 序列电压最好为加热电压的10-20%，M 序列的采样周期为数据采样周期的整数倍，因为实验时间有限，选择了2组数据，即M 序列信号为6V 、10V ，采样周期60s.
2.2 数据获取
高温老化试验温箱，以控制电压作为炉温控制系统的输入控制变量，即，设备的输入量是燃料供给量或电压、电流，而输出量是炉膛内腔的温度。
在热稳定工况的基础上，在电压稳定值上再附加一个辨识信号，即M 序列电压信号。加热炉热惯性大，升温过程较长，所以采样周期较长，M 序列的周期也较长。这样施加M 序列周期信号之后，记录几个周期的炉温试验数据。
2.3 离线辨识
利用处理过的数据，选择某种辨识方法；如RLS 、RELS 、RIV 或RML 等参数估计方法计算，以及F 检验方法或AIC 定阶法。离线估计出参数模型参数，并计算相应的模型静态增益，同时比较利用不同方法所得到的辨识结果。最后，模型验证。
三．数据处理
数据处理主要包括输入输出数据，模型结构确定，然后辨识参数。
3.1 数据预处理
在实验中采集了四组数据，仿真时选择M 序列为10的数据组。
其中采样/20s T ?=，为隔20点采样。在离线辨识时发现，如果数据全部运用则采样时间长度过短，则计算量大并容易产生病态方程，所以数据采用隔点取数据的方法。
程序见附录1。
图3.1 M 序列及其温度变化
3.2 结构辨识
残差校验估计方差最小定阶。利用最小二乘法求出θ值估值，计算残差。 不考虑相关误差
可以假设系统的传递函数如下：
10()(1)()()()na nb y k a y k a y k na b u k d b u k d nb ε
+-++-=-++--+L L
即 10()-(1)--()()()na nb y k a y k a y k na b u k d b u k d nb ε=--+-++--+L L
Y θε=Φ+
需要估计的参数：
1201[]T na nb a a a b b b θ=L
L 已知数据构成的向量：
()[(1)()()()]k y k y k na u k d u k nb d Φ=-------L L
根据最小二乘法
1?()T T Y θ-=ΦΦΦ
进行系统定阶时的误差平方和：
21(()())n N
k n J y k k θ+=+=
-Φ?∑
在F 检验法中 )
(22)()()()(12222121n n n N n J n J n J n n t --?-=
， 代入数据由matlab 计算得：
图3.2 残差校验定阶
有数据及图可以看出，当n=3时，AIC 最小，而且此时t<3，与工程实践结
实验报告 --实验1.基于matlab的4阶系统辨识实验 课程：系统辨识 题目：基于matlab的4阶系统辨识实验 作者： 专业：自动化 学号：11351014 目录 实验报告 (1) 1.引言 (2) 2.实验方法和步骤 (2) 3.实验数据和结果 (2) 4.实验分析 (4)
1、 引言 系统辨识是研究如何确定系统的数学模型及其参数的理论。而模型化是进行系统分析、仿真、设计、预测、控制和决策的前提和基础。 本次实验利用matlab 工具对一个简单的4阶系统进行辨识，以此熟悉系统辨识的基本步骤，和matlab 里的一些系统辨识常用工具箱和函数。 这次实验所采取的基本方法是对系统输入两个特定的激励信号，分别反映系统的动态特性和稳态特性。通过对输入和输出两个系统信号的比较，来验证系统的正确性。 2、 实验方法和步骤 2.1 实验方法 利用matlab 对一个系统进行辨识，选取的输入信号必须能够反映系统的动态和稳态两个方面的特性，才能更好地确定系统的参数。本次实验采取了两种输入信号，为反映动态特性，第一个选的是正弦扫频信号，由下面公式产生： 选定频率范围 ，w(t)是时间t 的线性函数，具有扫频性质，可以反映系统的动态特性。 为反映稳态特性，选的输入信号是阶跃信号。以上的到两组数据，利用matlab 的merge()函数，对两组数据融合，然后用matlab 系统辨识工具箱中的基于子空间方法的状态空间模型辨识函数n4sid()来对系统进行辨识 2.2 实验步骤 (1)建立一个4阶的线性系统，作为被辨识的系统，传递函数为 3243211548765 ()125410865 s s s G s s s s s -+-+=++++ (2)产生扫频信号u1和阶跃信号u2 (3)u1、u2作为输入对系统进行激励，分别产生输出y1和y2 (4)画出稳态测试输入信号u1-t 的曲线，和y1-t 的曲线 画出动态测试输入信号u2-t 的曲线，和y2-t 的曲线 (5)使用merge()函数对u1-y1数据和u2-y2数据进行融合，并使用n4sid()函数对系统进行辨识。 (6)画出原系统和辨识出的系统的零极点图，画出原系统和辨识出的系统的阶跃响应特性曲线，通过对比，验证辨识出的系统的准确性。 3、 实验数据和结果 （1） 分别以扫频正弦函数、阶跃函数作为系统的激励，得到的输出：
信号与系统仿真实验报告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)
在下面的题目中选做100分的题目，给出详略得当的答案。 一.通过举例简要说明数学建模的一般过程或步骤。(15分) 答：建立数学模型的方法大致有两种，一种是实验归纳的方法，即根据测试或计算数据，按照一定的数据，按照一定的数学方法，归纳出系统的数学模型；另一种是理论分析的方法，具体步骤有五步（以人口模型 为例）： 1、明确问题，提出合理简化的假设：首先要了解问题的实际背景，明确题目的要求，收集各种必要的信息 2、建立模型：据所做的假设以及事物之间的联系，构造各种量之间的关系。（查资料得出数学式子或算法）。 3、模型求解：利用数学方法来求解上一步所得到的数学问题，此时往往还要做出进一步的简化或假设。注意要尽量采用简单的数学公具。例如：马尔萨斯模型，洛杰斯蒂克模型 4、模型检验：根据预测与这些年来人口的调查得到的数目进行对比检验 5、模型的修正和最后应用：所建立的模型必须在实际应用中才能产生效益，根据预测模型，制定方针政策，以实现资源的合理利用和环境的保护。 二.把一张四条腿等长的正方形桌子放在稍微有些起伏的地面上，通常只有三只脚着地，然而 只需稍为转动一定角度，就可以使四只脚同时着地，即放稳了。(1) 请用数学模型来描述和证明这个实际问题; (2)讨论当桌子是长方形时，又该如何描述和证明？(15分) 答： 模型假设： 1．椅子四条腿一样长，椅脚与地面的接触部分相对椅子所占的地面面积可视为一个点。 2．地面凹突破面世连续变化的，沿任何方向都不会出现间断（没有向台阶那样的情况），即地面可看作数学上的连续曲面。 3．相对椅脚的间距和椅子腿的长度而言，地面是相对平坦的，即使椅子在任何位置至少有三条腿同时着地。4．椅子四脚连线所构成的四边形是圆内接四边形，即椅子四脚共圆。 5．挪动仅只是旋转。 我们将椅子这两对腿的交点作为坐标原点，建立坐标系，开始时AC、BD这两对腿都在坐标轴上。将AC和BD这两条腿逆时针旋转角度θ。记AC到地面的距离之和为f(θ)。记BD到 地面的距离之和为g(θ)。易得f(θ)，g(θ)至少有一个为零。
一、相关分析法 (1)实验原理 图1 实验原理图 本实验的原理图如图1。过程传递函数()G s 中12120,8.3, 6.2K T Sec T Sec ===；输入变量()u k ，输出变量()z k ，噪声服从2(0,)v N σ，0()g k 为过程的脉冲响应理论 值，?()g k 为过程脉冲响应估计值，()g k 为过程脉冲响应估计误差。 过程输入()u k 采用M 序列，其输出数据加白噪声()v k 得到输出数据()z k 。利 用相关分析法估计出过程的脉冲响应值?()g k ，并与过程脉冲响应理论值0()g k 比较，得到过程脉冲响应估计误差值()g k 。 M 序列阶次选择说明：首先粗略估计系统的过渡过程时间T S (通过简单阶跃响应)、截止频率f M (给系统施加不同周期的正弦信号或方波信号，观察输出)。本次为验证试验，已知系统模型，经计算Hz T T f M 14.01 2 1≈= ，s T S 30≈。根据式M f t 3 .0≤ ?及式S T t N ≥?-)1(，则t ?取值为1，此时31≥N ，由于t ?与N 选择时要求完全覆盖，则选择六阶M 移位寄存器，即N =63。
(2)编程说明 图2 程序流程图 (3)分步说明 ① 生成M 序列： M 序列的循环周期63126=-=N ，时钟节拍1t Sec ?=，幅度1a =，移位寄存器中第5、6位的内容按“模二相加”，反馈到第一位作为输入。其中初始数据设为{1,0,1,0,0,0}。程序如下：
② 生成白噪声序列： 程序如下： ③ 过程仿真得到输出数据： 如图2所示的过程传递函数串联，可以写成形如1212 11 ()1/1/K G s TT s T s T = ++， 其中112 K K TT = 。 图2 过程仿真方框图 程序如下： ④ 计算脉冲响应估计值：
中南大学 计算机仿真与建模 实验报告 题目：理发店的服务过程仿真 姓名：XXXX 班级：计科XXXX班 学号：0909XXXX 日期：2013XXXX
理发店的服务过程仿真 1 实验案例 (2) 1.1 案例：理发店系统研究 (2) 1.1.1 问题分析 (3) 1.1.2 模型假设 (3) 1.1.3 变量说明 (3) 1.1.4 模型建立 (3) 1.1.5 系统模拟 (4) 1.1.6 计算机模拟算法设计 (5) 1.1.7 计算机模拟程序 (6) 1实验案例 1.1 案例：理发店模拟 一个理发店有两位服务员A和B顾客随机地到达该理发店,每分钟有一个顾客到达和没有顾客到达的概率均是1/2 , 其中60%的顾客理发仅用5分钟,另外40%的顾客用8分钟. 试对前10分钟的情况进行仿真。 （“排队论”，“系统模拟”，“离散系统模拟”，“事件调度法”）
1.1.1 问题分析 理发店系统包含诸多随机因素，为了对其进行评判就是要研究其运行效率， 从理发店自身利益来说，要看服务员工作负荷是否合理，是否需要增加员工等考 虑。从顾客角度讲，还要看顾客的等待时间，顾客的等待队长，如等待时间过长 或者等待的人过多，则顾客会离开。理发店系统是一个典型的排队系统，可以用 排队论有关知识来研究。 1.1.2 模型假设 1． 60％的顾客只需剪发，40％的顾客既要剪发，又要洗发； 2． 每个服务员剪发需要的时间均为5分钟，既剪发又洗发则花8分钟； 3． 顾客的到达间隔时间服从指数分布； 4． 服务中服务员不休息。 1.1.3 变量说明 u ：剪发时间（单位：分钟），u=5m ； v: 既剪发又理发花的时间（单位：分钟），v=8m ； T ： 顾客到达的间隔时间，是随机变量，服从参数为λ的指数分布，（单位： 分钟） T 0：顾客到达的平均间隔时间（单位：秒），T 0＝λ 1； 1.1.4 模型建立 由于该系统包含诸多随机因素，很难给出解析的结果，因此可以借助计算机 模拟对该系统进行模拟。 考虑一般理发店的工作模式，一般是上午9：00开始营业，晚上10：00左 右结束，且一般是连续工作的，因此一般营业时间为13小时左右。 这里以每天运行12小时为例，进行模拟。 这里假定顾客到达的平均间隔时间T 0服从均值3分钟的指数分布， 则有 3小时到达人数约为603 603=?人， 6小时到达人数约为1203 606=?人， 10小时到达人数约为2003 6010=?人， 这里模拟顾客到达数为60人的情况。 （如何选择模拟的总人数或模拟总时间）
物流系统建模与仿真课程实验报告 实验名称：物流系统建模与仿真Flexsim实验 学院：吉林大学机械与航空航天工程学院 专业班： 141803 姓名：龙振坤 学号： 14180325 2019年5月19日
一、实验目的 用flexsim模拟仓库分拣系统。 二、仿真实验内容（简要阐述仿真模型） 将五种不同货物通过分拣传送带分拣到五条传送带上，再由叉车将这五种货物分别运送到不同的货架之上。 三、仿真模型建模步骤 1、打开软件flexsim，并新建文件。 2、拉出所需要的离散实体： 发生器、暂存区、分拣传送带、传送带（5个）、叉车（3个）、货架（5个）。（如图） 3、设置分拣传送带、传送带、货架参数，并调整位置 分拣传送带布局：第一段平直，长度为5；第二段弯曲，角度为90°，半径为5；第三段平直，长度为20。传送带布局：长度为10。 货架布局：10层10列。
4、连接各个离散实体 将发生器与暂存区用“A”连接；暂存区与分拣传送带用“A”连接； 分拣传送带与传送带1、2、3、4、5分别用“A”连接； 传送带1、2、3、4、5与货架1、2、3、4、5分别用“A”连接； 传送带1、2与叉车1用“S”连接；传送带3、4与叉车2用“S”连接；传送带5与叉车3用“S”连接 5、设置各个离散实体的参数 发生器： 分拣传送带： 传送带：在临时实体流处勾选使用运输工具
6、运行文件 运行结果 四、课程体会及建议 课程体会： 作为flexsim软件的初学者，一开始在完成各种实例，熟悉各种操作的命令时确实遇到了不少的问题，但由于老师的耐心解答、同学的帮助、以及自己通过网络所寻求到的帮助，最终能够逐布掌握flexsim的一些基本使用方法。以目前的眼光看来，flexsim是一个功能非常强大的管理类模拟软件，这是我作为一名机械专业的学生在今后的学习中很少有机会能够接触到的。对于我来说，物流系统建模与仿真这门课不仅让我了解并掌握了一种从新的软件、一种没有见过的工具，更重要的是他对于我的一种工程思想的培养。在使用flexsim的过程中，深感整体性思想的重要性，操作过程中，每一个功能的实现都离不开各个离散实体的配合，选择何种实体型，使用何种函数命令，构成怎样的连接，这些都是功能可以最终实现的关键。 课程建议： ①没有使用麦克，声音过小，中后排听课效率低；②投影设备老化，颜色浅，清晰度低，部分操作难以看清，尤其是在输入一些代码的时候；③版本存在差异，属性界面略有区别，在一开始学习的时候很难跟上老师的脚步，强烈建议以后将该课程改为在机房上课。
南京理工大学 电加热炉动态特性辨识实验报告 作者: 张志鹏（94）学号：813001010014 实验时间2013年11月24日 组员: 刘心刚（63）李昊（88）倪镭（90） 任课老师：郭毓教授 2013 年 11 月
1.熟悉对实际控制系统的辨识与参数估计，并利用所得模型进行控制仿真，进而控制实际系统。 2.掌握实际工程中常用的辨识方法，如LS，RLS，RLES等。 二、实验平台: 嵌入式温度控制系统主要由嵌入式温度控制器、立式RGL-9076A 型温箱、NETGEAR 无线路由器和24V 开关电源等组成。系统电气连接如图1 所示。系 统采用CS（客户端—服务器）模式实现了一对一的服务器、客户端的数据通信。 嵌入式控制系统软软硬件运行平台. 硬件：PC 机、嵌入式温度控制器、NETGEAR 无线路由器等。 软件：Windows XP、Microsoft Visual C++ 6.0、Matlab 2007a 等。 图1 实验硬件平台
1.设置硬件。根据实验手册上的连接方式，确认硬件连接是否正确。根据使用手册进行IP设置、系统参数设置，直至软件可以实时显示温度曲线。 2.达到稳态。我们首先采用81V的加热电压加热使系统尽快到达某一较稳定温度。使用3S的采样周期进行采样温度信号。当温箱实际温度达到135度左右时，温度变化曲线几乎持平，我们认定此时温箱系统处于稳态。 3.加入辨识信号。这里选选取M序列进行辨识，在试验阶段我们组做了一组数据：选取M序列幅值为+20，-20，，辨识信号的采样周期为40s。加入辨识信号后继续进行数据采集。 4.数据处理、辨识系统模型。 5.分析辨识结果得出结论。 四、辨识算法及过程 经过分析研究，确定使用计算残差平方和的RELS方法验证模型的阶次及延时并辨识系统模型参数。 1、确定系统的延迟d
重庆大学 学生实验报告 实验课程名称物流系统建模与仿真 开课实验室物流工程实验室 学院自动化年级12 专业班物流工程2班学生姓名段竞男学号20124912 开课时间2014 至2015 学年第二学期 自动化学院制
《物流系统建模与仿真》实验报告
(2)属性窗口（Properties Window） 右键单击对象，在弹出菜单中选择 Properties；用于编辑和查看所有对象都拥有的一般性信息。 (3)模型树视图（Model Tree View） 模型中的所有对象都在层级式树结构中列出；包含对象的底层数据结构；所有的信息都包含在此树结构中。 4)重置运行 (1)重置模型并运行 (2)控制仿真速度（不会影响仿真结果） (3)设置仿真结束时间 5)观察结果 (1)使用“Statistics”（统计）菜单中的Reports and Statistics（报告和统计）生成所需的 各项数据统计报告。 (2)其他报告功能包括：对象属性窗口的统计项；记录器对象；可视化工具对象；通过触发器 记录数据到全局表。
五、实验过程原始记录(数据、图表、计算等) 1、运行结果的平面视图： 2、运行结果的立体视图 3、运行结果的暂存区数据分析结果图：
第一个暂存区 第二个暂存区 由报表分析可知5次实验中，第一个暂存区的平均等待时间为11.46，而第二个暂存区的平均等待时间为13.02，略大于第一个暂存区，由此可见，第二个暂存区的工作效率基本上由第一个暂存区决定。 4、运行结果三个检测台的数据分析结果图，三个检测台的state饼图： （1）处理器一：
数学建模实验报告
一、实验目的 1、通过具体的题目实例，使学生理解数学建模的基本思想和方法，掌握 数学建模分析和解决的基本过程。 2、培养学生主动探索、努力进取的的学风，增强学生的应用意识和创新 能力，为今后从事科研工作打下初步的基础。 二、实验题目 （一）题目一 1、题目：电梯问题有r个人在一楼进入电梯，楼上有n层。设每个 乘客在任何一层楼出电梯的概率相同，试建立一个概率模型，求直 到电梯中的乘客下完时，电梯需停次数的数学期望。 2、问题分析 （1）由于每位乘客在任何一层楼出电梯的概率相同，且各种可能的情况众多且复杂，难于推导。所以选择采用计算机模拟的 方法，求得近似结果。 （2）通过增加试验次数，使近似解越来越接近真实情况。 3、模型建立 建立一个n*r的二维随机矩阵，该矩阵每列元素中只有一个为1，其余都为0，这代表每个乘客在对应的楼层下电梯（因为每 个乘客只会在某一层下，故没列只有一个1）。而每行中1的个数 代表在该楼层下的乘客的人数。 再建立一个有n个元素的一位数组，数组中只有0和1,其中1代表该层有人下，0代表该层没人下。 例如： 给定n=8;r=6（楼8层，乘了6个人）,则建立的二维随机矩阵及与之相关的应建立的一维数组为： m = 0 0 1 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 0 c = 1 1 0 1 0 1 1 1 4、解决方法（MATLAB程序代码）：
n=10;r=10;d=1000; a=0; for l=1:d m=full(sparse(randint(1,r,[1,n]),1:r,1,n,r)); c=zeros(n,1); for i=1:n for j=1:r if m(i,j)==1 c(j)=1; break; end continue; end end s=0; for x=1:n if c(x)==1 s=s+1; end continue; end a=a+s; end a/d 5、实验结果 ans = 6.5150 那么，当楼高11层，乘坐10人时，电梯需停次数的数学期望为6.5150。 （二）题目二 1、问题：某厂生产甲乙两种口味的饮料,每百箱甲饮料需用原料6 千克,工人10名,可获利10万元;每百箱乙饮料需用原料5千 克,工人20名,可获利9万元.今工厂共有原料60千克,工人 150名,又由于其他条件所限甲饮料产量不超过8百箱.问如何 安排生产计划,即两种饮料各生产多少使获利最大.进一步讨 论: 1)若投资0.8万元可增加原料1千克,问应否作这项投资. 2)若每百箱甲饮料获利可增加1万元,问应否改变生产计划. 2、问题分析 （1）题目中共有3个约束条件，分别来自原料量、工人数与甲饮料产量的限制。 （2）目标函数是求获利最大时的生产分配，应用MATLAB时要转换
广东外语外贸大学 物流系统仿真实验 通达企业立体仓库实验报告 指导教师：翟晓燕教授专业：物流管理1101 姓名：李春立 20110402088 吴可为 201104020117 陈诗涵 201104020119 丘汇峰 201104020115
目录 一、企业简介 (2) 二、通达企业立体仓库模型仿真 (2) 1................................ 模型描述：2 2................................ 模型数据：3 3.............................. 模型实体设计4 4.................................. 概念模型4 三、仿真模型内容——Flexsim模型 (4) 1.................................. 建模步骤4 2.............................. 定义对象参数5 四、模型运行状态及结果分析 (7) 1.................................. 模型运行7 2................................ 结果分析：7 五、报告收获 (9) 一、企业简介 二、通达企业立体仓库模型仿真 1. 模型描述： 仓储的整个模型分为入库和出库两部分，按作业性质将整个模型划分为暂存区、分拣区、
储存区以及发货区。 入库部分的操作流程是： ①.（1）四种产品A，B，C，D首先到达暂存区，然后被运输到分类输 送机上，根据设定的分拣系统将A，B，C，D分拣到1，2,3,4，端口； ②.在1,2,3,4，端口都有各自的分拣道到达处理器，处理器检验合格 的产品被放在暂存区，不合格的产品则直接吸收掉；每个操作工则将暂存 区的那些合格产品搬运到货架上；其中，A，C产品将被送到同一货架上， 而B，D则被送往另一货架； ③.再由两辆叉车从这两个货架上将A/B，C/D运输到两个暂存区上； 此时，在另一传送带上送来包装材料，当产品和包装材料都到达时，就可 以在合成器上进行对产品进行包装。 出库部分的操作流程是：包装完成后的产品将等待被发货。 2. 模型数据： ①.四种货物A，B，C，D各自独立到达高层的传送带入口端： A: normal(400,50) B: normal(400,50) C: uniform(500,100) D: uniform(500,100) ②.四种不同的货物沿一条传送带，根据品种的不同由分拣装置将其推 入到四个不同的分拣道口，经各自的分拣道到达操作台。 ③.每检验一件货物占用时间为60,20s。 ④.每种货物都可能有不合格产品。检验合格的产品放入检验器旁的暂 存区；不合格的吸收器直接吸收；A的合格率为95%，B为96%，C的合格 率为97%，D的合格率为98%。 ⑤.每个检验操作台需操作工一名，货物经检验合格后，将货物送至货 架。 ⑥.传送带叉车的传送速度采用默认速度（包装物生成时间为返回60 的常值），储存货物的容器容积各为1000单位，暂存区17,18，21容量为 10；
生产系统建模与及仿真 实验报告 实验一Witness仿真软件认识 一、实验目的 1、学习、掌握Witness仿真软件的主要功能与使用方法； 2、学习生产系统的建模与仿真方法。 二、实验内容 学习、掌握Witness仿真软件的主要功能与使用方法 三、实验报告要求 1、写出实验目的: 2、写出简要实验步骤； 四、主要仪器、设备 1、计算机（满足Witness仿真软件的配置要求） 2、Witness工业物流仿真软件。 五、实验计划与安排 计划学时4学时 六、实验方法及步骤 实验目的： 1、对Witness的简单操作进行了解、熟悉，能够做到基本的操作，并能够进行简单的基础建模。 2、进一步了解Witness的建模与仿真过程。 实验步骤： Witness仿真软件是由英国lanner公司推出的功能强大的仿真软件系统。它可以用于离散事件系统的仿真，同时又可以用于连续流体（如液压、化工、水力）系统的仿真。目前已成功运用于国际数千家知名企业的解决方案项目，有机场设施布局
优化、机场物流规划、电气公司的流程改善、化学公司的供应链物流系统规划、工厂布局优化和分销物流系统规划等。 ◆Witness的安装与启动： ?安装环境：推荐P4 1.5G以上、内存512MB及以上、独立显卡64M以上显存，Windows98、Windows2000、Windows NT以及Windows XP的操作系统支持。 ?安装步骤：⑴将Witness2004系统光盘放入CD-ROM中，启动安装程序； ⑵选择语言（English）；⑶选择Manufacturing或Service；⑷选择授权方式（如加密狗方式）。 ?启动：按一般程序启动方式就可启动Witness2004，启动过程中需要输入许可证号。 ◆Witness2004的用户界面： ?系统主界面：正常启动Witness系统后，进入的主界面如下图所示： 主界面中的标题栏、菜单栏、工具栏状态栏等的基本操作与一般可视化界面操作大体上一致。这里重点提示元素选择窗口、用户元素窗口以及系统布局区。 ?元素列表窗口：共有五项内容，分类显示模型中已经建立和可以定义的模型元素。Simulation中显示当前建立的模型中的所有元素列表；Designer中显示当前Designer Elements中的所有元素列表；System中显示系默认的特殊地点；Type中
学生学号实验课成绩 学生实验报告书 实验课程名称数据分析与建模 开课学院 指导教师姓名 学生姓名 学生专业班级 2015 —2016 学年第 1 学期
实验报告填写规范 1、实验是培养学生动手能力、分析解决问题能力的重要环节；实验报告是反映实验教学水 平与质量的重要依据。为加强实验过程管理，改革实验成绩考核方法，改善实验教学效果，提高学生质量，特制定本实验报告书写规范。 2、本规范适用于管理学院实验课程。 3、每门实验课程一般会包括许多实验项目，除非常简单的验证演示性实验项目可以不写实 验报告外，其他实验项目均应按本格式完成实验报告。在课程全部实验项目完成后，应按学生姓名将各实验项目实验报告装订成册，构成该实验课程总报告，并给出实验课程成绩。 4、学生必须依据实验指导书或老师的指导，提前预习实验目的、实验基本原理及方法，了 解实验内容及方法，在完成以上实验预习的前提下进行实验。教师将在实验过程中抽查学生预习情况。 5、学生应在做完实验后三天内完成实验报告，交指导教师评阅。 6、教师应及时评阅学生的实验报告并给出各实验项目成绩，同时要认真完整保存实验报 告。在完成所有实验项目后，教师应将批改好的各项目实验报告汇总、装订，交课程承担单位（实验中心或实验室）保管存档。
画出图形 由图x=4时，y最大等于1760000 (2)求关于所做的15%假设的灵敏性 粗分析： 假设C=1000 即给定r y=f(x)=(1500-100x)1000(1+rx)=-100000rx^2+1500000rx-100000x+1500000 求导，f’(x)=-200000rx+1500000r-100000,令f’(x)=0,可得相应x值，x=(15r-1)/2r Excel画出相应图形
系统辨识实验报告
实验一 最小二乘法 1 最小二乘算法 1.1 基本原理 系统模型 )()()()()(11k n k u z B k z z A +=-- a a n n z a z a z a z A ----++++= 221111)( b b n n z b z b z b z B ----+++= 22111)( 最小二乘格式 )()()(k n k h k z T +=θ [][] ?????=------=T n n T b a b a b b a a n k u k u n k z k z k h 11)()1()()1()(θ 对于L k ,,2,1 =，构成线性方程组 L L L n H z +=θ 式中， []T L L z z z z )()2()1( = []T L L n n n n )()2()1( = ? ????? ???? ??--------------= ??????????????=)()1()()1()2()1()2()1()1() 0() 1()0()()2()1(b a b a b a T T T L n L u L u n L z L z n u u n z z n u u n z z L h h h H 参数估计值为 ()L T L L T L LS z H H H 1 ?-=θ 1.2 Matlab 编程 % 基本最小二乘法LS clear;clc A=ones(5,1);B=ones(4,1);%A 为首1多项式，B 中体现时滞（d=1） na=length(A)-1;nb=length(B); load dryer2
实验1 Witness仿真软件认识 一、实验目的 熟悉Witness 的启动；熟悉Witness2006用户界面；熟悉Witness 建模元素；熟悉Witness 建模与仿真过程。 二、实验内容 1、运行witness软件，了解软件界面及组成； 2、以一个简单流水线实例进行操作。小部件（widget）要经过称重、冲洗、加工和检测等操作。执行完每一步操作后小部件通过充当运输工具和缓存器的传送带（conveyer）传送至下一个操作单元。小部件在经过最后一道工序“检测”以后，脱离本模型系统。 三、实验步骤 仿真实例操作: 模型元素说明：widget 为加工的小部件名称；weigh、wash、produce、inspect 为四种加工机器，每种机器只有一台；C1、C2、C3 为三条输送链；ship 是系统提供的特殊区域，表示本仿真系统之外的某个地方； 操作步骤： 1：将所需元素布置在界面：
2：更改各元素名称： 如; 3:编辑各个元素的输入输出规则：
4: 运行一周（5 天*8 小时*60 分钟=2400 分钟），得到统计结果。5:仿真结果及分析： Widget： 各机器工作状态统计表：
分析：第一台机器效率最高位100%，第二台机器效率次之为79%，第三台和第四台机器效率低下，且空闲时间较多，可考虑加快传送带C2、C3的传送速度以及提高第二台机器的工作效率，以此来提高第三台和第四台机器的工作效率。 6：实验小结： 通过本次实验，我对Witness的操作界面及基本操作有了一个初步的掌握，同学会了对于一个简单的流水线生产线进行建模仿真，总体而言，实验非常成功。
数学建模实验报告 实验一计算课本251页A矩阵的最大特征根和最大特征向量 1 实验目的 通过Wolfram Mathematica软件计算下列A矩阵的最大特征根和最大特征向量。 2 实验过程 本实验运用了Wolfram Mathematica软件计算，计算的代码如下：
3 实验结果分析 从代码的运行结果，可以得到最大特征根为5.07293，最大特征向量为 {{0.262281},{0.474395},{0.0544921},{0.0985336},{0.110298}},实验结果 与标准答案符合。
实验二求解食饵-捕食者模型方程的数值解 1实验目的 通过Wolfram Mathematica或MATLAB软件求解下列习题。 一个生物系统中有食饵和捕食者两种种群，设食饵的数量为x（t），捕食者为y（t），它们满足的方程组为x’(t)=(r-ay)x,y’(t)=-(d-bx)y,称该系统为食饵-捕食者模型。当r=1，d=0.5，a=0.1，b=0.02时，求满足初始条件x（0）=25，y（0）=2的方程的数值解。 2 实验过程 实验的代码如下 Wolfram Mathematica源代码： Clear[x,y] sol=NDSolve[{x'[t] (1-0.1y[t])x[t],y'[t] 0.02x[t]y[t]-0.5y[t],x[0 ] 25,y[0] 2},{x[t],y[t]},{t,0,100}] x[t_]=x[t]/.sol y[t_]=y[t]/.sol g1=Plot[x[t],{t,0,20},PlotStyle->RGBColor[1,0,0],PlotRange->{0,11 0}] g2=Plot[y[t],{t,0,20},PlotStyle->RGBColor[0,1,0],PlotRange->{0,40 }] g3=Plot[{x[t],y[t]},{t,0,20},PlotStyle→{RGBColor[1,0,0],RGBColor[ 0,1,0]},PlotRange->{0,110}] matlab源代码 function [ t,x ]=f ts=0:0.1:15; x0=[25,2]; [t,x]=ode45('shier',ts,x0); End function xdot=shier(t,x)
《物流系统仿真》 实验报告书 实验报告题目: 物流系统仿真学院名称: 专业: 班级: 姓名: 学号: 成绩: ２０15年5月
实验报告 一、实验名称 物流系统仿真 二、实验要求 ⑴根据模型描述与模型数据对配送中心进行建模; ⑵分析仿真实验结果，找出配送中心运作瓶颈,提出改进措施。 三、实验目得 1、掌握仿真软件Ｆleｘsｉｍ得操作与应用，熟悉通过软件进行物流仿真建模。 2、记录Ｆleｘsim软件仿真模拟得过程,得出仿真得结果。 3、总结Flexsim仿真软件学习过程中得感受与收获。 三、实验设备 PC机，Windｏwｓ XP，Flｅxsｉｍ教学版 四、实验步骤 1 货物得入库检验过程模型描述 三种货物以特定得批量在特定得时间送达仓库得暂存区，由两名操作员将它们搬运到相对应得检验台上去,检验时需要操作员对检验设备进行预置,并在完成检验时自动贴上相应得标签。货物经过检验后，通过不同得三个传输带传送到同一个位置。 构建模型布局 为验证Flexsｉm软件已被正确安装,双击桌面上得Fｌeｘsim图标打开应用程序。一旦软件安装好您应该瞧到Flexｓim菜单与工具条、实体库，与正投影模型视窗.
第1步:在模型中生成所需实体 从左边得实体库中拖动一个发生器到模型(建模）视窗中。具体操作就是,点击并按住实体库中得实体,然后将它拖动到模型中想要放置得位置,放开鼠标键。这将在模型中建立一个发生器实体,把其余实体按照同样得方法生成。如下图所示。一旦创建了实体,将会给它赋一个默认得名称,在以后定义得编辑过程中,可以对模型中得实体进行重新命名。 完成后,将瞧到上面这样得一个模型.模型中有1个发生器、1个暂存区、3个处理器、3个输送机、1个分配器、2名操作员与１个吸收器。 第2步:定义物流流程 (１）连接端口
实验报告 课程名称：现代控制工程与理论实验课题：最优控制 学号：12014001070 姓名：陈龙 授课老师：施心陵
最优控制 一、最优控制理论中心问题： 给定一个控制系统（已建立的被控对象的数学模型），选择一个容许的控制律，使被控对象按预定要求运行，并使给定的某一性能指标达到极小值（或极大值） 二、最优控制动态规划法 对离散型控制系统更为有效，而且得出的是综合控制函数。这种方法来源于多决策过程，并由贝尔曼首先提出，故称贝尔曼动态规划。 最优性原理：在一个多级决策问题中的最优决策具有这样的性质，不管初始级、初始状态和初始决策是什么，当把其中任何一级和状态做为初始级和初始状态时，余下的决策对此仍是最优决策 三、线性二次型性能指标的最优控制 用最大值原理求最优控制，求出的最优控制通常是时间的函数，这样的控制为开环控制当用开环控制时，在控制过程中不允许有任何干扰，这样才能使系统以最优状态运行。在实际问题中，干扰不可能没有，因此工程上总希望应用闭环控制，即控制函数表示成时间和状态的函数。 求解这样的问题一般来说是很困难的。但对一类线性的且指标是
二次型的动态系统，却得了完全的解决。不但理论比较完善，数学处理简单，而且在工际中又容易实现，因而在工程中有着广泛的应用。 一．实验目的 1.熟悉Matlab的仿真及运行环境； 2.掌握系统最优控制的设计方法； 3.验证最优控制的效果。 二．实验原理 对于一个给定的系统，实现系统的稳定有很多途径，所以我们需要一个评价的指标，使系统在该指标下达到最优。如果给定指标为线性二次型，那么我们就可以利用MATLAB快速的计算卡尔曼增益。 三．实验器材 PC机一台，Matlab仿真平台。 四．实验步骤 例题1 （P269）考虑液压激振系统简化后的传递函数方框图如下，其中K a为系统前馈增益，K f为系统反馈增益，w h为阻尼固有频率。（如图5-5所示） 将系统传递函数变为状态方程的形式如下： ,
《MATLAB与控制系统仿真》 实验报告 班级： 学号： 姓名： 时间：2013 年 6 月
目录实验一 MATLAB环境的熟悉与基本运算（一）实验二 MATLAB环境的熟悉与基本运算（二）实验三 MATLAB语言的程序设计 实验四 MATLAB的图形绘制 实验五基于SIMULINK的系统仿真 实验六控制系统的频域与时域分析 实验七控制系统PID校正器设计法 实验八线性方程组求解及函数求极值
实验一 MATLAB环境的熟悉与基本运算（一） 一、实验目的 1．熟悉MATLAB开发环境 2．掌握矩阵、变量、表达式的各种基本运算 二、实验基本原理 1.熟悉MATLAB环境: MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器、文件和搜索路径浏览器。 2.掌握MA TLAB常用命令 表1 MA TLAB常用命令 3.MATLAB变量与运算符 3．1变量命名规则 3．2 MATLAB的各种常用运算符 表3 MATLAB关系运算符 表4 MATLAB逻辑运算符
| Or 逻辑或 ~ Not 逻辑非 Xor 逻辑异或 符号功能说明示例符号功能说明示例 ：1:1:4;1:2:11 . ；分隔行.. ，分隔列… （）% 注释 [] 构成向量、矩阵！调用操作系统命令 {} 构成单元数组= 用于赋值 4.MATLAB的一维、二维数组的寻访 表6 子数组访问与赋值常用的相关指令格式 三、主要仪器设备及耗材 计算机 四．实验程序及结果 1、新建一个文件夹（自己的名字命名，在机器的最后一个盘符） 2、启动MATLAB，将该文件夹添加到MATLAB路径管理器中。 3、学习使用help命令。
专业仿真课程设计题目： 学院： 专业班级： 学号： 学生姓名： 指导教师： 设计时间：
专业仿真课程设计题目 主要研究内容: 从所拍摄的多个目标物中检测三角形物，给出三角形物几何中心、三个边长以及边长的方向、面积。 设计要求： （1）提交能够实现题目要求、并通过演示验收的可执行文件。 （2）提交课程设计报告（包括程序清单）。 （3）通过答辩，答辩成绩满分20分，其中个人设计部分10分，非个人设计部分10分。 （4）软件设计要求：有一个人机交互界面，模块化设计，在模块之间通过BMP文件或者文本文件传送数据，可以查看中间结果。 （5）5个人一组，组长协调分工，每个组员一定要有具体任务，以便考核。预期达到的目标： 1、能够通过相关文献查阅、文献综述和总结，给出问题求解的多种可行方案。 2、能够综合运用测控技术与仪器专业理论和技术手段，设计实验方案、分析实验结果，得出有效的结论。 3、能够借助MATLAB仿真软件，进一步掌握高等数学、复变函数与积分变换等相关数学和自然科学知识以及测控技术与仪器专业的基本理论知识，能够结合本专业“自动控制原理”、“数字信号处理”、“误差理论”等相关课程，采用MATLAB软件对复杂工程问题建立模型并进行预测与模拟； 4、能够与团队中其他学科成员合作开展工作，能够与其他队员很好地沟通和交流意见，能够通过口头或书面方式表达自己的设计思路，具有一定的表达能力和人际交往能力。
目录 第一章课程设计相关知识综述 1.1 MATLAB相关知识叙述 1.1.1 MATLAB基本知识介绍 1.1.2 MATLAB的优势特点 1.1.3 MATLAB的发展历程 1.2 MATLAB工具箱与函数 1.2.1 MATLAB图像处理工具箱 1.2.2 课程设计所用图像处理函数介绍第二章课程设计内容和要求 2.1 课程设计主要研究内容 2.2 课程设计要求 2.3 课程设计预期目标 第三章设计过程 3.1 设计方案 3.2 设计步骤及流程图 3.3 程序清单及相关注释 3.4 实验结果分析 3.5 结论 第四章团队情况 第五章总结 第六章参考文献
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”，参数设置如下图：