[数学]常微分方程的数值解.pptVIP

称中点公式，相当于数值积分的中矩形公式： 如取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./(1+x.^2); subplot(221) plot(x,y1,x,y1,b*) subplot(222) plot(x,y,x,y,o) subplot(223) pl