1樓:夜宵賭神
這個方程最簡單的方法是用shooting method.可wiki 之。
主要就假設一個y(a1), 然後看,y'(a2)是否等於b2。
可以用matlab寫兩個小函式,一個用在ode45(fun1),然後y(a1)的值可以用fzero(fun2).
fun1就是把二階微分化成一階方程組,fun2就是呼叫ode45並返回y'(a2)-b2. 這樣會自動求出初值的。
別外這個方程可能有解析解。可以利用格林函式來解,得到一個積分形式的解。
給你看一個我以前寫的例子吧。適當改下就行了。
function out=fdif(eta,f)
out=zeros(size(f));
out(1)=f(2);
out(2)=f(3);
out(3)=f(2)^2-1-f(1)*f(3);
function result=soldif(x)
[etaout,fout]=ode45('fdif',[0,10],[0,0,x]);
result=fout(end,2)-1;
2樓:哆嗒數學網
這個是線性非其次方程,解出來的通解一定是 齊次式通解加特解形式。
特解,可以用你的數值解法來找一個。
這樣就能得到原方程的通解。通解中有兩個待定常數。
把y'(a1)=b1,y'(a2)=b2代入進去,得到兩個方程。可以把兩個常數解出來
3樓:匿名使用者
使用matlab用龍格-庫塔解法求數值解,可以呼叫函式,也可以利用simulink做。
簡單的用matlab的符號運算可以得到解析解。
有限元中一般使用伽遼金法(餘量加權法的一種),簡單的例子可參考
解微分方程組對工程問題來說十分重要,正不斷學習中。
微分方程的通解求法,微分方程的通解怎麼求
二階常係數齊次線性微分方程解法 特徵根法是解常係數齊次線性微分方程的一種通用方法。設特徵方程r r p r q 0兩根為r1,r2。1 若實根r1不等於r2 y c1 e r1x c2 e r2x 2 若實根r1 r2 y c1 c2x e r1x 3 若有一對共軛復根 略 關於一階微分方程 齊次方...
什麼是微分方程初始條件,邊界條件,定解
微分方程在給出的bai時候是du有給初始的未知數的zhi值,這個叫做初始條dao件,微分方程時其通解都版包權含有未知常數c,c可以是任意常數。這些未知常數是由微分方程的定解條件確定的,也就是把初始條件帶入之後得出的c的值。微分方程的最後的解既滿足微分方程又滿足定解條件。微分方程的定解條件分為兩類 一...
關於偏微分方程的求解精通數值計算的進
雖然不知樓上從抄 抄的答 bai案,但是回答基本是正確 du的。初始條zhi 件和邊界條件構成dao了偏微分方程的定解條件。舉個例子說,連續方程和n s方程適用於描述一般的流動問題了,而邊界條件就決定你要解決的具體問題了,如果是非穩態問題就還要有初始條件。另外描述你問題的整個方程組中有多少個未知變數...