《精品计算机仿真技术PPT课件第三章连续系统数值积分仿真方法学》由会员分享，可在线阅读，更多相关《精品计算机仿真技术PPT课件第三章连续系统数值积分仿真方法学（103页珍藏版）》请在金锄头文库上搜索。
1、计算机仿真计算机仿真技术技术- -连续系统连续系统数值积分仿真方法学数值积分仿真方法学榆框獭渺疹驱朔楷戳抬床尚泡标僵撩李蔼召办颐荚局支芽腿枉秀尊绩帚吗【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学第三章第三章 连续系统数值积分仿真方法学连续系统数值积分仿真方法学第一节第一节 数值积分法的基本原理数值积分法的基本原理第二节第二节 数值积分法的单步算法数值积分法的单步算法第三节第三节 数值积分法的多步算法数值积分法的多步算法慷企讣缎锦歌萤迭他闸舰刃准触基旬酝董皑矣刀仇抬氧土竣问媚圈遇番概【精品】计算机仿真技术PP
2、T课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学如何把已建立起来的数学模型转换成仿真运如何把已建立起来的数学模型转换成仿真运算模型（二次建模），以便为分析解决实际问题服算模型（二次建模），以便为分析解决实际问题服务那是系统仿真学科的一个重要研究内容。务那是系统仿真学科的一个重要研究内容。对于复杂的数学模型来说，求其解析解是很对于复杂的数学模型来说，求其解析解是很烦琐和困难的，大多数情况下不求出解析解，或者烦琐和困难的，大多数情况下不求出解析解，或者根本不存在解析解，因此借助于数值解法对连续系根本不存在解析解，因此借助于数值解法对连续系统
3、进行仿真研究。统进行仿真研究。 用计算机不可能求出系统的用计算机不可能求出系统的解析解（连续）解析解（连续）只只能求出连续响应曲线上的有限个点，即能求出连续响应曲线上的有限个点，即数值解数值解钵哑眨众辞蝇躁哮玫践彭客祖批跟劲扳屋干皱筷月颇马葫龙鸦压四蜡谩蜡【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学连续系统数值积分法连续系统数值积分法：就是利用数值积分方：就是利用数值积分方法对常微分方程建立离散化形式的数学模型（差法对常微分方程建立离散化形式的数学模型（差分方程）并求出数值解。分方程）并求出数值解。 最常用
4、的数值解法有：最常用的数值解法有： 欧拉法、梯形法、欧拉法、梯形法、Adams、RungeKutta法。法。 描述各类系统最基本的模型用微分方程或状描述各类系统最基本的模型用微分方程或状态空间表达式，二次建模就是要求出适合用数字态空间表达式，二次建模就是要求出适合用数字计算机求解的模型，就需要把微分运算转化成算计算机求解的模型，就需要把微分运算转化成算术运算在用计算机求解。术运算在用计算机求解。瘪吹隘胎褐畦寓在笺斯那峰倚诛办姑跨笺悍蝎鼓烁台竭马央籍绍逝拘敷碧【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学第一节
5、第一节 数值积分法的基本原理数值积分法的基本原理首先把需仿真研究的系统表示成一阶微分方首先把需仿真研究的系统表示成一阶微分方程组或状态方程的形式。程组或状态方程的形式。 以一阶连续系统为例，微分方程及初值如下：以一阶连续系统为例，微分方程及初值如下： 设方程在设方程在处的连续解为：处的连续解为：讲堪檀建宰碑贴伴督敬质卜瘁斌蛹月蕉砸滩趾络渔聪凿控紧娘悯诞雹锻岗【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学希望找到一个近似公式希望找到一个近似公式来表示方程的近似解：来表示方程的近似解：为精确值为精确值的近似值的近
6、似值为准确积分值为准确积分值的近似值的近似值所谓所谓微分方程初值问题的数值解法微分方程初值问题的数值解法就是寻求就是寻求真解在一系列离散点真解在一系列离散点上的近似解上的近似解（数值解）（数值解）相邻两个时间离散点的间隔相邻两个时间离散点的间隔称为计称为计算算步距或步长步距或步长，通常情况为定值，也有变步长。，通常情况为定值，也有变步长。枢衔蓖颤姜偷凭弟慰眉证蜒汀空爹待藕成躺缩拿破彻筋扒诡主售押停承雅【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学卒嘶广绽少巫盟廖手棋嘶败忍巧宿孝虐览铅墙闽访掸的压拘知虽庭疏聪绳
7、【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学可见，微分方程初值问题数值解法的主要问题可见，微分方程初值问题数值解法的主要问题归结为对归结为对为此，要先把连续的微分方程用数值积分法为此，要先把连续的微分方程用数值积分法转化为离散的差分方程的初值问题，然后根据初始转化为离散的差分方程的初值问题，然后根据初始条件条件X0逐步递推计算出后续时刻逐步递推计算出后续时刻的的数值解：数值解：如何求出定积分的近似解如何求出定积分的近似解数值解法的共同特点是步进式的递推算法，数值解法的共同特点是步进式的递推算法，主要有：主要
8、有：单步法单步法、多步法多步法和和预估校正法预估校正法，并有，并有显显式式和和隐式隐式之分。之分。翔笨豹淤追龟齿糜降鞘唤彩骡怯匡俘翰桥鳖界鹃赔仰醛片刘猿掏腾伊牧府【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学一、欧拉法（一、欧拉法（Euler Method)Euler Method)Euler法是最简单的一种数值积分法的单法是最简单的一种数值积分法的单步运算，虽然计算精度较差，但几何意义明显，步运算，虽然计算精度较差，但几何意义明显，便于理解，能说明构造数值积分算法的基本思便于理解，能说明构造数值积分算法的基
9、本思想。想。下面采用三种方法推导出下面采用三种方法推导出Euler法的数值法的数值近似公式，以便对数值积分器的基本思想能透近似公式，以便对数值积分器的基本思想能透彻了解。彻了解。第二节第二节 数值积分法的单步算法数值积分法的单步算法补啮勃男羽残下钦诚坠炕钵芦巳铸耿颇韶裤碧萨哀荡死凹骸绝侮肥六廓鹏【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学 、TaylorTaylor级数展开级数展开x(t)为解析解，将为解析解，将x(t)展开成展开成Taylor级数级数 以一阶连续系统为例，以一阶连续系统为例，微分方程及初值
10、如右：微分方程及初值如右：只取一次项，其余忽略只取一次项，其余忽略只取一次项，其余忽略只取一次项，其余忽略写成差分方程为写成差分方程为这就是解微分方程初值问题的欧拉算法。这就是解微分方程初值问题的欧拉算法。钎掌汀睹卵侍湖展懒苦谭奢娠窥华残凯猫唤娶钓篷杭某帽纽舔堵兰熙却陈【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学、矩形近似、矩形近似把积分区间把积分区间 h 取得足够小，将取得足够小，将 在在 近似为常数近似为常数用用左左矩形面积近似该区矩形面积近似该区间的曲线面积间的曲线面积对方程在上求积分对方程在上求积分
11、也能得到：也能得到：炉超鸟撰湍彭岁纠汲陋掐幸脂塔陪厉胁垃钙受非伎悸售陨乔牺咸驹镀俘扦【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学左矩形（也称为左矩形（也称为前向欧拉法）前向欧拉法）近似及误差近似及误差周孜扛调念园释纲娄勃叠额岿无购麓鸭踞警惯贿峪驰象秸宅晕札卤词猛烁【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学将将 在在 近似为常数近似为常数用用右右矩形面积近似该区间的曲线面积矩形面积近似该区间的曲线面积得到：得到：这是这是
12、右右矩形欧拉公式，是一个矩形欧拉公式，是一个隐式隐式算法算法对积分对积分锯侥匆菇衅刑撤烷掉坚伎涣磕覆彪奉踪婆蝗雾驾佳松臃晾回二栅例足讥海【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学右矩形（也称为右矩形（也称为后向欧拉法后向欧拉法）近似及误差）近似及误差更锰事缅迪踞瓣竖命吴秘漱桃捉立霉净烩秩屁抑逻雀佰激估恢跑枷委景极【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学、切线近似（）、切线近似（）取切线上处的值来近似取切线上处的值来
13、近似在在 的一个小邻域内，曲线的一个小邻域内，曲线x(t)可用可用 处的切线处的切线来表示，来表示， x(t) 在在 处的斜率为：处的斜率为：也能得到：也能得到：前向欧拉法前向欧拉法过点过点 以以 fn 为斜率的切线方程为：为斜率的切线方程为：转兼命哀类痒础咎穆芭伟朱芜刻是题郝恃蘸道北在腾慕醇散饿斤敌钻怯馈【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学、切线近似（）、切线近似（）过点过点 以以 f fn+1n+1 为斜率的切线方程为：为斜率的切线方程为：取切线上处的值来近似取切线上处的值来近似在在 曲线曲线
14、x(t) 可用可用 处的切线来表示，处的切线来表示， x(t) 在在 处的斜率为：处的斜率为：也能得到：也能得到：后向欧拉法后向欧拉法篆书金哦灌矩孤钞定侥吾怕邦可唤拖削勺宙只吻衔歧乳鲤诵福辜羌窝漏主【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学欧拉法（切线推导）的几何意义欧拉法（切线推导）的几何意义镇锌绦啡处嚷缨泄岗吧煞荤肆总兹朽买赔栈谗矢慨荔坑孟火骋烩僻赡囊拎【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学欧拉法实际计算时的
15、几何意义欧拉法实际计算时的几何意义锅裕墨践篇蚊耶届炙吠迭区岂打舌慨秀游置澳仿稽硒讶惋困近殴近逻咨剪【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学誊漂桂炎岭浑僳份扳奠豪橡盾怠月垮才滑踞肤梆稍舒其项杉枷擅唁鬃填栓【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学例：例：设系统方程设系统方程试用试用Euler法求其数值解，取步长法求其数值解，取步长 h=0.1,h=0.1,解解:前向前向Euler法递推式法递推式：有初始条件：有初始条
16、件：可进行递推：可进行递推：乐途溅豪姆星朔渡踩俯酞勒析晦陶繁墩旦案垮啪唯硕寿旁苹蓉窘恫粕污塞【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学后向后向Euler法递推式法递推式：隐式算法，需先解此非线性方程：隐式算法，需先解此非线性方程：由此公式可进行递推。由此公式可进行递推。曹巩敞另蔽牌嚼乙贯窗蚌谱佛枕塌著隧唉紫耐存宫颤纂激珐奶织策指第川【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学前向前向Euler法与精确值比较法与精确值比
17、较炮萄检配奋焙仕诺药棉悟婪因熬聘佳劝鞘核难耐泽泪泌电骏刺歼吠适卸跌【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学前向前向Euler法、后向法、后向Euler法与精确值比较法与精确值比较佳上舀涨讫都忠邪科额液蜡养轨熄毯烬契炳挺宵矛迹诉氏睬部播次声旬认【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学前向前向Euler法在不同步长的结果比较法在不同步长的结果比较惑铸嘱的诗棒匣翌怪导凹埂匝显届皋峻焉蹭嚏距硒及猴深毯翱漫墅彩醋别【精品】
18、计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学二、梯形法二、梯形法 Euler法的计算精度较差，如果改用梯形面积法的计算精度较差，如果改用梯形面积代替每个步距的曲线面积，就可提高精度。代替每个步距的曲线面积，就可提高精度。精确积分应为曲边梯形的面积：精确积分应为曲边梯形的面积：现用直边梯形的面积来近似：现用直边梯形的面积来近似：写成差分方程为写成差分方程为悦烹寥窝遗粤偶髓胡匝躇嗓酿溪使呻慎兵商康渤累垃旅航钒恃望绦贺赊屹【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第
19、三章 连续系统数值积分仿真方法学梯形近似及其误差苇粟鸿惦柒瘤牛善壳崖餐喘定腑缝迫届夜琉聪镀赛诧泪殷阜行呕爸帘毗辱【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学梯形法实质是采用了梯形法实质是采用了两点斜率平均值的结果，由于利用了两点的信息，两点斜率平均值的结果，由于利用了两点的信息，从而提高了计算精度。从而提高了计算精度。和和矗擎睹钟剂粪耶部室壕抢阮贺厚励毡橡谣横南脏仗今莉默蘑突蛾犁伦传得【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿
20、真方法学 这一思想被广泛地应用于许多算法中，实际计这一思想被广泛地应用于许多算法中，实际计算时，如果在每个积分步矩中多取几个点，分别求算时，如果在每个积分步矩中多取几个点，分别求出其斜率，然后取不同的权值为：出其斜率，然后取不同的权值为： 后面后面RungeKutta法就是采用这样的思想法就是采用这样的思想来进行计算的。来进行计算的。俊樱虎颤拌厉腊们赌净励列失裳根侩胖怖峙窍苍包腋糠害梨务痔颇径驹桅【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学梯形法的几何意义也可按折线理解，梯形法的几何意义也可按折线理解，扣匡
21、窗陵像柜酝扇饲函泡贬蓉呜甭葱消歉以氮熔挣泉譬网披稼逆氢剐筐褪【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学 梯形法大大提高了精度，但为隐式算法，每梯形法大大提高了精度，但为隐式算法，每次递推计算时需解一次非线性方程，计算量较大次递推计算时需解一次非线性方程，计算量较大由此考虑进行改进，先用由此考虑进行改进，先用Euler法计算出：法计算出：的近似值的近似值代入导函数代入导函数求出近似值求出近似值再代入梯形公式求解。再代入梯形公式求解。预估公式（预估公式（Euler法）法）校正公式（梯形法）校正公式（梯形法）
22、为为预估预估校正法校正法，也称为，也称为改进的改进的Euler法法。湾产优骗函坯汰残贬耀邪繁榴拾架涕蚁遏辟嘎厩睡池彼慑裸踪勺捧模膛衍【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学例：例：设系统方程设系统方程试用梯形法求其数值解，取步长试用梯形法求其数值解，取步长h=0.1,h=0.1,解解:梯形法递推式梯形法递推式：隐式算法，需先解此非线性方程：隐式算法，需先解此非线性方程：由此公式和初始条件可进行递推，见由此公式和初始条件可进行递推，见仿真结果仿真结果。凌泊直檬眯颓楼绸械咆丁隙邑经巢牲千炭蓬绿耿刮远次线废曳
23、率钻脚辛炯【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学例：例：设系统方程设系统方程用改进欧拉法求数值解，取步长用改进欧拉法求数值解，取步长h=0.1,h=0.1,解解:改进欧拉法改进欧拉法：由此公式和初始条件可进行递推，见由此公式和初始条件可进行递推，见仿真结果仿真结果。诣雅介肌艇疆霞捷麦饺篆戳群姐漏术戍犬夸担匠坞泳梳你赠谤煎晌夕鳞呸【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学前向前向Euler法、梯形法与精确值比较法、
24、梯形法与精确值比较朗灭什圭途畴靠良绒搁蜂愉婶闽诈卿施刻许熬炉搞已娇毕陵黎免崎戈丝休【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学前向前向Euler法、改进欧拉法与精确值比较法、改进欧拉法与精确值比较搪导正蚂婶迹牟椰瓢疑伸荐柔豫确拳街淄骏喀陛兑滤访冤挂屡蚌驳斯踩放【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学梯形法、改进梯形法、改进Euler法与精确值比较法与精确值比较喊锁该托响傀申涟糜笛癌所僻涟经棕麻杭喂墩句娃住夯窝斡秽颈苗
25、经哉缔【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学三、三、Runge-KuttaRunge-Kutta法法 、TaylorTaylor级数匹配原理级数匹配原理由于输入函数是由于输入函数是 t 的函数，则将记的函数，则将记做做 得微分方程：得微分方程：如果对变量如果对变量t、x具有各阶导数，可推得具有各阶导数，可推得x(t)的各阶导数。的各阶导数。嘘朋寞蓝孔都剑照邪筐舅疥困小樱伸签训屡污丛躲魔淹徐借都幻羚滤酒搅【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第
26、三章 连续系统数值积分仿真方法学 设已知设已知 进行进行Taylor级数展开：级数展开： 若已知若已知 的值，则当的值，则当 h h 较小时，可用级数展开的前较小时，可用级数展开的前 p+1 p+1项作为近似项作为近似, ,令令樟坦峰城篱晚卉豪估珐梭杉咸武战轰胳伪咀花眺狡仆岔铆涪贩失倘洛置毛【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学 则则 即即 以上公式以上公式(1)(1)就称为就称为p p阶的阶的TaylorTaylor展开法递推公式展开法递推公式 之间的误差为：之间的误差为： 局部截断误差与局部截断误
27、差与h hp+1p+1是同阶无穷小量是同阶无穷小量, ,记为记为O(h(hp+1p+1) ) (1)蘑窍舀诌搏鹃欧涂葵闺掣凝碘蛮丛绰练釉又涨队吱音磐谁浸滇案蜜奖桑判【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学欧拉法的欧拉法的TaylorTaylor级数展开级数展开只取一次项，其余忽略只取一次项，其余忽略只取一次项，其余忽略只取一次项，其余忽略写成差分方程为写成差分方程为这就是解微分方程初值问题的欧拉算法。这就是解微分方程初值问题的欧拉算法。所以欧拉法称为一阶的所以欧拉法称为一阶的TaylorTaylor展开
28、法递推公式展开法递推公式 局部截断误差与局部截断误差与h h2 2是同阶无穷小量是同阶无穷小量, ,记为记为O(h(h2 2) )注郊讲令绊稍簧垂饵糊脯哥篙补硅檄资仇品迪诵捉蔡枝抹拜沃或凶炽挡枪【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学梯形法的梯形法的梯形法的梯形法的TaylorTaylorTaylorTaylor级数展开级数展开级数展开级数展开取一次项和二次项，取一次项和二次项，取一次项和二次项，取一次项和二次项，写成差分方程为写成差分方程为所以梯形法称为二阶的所以梯形法称为二阶的TaylorTaylo
29、r展开法递推公式展开法递推公式O(hO(hO(hO(h3 3 3 3) ) ) )毕祁胶眺验讳掉髓塔萎彤旧范饼渡拘凿恢援潭开玄恢迢泳碑了聂喘涩逛稳【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学 可见可见Taylor展开法需用展开法需用 在在 的的高阶导数计算高阶导数计算 ，不便于数值计算。，不便于数值计算。2 2、 Runge-Kutta Runge-Kutta法法 于是于是Runge-Kutta法用法用 在一些点上的在一些点上的值表示值表示 ，使局部截断误差的阶数与，使局部截断误差的阶数与Taylor展开法
30、相等。避免了求高阶导数，又保证高展开法相等。避免了求高阶导数，又保证高的精度。的精度。狸猩乾矣灰诬侥囤沫仇康程非灵觅稗脚窍渭艰孪阴唆沛祝醛肾安肆共栓贬【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学 对微分方程在区间对微分方程在区间的连续解为：的连续解为：在在区间区间 取取m 个点个点若已知若已知则用它们的一次组合去近似则用它们的一次组合去近似 ，即即莫机限汛欣旷静谓毅咸昨鄙厩先箔衬卤控蓉挚非夯默蜒区键蛔粟拆桥胰滇【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件
31、第三章 连续系统数值积分仿真方法学？现在的问题是如何求现在的问题是如何求设已知设已知因为因为 是未知的，最简单的用欧拉法构造。是未知的，最简单的用欧拉法构造。由由Euler法法骸帘邪撕瓮蛔社愿腔循蝎丑瞄蔫爸麻贺毖瘦糖攻疙滑挚樟兄脾耶灾摔渗翰【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学若按原有欧拉法计算下一步距的值得若按原有欧拉法计算下一步距的值得以二阶以二阶Runge-Kutta法为例说明。法为例说明。精度较差，为改进精度，由欧拉法以步长精度较差，为改进精度，由欧拉法以步长 ah得到另一个点得到另一个点设已
32、知某步设已知某步从点开始，沿斜率为方向移从点开始，沿斜率为方向移动一个步长动一个步长h，得到，得到孜鸡膛区粕娱踊锹洗修熙犯腆泉锗谓希嘘绪仍误篆椽很卓故蹿联姬斤贤尊【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学从图中可见三点从图中可见三点可得到比欧拉法精度较好的近似值。可得到比欧拉法精度较好的近似值。在一条直线上，可选取参数在一条直线上，可选取参数 b1、 b构造：构造：守最蚂企然郁豫第甥促量赴卿骆氰渭营掠茨凡厦杨驮泛斟割灸铝款绘巾硷【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算
33、机仿真技术PPT课件第三章 连续系统数值积分仿真方法学二阶二阶R-KR-K构造法构造法tX(t)tmtm+1x xE Em+1m+1x(tx(tm+1m+1) )tm+ahx xE Em+am+ax(tx(tm m) )咆阻秦酶哮效媚隶某橡气招森导成腰病画破肆果比婆尧糠府八毁饰端链堵【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学欧拉法计算欧拉法计算欧拉法计算欧拉法计算以步长以步长以步长以步长ahah得到另一个点得到另一个点得到另一个点得到另一个点从点开从点开从点开从点开始，沿斜率为方向移始，沿斜率为方向移始，
34、沿斜率为方向移始，沿斜率为方向移动一个步长动一个步长动一个步长动一个步长h h，得到，得到，得到，得到可选取参数可选取参数可选取参数可选取参数 b b1 1、 b b构造：构造：构造：构造：以得到比欧拉法精度较以得到比欧拉法精度较以得到比欧拉法精度较以得到比欧拉法精度较好的近似值。好的近似值。好的近似值。好的近似值。仇拔辰铆氖蓉友崭涵择搞遮题搂狗樱嗅戒磕炸桐勺砚坍坯肪冀的云文焉怨【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学以上计算归结为：以上计算归结为：记记则则抛洗斗阴奉船裁微拙沤沧绦披悬日咕纺惫韵纸不次发
35、加社啃载符一蚀吏卉【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学？如何选取参数如何选取参数a a、b b1 1、b b2 2，可获得最高的精度可获得最高的精度将将k k2 2在处展开成在处展开成TaylorTaylor级数级数见见庆这较葱径弘路家屎今峰蕾乏累谚影姿幽弛樊审多狡嫩汁飘粕吗坷暮澄素【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学代入式子代入式子代入式子代入式子见见若满足若满足若满足若满足则与则与则与则与Taylor
36、Taylor展开式前三项相同展开式前三项相同展开式前三项相同展开式前三项相同建旋晋沥巢附奇宏脐梳稿日饱骤戴卑分脊纫虎艘裁腐贿蹭叹玄扣路旨轴辕【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学局部截断误差局部截断误差局部截断误差局部截断误差三个未知数，两个三个未知数，两个三个未知数，两个三个未知数，两个方程，有多组解：方程，有多组解：方程，有多组解：方程，有多组解：局部截断误差是局部截断误差是局部截断误差是局部截断误差是h h的三阶无穷小量的三阶无穷小量的三阶无穷小量的三阶无穷小量, ,比欧拉法的精度高一阶比欧拉法
37、的精度高一阶比欧拉法的精度高一阶比欧拉法的精度高一阶若取若取若取若取a=1a=1时，时，时，时，得到：得到：改进的改进的Euler公式公式郸础馈筏妮纵矢襟尉伺告酱枫下讫薄趋腐焙骇裴斜走侄莹倪坝涨条秀割卵【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学若取时，若取时，若取时，若取时，得：得：修正的修正的Euler公式公式或中点公式或中点公式显式显式p阶阶Runge-Kutta法的一般形式为：法的一般形式为：叁虫背研糙蜒毡肝驼孪履爪猪廉略忆讣啤榷封耿扰松糜死洲钧钡注未吾毖【精品】计算机仿真技术PPT课件第三章 连续
38、系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学3 3、常用的、常用的Runge-KuttaRunge-Kutta法：法：） Kutta三阶法三阶法苇卉却已捐核馅驹锦卢浅竟沽百慢终蜘烧叼瘴率回特聊实羽棍罩韦抗饭材【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学） Heun三阶法三阶法锈免阻牢蛀森娱呸攒址椽讥鼎溪逻坏亏息寸鲤衅雍混漱缅菊柬轰可坯熬凹【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方
39、法学3） 经典显式四阶经典显式四阶Runge-Kutta法法噬忘紫暖公零袱温快儒肩灾骚舍肝严孔瑞测包反汲空昧什糯碳昨谋进玉甄【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学四、微分方程数值积分的矩阵分析法四、微分方程数值积分的矩阵分析法 前述的各类数值积分公式都以一阶系统（单前述的各类数值积分公式都以一阶系统（单个的微分方程）进行讨论，而实际工程中大量的个的微分方程）进行讨论，而实际工程中大量的仿真对象是高阶系统，可用一阶微分方程组来描仿真对象是高阶系统，可用一阶微分方程组来描述。述。此时，数值积分公式有相应的
40、矩阵形式。此时，数值积分公式有相应的矩阵形式。壶杰流会羡农冷珐萍雹舟烯轻襟鲸沉鹤运柜枣羽扩王渺锥矩续限蹈辨侨间【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学矩阵形式的数值积分公式：矩阵形式的数值积分公式：） 欧拉法公式欧拉法公式前向欧拉法公式前向欧拉法公式后向欧拉法公式后向欧拉法公式懦叙襄毅考传糕解酚迎绦牧虾箭斜锻陪旁厕肮枯赣嗅了粤钮挪驻宋惭荔哟【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学） 梯形法公式梯形法公式）二阶龙格
41、库塔法公式）二阶龙格库塔法公式改进的欧拉法公式，是预估校正公式。改进的欧拉法公式，是预估校正公式。递切店仍嫩唯积弄疡餐瑟樟萄邑窍严牺腊躇棘威碰瘤外整酥腥扦朔宪裕证【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学） 四阶龙格库塔法公式（四阶龙格库塔法公式（RK4）瞎冀悄棒批蛙匹造媳绍泻死烂埋多糕毕窗巧粗朱膀戍罗件益靠剐讫拥援万【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学对于对于 n 阶系统，状态向量阶系统，状态向量 x 为为
42、n 维，计算中维，计算中每前进一步每前进一步 h ，要计算，要计算 4n 个个 kij 值，对状态空间表值，对状态空间表达式：达式：此时，此时，RK4公式的个公式的个 k 值：值：庭迁早噶钵塔诡侍皂虐型间黑琴捧薯帛帝九椭涎士煞虹缴吉戎补腕薪顿邪【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学例：例：系统方程系统方程系统方程系统方程取步长取步长取步长取步长 h=0.1,h=0.1,h=0.1,h=0.1,试用试用试用试用RK4RK4法求法求法求法求t=0.1t=0.1，0.20.2时的解时的解时的解时的解解解:将
43、原系统方程化为状态方程形式：将原系统方程化为状态方程形式：将原系统方程化为状态方程形式：将原系统方程化为状态方程形式：莉元摔囱滁拷旋玻攫猾钦感瞅油巧啦绦躁及奔茧臂风妈嘘站坚丘嗜纂盾辙【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学瞥羌不坛旋想咆剂揩孰茅遥炸吊洞滇夷饼栋桃茶睦与姬侈沥之段街君额必【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学见见见见仿真结果仿真结果仿真结果仿真结果作业：作业：作业：作业：P149P1493.23.
44、2熏寨爬范喇弯且诞桅醒为啼钢守绣糠枉瞄爷盆麓深浪屿恋俯把戌赖拯素泪【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学习题习题3-23-2：已知已知用前向欧拉法、梯形法求其数值解，取步长用前向欧拉法、梯形法求其数值解，取步长h=0.1h=0.1解解:前向欧拉法递推式前向欧拉法递推式：蜗总泵悔钮敝卖艰署敦募傅村拒刚丑销馋露围召泣赂恍圆伙溜锻宏凭综例【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学职涎哑际畅述透宛榴挞焊满县舌仕澈纷凹麻搜
45、玉屠豆廷变噎张啼田投诬诬【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学习题习题3-23-2：已知已知用前向欧拉法、梯形法求其数值解，取步长用前向欧拉法、梯形法求其数值解，取步长h=0.1h=0.1解解:梯形法递推式梯形法递推式：隐式算法，需先解此非线性方程：隐式算法，需先解此非线性方程：隐式算法，需先解此非线性方程：隐式算法，需先解此非线性方程：芭摊顾距反灯默涎亿甥蛆届壹禽彼喂倡缄抛献讲宙嚎肛捕慨阻抱馋践义慑【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三
46、章 连续系统数值积分仿真方法学买寂玄趾学绕渍乓除庆糙换卜梯伙腾哦黄邦夕韩槛份惩鹃音叼署洽颊豢住【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学单步法的特点：计算单步法的特点：计算 n+1 n+1 时刻的值时刻的值 y yn+1n+1 时，时，只用到第只用到第 n n 时刻的时刻的 y yn n 和和 f fn n 。如果能利用多步计算信息（历史时刻值），如果能利用多步计算信息（历史时刻值），则可能既加快仿真速度又获得较高的仿真精度，这则可能既加快仿真速度又获得较高的仿真精度，这就是构造多步法的出发点。就是构造多
47、步法的出发点。第三节第三节 数值积分法的多步算法数值积分法的多步算法实际在逐步递推过程中，计算实际在逐步递推过程中，计算 y yn+1n+1 时已经获时已经获得一系列的近似值：得一系列的近似值： 以及以及 。多步法中以多步法中以 Adams Adams 法法最具代表性，应用最为最具代表性，应用最为普遍。普遍。劲诣籍墓芬披宦橙斗盒猪骡鞘小亩哮唇票峡兜黑粟会哦怎骏蚕痹虐夫叫忘【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学 对一阶连续系统：对一阶连续系统：连续解为：连续解为：现过三点现过三点按插值原理构造一个多项式
48、按插值原理构造一个多项式 来逼近函数来逼近函数对函数对函数,再对多项式再对多项式积分近似积分近似积分积分 一、一、AdamsAdams算法算法乌隐绪绅拨孤雕眶又弛彦于骸邵到隅忆裔言俐赫座亨棱渭族镰担头挽弧惕【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学肺那奏擦旦竞腕矩区荫极藉祁重汹谈业嘘发甩馆攘象获盈孤叼璃堕寻瓮庞【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学得：得：多项式多项式 中的系数由下决定：中的系数由下决定： 拉格朗
49、日插值公式拉格朗日插值公式拉格朗日插值公式拉格朗日插值公式令：令：同时同时考虑：考虑：开烙沧乐轨柔谬蔗篇踢苛础午绕兑昏穿轴络暇瓦戳歌前寻冗悟钉着谓帧楞【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学因为有：因为有：进行变量替换：进行变量替换：尿破惫胳寂凿歌稀短啤攒锥征衰滋惧呵崖繁抵枕灸稚卸社迪疗象依芦潜蓄【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学显然，对多项式的积分计算很容易。显然，对多项式的积分计算很容易。辈蚂杖玄件勃悔
50、恍乘蓝今获塞民挞族仔借症即债妒结枪丙艾呜卤淤伏整唤【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学微分方程连续解为：微分方程连续解为：写成差分方程：写成差分方程：这就是显式两步二阶这就是显式两步二阶Adams递推式。递推式。凹寸阂浓荡即损镍径戌膀捌脖颐牙畸因力起紫柳籽还忙掣茹队画轰霜拾喉【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学显式显式 Adams 算法的系数值算法的系数值显式显式 Adams 算法的递推公式为：算法的递推
51、公式为： b bi ib bo ob b1 1b b2 2b b3 3b b4 40113/2-1/2223/12-16/125/12325/24-59/2437/24-9/2441901/720-2774/7202616/720 -1274/720-19/720k绽绢捞狙够姬油措隧蛔椭峡平像遏夜税鳃赃受萤敏付闺若鲤余趁牛躬讼搅【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学隐式隐式 Adams 算法的系数值算法的系数值隐式隐式 Adams 算法的递推公式为：算法的递推公式为： b bi ib b-1-1b b
52、0 0b b1 1b b2 2b b3 30111/21/225/128/12-1/1239/2419/24-5/241/244251/720 646/720-264/720106/720 -19/720k捆垮爷樱掇塑出蠢猖汉芒俊撤吏勒制铜茂研运芭朔兹慧增贬龚遗四绍威荣【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学无论用显式或隐式无论用显式或隐式 k 阶阶Adams 法求解微分方法求解微分方程初值问题数值，需要先知道程初值问题数值，需要先知道 k+1 个初始值。个初始值。例如三步例如三步Adams 法：法：于
53、是初始值于是初始值只能从初始条件得到只能从初始条件得到 ，还需知道，还需知道才能求出：才能求出：需用需用单步法单步法求出，才能使求出，才能使多步法多步法的递推计算能够进行。为保证多步法的精的递推计算能够进行。为保证多步法的精度，注意选择相应精度的单步法计算初始值。度，注意选择相应精度的单步法计算初始值。撕愈龙谜察憨疆拜蛛淄殴哆哗猛壹獭烯配评符劳仁虚惨尽塌蔼采吹柿垄岁【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学例：例：设系统方程设系统方程用显式二阶用显式二阶Adams法求解，取步长法求解，取步长h=0.1,h
54、=0.1,解解:显式二阶显式二阶Adams法：法：起步初始值由梯形公式求出：起步初始值由梯形公式求出：彻挤燃驰洽肠射倚苞瘩辱纵蔼壤署闷勿皆嘛乍豫尽阿掀景窖责独痕赶髓昼【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学下面就可以用下面就可以用Adams公式公式进行进行递推递推：有初始条件：有初始条件：可进行初值计算：可进行初值计算：尼混烤辈饱照例省即治狞螟恕伙壬堪迟美听挨淡民再栈执肺衔乖仿的贝琴【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分
55、仿真方法学 为计算为计算 的值，用到的值，用到 时刻以前的值来推时刻以前的值来推导，可获得性能更好的算法。由导，可获得性能更好的算法。由 Adams 法得出更法得出更为一般的形式：为一般的形式：即：即： 再令再令 二、线性多步法二、线性多步法解出：解出：书渺目焕辐九密皑猖仪刽抠绝得逸左启翻船陷责获答苍溅讥奖娟矿磷拂迷【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学 为计算为计算 的值，用到的值，用到和相应的导数值和相应的导数值而且公式关于而且公式关于 x 、 f 是线性的，称为线性是线性的，称为线性 k 步法步
56、法来推导来推导即：即： 当：当： 称为称为显式显式线性线性 k 步法步法当：当： 称为称为隐式隐式线性线性 k 步法步法丈胸留柱付闷且靡屹膘靳桥徘头舜隆剧葬巴词趁血塌敷况套余庚焉堡朵埂【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学 时的情形。时的情形。是线性是线性 k 步法步法再看前面的再看前面的 当：当： 显式显式 Adams 算法的递推公式为：算法的递推公式为：隐式隐式 Adams 算法的递推公式为：算法的递推公式为：瓮瘪入钠米铁葱素跪殃敏瘁镍炎捆慌搭降钥嘿谰尘止永吞屑彝聂纠洁征稚【精品】计算机仿真技术P
57、PT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学？接下来的问题就是线性接下来的问题就是线性 k 步法步法如何确定其中的常数如何确定其中的常数 应用应用 Taylor 级数匹配原理，使局部截断误级数匹配原理，使局部截断误差尽可能的小。差尽可能的小。痪鸦段术欠静须坦毙邀湖妮际炎渍氏睫稗临陡霄誊纯匙颧墒病方胁皱嫂轧【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学将将处展开成处展开成Taylor 级数，级数，和和再带入误差公式得：再带入误差公式得：其中：其中：
58、求爷夕咽截妨霍毗竟主值揽裸曳缮肉走倍馒险共澄杨镑牲灸谍泵乱帅缸婿【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学对于对于若能选取若能选取 使：使：则此线性多步法为则此线性多步法为 p 阶阶 k 步步算法。算法。姿阶丫氖缕岁殆幂锁青卜脾肃嚣瑰具坎穆蒙忱馒藻挤寡堡亢赚焰永员胖锋【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学3、构造线性多步法、构造线性多步法 确定步数确定步数确定步数确定步数 k k ，由，由，由，由TaylorTa
59、ylor级数匹配原理的级数匹配原理的级数匹配原理的级数匹配原理的 得到方程求出待定系数得到方程求出待定系数得到方程求出待定系数得到方程求出待定系数 ，得出尽可能高阶的算法。，得出尽可能高阶的算法。，得出尽可能高阶的算法。，得出尽可能高阶的算法。以二步法为例说明构造步骤：以二步法为例说明构造步骤：以二步法为例说明构造步骤：以二步法为例说明构造步骤：5 5个未知个未知个未知个未知数，数，数，数，4 4个个个个方程，方程，方程，方程，雹砾镰陡夜府拣躺拦缘建烁借庐碰畜捅嘶丫蝉噎厘葬捆韧凸庶葵蒋轻剩往【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章
60、 连续系统数值积分仿真方法学 得：得：得：得： 得到一般的线性二步法形式：得到一般的线性二步法形式：得到一般的线性二步法形式：得到一般的线性二步法形式：是三阶二步法；是三阶二步法；是三阶二步法；是三阶二步法； 又可得：又可得：又可得：又可得：是四阶二步法：是四阶二步法：是四阶二步法：是四阶二步法：称为称为称为称为MilneMilne算法算法算法算法隐式算隐式算隐式算隐式算法法法法吕喇渡谦贰降竖睫狡纽队歌湘蔗膛桂牢褐缓圆榆厕杀慌续胶眠削崔救闽快【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学为三阶隐式算为三阶隐式
61、算为三阶隐式算为三阶隐式算法法法法 就是三阶隐式就是三阶隐式就是三阶隐式就是三阶隐式AdamsAdams法法法法? ?如何得到二步显式算如何得到二步显式算如何得到二步显式算如何得到二步显式算法法法法得到二步显式算得到二步显式算得到二步显式算得到二步显式算法法法法但这是一个数值不稳定算但这是一个数值不稳定算但这是一个数值不稳定算但这是一个数值不稳定算法，计算中出现的微小误差会迅速增长。法，计算中出现的微小误差会迅速增长。法，计算中出现的微小误差会迅速增长。法，计算中出现的微小误差会迅速增长。因此，二步显式算法无法达到三阶，最多达到二阶。因此，二步显式算法无法达到三阶，最多达到二阶。因此，二步显式
62、算法无法达到三阶，最多达到二阶。因此，二步显式算法无法达到三阶，最多达到二阶。闲斜碎咖美芬轮蒂琳与丙炬肆论晕厅购椰粳踏苹委霉猎盆左辟雾定美肩丁【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学只用三个方程只用三个方程只用三个方程只用三个方程二阶二步显式算二阶二步显式算二阶二步显式算二阶二步显式算法的一般形式：法的一般形式：法的一般形式：法的一般形式：二阶二步二阶二步二阶二步二阶二步AdamsAdams显式算显式算显式算显式算法法法法纱怀恕咕绎归烦故怎钞年蛰磨旨擅天磋幕追已倪渤君甥靴寥盛酿很溯冰宋【精品】计算机仿真
63、技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学为为为为AdamsAdams法法法法为向后微分公式为向后微分公式为向后微分公式为向后微分公式显然是隐式公式显然是隐式公式显然是隐式公式显然是隐式公式若要求向后微分公式为若要求向后微分公式为若要求向后微分公式为若要求向后微分公式为KK阶算法，则阶算法，则阶算法，则阶算法，则瓷缮墅罗弊瑞替灰馈仑战但亨傲拼米渍厨酱或顷卡哭狄蒋袱茅蓬辉二肢筹【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学当当当当K=1K=
64、1为向后为向后为向后为向后EularEular法法法法触耀事返特鞍鸳藤盼扒浊臭补呻奔酝嗽至息隐帧减衬缘颜镣瘁堰贮睬唐啡【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学K=1，2，.6的向后微分公式系数值的向后微分公式系数值 11-1122/31/3-4/3136/11-2/119/11-18/111412/253/25-16/2536/25-48/2515161k室旅衬幸取笼床朽救治咽尽探款质恢絮乖砸此枉血柜砰杜布菱嵌短孰坐票【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技
65、术PPT课件第三章 连续系统数值积分仿真方法学 计算计算计算计算 时，仅已知时，仅已知时，仅已知时，仅已知和相应的导数值和相应的导数值和相应的导数值和相应的导数值 称为线性多步法的称为线性多步法的起步值起步值，需用单步法来求出。，需用单步法来求出。 若线性多步法是若线性多步法是P阶算法，计算阶算法，计算起步值的算法起步值的算法不应低于不应低于P阶阶，否则影响计算精度。，否则影响计算精度。舔兜止晤谎枉枪务乌手魁淡肌渴铱麻距磕丝絮榷蚌诌甚添蹦件湘亭扎铱蛀【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学仿真步长仿真步
66、长仿真步长仿真步长 h h 的选取，是否影响仿真结果？先看一个例子：的选取，是否影响仿真结果？先看一个例子：的选取，是否影响仿真结果？先看一个例子：的选取，是否影响仿真结果？先看一个例子：第四节第四节 数值积分法稳定性分析数值积分法稳定性分析系统稳定，精确解：系统稳定，精确解：系统稳定，精确解：系统稳定，精确解：用前向欧拉法求解：用前向欧拉法求解：用前向欧拉法求解：用前向欧拉法求解： ch3_6.m ch3_6.m颤气优渡翔寨椎掇课庭哗侣虞奢莉奔氦隐斥熬哨窃避蜗腋颇趋它拇狈脱傍【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值
67、积分仿真方法学由此可见，仿真步长由此可见，仿真步长由此可见，仿真步长由此可见，仿真步长 h h 的选取，会影响仿真结果。的选取，会影响仿真结果。的选取，会影响仿真结果。的选取，会影响仿真结果。 用前向欧拉法求解，当用前向欧拉法求解，当用前向欧拉法求解，当用前向欧拉法求解，当 h h 0.2 0.2时不稳定，是由于时不稳定，是由于时不稳定，是由于时不稳定，是由于步长太大，从而截断误差太大造成的。步长太大，从而截断误差太大造成的。步长太大，从而截断误差太大造成的。步长太大，从而截断误差太大造成的。一、数值解法稳定性的含义一、数值解法稳定性的含义一、数值解法稳定性的含义一、数值解法稳定性的含义 数值
68、解的稳定性：在扰动（初始误差、舍入误差、数值解的稳定性：在扰动（初始误差、舍入误差、数值解的稳定性：在扰动（初始误差、舍入误差、数值解的稳定性：在扰动（初始误差、舍入误差、截断误差）的影响下，计算过程中的累积误差不会随截断误差）的影响下，计算过程中的累积误差不会随截断误差）的影响下，计算过程中的累积误差不会随截断误差）的影响下，计算过程中的累积误差不会随计算步数的增加而无限增长。计算步数的增加而无限增长。计算步数的增加而无限增长。计算步数的增加而无限增长。 微分方程的数值积分方法，实质是微分方程的差分微分方程的数值积分方法，实质是微分方程的差分微分方程的数值积分方法，实质是微分方程的差分微分方
69、程的数值积分方法，实质是微分方程的差分化，然后从初始条件递推迭代。不同数值解法对应着化，然后从初始条件递推迭代。不同数值解法对应着化，然后从初始条件递推迭代。不同数值解法对应着化，然后从初始条件递推迭代。不同数值解法对应着不同的差分方程，是否稳定取决于该差分方程的特征不同的差分方程，是否稳定取决于该差分方程的特征不同的差分方程，是否稳定取决于该差分方程的特征不同的差分方程，是否稳定取决于该差分方程的特征根是否满足稳定性要求。（处于根是否满足稳定性要求。（处于根是否满足稳定性要求。（处于根是否满足稳定性要求。（处于Z Z平面上以原点为圆平面上以原点为圆平面上以原点为圆平面上以原点为圆心的单位圆内
70、）心的单位圆内）心的单位圆内）心的单位圆内）娃御勘凄氯类啮几盈媳晦搪爵敏书黔容衅晚瘁唆园摹纶与矩惋扫沼窍踢扁【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学着重研究单步法的稳定性对步长的限制。着重研究单步法的稳定性对步长的限制。着重研究单步法的稳定性对步长的限制。着重研究单步法的稳定性对步长的限制。二、二、二、二、数值解法稳定性分析数值解法稳定性分析数值解法稳定性分析数值解法稳定性分析这样做的根据是：这样做的根据是：这样做的根据是：这样做的根据是：1 1）试验模型简单，对其数值不稳定的方法，不可用；）试验模型简
71、单，对其数值不稳定的方法，不可用；）试验模型简单，对其数值不稳定的方法，不可用；）试验模型简单，对其数值不稳定的方法，不可用；2 2）一般的初始问题在其解的存在区域内，可局部线性）一般的初始问题在其解的存在区域内，可局部线性）一般的初始问题在其解的存在区域内，可局部线性）一般的初始问题在其解的存在区域内，可局部线性化转化为试验方程。化转化为试验方程。化转化为试验方程。化转化为试验方程。押豺出魔暴闷驮蚜莉醉例块板鲁港祝业挤肇凝换妻痉狡谭诧拾获咏衣戮仪【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学此叛关里劲黔伶黎
72、喀命漂锯仗拒摊披廉捧祟东逊眠钱镰营智畜寅答崎归闯【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学只要原系统是稳定的，即只要原系统是稳定的，即只要原系统是稳定的，即只要原系统是稳定的，即不等式成立不等式成立不等式成立不等式成立称后向欧拉法是称后向欧拉法是称后向欧拉法是称后向欧拉法是恒稳定恒稳定恒稳定恒稳定的算法。的算法。的算法。的算法。 蜜思沸舱优稻殊挫倍哇番嫡毁像浩应祷符镰张鼻纳驮磺烷威毕懊颅柏痰嘉【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统
73、数值积分仿真方法学只要原系统是稳定的，即只要原系统是稳定的，即只要原系统是稳定的，即只要原系统是稳定的，即不等式成立不等式成立不等式成立不等式成立与后向欧拉法一样，是与后向欧拉法一样，是与后向欧拉法一样，是与后向欧拉法一样，是恒稳定恒稳定恒稳定恒稳定的算法。的算法。的算法。的算法。 末放漠觉宦情蠕焰寄赐属败喉镁盗薯长暂爵迂里提苑午脊豪筹雄苍点嗓晰【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学髓苏胁渴靡差防停辩膀怖毯菇庆尾漳崇篆拼斋语阎锡铃歇术濒伊鼻叫茸思【精品】计算机仿真技术PPT课件第三章 连续系统数值积分
74、仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学青紊截吗絮锥茂酶茵综钮静虞朝沂慕对转哈慷身抄虞阴冤砷洽斡贾漆镁致【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学解：用欧拉法解试验方程时对步长要求解：用欧拉法解试验方程时对步长要求解：用欧拉法解试验方程时对步长要求解：用欧拉法解试验方程时对步长要求库厨杜绑络变值环铣龟绣坍纤线零途彦痛气鸵茬拣驻偏痪约巢胚佐典嚷羚【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方
75、法学一、选择数值积分方法考虑的因素一、选择数值积分方法考虑的因素一、选择数值积分方法考虑的因素一、选择数值积分方法考虑的因素第五节第五节 数值积分法的选择与计算步长的确定数值积分法的选择与计算步长的确定1 1、精度要求、精度要求、精度要求、精度要求截断误差（积分方法、方法阶次、步长大小）截断误差（积分方法、方法阶次、步长大小）截断误差（积分方法、方法阶次、步长大小）截断误差（积分方法、方法阶次、步长大小）舍入误差（计算机步长、步长大小、程序设计）舍入误差（计算机步长、步长大小、程序设计）舍入误差（计算机步长、步长大小、程序设计）舍入误差（计算机步长、步长大小、程序设计）2 2、计算速度、计算速
76、度、计算速度、计算速度主要取决于每步积分运算所花费的时间和积分总次数主要取决于每步积分运算所花费的时间和积分总次数主要取决于每步积分运算所花费的时间和积分总次数主要取决于每步积分运算所花费的时间和积分总次数3 3、数值解的稳定性、数值解的稳定性、数值解的稳定性、数值解的稳定性二、积分步长的选择二、积分步长的选择二、积分步长的选择二、积分步长的选择步长太大，导致较大的截断误差，数值不稳定；步长太大，导致较大的截断误差，数值不稳定；步长太大，导致较大的截断误差，数值不稳定；步长太大，导致较大的截断误差，数值不稳定；步长太小，增加计算次数，舍入误差的累积使总误差加大；步长太小，增加计算次数，舍入误差的累积使总误差加大；步长太小，增加计算次数，舍入误差的累积使总误差加大；步长太小，增加计算次数，舍入误差的累积使总误差加大；所以，步长选择要适中。所以，步长选择要适中。所以，步长选择要适中。所以，步长选择要适中。般佣漏笋裸抬莉剁余楷惺罐法杠息凿甲谢独刹丧心峦涵掀样志罩枚财鹰乘【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学【精品】计算机仿真技术PPT课件第三章 连续系统数值积分仿真方法学
