2019年12月30日
摘要：
这段代码可以给出常用的4类正交多项式的具体表达式，后续将给出可自定义的任意正交多项式代码 %%正交多项式 %%此函数包括勒让德正交多项式，切比雪夫正交多项式（两类），拉盖尔正交多项式，埃尔米特正交多项式，输入项数应从1开始 %%n是多项式的项数，n>=0,type是类型，分为Legendre、Che 阅读全文
摘要：
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 阅读全文
摘要：
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 阅读全文
摘要：
1.代码 %%雅可比迭代法（此迭代法对于病态矩阵的解不理想） %%线性方程组M*X = b，M是方阵，X0是初始解向量,epsilon是控制精度 function JIM = Jacobian_iteration_method(M,b,X0,epsilon) [m,n] = size(M); d = 阅读全文
摘要：
1.代码 %%高斯-塞得勒迭代法 %%线性方程组M*X = b，M是方阵，X0是初始解向量，epsilon是控制精度 function GSIM = Gauss_Seidel_iterative_method(M,b,X0,epsilon) [m,n] = size(M); d = diag(M); 阅读全文
摘要：
1.代码 %%超松弛迭代法（此方法适用于大型稀疏矩阵但不适合与病态方程的解 %%线性方程组M*X = b，M是方阵，X0是初始解向量，epsilon是控制精度，omiga是松弛因子 function OIM = Overrelaxation_iterative_method(M,b,X0,epsil 阅读全文
摘要：
1.代码 %%共轭梯度法（用于求解正定对称方程组） %%线性方程组M*X = b，M是方阵，X0是初始解向量,epsilon是控制精度 function CGM = Conjugate_gradient_method(M,b,X0,epsilon) m = size(M);up = 1000;e = 阅读全文
摘要：
1.代码 %%最速下降法（用于求解正定对称方程组） %%线性方程组M*X = b，M是方阵，X0是初始解向量,epsilon是控制精度 function TSDM = The_steepest_descent_method(M,b,X0,epsilon) m = size(M);up = 1000; 阅读全文
摘要：
这段代码实现了埃特金加速迭代法和斯特芬森加速迭代法，我们以斯特粉森迭代为例 1.代码 %%注意，这里的fei不再是形如f(x)=0的形式而是x=fei(x)的形式，有些fei(x)不收敛，需要寻找，X0是初始值，method取值0和1代表上述两种方法 function AIM = Accelerat 阅读全文
摘要：
二分法求根主要应用了区间套定理，这一算法实现简单且结果也迭代的较好，但对于复杂函数其结果不理想 1.代码 %%二分法求根 %%f为函数表达式，interval0为初始区间，epsilon为控制精度 function RD = Roots_dichotomy(f,interval0,epsilon) 阅读全文
摘要：
这段代码实现了牛顿切线法、简化牛顿法和牛顿下山法这三种方程求解法，由于输出结果较长，只以牛顿下山法为例写一段例题 1.代码 %%牛顿迭代法 %%method为-1时为牛顿切线法，method为0时为简化牛顿法，method为1时为牛顿下山法 %%f是表达式f(x) = 0，X0是初值，epsilon 阅读全文
