bvp4c针对两点边界值问题进行数值求解,实现对以下n个一阶常微分方程的数值解:x的取值范围为a≤x≤b区间内,初始条件bvp4c调用格式functiondydx=FunctionName(x,y,p1,p2,...)其中x是一标量,y为yj列向量,p1,p2,...等为fj的参数,参数值已知,输出dydx是列向量,对应于fj列向量functionRes=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的输出。如有侵权请联系告知删除,感谢你们的配合!
第四讲 bvp4c函数使用说明 来自淘豆网m.daumloan.com转载请标明出处.
