PRECONDITIONING TECHNIQUES 预条件子技术

最新推荐文章于 2025-04-17 14:52:05 发布
最新推荐文章于 2025-04-17 14:52:05 发布
阅读量1.6k 收藏 2
点赞数
分类专栏： 读书笔记 文章标签： 稀疏线性系统
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/lifestxx/article/details/102534007
版权

寻找一个好的预条件子来解决给定的稀疏线性系统通常被认为是艺术与科学的结合。

预条件子基本上可以简单认为是通过变换求解一些条件数比较小的线性方程组来求原始问题。

1 Jacobi, SOR, SSOR预条件子

Ax=b 可视为

x_{k+1}=M^{-1}Nx_{k}+M^{-1}b

其中A=M-N

x_{k+1}=Gx_{k}+f 

则：f=M^{-1}b, G=M^{-1}N=M^{-1}(M-A)=I-M^{-1}A

 

 

预条件共轭梯度下降法PCG浅谈
m0_47096428的博客
04-26 2300
共轭梯度法是一种求解SPD系统线性方程组的迭代方法。它本来是一种直接法，但是通过迭代法求解后，配合复杂的预条件方法，反而更受欢迎。
用SSOR迭代法和预优共轭梯度法求解Ax=b方程组(附MATLAB代码)
waitingwinter的博客
01-07 4570
一、 SSOR 迭代法 SSOR迭代法的迭代格式为 对应的MATLAB代码为 function [x,k,error,time]=myssor(A,b,x,w,tol,max_it) %% SSOR迭代算法 %% To solve the equation Ax=b; %% Iterative formula : (D-wL)x_(k+1)=((1-w)D+wU)x_k+wb; (D-...
【数值分析×机器学习】使用CNN进行雅可比预条件子的生成（烦）
梦想破三的奔三狗
05-05 1万+
英文标题：Machine Learning-Aided Numerical Linear Algebra: Convolutional Neural Networks for the Efficient Preconditioner Generation 中文标题：机器学习辅助数值线性代数：用于高效预处理器生成的卷积神经网络 论文下载链接：ResearchGate 论文项目代码：GitHub@block-prediction Prescipt： 最近有点伤，不过还是及时调整过来，赶在五一结束把摘要赶..
稀疏矩阵的Krylov子空间方法（KSP）迭代求解中，使用Chebyshev多项式作为预条件子
最新发布
weixin_42849849的博客
04-17 828
在稀疏矩阵的Krylov子空间方法（KSP）迭代求解中，使用Chebyshev多项式作为预条件子（Chebyshev加速）可以显著加速收敛。，即能在所有多项式中最均匀地振荡。通过将矩阵的特征值区间映射到 ([-1, 1])，Chebyshev多项式可以最大化地抑制该区间外的分量，从而加速收敛。\alpha_k = \frac{2}{\lambda_{\max} + \lambda_{\min}}, \quad \text{或动态调整}此时，( \tilde{\lambda} \in [-1, 1] )。
预条件共轭梯度法求解线性方程组2
chd_lkl的博客
12-16 1201
!// 在病态矩阵问题上，共轭梯度法比部分主元的高斯消去法还要差 !// 这种情况下，可通过预条件得到缓解，主要是将问题转化为良态矩阵系统，再实施CG方法 !// 此代码的预条件因子M= ( D + wL ) * Inv(D) * ( D + wU )，为高斯-赛德尔预条件因子 !// A = L + D + U，L为A的下三角，U为A的上三角，D为A的对角线元素 Program Conjuga...
论文（梳理自用）
weixin_48408072的博客
08-27 330
来自图分析（graph analytics），对未知数重排序（reordering of the unknowns），以及优先级块（priority blocking）的诸多聚合技术（clustering techniques）通常可以用于生成高效的预条件子（参考文献[12]），但是高昂的计算成本与较差的可扩展性使得这些算法并不常用，结果就是为块。基于图分析的卷积模式生成工具不太可行。备注：邻接变量应该就是系统矩阵中相邻的两行或者两列，雅可比矩阵好像就是对角线和对角线邻近的两个副对角线上存在元素的矩阵。
Preconditioning:前提条件 A，ta Frob(IA*W)。-matlab开发
05-31
在MATLAB环境中，"Preconditioning"（预处理）是一种优化数值计算方法的技术，特别是用于线性代数问题，如求解线性方程组或矩阵特征值问题。这个技术的核心在于通过改变原始问题的形式来改善算法的效率和稳定性。在...
预条件迭代算法在时域有限元法分析波导中的应用
本文分析了几种常见的预条件技术，如 Jacobi 预条件、SSOR（Symmetric Successive Over-Relaxation）预条件和 ILU0（Incomplete LU factorization 的零填充版本）预条件。 预条件的主要目的是改善线性系统的条件数...
改进的SSOR预条件法加速电磁散射问题求解
本文《Shifted SSOR预条件技术在电磁波散射问题中的论文研究》由陈嘉琪和陈如山两位作者，来自南京理工大学通信工程系，探讨了在解决电磁散射问题中大型密集线性系统的高效算法。他们聚焦于利用电场积分方程...
Hypre简介3
weixin_46738358的博客
07-20 714
需要注意的是，预条件子和迭代法的结合使用可以在不同的应用和问题中具有不同的效果。选择适当的预条件子：根据问题的性质和线性方程组的特点，选择合适的预条件子。常见的预条件子包括对角预条件子、块对角预条件子、代数多重网格（AMG）预条件子等。这涉及对不同组合的预条件子和迭代法进行测试，比较收敛速度、迭代步数、解的精度等指标，并根据问题的特点选择最佳的组合。预条件子可以有多种类型，如对角预条件子、块对角预条件子、代数多重网格（AMG）预条件子等。这取决于线性方程组的性质、矩阵的结构以及预条件子和迭代法的选择。
线性方程组求解——预处理Preconditioning介绍
保持写作习惯，完成知识沉淀
09-11 1562
通俗地说, 预处理就是将难以求解的问题转化成等价的容易求解的新问题对于线性方程组而言, 预处理就是对(病态) 系数矩阵进行适当的线性变换,转换为一个(良态) 新矩阵, 从而达到改善迭代法收敛性的目的.预处理子选取基本准则：一个好的预处理子P 通常需满足下面两个要求:(1)具有更小的条件数和(或) 更好的特征值分布;P是A的一个很好的近似。(2) 线性方程组Pz = r 容易求解, 即预处理子P 的使用成本低廉.▶ 第一条是为了确保预处理后的线性方程组更容易求解, 即预处理子有效。
SSOR矩阵分解法的流程计算
abner361的博客
12-15 704
基本思想是将系数矩阵(有限元和有限差分中都是稀疏矩阵)A,进行分解为。可见SSOR的收敛性相较于CG较好。1：给定初始x,w,计算。
MATLAB中ilu函数用法
jk_101的博客
09-22 1071
ilu函数的功能是对矩阵进行不完全 LU 分解。
SSOR-PCG FORTRAN版CSR压缩存储
摄心神，致虚极，守静笃
06-08 1674
SSOR-PCG FORTRAN版CSR压缩存储 http://blog.sina.com.cn/s/articlelist_1493956393_0_1.html   水滴石穿的博客   最近因为项目需要，自己写了个迭代的SSOR-PCG程序，基于上三角CSR压缩存储（或者是下三角CSC），我认为这些东西应该没有必要让大家重复开发，特发表在次，希望
预处理共轭梯度法(PCG)
热门推荐
OORRANNGGE的博客
09-19 2万+
开始读研了，研究方向是SLAM。SLAM中很重要的一部分是位姿估计，需要用到一些优化算法对位姿进行优化。预处理共轭梯度法PCG(Preconditioned Conjugate Gradient)是一种SLAM中常用的位姿优化算法，从共轭梯度法（Conjugate Gradient）衍生而来，用于快速计算最优值。 直接法 假设有一个线性方程： Ax=bAx=b\textbf{A}\t...
MSL（湿气等级）2020-09-29
qq_42878026的博客
09-29 3185
MSL（湿气等级） 编辑 本词条缺少概述图，补充相关内容使词条更完整，还能快速升级，赶紧来编辑吧！ MSL是Moisture Sensitivity Level的缩写，是湿气敏感性等级的意思。 MSL的提出就是为了给湿度敏感性SMD元件的封装提供一种分类标准。 中文名 湿气敏感性等级 外文名 Moisture Sensitivity Level 类型 湿气等级 简称 MSL MSL：从而使不同类型的元件能够得到正确的封装、储藏和处理，避免在装配或修理过程中出现事故。...
matlab ichol,求解线性方程组 - 预条件共轭梯度法
weixin_29310605的博客
03-17 1032
通过为 pcg 提供用来计算 A*x 的函数句柄(而非系数矩阵 A)来求解线性方程组。使用 gallery 生成一个 20×20 正定三对角矩阵。上对角线和下对角线上的元素都是 1，而主对角线上的元素是从 20 递减到 1。预览该矩阵。n = 20;A = gallery('tridiag',ones(n-1,1),n:-1:1,ones(n-1,1));full(A)ans = 20×2020 ...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值