[問題] dsolve 與 ode 解法器結果不同 已解決已回收
各位前輩大家好,
目前在解一階微分方程遇到問題,
用dsolve算出來後再帶數值,
與用ode解法器算出來差異很大,
不知道是哪邊出錯?
還是說原本就會這樣?
附上程式碼:
dsolve:
T=dsolve('Dy=0.4*y*2010*sin(t)/(402+2010*(1+cos(t)))','y(0)=300');
ezplot(T)
ode解法器:
function CaseII20100318
global a
a=0.4;
[t x] = ode15s(@CaseII2010,[0 500],300);
plot(t,x)
低階的以ode23.ode45解即可
function f = CaseII2010(t,x)
global a
f = a.*x.*2010*sin(t)/(402+2010*(1+cos(t)));
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.120.25.178
※ 編輯: YangJimmy 來自: 140.120.25.178 (03/18 16:37)