% [z, info] = BulirschStoer(dynFun,t,z0,tol) % % 使用 Bulirsch-Stoer 方法求解初始值问题。 这% 方法非常适合用于平滑初始值的高精度解％ 问题。 % % 计算 z(t) 使得 dz/dt = dynFun(t,z)，从初始值开始% 状态 z0。 网格点处的解将精确到 tol 之内。 % % 如果提供的网格不足，此功能会自动% 引入中间网格点以达到所需的精度。 % % 输入： % dynFun = 系统动力学的函数句柄％dz = dynFun（t，z） % t = 标量时间% z = [nz,1] = 状态为列向量% dz = [nz,1] = 状态的导数作为列向量% t = [1,nt] = 时间网格点向量% z0 = [nz,1] = 初始状态向量% tol = [nz,1] = 每个维度的误差容限。 如果 tol 是一个%
2022-02-13 21:42:21
4KB
1