Re: [問題] 解線性系統已回收
※ 引述《fsm (新的生活新的開始)》之銘言：
: 我在解 u''(x) = f(x)
: u'(0) = α, u'(1) = β 的時候
: 利用second difference method
: ┌ ┐
: 得到 │ -1 1 │
: │ 1 -2 1 │
: │ 1 -2 1 │x = B
: │ ....... │
: │ 1 -2 1 │
: │ 1 -1 │
: └ ┘
把原 GE 的解看成 (GE2) C．du/dt = u"(x)-f(x) 在時間無限大，
或者是收斂狀態下(steady state)的解。
請把方程式(GE2)離散化之後，小心排列並檢查一下數值，
讓主對角線上的數值變成正（最好在同一個 column 中是最大的數值 ）
: 但是 A 是 singular ， 也沒有辦法用A\B的方式去做
: 請問要怎麼解呢....
: 我試過利用(第 k 列) +(第 k+1 列) 的方式去做，最後一列會變成[0 .... 0]
: 所以會有無窮多組解，
: 但是如果不丟值讓他去跑結果出來的話，就不知道怎麼去求解了....
: 請問有人可以幫助我嗎..卡兩三天了.....
: 謝謝...
你可以上網查詢一下 tri-diagonal matrix 的相關資訊
應該可以找到簡單的線性代數求解方法，並不用作反矩陣的計算。
以上供你參考一下
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.123.122.179
※ 編輯: hilosi 來自: 140.123.122.179 (10/15 09:54)
※ 編輯: hilosi 來自: 140.123.122.179 (10/15 10:03)
推
10/15 13:13, , 1F
10/15 13:13, 1F
討論串 (同標題文章)
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
