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

MATLAB常微分方程数值解 作者：凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 1.一阶常微分方程初值问题 2.欧拉法 3.改进的欧拉法 4.四阶龙格库塔方法 5.例题 用欧拉法，改进的欧拉法及4阶经典 ...

举例：分别用欧拉法和龙哥库塔法求解下面的微分方程 我们知道的欧拉法（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} \] 如何使用四阶龙格-库塔法求解该微分方程？ 一阶 ...

非刚性常微分方程的数值解法通常会用四阶龙格库塔算法，其matlab函数对应ode45。 对于dy/dx = f(x,y)，y(0)=y0。 其四阶龙格库塔公式如下： 对于通常计算，四阶已经够用，四阶以上函数f(x,y)计算工作量大大增加而精度提高较慢。 下面以龙格库塔法解洛伦兹方程为例 ...

[转] http://blog.sina.com.cn/s/blog_46e9b2010100tsqv.html 用matlab时间也不短了，可是一直没有接触过微分方程。这次看看书，学习学习，记点儿笔记。 1.可以解析求解的微分方程。 dsolve() 调用格式 ...

一、常微分方程的求解 例1、 例2、 例3、 通常我们使用syms 和dsolve来求解； first： second：表示 third：如果有必要 功能函数 diff可以完成 一元或多元函数任意阶数的微分： （对于自变量的个数多于一个的符号矩阵 ...

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