有時候我們要數值求解超越方程的多個根，但是數值方法都要給定一個初值。 matlab有內建函數fsolve求解非線性方程（組），但是只能求一組給定初值的解，如果要求多個根（如頻率方程），可以先用mathematica畫圖，然后觀察得到零點個數和大概位置，然后調用fsolve求解。 這里提供了一個 ...

詳細實驗指導見上一篇，此處只寫內容啦 實驗內容： 求解如下4元線性方程組的近似解。 • Jacobi迭代過程 運行結果： •Gauss-Seidel迭代 運行 ...

Doolittle分解復習 LDU分解 定理：矩陣A的各階順序主子式≠0，A可唯一分解為：A=LDU 其中，L單位下三角矩陣，U單位上三角矩陣，D非奇 ...

1.追趕法--先對矩陣A作Crount分解A=LU 例：用追趕法求解Ax=b 2.向量的范數（絕對者的和，絕對值的最大值，平方和開方） 1.絕對值的和： 2.絕對值的最大值： 3.平方和開方： 3.矩陣的范數（最大列 ...

一、由於任何偏微分方程初邊值問題通過半離散之后都會轉化為兩點邊值問題，以下僅以兩點邊值問題開始簡要討論Galerkin-Legendre 譜方法的數值實現過程 考慮如下分數階兩點邊值問題 \begin{align}&-(-\Delta)^{\frac{\alpha}{2}}u(x)=f ...

　　在乘法逆元里我們對於僅滿足b，m互質的情況，我們需要求解的是一個同余方程：b*x≡1（mod m），那么接下來我們就討論一下類似的線性同余方程的求解。 線性同余方程： 　　給定整數a，b，m，求一個整數滿足：a*x≡b（mod m），或給出無解。 　　因為未知數的次數為1，所以我 ...

簡介 van der Pol 方程 code ...

1）求解線性不定方程 　　ax + by = c 　　先求出一組解， 然后考慮如何表示通解， 設d = gcd(a, b), 假設c不是d的倍數， 則左邊是d的倍數而右邊不是， 則方程無解， 所以方程有解當且僅當d | c. 　　設c = c' * d, 我們先考慮方程 ax ...