作者weili419 (fem)
看板MATLAB
标题Re: [问题] 请问联立解两个二阶常微分方程式
时间Tue May 23 00:54:25 2006
※ 引述《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)