※ 引述《chickenribs (雞排)》之銘言：
: d2x/dt2 = a*(b-dx/dt)*sqrt((b-dx/dt)^2+(-dy/dt)^2)
: d2y/dt2 = a*(-dy/dt)*sqrt((b-dx/dt)^2+(-dy/dt)^2)+c
: a為y和dx/dt和dy/dt的函數
: b為y的函數
: c為常數
: 借問這樣的形式ode的指令要怎麼寫呢.....
: 我是新手>"<只會寫單一一階ode
: 感謝大家了~!
你需要設立新的state variable因為matlab只能解一階的ode
let new state variable p=dx/dt, q=dy/dt
so dp/dt= d2x/dt2= a(b-p)*sqrt((b-p)^2+(-q))
dq/dt= d2y/dt2= a(-q)*sqrt((b-p)+(-q)^2)+c
dx/dt=p
dy/dt=q
solve this 4-order state space equation!!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.167.194.197
※ 編輯: weili419 來自: 218.167.194.197 (05/23 01:06)