[問題]矩陣 求出行列式值已回收

看板MATLAB作者 (icywind)時間17年前 (2008/12/15 23:23), 編輯推噓4(4014)
留言18則, 4人參與, 最新討論串1/1
有一個如下的矩陣 Y= [ 24*EI/h^3-166/35*w^2*m*h, 6*EI/h^2-11/210*w^2*m*h^2, 6*EI/h^2-11/210*w^2*m*h^2] [ 6*EI/h^2-11/210*w^2*m*h^2, 5*EI/h-17/105*w^2*m*h^2, 1/2*EI/h+4/35*w^2*m*h^3] [ 6*EI/h^2-11/210*w^2*m*h^2, 1/2*EI/h+4/35*w^2*m*h^3, 5*EI/h-17/105*w^2*m*h^3] 小弟我已經利用 syms EI h w m 現在我想計算這個矩陣的行列式值 並令其行列式為0 想要解出w=.... 這個w 希望是以 EI h m來表示 最後的答案有一個解是 w=1.5354*(EI/h^2)^(1/2) 我指令 [w]=solve(det(Y)) 可是 出來的答案 卻很長一串 跟我的解答有很大的出入 還是說有其他的方法可以解出跟我一樣的答案 不知道有沒有大大可以幫我排解一下這怪問題= = 小弟感激不盡! -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 61.228.246.63

12/15 23:32, , 1F
你如果輸入 w=1.5354*(EI/h^2)^(1/2) 好像數字會變成浮
12/15 23:32, 1F

12/15 23:32, , 2F
點 建議你 w=w(第幾個) 一次跑完 不要用手動的
12/15 23:32, 2F

12/15 23:36, , 3F
你矩陣一開始輸入可能就有問題~~給值要sym()
12/15 23:36, 3F

12/15 23:59, , 4F
請不要太依靠symbolic去解東西,簡單的還可以,複查的通
12/15 23:59, 4F

12/15 23:59, , 5F
常不是人能夠看的東西,就算去化簡最後還是要人去做處理
12/15 23:59, 5F

12/16 00:50, , 6F
alam大 那個w是正的值 其實小弟我不大懂你的意思
12/16 00:50, 6F

12/16 00:51, , 7F
我比較笨= = 抱歉 可以請你用其他方式解釋嗎 感恩
12/16 00:51, 7F

12/16 00:52, , 8F
Gwa大 那個結果真的不好整理 不知道你有沒有好的方
12/16 00:52, 8F

12/16 00:52, , 9F
法 可以供小弟參考 感謝您
12/16 00:52, 9F

12/16 01:01, , 10F
symbolic解完的可以試試再用pretty 會好看一點
12/16 01:01, 10F

12/16 01:09, , 11F
請問pretty的用法跟功用 拜託你解釋一下 感謝
12/16 01:09, 11F

12/16 02:05, , 12F
就跟字面一樣啊 "好看"一點...help有啊
12/16 02:05, 12F

12/16 22:05, , 13F
Y=sym([你的矩陣])
12/16 22:05, 13F

12/16 22:06, , 14F
w=solve(det(Y)) 如果第一個解釋你要的w1=w(1)
12/16 22:06, 14F

12/16 22:07, , 15F
建議你矩陣要先整理一下~~同一列要消掉一些東西
12/16 22:07, 15F

12/16 22:07, , 16F
我覺得你矩陣的因次怪怪的 沒有一致是正確的嗎?
12/16 22:07, 16F

12/17 02:23, , 17F
修正一下w=solve(det(Y),'w')
12/17 02:23, 17F

12/17 22:25, , 18F
矩陣的確是3*3 感恩 我去試試看 謝謝你
12/17 22:25, 18F
文章代碼(AID): #19HdQ0J1 (MATLAB)
文章代碼(AID): #19HdQ0J1 (MATLAB)