這是我在學習飛行器制導與控制時的課程作業。用四階龍格庫塔法解微分方程組。我一開始的想法是分別利用龍格庫塔法解每一個微分方程，但變量很多，算法會比較復雜。后來明白可以把多變量看作是一個變量，利用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 ...