作者birdhackor (夜残狼)
看板Grad-ProbAsk
标题[理工] [数学or程设]Mathematica
时间Mon Dec 7 04:04:29 2009
这个我真不知道该PO哪边...
因为有MatLab版但是没有mathematica版...
PO到这的话...标题分类又不知道该怎麽打....
如果不能PO在这的话麻烦站内信跟我讲一下我会自D
这是我在用电脑解力学时遇到的问题(因为微方一长串实在不想用手解)
解的是阻尼震荡
我想问的是
解2阶常系数ODE的时候
当特徵方程式为负数根a+-bi时一般表示为exp(ax)*(C1*cos[bx]+C2*sin[bx])
可是如果是把有符号的微方拿去解的话
mathematica只会解出通解
也就是比如说输入
DSolve[a*x''[t] + b*x'[t] + c*x[t] == 0, x[t], t]
输出会是
{x[t] -> E^(((-b - Sqrt[b^2 - 4 a c]) t)/(2 a)) C[1] +
E^(((-b + Sqrt[b^2 - 4 a c]) t)/(2 a)) C[2]}
就是他会直接把特徵方程式带入一元二次方程式公式解
然後写通解输出
可是假设我知道b^2 - 4 a c < 0 希望他输出成exp(ax)*(C1*cos[bx]+C2*sin[bx])
这样的形式
有办法吗?我一直找不到有没有办法可以输入常数之间的条件
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 122.121.204.117
※ 编辑: birdhackor 来自: 122.121.204.117 (12/07 04:06)
※ 编辑: birdhackor 来自: 122.121.204.117 (12/07 04:36)