边界未知的Bvp4c

Question

我有一个描述悬链垂垂的BVP，其中右边界是未知的（8 DE和9 BC）。 有人知道如何使用bvp4c函数在Matlab中对此建模吗？ 谢谢。

Answer 1

采用与链接文章中相同的方法，如果您具有函数right_side(t,x)和boundary_conditions(xa,xb)，则需要包装器函数来处理间隔长度的其他分量，

function doty = wrapper_rhs(s,y)
  T = y(end)
  doty = T*right_side(t0+s*T, y(1:end-1))
  doty(end+1)=0
end

和

function bc = wrapper_bc(ya,yb)
  bc = boundary_conditions(ya(1:end-1), yb(1:end-1))
end

这些应该允许在固定间隔“ [0,1]”内调用BVP求解器。