|
In numerical analysis, the Runge–Kutta methods are an important family of implicit and explicit iterative methods, which are used in temporal discretization for the approximation of solutions of ordinary differential equations. These techniques were developed around 1900 by the German mathematicians C. Runge and M. W. Kutta. See the article on numerical methods for ordinary differential equations for more background and other methods. See also List of Runge–Kutta methods. ==''The'' Runge–Kutta method== One member of the family of Runge–Kutta methods is often referred to as "RK4", "classical Runge–Kutta method" or simply as "''the'' Runge–Kutta method". Let an initial value problem be specified as follows. : Here, ''y'' is an unknown function (scalar or vector) of time ''t'' which we would like to approximate; we are told that , the rate at which ''y'' changes, is a function of ''t'' and of ''y'' itself. At the initial time the corresponding ''y''-value is . The function ''f'' and the data , are given. Now pick a step-size ''h''>0 and define : for ''n'' = 0, 1, 2, 3, . . . , using : 〔; 〕 :''(Note: the above equations have different but equivalent definitions in different texts).''〔, , and leave out the factor ''h'' in the definition of the stages. , and use the ''y''-values as stages.〕 Here is the RK4 approximation of , and the next value () is determined by the present value () plus the weighted average of four increments, where each increment is the product of the size of the interval, ''h'', and an estimated slope specified by function ''f'' on the right-hand side of the differential equation. * is the increment based on the slope at the beginning of the interval, using , (Euler's method) ; * is the increment based on the slope at the midpoint of the interval, using ; * is again the increment based on the slope at the midpoint, but now using ; * is the increment based on the slope at the end of the interval, using . In averaging the four increments, greater weight is given to the increments at the midpoint. If is independent of , so that the differential equation is equivalent to a simple integral, then RK4 is Simpson's rule. The RK4 method is a fourth-order method, meaning that the local truncation error is on the order of , while the total accumulated error is order . 抄文引用元・出典: フリー百科事典『 ウィキペディア(Wikipedia)』 ■ウィキペディアで「Runge–Kutta methods」の詳細全文を読む スポンサード リンク
|