搜索
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的输出。(,1,2,,,)j=1,2,3njjndyfxyyyqdx()(),1,2,,jjjjjjyaaybbjnab及和为常数124(@,@BC,,,,,)solbvpcFunctionNameFunctionsolinitoptionspp