sxjm53微分方程数值解.pptVIP

即 c1 = 1-a ， ?2 = ?21 =1/(2a) 要使上式p=2阶，则需 方程组解不唯一，可令c2=a ? 0 ，则 满足上述条件的公式都为2阶R-K公式。 称中点公式，相当于数值积分的中矩形公式： 如取a= ? ，则c1= c2= ?， ?2=?21 =1，即为改进Euler公式。 若取a= 1，则c1= 0，c2= 1， ?2=?21 = ? ，得 例:蛇形曲线的初值问题 令f(x,y)=y/x –2y2, 取f(0,0)=1, h=0.2, xn=nh , ( n = 1,2,…,15) 2阶龙格-库塔公式计算格式： k1=yn/xn – 2yn2, k2 = (yn+hk1)/(xn+h) – 2(yn+hk1)2 yn+1=yn + 0.5h[ k1 + k2] x0=0;y0=0;h=.2; x=.2:h:3; k1=1; k2=(y0+h*k1)/x(1)-2*(y0+h*k1)^2; y(1)=y0+.5*h*(k1+k2); for n=1:14 k1=y(n)/x(n)-2*y(n)^2; k2=(y(n)+h*k1)/x(n+1)-2*(y(n)+h*k1)^2; y(n+1)=y(n)+0.5*h*(k1+k2); end y1=x./