第四讲bvp4c函数应用重点分析.ppt

第四讲 bvp4c函数的应用 bvp4c针对两点边界值问题进行数值求解，实现对以下n个一阶常微分方程的数值解： x的取值范围为a≤x≤b区间内，初始条件 bvp4c调用格式 FunctionName函数的接口 function dydx=FunctionName(x,y,p1,p2,...) 其中x是一标量,y为 yj列向量, p1,p2,...等为fj的参数，参数值已知，输出dydx是列向量，对应于fj列向量 BCFunction函数的创建 function Res=BCFunction (ya,yb,p1,p2,...) ya是表示 yj(a)的列向量,yb是表示 yj(b)的列向量，即使不要求边界条件，已知参数p1,p2等也必须出现在接口定义语句中，输出Res为列向量。 调用bvp4c函数所涉及到的另外两个函数 变量solinit为一结构？，可由函数bvpinit得到：solinit=bvpinit(x,y)，x向量为初始网格点的估计值，向量y为每个yj的估计值。向量x和y的长度互不相关。 bvp4c的输出sol是一结构，为特定数量点对应的解。为使曲线变得更光滑，需要在中间插入一些点，使用：sxint=deval(sol,xint),xint为点向量，函数deval根据这些点向量求解。sol为函数bvp4c的输