Re: [問題] 設有未知數的矩陣Det=0,代回驗算不為0?
看板Mathematica作者chungyuandye (養花種魚數月亮賞星星)時間12年前 (2012/11/03 11:34)推噓0(0推 0噓 0→)留言0則, 0人參與討論串3/3 (看更多)
※ 引述《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)
討論串 (同標題文章)
本文引述了以下文章的的內容:
完整討論串 (本文為第 3 之 3 篇):
Mathematica 近期熱門文章
PTT數位生活區 即時熱門文章