应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 【讨论】预调件共轭梯度法（PCG）
54/1返回列表
查看: 3632  |  回复: 50
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 数学EPI ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xjw0413

铜虫 (初入文坛)

[交流] 【讨论】预调件共轭梯度法（PCG） 已有4人参与

有限元计算经常碰到大型稀疏矩阵，由于此类线性方程组通常条件数是比较大的，方程组的性态不好，所以最好用迭代方法求解，比方说是预调件共轭梯度法，但此方法在选择预调件矩阵时似乎没有一个同一的标准，大多推荐的是采用incomplete LU decomposition做为预调件矩阵。incomplete LU decomposition的计算方法似乎又有很多种。
1. incomplete LU decomposition 的计算时间应该比 LU decomposition要快速的多吧，不然直接用LU decomposition不就解出来了吗，又何必再来PCG迭代呢？
2. 采用PCG方法的前提应该是系数矩阵对称、正定吧，因为其原理是一个相当于势函数的东西取极小值。那对于非正定的系数矩阵能求解吗，我构造了几个非正定的，有的似乎是能够收敛到正确结果的。

希望各位虫用解答和讨论。
回复此楼

» 猜你喜欢
The life I want, there's no short cut.
1楼 2010-09-06 22:35:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

saladin983

铁杆木虫 (正式写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小雨萌萌: 金币+2, 谢谢关注~ 2012-09-03 12:58:34
引用回帖:
2405247楼: Originally posted by mastergxm at 2012-08-03 08:23:46
对于PCG方法有个疑问，不知道是我没有弄情楚，还是实际上就是如此。
81/2a/905677_1343974829_332.jpg
这是PCG的计算步骤。其中用到的P就是preconditioner。书中介绍说把A作ILU分解后，用LU的乘积矩阵作为这 ...

ILU分解得到的是上三角矩阵和下三角矩阵，应用P的时候，实际上是求解相关的两个方程组，由于是三角矩阵，计算量并不大。事实上，在应用中，几乎不会出现使用P的逆矩阵的情况，而是以求解相关线性方程组，亦即Py=r的方式实现。
赞一下(1人) 回复此楼
38楼2012-08-05 04:03:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 51 个回答

saladin983

铁杆木虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lovibond(金币+1): 鼓励交流 2011-06-18 10:34:45
1、还有改良条件数、数值稳定性的考虑
2、系数矩阵非正定的话，迭代过程中应该出现残量上升的情况，能够求解只能算是运气，个人以为实在迭代到负方向之前程序就已经终止。非正定的情况下通常用另一种Krylov子空间迭代法——GMRES求解。
赞一下(2人) 回复此楼
2楼2010-09-06 22:58:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xjw0413

铜虫 (初入文坛)
引用回帖:
Originally posted by saladin983 at 2010-09-06 22:58:01:
1、还有改良条件数、数值稳定性的考虑
2、系数矩阵非正定的话，迭代过程中应该出现残量上升的情况，能够求解只能算是运气，个人以为实在迭代到负方向之前程序就已经终止。非正定的情况下通常用另一种Krylov子空间 ...

“个人以为实在迭代到负方向之前程序就已经终止”，你的意思是说，设定判断标准，如果残量变大就终止程序并报错吗？
我试一试你推荐的那个方法。
谢谢你了！
赞一下 回复此楼
The life I want, there's no short cut.
3楼2010-09-07 14:04:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

saladin983

铁杆木虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lovibond(金币+1): 鼓励耐心交流 2011-06-18 10:35:05
引用回帖:
Originally posted by xjw0413 at 2010-09-07 08:04:33:

“个人以为实在迭代到负方向之前程序就已经终止”，你的意思是说，设定判断标准，如果残量变大就终止程序并报错吗？
我试一试你推荐的那个方法。
谢谢你了！

在特定精度下终止程序实际上就是在一定的子空间内取得了解的近似。我的意思是你得到看似正确的解可能是因为这个子空间还没有包含任何系数矩阵的负曲率方向，这种情况下能得到满足精度要求的一个近似解似乎也很自然。一旦侦测到一个负的方向，迭代序列很可能因此发散。具体的表现是什么，还是需要比较严谨的推导。总之这种情况下用CG得到的结果是没有说服力的，这个不是改变终止准则就能修正的。
赞一下(2人) 回复此楼
4楼2010-09-07 16:06:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 51 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 2rKfVQg2jQww 2026-11-06 5/250 2026-11-07 17:34 by Xc2FxHsHkBNc
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 5sJFDKp2kwuI 2026-11-06 4/200 2026-11-07 16:50 by Xc2FxHsHkBNc
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 4/200 2026-11-07 16:28 by Xc2FxHsHkBNc
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 16:17 by Xc2FxHsHkBNc
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +10 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 17:18 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 7/350 2026-11-06 16:43 by 5sJFDKp2kwuI
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 13:47 by BkA9Bi5cMhkD
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +8 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:36 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 14/700 2026-11-06 05:11 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
信息提示
关闭
请填处理意见
关闭