这是我在学习飞行器制导与控制时的课程作业。用四阶龙格库塔法解微分方程组。我一开始的想法是分别利用龙格库塔法解每一个微分方程，但变量很多，算法会比较复杂。后来明白可以把多变量看作是一个变量，利用matlab的feval函数进行代入变量的函数运算。 matlab中feval函数的作用：feval(f ...

举例：分别用欧拉法和龙哥库塔法求解下面的微分方程 我们知道的欧拉法（Euler）"思想是用先前的差商近似代替倒数"，直白一些的编程说法即：f(i+1)=f(i)+h*f(x,y)其中h是设定的迭代步长，若精度要求不高，一般可取0.01。在定义区间内迭代求解即可。龙哥库塔法一般用于高精度 ...

引言 考虑存在以下二阶偏微分方程 \[\begin{align} f_2 \cdot \ddot{X(t)}+f_1 \cdot \dot{X(t)} +f_0 \cdot {X(t)} =F(t) \end{align} \] 如何使用四阶龙格-库塔法求解该微分方程？ 一阶 ...

欧拉法解微分方程 本文介绍如何使用简单的欧拉法求解微分方程，大部分内容出自吴一东老师在他的B站个人空间发布的课程 方法介绍 对于一个一般的微分方程： \[\begin{cases} \begin{aligned} \frac{\mathrm{d} y}{\mathrm{d ...

欧拉法、改进的欧拉法、龙格-库塔法求解初值问题 简介 通过求解简单的初值问题： \[\begin{cases}\dfrac{du}{dx}=f(x,u)&&&&&&(1)\\u(x_0)=u_0&&&& ...

实例： u'=-3u+6x+5 u(0)=3 解析解：u=2e^(-3x)+2x+1 欧拉法 改进的欧拉法 ode45求解 总体 ...

任何其他的语言实现。欧拉（Euler）和中心差分逼近，是最朴素的想法，可惜代数精度太低了，而龙格库塔的稳 ...

老师说系统给的ode45好多都解决不了。 1.lorenz系统 test.m runge_kutta.m test_fun.m ...