[問題] 方程式的求解?已回收

看板MATLAB作者 (小熊)時間15年前 (2010/10/03 12:01), 編輯推噓2(209)
留言11則, 2人參與, 最新討論串1/1
cha = (M*s^2 + c1*s + k1)*(m*s^2 + (c1 + c2)*s + k1 + k2) - (k1 + c1*s)^2 expands = (M*m)*s^4 + (c1*m + M*(c1 + c2))*s^3 + (M*k1 + M*k2 + c1*c2 + k1*m)*s^2 + (c1*k2 + c2*k1)*s + k1*k2 >> solve(cha,s) ans = RootOf(M*X1^4*m + X1^3*c1*m + M*X1^3*c2 + M*X1^3*c1 + X1^2*k1*m + X1^2*c1*c2 + M*X1^2*k2 + M*X1^2*k1 + X1*c2*k1 + X1*c1*k2 + k1*k2, X1) cha是我的方程式 展開後可以得到 以s為變數的四次方的多項式 想求解多項式的根 照理講會有4個解 但是我的ans 卻怪怪的 是...哪邊出問題了?? -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 114.37.102.67

10/03 12:40, , 1F
四次的解本來就很難看了@@
10/03 12:40, 1F

10/03 12:42, , 2F
先用一個簡單方程式x4+ax3+bx2+cx+d=0試試看
10/03 12:42, 2F

10/03 12:56, , 3F
他的ans的結果就是我的變數s變X1 感覺很奇怪
10/03 12:56, 3F

10/03 12:57, , 4F
應該沒有說變數不能用s吧= =?
10/03 12:57, 4F

10/03 13:55, , 5F
結果是
10/03 13:55, 5F

10/03 13:55, , 6F
RootOf(X3^4*a + X3^3*b + X3^2*c + X3*d + e, X3)
10/03 13:55, 6F

10/03 23:06, , 7F
你是不是沒有加 syms X a b c d e
10/03 23:06, 7F

10/03 23:07, , 8F
要先命令他是符號!
10/03 23:07, 8F

10/03 23:45, , 9F
有呀~我有事過 ax^2+bx+c=0沒問題(只是沒把完整的列出)
10/03 23:45, 9F

10/04 00:04, , 10F
a大有試過嗎? 我2次跟3次結果都跑得出來
10/04 00:04, 10F

10/04 00:05, , 11F
但是四次就怪怪的。
10/04 00:05, 11F
文章代碼(AID): #1Cg00hG5 (MATLAB)