第1章建模与仿真的基本概念
参照P8例子,列举一个你相对熟悉的简单实际系统为例,采用非形式描述出来。
第2章建模方法论
1、什么是数学建模形式化的表示?试列举一例说明形式化表示与非形式化表示的区别。
模型的非形式描述是说明实际系统的本质,但不是详尽描述。是对模型进行深入研究的基础。主要由模型的实体、包括参变量的描述变量、实体间的相互关系及有必要阐述的假设组成。模型的非形式描述主要说明实体、描述变量、实体间的相互关系及假设等。
例子:环形罗宾服务模型的非形式描述:
实体
CPU,USR1,…,USR5
描述变量
CPU:Who,Now(现在是谁)----范围{1,2,…,5}; Who.Now=i表示USRi由CPU服务。
USR:Completion.State(完成情况)----范围[0,1];它表示USR完成整个程序任务的比例。参变量
X-----范围[0,1];它表示USRi每次完成程序的比率。
i
实体相互关系
(1)CPU 以固定速度依次为用户服务,即Who.Now为1,2,3,4,5,1,2…..循环运行。
X工作。假设:CPU对USR的服务时间固定,不(2)当Who.Now=I,CPU完成USRi余下的
i
X决定。
依赖于USR的程序;USRi的进程是由各自的参变量
i
2、何谓“黑盒”“白盒”“灰盒”系统?
“黑盒”系统是指系统内部结构和特性不清楚的系统。对于“黑盒”系统,如果允许直接进行实验测量并通过实验对假设模型加以验证和修正。对属于黑盒但又不允许直接实验观测的系统,则采用数据收集和统计归纳的方法来假设模型。
对于内部结构和特性清楚的系统,即白盒系统,可以利用已知的一些基本定律,经过分析和演绎导出系统模型。
3、模型有效性和模型可信性相同吗?有何不同?
模型的有效性可用实际系统数据和模型产生的数据之间的符合程度来度量。它分三个不同级别的模型有效:复制有效、预测有效和结构有效。不同级别的模型有效,存在不同的行为水平、状态结构水平和分解结构水平的系统描述。
模型的可信度指模型的真实程度。一个模型的可信度可分为:
在行为水平上的可信性,即模型是否重现真实系统的行为。
在状态结构水平上可信性,即模型能否与真实系统在状态上互相对应,通过这样的模型可以对未来的行为进行唯一的预测。
在分解结构水平上的可信性,即模型能否表示出真实系统内部的工作情况,而且是惟一表示出来。
不论对于哪一个可信性水平,可信性的考虑贯穿在整个建模阶段及以后各阶段,必须考虑以下几个方面:
1在演绎中的可信性。2在归纳中的可信性。3在目的方面的可信性。
4、基于计算机建模方法论与一般建模方法论有何不同?(P32)
经典的建模与仿真的主要研究思路,首先界定研究对象-实际系统的边界和建模目标,利用已有的数学建模工具和成果,建立相应的数学模型,并用计算装置进行仿真。这种经典的建
模与仿真方法,对研究和解决物理系统的有关问题是非常有效。但对于复杂系统如社会系统、经济系统、水资源环境系统等,它是无能为力的。
研究基于计算机的建模理论目的:
在交互式计算机的辅助下,使系统建模过程进行得更加有效,现实和高速。
5、何谓结构模型?试举一例说明之。(P47例子)
结构模型就是描述系统各实体间的关系,以表示一个作为实体集合的系统模型,其中用集合S={s1,s2,…,sn}表示实体集合,si表示实体结合中元素,R={
第3章模型的简化
1、模型描述变量化简与动态系统的模型化简二者间的区别;
有关模型描述变量的简化方法,它适应于各类系统模型的简化,模型描述变量是系统建模的基础,它们选取的主要依据是建模的目标,而它们的选取则决定了模型的复杂程度。
建模过程中,在能满足建模的前提下,系统的描述变量应是愈简单愈好。模型描述变量一般有以下四种方法:
(1)、淘汰一个或多个实体、描述变量或相互关系规则;
(2)、随机变量取代确定性变量;
(3)、粗化描述变量;
(4)、粗化描述变量和归组实体及聚焦变量。
有关动态系统的模型简化的时域方法,它包括“集结法”和“摄动法”。
这两种方法分别是从经济理论与数学中引进来的。系统的集结法是指用一组“较粗略的”
状态变量来描述系统的模型,但应使这个系统的关键性不变。摄动法的基本概念是略去模型内部的某些相互作用,从而用一个低阶模型来逼近系统的结构,它是一种近似集结法,包括弱耦合和强耦合模型,也叫做非奇异摄动法和奇异摄动法。
2、模型描述变量化简的四种方法比较;
建模过程中,在能满足建模的前提下,系统的描述变量应是愈简单愈好。模型描述变量一般有以下四种方法:
(1)、淘汰一个或多个实体、描述变量或相互关系规则;建模者决定淘汰那些次要因素,只要忽略的因素不会显著地改变整个模型行为,相反却使不必要的复杂了。
淘汰一个实体可能要淘汰或修改其他实体:淘汰一个实体,需要淘汰所有涉及这个实体的描述变量;淘汰一个描述变量,需要淘汰或修改涉及该变量的相互关系。
(2)、随机变量取代确定性变量;在一个确定性模型中,相互关系的规则控制着整个描述变量的值。有些随机值也是由相互关系的规则确定,为了使模型相对简化,可利用概率原理,用随机变量来取代某些变量的相互关系规则,从而将影响变量转换成随机变量。
(3)、粗化描述变量;描述变量是描述模型实体条件的一种方法,变量可能出现的值表示在某一时间可找到这个实体的一种可能条件,其变量的范围集是变量可能出现的所有值的集合。描述变量的范围粗化也是一个简化过程。
粗化有以下2种方法:①舍入.根据需要,将描述变量的范围进行一定的缩小。例如,记账常用元角分,简化后只有元,角和分舍入。②归类和非一致粗化。对于归类和非一致粗化,简化前后的描述变量虽然还是一一对应,但是它们所代表的物理意义已经不同。
(4)、粗化描述变量和归组实体及聚焦变量。把具有相同性质的实体或描述变量聚焦起
来,合并并成一个实体或描述变量,这称为实体的归组和聚焦。特点:在聚焦过程中信息不受损失,且合成变量的范围粗化。
3、 什么是集结法?(P55)
4、 在集结法中,小结以下几种方法的针对对象、原理等方面异同点:精确集结法、模态集
6、 摄动法是从什么方面入手对模型进行化简的?如何实现的?
摄动法的基本概念是略去模型内部的某些相互作用,从而用一个低阶模型来逼近系统的结构,它是一种近似集结法,包括弱耦合和强耦合模型,也叫做非奇异摄动法和奇异摄动法。
奇异摄动法是针对弱耦合系统,系统可解耦为两个独立的系统,原有系统也可集结为以上两个模型。这样整个系统就可按两个分散系统来进行设计、仿真,可明显减少计算量。
非奇异摄动法针对强耦合系统,强耦合模型有一类系统,它的状态变量之间耦合较强不允许被忽略,但特征值的模之间数值差别较大,也就是过渡过程速度上差别较大。过渡过程较快的模所对应的元素(ε/21A ,ε/2A )比上半部分大一个数量级,当ε=0,矩阵A 的元素为无界,ε从ε=0到ε>0的变化称为奇异摄动。
第4章 建模的一般系统理论
1、 何谓轨迹及其意义?
令Z 表示一个集合,比如一个模型输入、状态或者输出集。基于Z 和T 的一个分段或者轨迹指的是一个从T 的区间到Z 的映射w :对某区间><10,t t ,有Z t t w >→<10,:.有时候用定义域><10,t t 分段。
其意义:轨迹Z t t w >→<10,:描述了从t0开始到t1结束的集的运动><10,t t ,在每个中,)(t w 描述了时间t 的Z 值。因此,),(T Z w ∈的含义为:存在T 的一个区间><10,t t ,使得Z t t w >→<10,:。
有状态轨迹和输出轨迹,状态轨迹为系统从初始状态运动到最终状态所通过的状态的时间纪录。状态本身可能不是可观测的。输出轨迹来纪录观测输出值的轨迹。
2、 什么是系统规范的分层描述?主要内容包含那些?
系统规范的分层描述是指对一个系统规范描述时从多个层次进行描述,而这些层次之间又是相互关联的,更高一层的描述以前一层为基础,其结构描述越来越复杂详细。系统规范的层次将不同级别的系统接口、内部结构和递归分解性的概念形式化,它的每一层将引入更多的具体化形式到内部结构的描述中.其主要内容有:观测框架,I/O 关系的观测,I/O 函数的观测,I/O 系统。
3、 I/O 关系的观测与I/O 函数的观测区别?(P79)
关系的观测概括了关于系统做为外部观察到黑箱的知识。由此,对它们就产生了两个方面的问题。第一个是从结构到行为的问题;如果知道黑箱内部的内容,用这样或那样的方法应该能够描述外部观察到的行为。第二个问题为相反的 情况:从行为到结构的问题,亦即试图通过外部观察来推断黑箱内部的结构问题。
I/O 函数观测的系统规范允许观测系统的初始状态(在进行试验之前),但不能够观测系统的最终状态(在试验之后)。而下一个高层次的系统规范将允许同时观测初始状态和最终状态,并用这种方法去直接了解转移函数和输出函数。这就是输入—输出系统规范。
4、 以I/O 系统为例说明从结构到行为的过程。(P87)
5、 系统观测同态包含了那些层次的同态?与I/O 系统同态有什么异同?
对应于系统可能是已知的、描述了的或者说明了的各个层次中的任何一个层次,都有一个适合于此层次中说明的成对系统的关系,这种关系称为同态,因为它建立起各系统之间的对应关系,借以把—个系统的特征保存在另—个系统中。 系统观测同态包含观测框架同态,I/O 关系的观测同态,I/O 函数的观测同态,I/O 系统同态。当s ’是S 的同态像,两系统证明有相同I/O 函数行为(见定理4.3),但差异是'S 的状态空间可能比S 的状态空间“小”很多,即S 的所有状态空间Q 被h 集总到'S 的所有状态空间'Q 。
6、 解释行为同态和结构同态。(P94)
第5章 连续系统的数学仿真
1已知y x dx dy +=,x=0时,y=1,取计算步长h=0.1,试用欧拉法、梯形法和4阶龙格—库塔
法求x=2h 时的y 值,并将求得的y 值与精确解x e x y x --=12)(比较,说明差异原因。
(1).欧拉法
由 计算可得
Y 0=1;
Y1=Y0+0.1(X0+Y0)=1+0.1=1.1
Y2=Y1+0.1(X1+Y1)=1.1+0.1(0.1+1.1)=1.22
精确解为:x=2h=0.2 y2=2428.12.0122.0=--e 相差:0.0228
由于欧拉法将曲线包围的面积,近似为矩形面积造成的误差相对较大。
(2).梯形法 由
可计算得
Y0=1
Y1=1.11
Y2=1.2426
相差:0.002
(3).龙格一库塔法
X0=0,y0=1;
X1=0.1, k1=1,k2=1.1,k3=1.105,k4=1.2105; y1=1.110342;
X2=0.2, k1=1.210342,k2=1.3208591,k3=1.326383,k4=1.442978; y2=1.2447;
相差:0.0019
第6章 模型的简化离散时间的系统仿真
排队系统(P141)
10、简单说明一般排队系统的三个基本组成部分。
一般的排队系统都有三个基本组成部分。
(1)到达模式:指临时实体按怎样的规律到达,一般用到达时间间隔的统计特征来描述。
(2)服务机构:指同一时刻有多少服务台(永久实体)可以接纳』临时实体,它们的服务要多少时间,它也具有一定的分布特征。
(3)排队规则:即服务台完成当前的服务后,从队列中选择下一个实体服务原则。
排队系统的基本结构
第1章建模与仿真的基本概念 参照P8例子,列举一个你相对熟悉的简单实际系统为例,采用非形式描述出来。 第2章建模方法论 1、什么是数学建模形式化的表示?试列举一例说明形式化表示与非形式化表示的区别。 模型的非形式描述是说明实际系统的本质,但不是详尽描述。是对模型进行深入研究的基础。主要由模型的实体、包括参变量的描述变量、实体间的相互关系及有必要阐述的假设组成。模型的非形式描述主要说明实体、描述变量、实体间的相互关系及假设等。 例子:环形罗宾服务模型的非形式描述: 实体 CPU,USR1,…,USR5 描述变量 CPU:Who,Now(现在是谁)----范围{1,2,…,5}; Who.Now=i表示USRi由CPU服务。 USR:Completion.State(完成情况)----范围[0,1];它表示USR完成整个程序任务的比例。参变量 X-----范围[0,1];它表示USRi每次完成程序的比率。 i 实体相互关系 (1)CPU 以固定速度依次为用户服务,即Who.Now为1,2,3,4,5,1,2…..循环运行。 X工作。假设:CPU对USR的服务时间固定,不(2)当Who.Now=I,CPU完成USRi余下的 i X决定。 依赖于USR的程序;USRi的进程是由各自的参变量 i 2、何谓“黑盒”“白盒”“灰盒”系统? “黑盒”系统是指系统内部结构和特性不清楚的系统。对于“黑盒”系统,如果允许直接进行实验测量并通过实验对假设模型加以验证和修正。对属于黑盒但又不允许直接实验观测的系统,则采用数据收集和统计归纳的方法来假设模型。 对于内部结构和特性清楚的系统,即白盒系统,可以利用已知的一些基本定律,经过分析和演绎导出系统模型。 3、模型有效性和模型可信性相同吗?有何不同? 模型的有效性可用实际系统数据和模型产生的数据之间的符合程度来度量。它分三个不同级别的模型有效:复制有效、预测有效和结构有效。不同级别的模型有效,存在不同的行为水平、状态结构水平和分解结构水平的系统描述。 模型的可信度指模型的真实程度。一个模型的可信度可分为: 在行为水平上的可信性,即模型是否重现真实系统的行为。 在状态结构水平上可信性,即模型能否与真实系统在状态上互相对应,通过这样的模型可以对未来的行为进行唯一的预测。 在分解结构水平上的可信性,即模型能否表示出真实系统内部的工作情况,而且是惟一表示出来。 不论对于哪一个可信性水平,可信性的考虑贯穿在整个建模阶段及以后各阶段,必须考虑以下几个方面: 1在演绎中的可信性。2在归纳中的可信性。3在目的方面的可信性。 4、基于计算机建模方法论与一般建模方法论有何不同?(P32) 经典的建模与仿真的主要研究思路,首先界定研究对象-实际系统的边界和建模目标,利用已有的数学建模工具和成果,建立相应的数学模型,并用计算装置进行仿真。这种经典的建
课程作业报告 课程名称 MATLAB计算及仿真院部名称机电工程 专业电气工程及其自动化班级) 学生姓名 学号 指导教师应明峰 金陵科技学院教务处制
第2章MATLAB概论 1、与其他计算机语言相比较,MA TLAB语言突出的特点是什么? MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 2、MATLAB系统由那些部分组成? MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 3、安装MA TLAB时,在选择组件窗口中哪些部分必须勾选,没有勾选的部分以后如何补安装? 在安装MATLAB时,安装内容由选择组件窗口中个复选框是否被勾选来决定,可以根据自己的需要选择安装内容,但基本平台(即MATLAB选项)必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行,只是在选择组件时只勾选要补装的组件或工具箱即可。 4、MATLAB操作桌面有几个窗口?如何使某个窗口脱离桌面成为独立窗口?又如何将脱离出去的窗口重新放置到桌面上? 在MATLAB操作桌面上有五个窗口,在每个窗口的右上角有两个小按钮,一个是关闭窗口的Close按钮,一个是可以使窗口成为独立窗口的Undock按钮,点击Undock按钮就可以使该窗口脱离桌面成为独立窗口,在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 5、如何启动M文件编辑/调试器? 在操作桌面上选择“建立新文件”或“打开文件”操作时,M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 6、存储在工作空间中的数组能编辑吗?如何操作? 存储在工作空间的数组可以通过数组编辑器进行编辑:在工作空间浏览器中双击要编辑的数组名打开数组编辑器,再选中要修改的数据单元,输入修改内容即可。 7、命令历史窗口除了可以观察前面键入的命令外,还有什么用途? 命令历史窗口除了用于查询以前键入的命令外,还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 8、如何设置当前目录和搜索路径,在当前目录上的文件和在搜索路径上的文件有什么区别? 当前目录可以在当前目录浏览器窗口左上方的输入栏中设置,搜索路径可以通过选择操作桌面的file菜单中的Set Path菜单项来完成。在没有特别说明的情况下,只有当前目录
动画建模与仿真技术(541032)教学大纲 01.教学单位:软件学院 02.课程编号:541032 03.课程名称:动画建模与仿真技术 04.课程英文名称:Technology of animation modeling and simulation 05.课程学时: 32学时,其中含实验0学时 06.课程学分: 2学分 07.课程类别:专业教育课 08.课程性质:选修 09.开课学期:第6学期 10.面向专业:软件工程 11.选用教材 1、侯鹏志等,《3ds Max 2010中文版从入门到精通》,电子工业出版社,2010年 2、杨丽等,《城市仿真建模工具:Creator软件教程》,同济大学出版社,2007年 3、王孝平等,《Vega Prime实时三维虚拟现实开发技术》,西南交通大学出版社,2012年 12.主要参考书 1、鲍虎军等,《计算机动画的算法基础》,浙江大学出版社,2000年。 2、Robert Bridson,《Fluid Simulation for Computer Graphics》,A K Peters,2007. 13.课程教学目的与任务 课程目的:主要讲授利用三维建模软件3ds Max 2010和Multigen Creator 3.0进行三维物体和场景建模的基本方法,利用实时视景仿真软件Multigen Vega Prime提供的Lynx Prime软件和Vega API进行视景仿真的基本方法。 课程任务:一、掌握利用3ds Max创建三维模型的方法 二、掌握利用Multigen Creator 创建三维模型及生成地形的方法 三、掌握利用Multigen Vega Prime进行实时视景仿真的方法。
利用Matlab 实现Romberg 数值积分算法 一、内容摘要 针对于某些多项式积分,利用Newton —Leibniz 积分公式求解时有困难,可以采用数值积分的方法,求解指定精度的近似解,本文利用Matlab 中的.m 文件编写了复化梯形公式与Romberg 的数值积分算法的程序,求解多项式的数值积分,比较两者的收敛速度。 二、数值积分公式 1.复化梯形公式求解数值积分的基础是将区间一等分时的Newton —Cotes 求积公式: I =(x)[f(a)f(b)]2 b a b a f dx -≈ +? 其几何意义是,利用区间端点的函数值、与端点构成的梯形面积来近似(x)f 在区间[a,b]上的积分值,截断误差为: 3" (b a)()12 f η-- (a,b)η∈ 具有一次的代数精度,很明显,这样的近似求解精度很难满足计算的要求,因而,可以采用将积分区间不停地对分,当区间足够小的时候,利用梯形公式求解每一个小区间的积分近似值,然后将所有的区间加起来,作为被求函数的积分,可以根据计算精度的要求,划分对分的区间个数,得到复化梯形公式: I =1 1 (b a)(b a) (x)dx [f(a)f(b)2(a )]2n b a k k f f n n -=--≈+++∑? 其截断误差为:
2" (b a)h ()12 R f η--= (a,b)η∈ 2.Romberg 数值积分算法 使用复化的梯形公式计算的数值积分,其收敛速度比减慢,为此,采用Romberg 数值积分。其思想主要是,根据I 的近似值2n T 加上I 与2n T 的近似误差,作为新的I 的近视,反复迭代,求出满足计算精度的近似解。 用2n T 近似I 所产生的误差可用下式进行估算: 12221 ()3 n n n I T T T -?=-=- 新的I 的近似值: 122 n n j T T -=?+ j =(0 1 2 ….) Romberg 数值积分算法计算顺序 i=0 (1) 002T i=1 (2) 102T (3) 012T i=2 (4) 202T (5) 112T (6) 022T i=3 (7) 302T (8) 212T (9) 122T (10) 032T i=4 (11) 402T (12) 312T (13) 222T (14) 132T … … … … 其中,第一列是二阶收敛的,第二列是四阶收敛的,第三列是六阶收敛的,第四列是八阶收敛的,即Romberg 序列。
三维建模及运动仿真 Pro/Engineer 软件集产品的三维造型设计、加工、分析、仿真及绘图等功能于一体,是一套使用方便、参数化造型精确的软件,其强大的造型功能及仿真分析功能受到众多工程人员的青睐。本节将采用Pro/E 软件,完成少齿数齿轮传动机构中所有零件的参数化建模,并对少齿数齿轮减速器进行虚拟装配,在此基础上,对传动机构进行运动仿真。 3.1 齿轮的参数化建模 3.1.1 零件分析 齿轮建模的操作步骤如下: (1)添加齿轮设计参数 (2)添加齿轮关系式 (3)创建齿轮的齿廓曲线 (4)创建螺旋线方程 (5)实体生成: 1)创建螺旋线线方程 2))拉伸 3))阵列 3.1.2 绘制齿轮 (1)新建文件: 启动PROE Wildfire4.0,单击工具栏新建工具,或单击菜单“文件/新建”。出现如图3.1所示对话框。选择系统默认“零件”,子类型“实体”方式,“名称”栏中输入“canshuhuachilun ”,同时注意关闭“使用缺省模板”。选择公制模板mmns-part-solid ,如图3.2所示,然后单击“确定”。 (2)创建齿轮程序。 选择菜单栏“工具/程序”命令,出现如图3.3所示对话框。单击“编辑设计”, 依次添加齿轮设计参数及初始值,添加完毕单击“确定”。选择工具菜单“工具/程序”命令,出现如图3.4信息窗口,在其中输入程序如下: Y0=(1/4)*PI*MT+XT*MT*TAN(α t) Xc=(HANX+CNX-XN)*MN-ρ
Yc=(1/4)*PI*MT+HANX*MN*TAN(αt)+ρ*COS(αt) (3)添加齿轮四个圆的关系式。 1)选择“插入/模型基准/ 草绘”特征工具,或单击工具栏 草绘命令,出现如图3.5所示对话框。单击“草绘”确认,进入二维草绘模式如图3.6所示。