（李荣荣）129基于优化方法的常微分方程边值问题数值解.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 dydx = twoode(x,y) dydx = [ y(2) -abs(y(1))]; function res = twobc(ya,yb) res = [ ya(1) yb(1)+ 2]; 图12.18 例12.1 数值解曲线 12.9.2 求解两点边值问题的打靶法 打靶法采用求初值的方法求解边值问题。用求解初值的方法求解式（12-2）和式（12-3）描述的问题时，缺少处的条件。如果能找到合适的值，通过求解恰好能使计算出的等于给定的边界值，则问题就得到解决。显然，寻找的过程是一个反复迭代的过程。就好比打靶，子弹以合适的角度和初速度射出才能命中目标。这就是打靶法的基本思想。 设一两点边值问题为： ， (12-4) 对应的初值问题为： ， （12-5） 由该初值问题求解得出的另一端边界值为：； 则应满足的方程为： 。 因此，应用求初值问题的方法求边值问题转换成求上述方程的根的问题。求一元函数根的方法很多，下面是用插值修正的方法求初始值，其迭代格式为： 【例122】 求二阶常微分方程的解。 初始条件为：（猜测值），计算程序如下 function shoot1 clc; ii=0; tol=1e-9; xspan=[0,2]; guess=1; g1=guess; target=1; y0=[0;g1]; [x,y]=ode45(@dEqs,xspan,y0); t1=y(end,1) subplot(2,1,1); plot(x,y(:,1)) g2=1.1*g1; y0=[0;g2]; [x,y]=ode45(@dEqs,xspan,y0); t2=y(end,1) subplot(2,1,2); plot(x,y(:,1)) while abs(t2-target)tol g=g2+(g2-g1)*(target-t2)/(t2-t1); t1=t2; y0=[0,g]; [x,y]=ode45(@dEqs,xspan,y0); t2=y(end,1); g1=g2; g2=g; g,g1,g2,t1,t2 ii=ii+1; if ii200 break end end figure(2) plot(x,y(:,1)) x,y function F=dEqs(x,y) % First-order differential F=[y(2);-3*y(1)*y(2)]; % equations. 计算结果为： y(0)=0,1.5145,y(2)=1, 0.0145。 （12-6） 动量方程为 （12-7） 在边界层外缘 （12-8） 能量方程为 ， 无量纲温度 （12-9） 对于自然对流，动量方程为 （12-10） 在边界层外缘 （12-11） 能量方程的形式与（12-）。