原创 解線性方程組的直接方法 ------ 三角分解法

矩陣分解 三角分解(LU分解) - billbliss的專欄 - CSDN博客  https://blog.csdn.net/billbliss/article/details/78559289 矩陣的三角分解其基本原理可見上面的博客鏈接

原创 解線性方程組的直接方法 --- 選主元的三角分解法(C語言實現)

  從直接三角分解法公式可以看出,當U矩陣中的Urr=0時,計算將中斷。或者當Urr絕對值很小時,按分解公式計算可能引起舍入誤差的累積。 當A滿秩的時候(非奇異), 我們可以通過交換A的行實現矩陣PA的LU分解,因此可以採用與列主元消去

原创 Gauss-Jordan 列主元消去法來實現矩陣的求逆(Python)-非增廣

如上所示,通常情況下的矩陣求逆,採用的是增廣矩陣的方式,然後通過不斷進行行變換來使得前半部分變成單位陣,後半部分將會是矩陣的逆。 爲了降低計算機編程過程中的內存單元的消耗,不必將單位矩陣存放起來,其具體原理可見於: gauss_jord

原创 c 採用列主元的高斯消去法來實現方程組的求解

因爲工作的需要,未來的學習和使用將多數集中在c的編寫上, Python的使用將會較少的使用。 本次主要是來練習C語言的編寫,並實現了列主元的高斯消去法來實現方程組的求解。   支持龍豬一分錢。

原创 解線性方程組的直接方法(c語言實現)

在自然科學和工程技術中很多問題的解決常常歸結爲解線性代數方程組,例如採用差分法或者有限元方法解常微分方程,偏微分方程編制問題等都會導致求解線性代數方程組。 關於線性方程組的數值解法一般分爲兩類:直接法和迭代法。 本博客主要來介紹直接法 直

原创 Python進行有限元編程-平面應力問題(三節點三角形單元)

參考書籍是:《有限元方法基礎教程》(國際單位制)(第五版) 章節爲:第6章 建立平面應力和平面應變剛度方程 重要的事情繼續強調(有限元的基本計算流程): Step 1: 選擇單元類型。 Step 2: 選擇位移函數。 Step 3: 定義

原创 Python 進行有限元分析 -- 有限元方法基礎教程 (第五版)(例題3.6)

在《有限元方法基礎教程》(第五版)中,第三章是桁架問題。 這裏再重新強調一遍有限元分析的基本步驟(重要的事情我們要經常說): step 1: 離散和選擇單元類型。 step 2: 選擇位移函數。(位移函數一般都是多項式函數,一般都選線性