假设你知道bvp4c的用法，仅说系数h处理方法的关键点。
公式解析：y1=u;y2=u'，则
设y1'=y2,则
y2'=(y1^2-h*y1)/λ
因此，在建子立函数时这样处理
function dydx = twoode(x,y)
H=[]; %H为已知的h序列，自己赋值
X=[]; %X为y=f(x)的定义域，即x的取值范围，注意与H同维，自己赋值
h=interp1(X,H,x,'cubic'); %即每个迭代中x对应的h值，用一维插值来做,插值方法可自己给定
dydx = [ y(2) ; (y(1)^2-h*y(1))/lamda];
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024