[問題] 高斯消去法的誤差問題已回收

看板MATLAB作者 (維尼)時間15年前 (2010/09/28 02:37), 編輯推噓1(100)
留言1則, 1人參與, 最新討論串1/1
我有一個矩陣A=[0.4096 0.1234 0.3678 0.2943 0.4043; 0.2246 0.3872 0.4015 0.1129 0.1550; 0.3645 0.1920 0.3781 0.0643 0.4240; 0.1784 0.4002 0.2786 0.3927 0.2557] 因為要寫個使用pivot的程式去做高斯消去 我找定了第一列作為主軸但問題是當我用0.2246/0.4096再乘與0.4096發現結果變成 0.2246000000000002 於是乎減0.2246後結果近乎為零,但不為零,再算下去代只會很大條 這應該是numerical error吧,但要怎麼讓消去後的矩陣某些元素結果為0 請熟悉的高手能解答一下,謝謝 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 64.121.116.115 ※ 編輯: jazz123457 來自: 64.121.116.115 (09/28 02:43)

09/29 16:36, , 1F
預設高斯消去法rref(X,Tol),可參考預設是怎麼寫的
09/29 16:36, 1F
文章代碼(AID): #1CeEHwFP (MATLAB)
文章代碼(AID): #1CeEHwFP (MATLAB)