试题二(龙格库塔方法_区间估计_约束优化).docVIP

计算方法（数学实验）试题（第2组） 2000.6.22 班级 姓名 学号 说明：（1）1，2题必做，答案直接填在试题纸上； （2）3，4题任选1题，将简要解题过程和结果写在试题纸上； （3）解题程序以网络作业形式提交，文件名用英文字母。 1. 设用数值解法算出 y(1)= 1.1635 ,你用的方法是 龙格-库塔方法 ，调用的 Matlab命令是： ts=0:0.1:2; y0=[1,0]; [x,y]=ode45(@cwf, ts,y0); ode45(@cwf, ts,y0) 算法精度为 4阶 。 %待解常微分方程组函数M文件源程序： function dy=cwf(x,y) dy=[y(2);y(1)*sin(x)]; %应用欧拉方法和龙格-库塔方法求解该常微分方程： ts=:0.1:2; y0=[1,0]; [x,y]=ode45(@cwf, ts,y0); %龙格-库塔方法求数值解 [, y(:,1)] 输出结果： 1.000000000000000 1.163536347595370 注