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

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

-*-------------------------- R-K法一般形式： ...

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

一阶线性微分方程求特解（附图）. ^letu= (x^3+1)ydu/dx = (x^3+1) dy/dx + 3x^2. y//y' +3x^2.y/(x^3+1) = y^2.(x^3+1). sinx(x^3+1)y' +3x^2.y = y^2.(x^3+1)^2. ...

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

1、如下公式： 其中x、y、m、n为参数，a、b为未知数，利用MATLAB求解方程； syms x y m n a b; [a,b]=solve('x=m*cos(a)+n*cos(a+b)','y=m*sin(a)+n*sin(a+b)','a','b'); 会得到四组解 ...

更新：9 APR 2016 ========方法======== 对于任意的二元二阶齐次线性偏微分方程， \(a_{11}\dfrac{\partial^2u}{\partial x^2}+2a_{12}\dfrac{\partial^2 u}{\partial x\partial y ...