💥💥💞💞欢迎来到本博客❤️❤️💥💥
🏆博主优势:🌞🌞🌞博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。
⛳️座右铭:行百里者,半于九十。
📋📋📋本文目录如下:🎁🎁🎁
目录
💥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.