（李荣荣）12.9基于优化方法的常微分方程边值问题数值解..docVIP

12.9 基于优化方法的常微分方程边值问题数值解 常微分方程边值问题的一般形式为： （12-2） （12-3） 式中：表示边界条件所满足的函数关系。 12.9.1 基于Matlab函数的求解方法 Matlab求解边值问题的函数为bvp4c，它采用有限差分法求解，其基本格式为： solinit = bvpinit(x, yinit, params) y = bvp4c(odefun,bcfun,solinit) 函数bvpinit输入参数依次为自变量的区间[],函数的一个猜测值。 函数bvp4c的输入参数依次为一阶微分方程或一阶微分方程组，用函数odefun定义，边界条件用函数bcfun定义，这两个函数名用户自行定义。 【例12-1】 求解二阶常微分方程 解：令 的取值区间为[]，和的猜测值分别为0和2。计算程序为： function ode_bvp1 clc;clf;clear all; solinit = bvpinit(linspace(0,4,5),[1 0]); sol = bvp4c(@twoode,@twobc,solinit); x = linspace(0,4); y = deval(sol,x); plot(x,y(1,:));%图12.18 function dy