资源描述：
《Matlab求解边值问题方法+例题.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、把待解的问题转化为标准边值问题因为边值问题可以多解，所以需要为期望解指定一个初始猜测解。该猜测解网（Mesh）包括区间[a,b]内的一组网点（Meshpoints）和网点上的解S(x)根据原微分方程构造残差函数利用原微分方程和边界条件，借助迭代不断产生新的S(x)，使残差不断减小，从而获得满足精度要求的解Matlab求解边值问题方法：bvp4c函数solinit=bvpinit(x,v,parameters)生成bvp4c调用指令所必须的“解猜测网”sol=bvp4c(odefun,bcfun,soli
2、nit,options,p1,p2,…)给出微分方程边值问题的近似解sxint=deval(sol，xint)计算微分方程积分区间内任何一点的解值Matlab求解边值问题的基本指令solinit=bvpinit(x,v,parameters)x指定边界区间[a,b]上的初始网络，通常是等距排列的（1×M）一维数组。注意：使x(1)=a，x(end)=b；格点要单调排列。v是对解的初始猜测solinit（可以取别的任意名）是“解猜测网（Mesh）”。它是一个结构体，带如下两个域：solinit.x是表示初
3、始网格有序节点的（1×M）一维数组，并且solinit.x(1)一定是a，solinit.x(end)一定是b。M不宜取得太大，10数量级左右即可。solinit.y是表示网点上微分方程解的猜测值的（N×M）二维数组。solinit.y(:,i)表示节点solinit.x(i)处的解的猜测值。初始解生成函数：bvpinit()sol=bvp4c(odefun,bcfun,solinit,options,p1,p2,…)输入参数：odefun是计算导数的m函数文件。该函数的基本形式为：dydx=odefu
4、n(x,y,parameters,p1,p2,…)，在此，自变量x是标量，y，dydx是列向量。bcfun是计算边界条件下残数的m函数文件。其基本形式为：res=bcfun(ya,yb,parameters,p1,p2,…)，文件输入宗量ya,yb是边界条件列向量，分别代表y在a和b处的值。res是边界条件满足时的残数列向量。注意：例如odefun函数的输入宗量中包含若干“未知”和“已知”参数，那么不管在边界条件计算中是否用到，它们都应作为bcfun的输入宗量。输入宗量options是用来改变bvp4c
5、算法的控制参数的。在最基本用法中，它可以缺省，此时一般可以获得比较满意的边值问题解。如需更改可采用bvpset函数，使用方法同odeset函数。输入宗量p1,p2等表示希望向被解微分方程传递的已知参数。如果无须向微分方程传递参数，它们可以缺省。边值问题求解指令：bvp4c()输出参数：输出变量sol是一个结构体sol.x是指令bvp4c所采用的网格节点；sol.y是y(x)在sol.x网点上的近似解值；sol.yp是y'(x)在sol.x网点上的近似解值；sol.parameters是微分方程所包含的未
6、知参数的近似解值。当被解微分方程包含未知参数时，该域存在。边值问题求解指令：bvp4c()原方程组等价于以下标准形式的方程组：solinit=bvpinit(linspace(0,1,10),[10]);sol=bvp4c(@ODEfun,@BCfun,solinit);x=[0:0.05:0.5];y=deval(sol,x);xP=[0:0.1:0.5];yP=[0-0.41286057-0.729740656...-0.95385538-1.08743325-1.13181116];plot(xP
7、,yP,'o',x,y(1,:),'r-')legend('AnalyticalSolution','NumericalSolution')%定义ODEfun函数functiondydx=ODEfun(x,y)dydx=[y(2);y(1)+10];%定义BCfun函数functionbc=BCfun(ya,yb)bc=[ya(1);yb(1)];求解两点边值问题：令：边界条件为：边值问题的求解原方程组等价于以下标准形式的方程组：solinit=bvpinit(linspace(0,1,10),[01]
8、);sol=bvp4c(@ODEfun,@BCfun,solinit);x=[0:0.1:1];y=deval(sol,x);xP=[0:0.1:1.0];yP=[11.07431.16951.28691.4284...1.59651.79472.02742.30042.62143];plot(xP,yP,'o',x,y(1,:),'r-')legend('AnalyticalSolution','NumericalSolution',..
