[問題] 微分方程組(ode45)已回收

看板MATLAB作者 (I'm OK...)時間15年前 (2010/07/09 10:35), 編輯推噓0(002)
留言2則, 2人參與, 最新討論串1/1
小弟目前正在使用ode45求解一組微分方程 遇到了一個小問題如下: dxdt=(x-a)*(afa/(afa*t-1)) dydt=(y-b)*(afa/(afa*t-1)) a, b為常數 afa=s/v,v也是常數 我有兩個case要求解 1. s=4.7713e-008 t=[0 261982.83 471569.1 785948.5 1047931.33 1205121.03]; 2. s=1.2479e-004 t=[0 100.169 180.305 300.508 400.678 460.779]; 我的問題是,雖然兩個case的值代入求解後,分母的afa*t是一樣的 但是分子的afa中的s分別是s=4.7713e-008與s=1.2479e-004 可是為什麼我求解出來的答案卻是相同的 單從公式上來看,應該是要解出不一樣的答案才對吧??? 麻煩matlab大大可以幫我想一下 謝謝!! 如有不清楚的地方我再補充,感恩~ -- Ben~ -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.138.138.238

07/11 20:52, , 1F
我比較好奇你怎麼從0開始算?? 分母不會暴掉?
07/11 20:52, 1F

07/12 08:41, , 2F
不好意思~我方程式打錯了~ 已改正!! 分母為afa*t-1
07/12 08:41, 2F
※ 編輯: ben1024 來自: 140.138.138.238 (07/12 08:42)
文章代碼(AID): #1CDehYj4 (MATLAB)