第六章控制糸统仿真 我们考虑计算不同点上f(x,y)的值,由欧拉公式: y+hf(rn,yn) 用高阶差分代替各级导数这一思想得到各阶 Runge- Kutta法, k1=hf(xn,yn)计算不同点上的函数值的线性组合 定义运算格式为: 机f(xn+ahy+∑bk 为待定参数 由式(1)用次重新整理得ym=yn+∑ck 通过台劳展开得 yO x, +h)=y(n)+hy(xn)+hy(xn)+ 2 相比较系数,尽可能重合到较高幂次,而求得待定参数。我们考虑计算不同点上 f x y ( , ) 的值，由欧拉公式： 用高阶差分代替各级导数这一思想得到各阶Runge-Kutta法。 定义运算格式为： 1 ( , ) n n k hf x y = 1 1 ( , ), 1,2, i i n i n ij j j k hf x a h y b k i n − = = + + =  计算不同点上的函数值的线性组合 , , i ij i a b c 为待定参数 由式（1）用h幂次重新整理得 1 1 n n n i i i y y c k + = = + 通过台劳展开得 1 ( , ) n n n n y y hf x y + = + （1） ' 2 '' 1 ( ) ( ) ( ) ( ) 2 ! n n n n y x h y x hy x h y x + = + + + 相比较系数，尽可能重合到较高幂次，而求得待定参数