- 1、本文档共47页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
第三章 连续系统的数字仿真; 用数字计算机来仿真或模拟一个连续控制系统的目的就是求解系统的数学模型。由控制理论知,一个n阶连续系统可以被描述成由n个积分器组成的模拟结构图。因此利用数字计算机来进行连续系统的仿真,从本质上讲就是要在数字计算机上构造出n个数字积分器,也就是让数字计算机进行n次数值积分运算。可见,连续系统数字仿真中的最基本的算法是数值积分算法。;3.1 数值积分法; 当t>t0时,x(t)是未知的,因此式(3-2)右端的积分是求不出的。为了解决这个问题,??积分间隔取得足够小,使得在tk与tk+1之间的f(t,x(t))可以近似看作常数f(tk,x(tk)),这样便得到用矩形公式积分的近似公式
或简化为
这就是欧拉公式。; 以x(t0)=x0作为初始值,应用欧拉公式,就可以一步步求出每一时刻tk的xk值,;3.1.2 梯形法 由上可知欧拉公式中的积分是用矩形面积f(tk,xk)h 来近似的。图3-2为梯形法的几何意义。; 由于上式右边包含未知量xk+1,所以每一步都必须通过迭代求解,每一步迭代的初值xk+1(0)通常采用欧拉公式来计算,因此梯形法每一步迭代公式为
(3-3)
式中 迭代次数R=0,1,2,…
;3.1.3 预估-校正法
虽然梯形法比欧拉法精确,但是由于每一步都要进行多次叠代,计算量大,为了简化计算,有时只对式(3-3)进行一次叠代就可以了,因此可得
通常称这类方法为预估-校正方法。它首先根据欧拉公式计算出xk+1的预估值xk+1(0) ,然后再对它进行校正,以得到更准确的近似值xk+1(1)。;3.1.4 龙格-库塔法
根据泰勒级数将式(3-1)在tk+1=tk+h时刻的解xk+1=x(tk+h) 在tk附近展开,有
(3-5)
可以看出,提高截断误差的阶次,便可提高其精度,但是由于计算各阶导数相当麻烦,所以直接采用泰勒级数公式是不适用的,为了解决提高精度问题,龙格和库塔两人先后提出了间接使用泰勒级数公式的方法,即用函数值f (t,x)的线性组合来代替f (t,x)的导数,然后按泰勒公式确定其中的系数, 这样既能避免计算f (t,x)的导数,又可以提高数值计算精度,其方法如下。;因
故式(3-5)可写成
(3-6)
为了避免计算式(3-6)中的各阶导数项,可令xk+1由以下多项式表示。
(3-7);式中 am为待定因子,v为使用f函数值的个数,km满足下列方程
(3-8)
即:
将式(3-7)展开成h的幂级数并与微分方程式(3-1)精确解式(3-6)逐项比较,便可求得式(3-7)和式(3-8)中的系数am ,bmj和cm等。;现以v=2为例,来说明这些参数的确定方法。设v=2,则有
(3-9)
将k1和k2在同一点(tk ,xk)上用二元函数展开为;将k1和k2代入式(3-9)整理后可得
(3-10)
将上式与式(3-6)逐项进行比较,可得以下关系式
若取
则
;于是由式3-9可得
(3-11)
由于式(3-11)只取到泰勒级数展开式的h2项,故称这种方法为两阶龙格-库塔法,其截断误差为0(h3)。; 同理当v=4时,仿照上述方法可得如下四阶龙格-库塔公式
; 通过上述龙格-库塔法的介绍,可以把以上介绍的几种数值积分法统一起来,它们都是基于在初值附近展开成泰勒级数的原理,所不同的是取泰勒级数多少项。欧拉公式仅取到h项,梯形法与二阶龙格库塔法相同,均取到h2项,四阶龙格库塔法取到h4项。从理论上讲,取得的项数愈多,计算精度愈高,但计算量愈大,愈复杂,计算误差也将增加,因此要适当的选择。目前在数字仿真中,最常用的是四阶龙格库塔法,其截断误差为0(h5), 已能满足仿真精度的要求。;3.1.5 关于仿
您可能关注的文档
- 《拆掉思维里的墙》阅读版推荐培训讲学.ppt
- 《收入与分配》课件幻灯片资料.ppt
- 《政治生活》第二单元第三课:政府的职能和责任教学提纲.ppt
- 《敬畏自然》教学教案.ppt
- 《数值分析》第五讲:常微分方程数值解培训教材.ppt
- 《数学分析》第三章函数极限4教材课程.ppt
- 《数学分析》第十六章多元函数的极限与连续1演示教学.ppt
- 《数据库性能优化》复习课程.ppt
- 《新材料作文的审题立意精教版》高峰教学教案.ppt
- 《时尚女魔头》管理学分析教材课程.ppt
- 2020年吉林省中级会计职称考试实行无纸化考试.pdf
- 2020建筑工程劳务合同书(示范合同).pdf
- 2023年卫生招聘考试之卫生招聘(文员)强化训练试卷A卷附答案.docx
- 2022-2023年教师招聘之幼儿教师招聘通关提分题库(考点梳理).docx
- 2023年军队文职人员招聘之军队文职法学提升训练试卷A卷附答案.docx
- 2022-2023年企业人力资源管理师之四级人力资源管理师真题附答案.docx
- 2022-2023年企业人力资源管理师之二级人力资源管理师高分通关题型题库附解析答案.docx
- 2023年企业人力资源管理师之三级人力资源管理师通关试题库(有答案).docx
- 广西桂林市七年级下学期英语期中试卷(含答案).docx
- 2015年人美版小学三年级下册美术教案.pdf
文档评论(0)