【使用积分器来模拟物体的运动】使用Matlab的ODE45积分器和标准的Runge-Kutta 4积分器研究（Matlab代码实现）

 💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

本文将演示如何利用积分器来模拟物体的运动。我们将重点介绍Matlab中的两种积分器：ODE45积分器和标准的Runge-Kutta 4积分器。通过本文的指导，您将学会如何在Matlab中使用这两种积分器来模拟物体的运动，为您的研究和实践提供有力的支持。

首先，我们将深入探讨ODE45积分器的使用方法。ODE45是Matlab中常用的一种积分器，它基于龙格-库塔方法，能够高效地求解常微分方程组。我们将演示如何利用ODE45积分器来模拟物体的运动轨迹，并详细介绍其参数设置和使用技巧。

其次，我们将介绍标准的Runge-Kutta 4积分器。这是一种经典的数值积分方法，通过迭代计算来逼近微分方程的解。我们将讨论如何在Matlab中实现Runge-Kutta 4积分器，并与ODE45积分器进行比较，以便更好地理解它们各自的特点和适用范围。

通过本文的学习，将掌握在Matlab中使用ODE45积分器和标准的Runge-Kutta 4积分器来模拟物体运动的技能，更加熟练地运用积分器来研究物体的运动。

📚2 运行结果

部分代码：

% Set up data for plotting - then plot based on options set.
dataP           = dataStore(dataP,0);

% Set plotting data structure with ODE45 data.
dataF.intType       = intType;
dataF.varTimeStep   = varTimeStep;
if intType == 1
    dataF.t         = t;
    dataF.stateVec  = stateVec;
% Set plotting data structure with Runga-Kutta 4th order data.
else
    dataF.t         = timeSpan';
    dataF.stateVec  = accumStateVec;
end

% Plot the entire trajectory set.
plotData(dataInit,dataF,dataP,appType);
%--------------------------------------------------------------------------

%--------------------------------------------------------------------------
% Save the data for comparison of the ODE45 and RK4 results.  The function
% compPerf.m is subsequently used to process the .mat files to overlay the 
% results for comparison purpose.

% Matlab ODE45 integration results
if intType == 1
    dataODE45.t         = t;
    dataODE45.stateVec  = stateVec;
    save('matode45.mat','dataODE45','appType');
% Runga-Kutta 4th order integration results
elseif intType == 2
    dataRK4.t           = timeSpan';
    dataRK4.stateVec    = accumStateVec;
    save('matrk4.mat','dataRK4','appType');
end

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

[1]叶远波,黄太贵,吴保文,等.基于对角隐式Runge-Kutta法的新型数字积分器研究[J].电测与仪表, 2019, 56(8):7.DOI:10.19753/j.issn1001-1390.2019.08.004.

[2] Calvo M , Montijano J I , Randez L .Algorithm 968: DISODE45: A Matlab Runge-Kutta Solver for Piecewise Smooth IVPs of Filippov Type[J].ACM transactions on mathematical software, 2017, 43(3):25.1-25.14.

🌈4 Matlab代码实现