这个函数实现了固定步龙格-库塔求解程序显式和隐式方法(和可选的自适应步长控制)。
函数支持显式和隐式方法和嵌万博1manbetx入方法。可以添加任何龙格-库塔方法只是通过指定屠夫生动的场面。算法本身是通用的,相对紧凑。约34目前实现的方法。
MATLAB解决的颂歌都变步,甚至不提供一个选项来运行与固定步长。这是因为自适应步长可以解算器相比更快和更精确的固定步长。然而,有时人们有充分的理由选择一个固定的解决:
参数研究(比较不同的模型参数的模拟结果)
——有限差分计算雅克比的仿真结果(自适应步长控制可能会引入大量噪声)
——执行逐点计算,解算器输出和测量数据必须引用同一个时间向量
——有预分配的数组为仿真结果和固定计算时间
接口和选项解释评论。有两个例子:
例二解决了阻尼和驱动谐振子与不同的方法和步骤大小比较的性能。
example2解决了杜芬方程有两个不同的隐式动力学。
例子(谐振子):
t = linspace (0, 10);
x y = rkfs (@ (~) (x (2) - x (1)), t, [1 0]);
情节(t ' y ');
引用作为
Tillmann Stubler (2023)。龙格-库塔固定步骤解决(//www.tianjin-qmedu.com/matlabcentral/fileexchange/69951-runge-kutta-fixed-step-solvers), MATLAB中央文件交换。检索。