《龙格库塔方法ppt课件.ppt》由会员分享，可在线阅读，更多相关《龙格库塔方法ppt课件.ppt（19页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分1变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分2龙格库塔法（Runge-Kutta）数值分析中，龙格库塔法（Runge-Kutta）是用于模拟常微分方程的解的重要的一类隐式或显式迭代法。这些技术由数学家卡尔龙格和马丁威尔海姆库塔于1900年左右发明。经典四阶龙格库塔法龙格库塔法的家族中的一个成员如此常用，以至于经常被称为“RK4”或者就是“龙格库塔法”。变电站电气主接线
2、是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分3变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分4112341211231224130022622/( ,.,.初始条件：对于问题方程如下：），nnnnnnnnnndydxf t yy tyhyyKKKKKfxyhKfxyKhKfRuxyKKfxyhngeKKutta变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一
3、个重要组成部分5这样，下一个值(yn+1)由现在的值(yn)加上时间间隔(h)和一个估算的斜率的乘积决定。该斜率是以下斜率的加权平均：k1是时间段开始时的斜率；k2是时间段中点的斜率，通过欧拉法采用斜率 k1 来决定 y在点 tn + h/2的值；k3也是中点的斜率，但是这次采用斜率 k2决定 y值；k4是时间段终点的斜率，其 y值用 k3 决定。当四个斜率取平均时，中点的斜率有更大的权值：变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分6误差分析：误差分析：注意上述公式对于标量或者向量函数(y可以是向量)
4、都适用。 四阶R-K方法的每一步需要计算四次函数值f，可以证明其局部截断误差为O(h5).变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分7R-K(R-K(高阶高阶) )方法不唯一方法不唯一, ,选择不同的选择不同的参数能得到参数能得到不同的不同的R-KR-K公式公式注意的问题注意的问题R-KR-K方法的推导是基于方法的推导是基于TaylorTaylor展开展开法，因而要求法，因而要求解具有较好的光滑性，如果光滑性解具有较好的光滑性，如果光滑性较差精度可较差精度可能不如改进能不如改进EulerEuler方法
5、方法, ,最好采用最好采用低低阶算法阶算法而将步长而将步长h 取小取小。Runge-Kutta法的主要运算在于计算法的主要运算在于计算 Ki 的值，即计算的值，即计算 f 的值。计算量与可达到的最高精度阶数的关系：的值。计算量与可达到的最高精度阶数的关系：753可达到的最高可达到的最高精度精度642每步须算每步须算Ki 的的个数个数)(2hO)(3hO)(4hO)(5hO)(6hO)(4hO)(2nhO8n变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分8四阶四阶Runge-KuttaRunge-Kutta
6、方法的方法的MATLABMATLAB实现原理：实现原理：变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分9开始开始0 0, ,x y hN 输输 入入10100200100200312342232210226;( , ),(,)(,), (,)()xxhhhKf x y Kf xyKhhKf xyKf xh yhKhyyKKKK输出输出x1,y1?n N10101 ;,n Nx x y y 结束结束YN变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力
7、系统接线组成中一个重要组成部分10function ff=rk(yy,x0,y0,h,a,b)%yy为y的导函数,x0,y0,为初值,h为步长,a,b为区间c=(b-a)/h+1;i1=1; %c为迭代步数；i1为迭代步数累加值y=y0;z=zeros(c,6); %z生成c行，6列的零矩阵存放结果；%每行存放c次迭代结果，每列分别存放k1k4及y的结果变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分11不断迭代运算：for x=a:h:b if i1=c k1=feval(yy,x,y); k2=feva
8、l(yy,x+h/2,y+(h*k1)/2); k3=feval(yy,x+h/2,y+(h*k2)/2); k4=feval(yy,x+h,y+h*k3); y=y+(h/6)*(k1+2*k2+2*k3+k4); z(i1,1)=x;z(i1,2)=k1;z(i1,3)=k2;z(i1,4)=k3;z(i1,5)=k4;z(i1,6)=y; i1=i1+1; endend变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分12 1) 0 () 10 (2) 1 . 0()43 (yxyxyyhKR步长步长解
9、初值问题解初值问题公式公式阶阶阶、阶、求求例例4解解yxyyxf/2),( .2,2,2,462131213211hKhKyhxfKKhyhxfKyxfKKKKhyynnnnnnnn例例 题题 4变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分13xnYn|yn-y(xn)|R-K3误差误差y(xn)0.11.0959 0.00051.095440.45e-41.09540.21.1841 0.0009 1.183220.17e-41.18320.31.2662 0.0013 1.264910.15e-41
10、.26490.41.3434 0.0018 1.341650.48e-41.34160.51.4164 0.0022 1.414220.25e-41.41420.61.4860 0.0028 1.483260.55e-41.48320.71.5525 0.0033 1.549210.14e-41.5492变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分14改进改进EulerEuler法一步需要计算两个函数值法一步需要计算两个函数值(h=0.1)(h=0.1)四阶四阶Runge-KuttaRunge-Kutt
11、a方法一步需要计算四方法一步需要计算四个函数值（个函数值（h=0.2h=0.2）总计算量大致相当，但四阶总计算量大致相当，但四阶Runge-Runge-KuttaKutta方法精度更高方法精度更高变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分15 从每一步看，步长越小，截断误差从每一步看，步长越小，截断误差越小；但随着步长的缩小，在一定越小；但随着步长的缩小，在一定求解范围内所要完成的步数就会增求解范围内所要完成的步数就会增加，步数的增加不但引起计算量的加，步数的增加不但引起计算量的增大，而且可能导致舍入
12、误差的严增大，而且可能导致舍入误差的严重积累，因此需要选择步长重积累，因此需要选择步长如何衡量和检验计算结果的精度如何衡量和检验计算结果的精度如何依据所判定的精度来处理步长如何依据所判定的精度来处理步长变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分16 以经典四阶以经典四阶Runge-KuttaRunge-Kutta方法为例方法为例155151( )( )( )( ) ,hnnhnnnxhyC hC y xx xn n+ +1 1设设 从从 节节 点点出出 发发 ， 先先 以以为为 步步 长长 求求 出出
13、 一一 个个 近近 似似 值值 ， 记记 为为因因 为为 经经 典典 格格 式式 的的 局局 部部 截截 断断 误误 差差 为为 O O( (h h) ), ,因因 此此 有有y y( (x x ) )- -y y其其 中中与与在在内内 的的 值值 有有 关关2112( )+hnnnhxxy将将 步步 长长 折折 半半 ， 即即 取取为为 步步 长长 从从跨跨 两两 步步 到到 ， 求求 得得 一一 个个 近近 似似 值值5521222( )( ),( )hnhhCCn n+ +1 1每每 跨跨 一一 步步 截截 断断 误误 差差 是是有有 y y( (x x ) )- -y y变电站电气主接
14、线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分1721111116116( )( )()()hnnhnny xyy xy步步 长长 折折 半半 后后 ， 误误 差差 大大 约约 减减 少少 为为 原原 来来 的的 ， 即即 有有221111115( )( )( )()hhhnnnny xyyy事事 后后 估估 计计 式式 为为变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分18 可以通过检查步长折半前后两次计可以通过检查步长折半前后
15、两次计算结果的偏差算结果的偏差211( )( )=hhnnyy 来判断所选取的步长是否合适来判断所选取的步长是否合适, （ 1 1） 对对 于于 给给 定定 的的 精精 度度， 如如 果果则则 反反 复复 将将 步步 长长 折折 半半进进 行行 计计 算算 直直 至至为为 止止 ， 这这 时时 取取 折折 半半 以以 后后 的的 “ 新新 值值 ”作作 为为 结结 果果 ；（ 2 2） 如如 果果， 则则 反反 复复 将将 步步 长长 加加 倍倍 ， 直直 至至为为 止止 ，这这 时时 取取 步步 长长 加加 倍倍 前前 的的 “ 老老 值值 ” 作作 为为 结结 果果变步长方变步长方法法变电站电气主接线是指变电站的变压器、输电线路怎样与电力系统相连接，从而完成输配电任务。变电站的主接线是电力系统接线组成中一个重要组成部分19Thanks！