這個我真不知道該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)