边界未知的 Bvp4c

Question

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

Answer 1

与链接 post 中相同的方法适用，如果您有一个函数 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 求解器。