常微分方程数值解实验报告3.docVIP

实验目的： 1.应用前面实验程序求解常微分方程初值问题，通过改变初值及步长，重点了解Euler方法及其改进方法、Runge-kutta方法的稳定性及收敛性。 2．熟练应用Matlab编写程序。 实验内容： 运用Euler方法及其改进方法、Runge-kutta方法分别求解无法获得精确解的常微分方程的初值问题。利用所设计的程序，改变初值及步长，观察计算结果的变化。 实验题目：，。 实验原理： 1）前述几种数值计算方法对于初值及步长的选取都有依赖性。 2）对于无法获得精确解的常微分方程初值问题，更需要分析计算结果的稳定性与收敛性。 3）程序： %fun为目标函数字符串 %x0为自变量初始值。 %y0为fun x0 ; %bou [a,b]自变量区间 %h为步长 fun y/x+sin y/x +1; bou [50,200]; a bou 1 ; b bou 2 ; x0 50; y0 5; h 10; n ceil b-a /h ; xx linspace a,b,n+1 ; yy1 zeros 1,n+1 ;yy2 zeros 1,n+1 ;yy3 zeros 1,n+1 ; xx 1 x0;yy1 1 y0;yy2 1 y0;yy3 1 y0; %%用欧拉法计算 for i 2:n+1 x xx i-1 ;y yy1 i-1 ; k eval fun ; yy1