第九章第一题作业.docVIP

第九章第一题 算法一：Runge-Kutta法 clear %从工作空间清除项，释放系统内存 clc %清除命令窗口 disp(算法一：Runge-Kutta法) %显示算法名称 %输入参数 a=0; %区间左端点 b=pi; %区间右端点 h=0.005; %步长 n=1; %记录迭代次数 x0=a; %自变量迭代初始值 y0=1; %已知初值 disp(迭代次数 数值解 精确解 误差); %输出迭代结果的标题栏 for i=a:h:b %开始迭代 x=x0+h; K1=-y0+2*cos(x0); K2=-(y0+h/2*K1)+2*cos(x0+h/2); K3=-(y0+h/2*K2)+2*cos(x0+h/2); K4=-(y0+h/2*K3)+2*cos(x0+h); y=y0+h/6*(K1+2*K2+2*K3+K4); plot(x,y,r) %绘制出由数值积分方法求得的函数图形（红色） hold on %保留该图形 exact=cos(x)+sin(x); %求出精确解 plot(x,exact,b) %绘制出精确解的函数图形（蓝色） error=exact-y; %数值解的误差 disp(sprintf(%d %15d %15d %16e\n,n,y,exact,error)) %显示迭代结果 n=n+1; %记录迭代次数 x0=x; %进行迭代 y0=y; end title(Runge-Kutta法数值解与精确解函数图) %添加图名 xlabel(x),ylabel(y) %添加坐标名 legend(数值解,精确解,2) %添加图例 算法二：四阶Adams预测——校正算法 clear %从工作空间清除项，释放系统内存 disp(算法二：四阶Adams预测——校正算法) %显示算法名称 figure; %打开另一个作图窗口 %输入参数 a=0; %区间左端点 b=pi; %区间右端点 h=0.005; %步长 n=1; %记录迭代次数 x0=a; %自变量迭代初始值 y0=1;