问答文章1 问答文章501 问答文章1001 问答文章1501 问答文章2001 问答文章2501 问答文章3001 问答文章3501 问答文章4001 问答文章4501 问答文章5001 问答文章5501 问答文章6001 问答文章6501 问答文章7001 问答文章7501 问答文章8001 问答文章8501 问答文章9001 问答文章9501
你好,欢迎来到懂视!登录注册
当前位置: 首页 - 正文

matlab 解大型稀疏线性方程组

发布网友 发布时间:2022-05-06 00:20

我来回答

1个回答

热心网友 时间:2022-06-29 13:53

这个问题不好直接回答,因为N=30000和稀疏性其实基本上不携带任何有价值的信息。
根据需求应该先搞清楚几件事,然后才能选择算法
1.是否具有对称性,如果有的话是否还有正定性
2.需要解一次还是解很多次(相同的A不同的b),或是矩阵还会变(不同的A)
3.对A的条件数或是部分特征值是否有已知信息或廉价的近似估计
4.A的非零元的分布特征
第1条通常决定了是否有相对廉价且高效的算法,2,3也是比较重要的参考信息,4则决定了稀疏直接法或者预条件的难度。

A\b采用的是稀疏直接法,是否可接受取决于第2条。
gmres和bicgstab都是针对非对称矩阵的Krylov子空间方法,默认参数可能不足以收敛,要用复杂的调用方式,自己指定参数,并且还要检查输出信息。这个你自己help一下用法。
一般来讲Krylov子空间迭代法需要预条件作为辅助,如果没有太多信息可以考虑用luinc来产生预条件。

楼上列的那篇paper不用看,基本没用。追答如果只解一次的话半分钟也算不了什么。已知条件只能用于估计A的非零元分布,没有过多信息。
GMRES的help看不懂说明你不知道GMRES的原理,如果经常要用的话最好先去学一下。
给你一个比较简单的例子
[L,U] = luinc(A,1e-3); %luinc(A,'0')也可以试一下,是一种完全不同的ILU
tol=1e-10; %残量的精度要求
restart=30; % 30-50之间吧,不要过小
maxit=100; %看情况,如果不收敛就适当调大
[x,flag]=gmres(A,b,restart,tol,maxit,L,U);
最后检查一下flag。

声明:本网页内容为用户发布,旨在传播知识,不代表本网认同其观点,若有侵权等问题请及时与本网联系,我们将在第一时间删除处理。
E-MAIL:11247931@qq.com
无水奶油对健康的影响是什么? 《米修斯之印》预约奖励介绍 蓝色星原旅谣在哪预约-游戏预约地址 阴阳师百闻牌taptap礼包怎么领_taptap礼包兑换码分享 猫之城抽卡怎么退出 猫之城测试资格怎么获得 内测名额申请教学 猫之城预约奖励有什么 预约奖励详细攻略 无人机怕什么干扰?解析无人机干扰的类型及应对措施 什么叫干扰三要素 观察动物日记越多越好 Krylov子空间迭代法二十世界比较经典的算法,但是还是不了解 管理出租车的是什么单位 Krylov子空间迭代法 巡游出租汽车什么情况是擅自转让 薛凯琪出了多少专辑 道路运输条例和巡游出租车管理规定可能并用吗 请问薛凯琪有哪些歌啊??? 巡游出租汽车经营者应当合理确定 方大同给薛凯琪写的所有歌 什么是出租车规费 进公司推出的基于微信生态环境下的新媒体平台有哪些? 薛凯琪有哪些歌好听, 求薛凯琪的复刻回忆mp3 薛凯琪有哪些歌比较好听?比较多人听的 巡游出租汽车车辆经营权的经营协议应当包括 推荐薛凯琪好听的歌 国家鼓励巡游出租汽车实行什么 薛凯琪的专辑有哪些? 微信小程序是一个全新的生态吗? 薛凯琪唱的好听的歌都有哪些? barzilai61borwein步长 由什么决定 如何解磁 Lanczos算法的概述 矩阵特征值问题。 系数矩阵的主对角元素为0,可以用迭代法求方程组吗 如题 证明:如果系数矩阵A至多有L个互不相同的特征值,则共轭梯度法至多L步就可得到方程组Ax=b的精确解 系数矩阵的主对角元素为0,可以用迭代法求方程组吗如题 严格对角占优矩阵换行后能构造迭代法么 教资考试小学语文学科教案怎么写 2019教师资格证小学语文《荷花》教案 教师资格证考试教学设计题怎么写 小学语文教师资格证面试的教案怎么写 小学教师资格证考试语文教案怎么写 教师资格证考试小学数学方向的课做语文的教学设计? 互动营销是什么 致即将步入大二的自己 作文800字 致自己成年的一篇完整的文章 什么是网络互动营销? 什么是互动营销?外部营销?内部营销? 致自己十八岁成年说说有哪些?
  • 焦点

最新推荐

猜你喜欢

热门推荐