数值计算整理
python练习系列
菜伙子
这个作者很懒,什么都没留下…
展开
-
【数值计算】python实现SOR迭代法
松弛法介绍 在之前我们所介绍的Jacobi迭代法和G-S迭代法中,我们发现,其在参数选取上都是固定的,没有选择余地,因此无法对于收敛速度进行改进,由此,引出了松弛法的提出。松弛法的关键在于加上松弛因子www。下面我们给出对于一般方程组的迭代法表示:x(k+1)=x(k)−w(A(k)−b)x^{(k+1)}=x^{(k)}-w(A^{(k)}-b)x(k+1)=x(k)−w(A(k)−b) 这种迭代法称之为Richardson迭代法,其中w∈Rw\in\mathbb{R}w∈R称之为松弛因子。对于松弛原创 2020-12-30 17:36:35 · 2862 阅读 · 1 评论 -
【数值计算】python实现Gauss-Seidel迭代法
G-S迭代法简介 G-S迭代法是在Jacobi迭代法的基础上采用了Seidel加速技巧来实现的,于是在这里我们首先对于Seidel迭代法做一个简要的介绍。简单迭代法和Seidel加速技巧 上一篇Jacobi迭代法中其实是简单迭代法的变式,简单迭代法是所有迭代法的基础,在这里我们只简单介绍推理过程,具体证明过程可以参考1中的证明。 依旧考虑求解如下方程组AX=b(1.1)AX=b\tag{1.1}AX=b(1.1) 将其化为等价形式:x=Hx+gx=Hx+gx=Hx+g 其中H=(hij)原创 2020-12-23 17:07:21 · 1749 阅读 · 0 评论 -
【数值计算】python实现Jacobi迭代法
【数值计算】python实现Jacobi迭代法上学期刚刚学完数值计算,写一篇博客就当是练习markdown和用latex了。Jacobi算法介绍 考虑方程组 AX=b(1.1)AX=b\tag{1.1}AX=b(1.1) 其中A∈Rn×nA\in \mathbb{R}^{n \times n}A∈Rn×n非奇异,将矩阵AAA写成如下形式:A=D−L−U(1.2)A = D-L-U\tag{1.2}A=D−L−U(1.2) 其中D=diag{a11,......ann}D=diag\{{a_原创 2020-12-21 17:27:07 · 2494 阅读 · 1 评论