Re: [問題] 解線性系統已回收

看板MATLAB作者hilosi時間16年前 (2008/10/15 17:57)推噓0(0推 0噓 0→)

留言0則, 0人參與討論串4/4 (看更多)

※ 引述《fsm (新的生活新的開始)》之銘言： : ※ 引述《hilosi ()》之銘言： : : 把原 GE 的解看成 (GE2) C．du/dt = u"(x)-f(x) 在時間無限大， : ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ : 不好意思，這邊我不太懂.. : : 或者是收斂狀態下(steady state)的解。 : : 請把方程式(GE2)離散化之後，小心排列並檢查一下數值， : : 讓主對角線上的數值變成正（最好在同一個 column 中是最大的數值） : : 你可以上網查詢一下 tri-diagonal matrix 的相關資訊 : 我有大致查了一下網路的資料 : 但是不是解不出來，就是爆掉(如果是高斯消去法的方式) 你會暴掉是因為 ┌ ┐ │ -1 1 │ │ 1 -2 1 │ │ 1 -2 1 │ 這個矩陣當中的主對角線上元素的數值的影響 │ ....... │ │ 1 -2 1 │ │ 1 -1 │ └ ┘ 你用的數值方法將 u''(x) = f(x) 拆解的話，應該怎麼拆都會得到你先前文章內容的形式你可以將原本要解的方程式 u''(x) = f(x) 等效為 C．du/dt=u"(x)-f(x) 此方程式在穩態時的解，那麼你在離散化之後得到的矩陣形式將變成很穩定的形式，求解就不會爆掉了你可以找幾本有關數值方法的書看看，就可以得到相當詳細的說明囉。 : : 應該可以找到簡單的線性代數求解方法，並不用作反矩陣的計算。 : : 以上供你參考一下 : 我有試過 x = pinv(A)*B的方式，雖然誤差變小許多 : 但是以數值解來看還是差了很多 : 請問這類的問題，要如何去解呢... google 找一下 TDMA 的相關資訊.... -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.123.122.179