微分方程數值解法,知道導數的邊界值,而不是函式的邊界值,怎麼做

時間 2021-08-11 17:04:08

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的符號運算可以得到解析解。

有限元中一般使用伽遼金法(餘量加權法的一種),簡單的例子可參考

解微分方程組對工程問題來說十分重要,正不斷學習中。

微分方程已知特解求通解,知道非其次微分方程的兩個特解怎麼求通解

墨汁諾 非齊次線性微分方程的解,等於一個特解加上對應齊次方程的通解。y 3 就是那個特解。x n a1x n 1 a2x n 2 a n 1 x an 0 這就是線性方程。右端等於0,說明它是齊次方程 右端不等於0,說明它是非齊次方程。這是針對齊次方程 非齊次方程來說的。那麼微分方程類似,無非是左端...

下列微分方程中,不是全微分方程的是()

我的名字不太冷 這是高數下冊打星號的部分,很少有人知道,你仔細找一下吧,我找不到高數的書了,不能告訴你具體頁碼 吉祿學閣 若p x,y dx q x,y dy du x,y 則稱pdx qdy 0為全微分方程,顯然,這時該方程通解為u x,y c c是任意常數 根據二元函式的全微分求積定理 設開區域...

考研高等數學微分方程問題請問微分方程的

為什麼不可以?c只是代表常數,具體是正還是負需要由實際情況求出,c和c是一個意思,c也不能說明是負的。是的,因為c就是一個常數,無所謂正負。望採納 c就像是一元一次方程裡面的x一樣,常數而已 高等數學 微分方程 做有關微分方程的題 有時候後面加c又有時候加lnc1到底怎麼加,還是都可 如果解中是 l...