※ 引述《candy88257 (阿泰斯)》之銘言：
: ※ 引述《candy88257 (阿泰斯)》之銘言：
: : Solve[Det[A]==0,Z]
: : 把求到的Z代回A算Det怎會不為0?
: : 有時候Det出來的值會很小，這樣就無所謂
: : 但是...
: : 有時候代回算出來的值很大!!!
: : 怎會這樣?
: 補上矩陣照片:(巴哈圖片上傳空間)
: truth.bahamut.com.tw/s01/201211/497a65dcb4f762be0adc8a12e4d5cd6f.JPG
: 懇求大師幫忙看一下原文問題出在哪!
: 感謝!
(* 沒力打圖片的矩陣資料，產生一下隨機矩陣 *)
A = RandomReal[{-1, 1}, {8, 8}] +
RandomReal[{-1, 1}, {8, 8}]*RandomChoice[{0, 1, h, h^2}, {8, 8}];
(* 求解 *)
ans = h /. Solve[Det@A == 0, h];
(* 驗算 *)
Det[A /. h -> #] & /@ Cases[ans, _Real]
只是多項式，不大可能有太大誤差。如果不想有計算上的誤差，請改用分數計算。
--
養花種魚數月亮賞星星
http://chungyuandye.twbbs.org
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.173.134.247
※ 編輯: chungyuandye 來自: 218.173.134.247 (11/03 11:36)