Matlab求解边值问题方法例题.pptVIP

把待解的问题转化为标准边值问题 因为边值问题可以多解，所以需要为期望解指定一个初始猜测解。该猜测解网（Mesh）包括区间[a, b]内的一组网点（Mesh points）和网点上的解S(x) 根据原微分方程构造残差函数 利用原微分方程和边界条件，借助迭代不断产生新的S(x)，使残差不断减小，从而获得满足精度要求的解 Matlab求解边值问题方法：bvp4c函数 solinit=bvpinit(x,v,parameters) 生成bvp4c调用指令所必须的“解猜测网” sol=bvp4c(odefun,bcfun,solinit,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是表示初始网格有序节点的（1×M）一维数组，并且solinit.x(1)一定是a，solinit.x(end)一定是b。M不宜取得