1.代码 %%高斯-塞得勒迭代法 %%线性方程组M*X = b，M是方阵，X0是初始解向量，epsilon是控制精度 function GSIM = Gauss_Seidel_iterative_method(M,b,X0,epsilon) [m,n] = size(M); d = diag(M); ...
分类：
其他好文 时间：
2019-12-30 14:51:22
阅读次数：
85
1.代码 %%列主元消去法 function ECPE = Elimination_of_column_pivot_entries(M,b) global n; [n,n] = size(M); B =[M,b]; R_A = rank(M);R_B = rank(B); if R_A ~= R_B ...
分类：
其他好文 时间：
2019-12-30 14:43:37
阅读次数：
95
1.代码 %%LU分解法 function LUDM = LU_Decomposition_method(A,b) global n;global B;global U;global L;global M; [n,n] = size(A); B = [A,b]; R_A = rank(A);R_B ...
分类：
其他好文 时间：
2019-12-30 14:38:05
阅读次数：
68
1.代码 %%最速下降法（用于求解正定对称方程组） %%线性方程组M*X = b，M是方阵，X0是初始解向量,epsilon是控制精度 function TSDM = The_steepest_descent_method(M,b,X0,epsilon) m = size(M);up = 1000; ...
分类：
其他好文 时间：
2019-12-30 14:35:44
阅读次数：
95
1.代码 %%超松弛迭代法（此方法适用于大型稀疏矩阵但不适合与病态方程的解 %%线性方程组M*X = b，M是方阵，X0是初始解向量，epsilon是控制精度，omiga是松弛因子 function OIM = Overrelaxation_iterative_method(M,b,X0,epsil ...
分类：
其他好文 时间：
2019-12-30 14:34:32
阅读次数：
103
1.代码 %%雅可比迭代法（此迭代法对于病态矩阵的解不理想） %%线性方程组M*X = b，M是方阵，X0是初始解向量,epsilon是控制精度 function JIM = Jacobian_iteration_method(M,b,X0,epsilon) [m,n] = size(M); d = ...
分类：
其他好文 时间：
2019-12-30 14:33:53
阅读次数：
73
1.代码 %%共轭梯度法（用于求解正定对称方程组） %%线性方程组M*X = b，M是方阵，X0是初始解向量,epsilon是控制精度 function CGM = Conjugate_gradient_method(M,b,X0,epsilon) m = size(M);up = 1000;e = ...
分类：
其他好文 时间：
2019-12-30 14:18:32
阅读次数：
142
Scipy的核心计算部分是一些Fortran数值计算库： 线性代数使用LAPACK库 快速傅立叶变换使用FFTPACK库 常微分方程求解使用ODEPACK库 非线性方程组求解以及最小值求解使用MINPACK 库 1. constants 模块 scipy的constants模块包含了众多的物理常数： ...
分类：
其他好文 时间：
2019-12-27 18:17:17
阅读次数：
149
1.numpy和scipy都提供了线性代数函数库linalg。但是scipy的线性代数库比numpy更加全面。 2.numpy中的求解线性方程组：numpy.linalg.solve(a, b)。而scipy中的求解线性方程组： scipy.linalg.solve(a, b, sym_pos=Fa ...
分类：
其他好文 时间：
2019-12-27 18:15:13
阅读次数：
88
1.scipy的optimize模块提供了许多数值优化算法。 2.求解非线性方程组： scipy.optimize.fsolve(func, x0, args=(), fprime=None, full_output=0, col_deriv=0, xtol=1.49012e-08, maxfev= ...
分类：
其他好文 时间：
2019-12-27 17:55:44
阅读次数：
80
