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

这里以三元二次常微分方程组做一个例子，更多元更高次的都类似。 比如下列方程组： x'' = x' - x + y' -z' y'' = y' - y - x' z'' = z' - z + x' matlab代码如下： main.m： testfun.m ...

原文地址：MATLAB解常微分方程组的解法（好东西要共享）作者： 流若星辰 1：问题 常微分方程的初值问题的标准数学表述为：y'=f(t,y),a<=t<=b,y(a)=y(0) ；我们要求解的任何高阶常微分方程都可以用替换法化为上式所示的一阶形式，其中y ...

和线性常微分方程组参数拟合类似，我们要用差分代替微分，然后进行插值处理，然后构造最小化函数。 最后用最优化方法处理该函数即可。 这里举个例子，先随便设一个非线性微分方程组，并给定初值： 然后定义最小化函数： 最后用之前介绍的非线性最优化方法解决。 matlab代码 ...

比如我们已经有了微分方程模型和相关数据，如何求模型的参数。 这里以SEIR模型为例子，SEIR模型可以参考之前的文章。 一般的线性方程我们可以用最小二乘来解，一般的非线性方程我们可以用LM来解。 这里是线性微分方程组，所以我们采用最小二乘来解。 关键是构造出最小二乘形式，微分 ...

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

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

@ 目录 前言 一、常微分方程 二、常微分方程组 1.普通常微分方程组 2.线性常微分方程组 参考书目 前言 本文将介绍如何用matlab求解一阶常微分方程（组）的特解，通解。 如果你对微分方程的常见解法感兴趣 ...