2019年12月30日
摘要: 1.代码 %%超松弛迭代法(此方法适用于大型稀疏矩阵但不适合与病态方程的解 %%线性方程组M*X = b,M是方阵,X0是初始解向量,epsilon是控制精度,omiga是松弛因子 function OIM = Overrelaxation_iterative_method(M,b,X0,epsil 阅读全文
posted @ 2019-12-30 13:40 谷梁天 阅读(2493) 评论(0) 推荐(0) 编辑
摘要: 1.代码 %%共轭梯度法(用于求解正定对称方程组) %%线性方程组M*X = b,M是方阵,X0是初始解向量,epsilon是控制精度 function CGM = Conjugate_gradient_method(M,b,X0,epsilon) m = size(M);up = 1000;e = 阅读全文
posted @ 2019-12-30 13:36 谷梁天 阅读(1603) 评论(0) 推荐(0) 编辑
摘要: 1.代码 %%最速下降法(用于求解正定对称方程组) %%线性方程组M*X = b,M是方阵,X0是初始解向量,epsilon是控制精度 function TSDM = The_steepest_descent_method(M,b,X0,epsilon) m = size(M);up = 1000; 阅读全文
posted @ 2019-12-30 13:31 谷梁天 阅读(1491) 评论(0) 推荐(0) 编辑
摘要: 这段代码实现了埃特金加速迭代法和斯特芬森加速迭代法,我们以斯特粉森迭代为例 1.代码 %%注意,这里的fei不再是形如f(x)=0的形式而是x=fei(x)的形式,有些fei(x)不收敛,需要寻找,X0是初始值,method取值0和1代表上述两种方法 function AIM = Accelerat 阅读全文
posted @ 2019-12-30 13:14 谷梁天 阅读(1423) 评论(0) 推荐(0) 编辑
摘要: 二分法求根主要应用了区间套定理,这一算法实现简单且结果也迭代的较好,但对于复杂函数其结果不理想 1.代码 %%二分法求根 %%f为函数表达式,interval0为初始区间,epsilon为控制精度 function RD = Roots_dichotomy(f,interval0,epsilon) 阅读全文
posted @ 2019-12-30 13:03 谷梁天 阅读(811) 评论(0) 推荐(0) 编辑
摘要: 这段代码实现了牛顿切线法、简化牛顿法和牛顿下山法这三种方程求解法,由于输出结果较长,只以牛顿下山法为例写一段例题 1.代码 %%牛顿迭代法 %%method为-1时为牛顿切线法,method为0时为简化牛顿法,method为1时为牛顿下山法 %%f是表达式f(x) = 0,X0是初值,epsilon 阅读全文
posted @ 2019-12-30 12:54 谷梁天 阅读(1756) 评论(0) 推荐(0) 编辑
2019年12月28日
摘要: 1.插值函数 %%分段线性插值 function PLI = Piecewise_linear_interpolation(X,f,precision) [m,n] = size(X);a = min(X);b = max(X); X = sort(X); F = subs(f,X); for k 阅读全文
posted @ 2019-12-28 19:43 谷梁天 阅读(3145) 评论(0) 推荐(0) 编辑
摘要: 1.插值函数 %%n次插值多项式 %%X是插值节点,n是插值多项式次数,若已知函数表达式则attribute为0,未知函数表达式但已知函数值时为1 function IPn = Interpolation_polynomials_of_degree_n(X,Y,precision,attribute 阅读全文
posted @ 2019-12-28 19:33 谷梁天 阅读(889) 评论(0) 推荐(0) 编辑
摘要: 1.多项式插值函数 %%多项式插值 %%说明:precision为精度,越大则图像越精细,attribute是属性值,当未知函数表达式但已知函数值时为1,否则为0 function PI = Polynomial_interpolation(f,X,precision,attribute) X = 阅读全文
posted @ 2019-12-28 19:24 谷梁天 阅读(2089) 评论(0) 推荐(0) 编辑
摘要: 1.三次样条插值函数 %%三次样条插值 %%bc为boundary conditions(边界条件),当已知两端点的一阶导数值时为-1,当已知两端的二阶导数时为0,当函数为周期函数时为1 %%X为节点值,Y为函数表达式(attribute=0)或者具体值(attribute=1) function 阅读全文
posted @ 2019-12-28 19:08 谷梁天 阅读(2021) 评论(0) 推荐(0) 编辑