| 查看: 3056 | 回复: 10 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
wangxn06金虫 (正式写手)
|
[求助]
解线性代数方程组
|
||
|
问题描述：有一个线性方程组AX=b（A是方阵）（系数矩阵A的大小为200*200以上） 系数矩阵A 是对称矩阵，对角不占优 不同的自由项b所计算得到的剩余残差不同，即的到r=b-AX'的量级上有很大差别。 问题：1：同样的系数矩阵，不同的自由项，为什么得到的残差量级上会有很大的差别，这是不是由于系数矩阵的某些特点引起的——比如说矩阵A的条件数过大？ 2.是否有合适的方法可以使得残差达到比较合适的量级。 我用过全主元高斯消去法、列主元高斯消去法、并且在全主元高斯消去法的基础上，应用了解的迭代改进法、解对称正定线性方程组的楚列斯基法，但是这些都没有提高解的精确度（这些方法对某些自由向量b，解的结果精度都比较高，但是对某些自由向量b解的精度都不足，而这些自由向量b之间到底有什么差别也不好说）。Sample Text |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【分享】2011考研数学二大纲
已经有9人回复
- 【求助】从AB=I怎么推导出BA=I，这里A、B和I都是同阶方阵，并且I是单位阵
已经有58人回复
- 【求助】方程组最优解
已经有18人回复
- 【转帖】理解矩阵
已经有51人回复
saladin983
铁杆木虫 (正式写手)
- 数学EPI: 4
- 应助: 6 (幼儿园)
- 贵宾: 0.2
- 金币: 9197.9
- 红花: 3
- 帖子: 687
- 在线: 115.7小时
- 虫号: 448168
- 注册: 2007-11-01
- 专业: 计算数学与科学工程计算
|
如果使用迭代法的话，可以利用preconditioning来改善条件数，具体的需要从系数矩阵的特性来考虑，常用的比如iLU。 关于后面一个问题，我想我们可以打个比方考虑这样的情况：系数矩阵有很小的特征值lambda（接近于0），而右端向量非常接近于这个特征值对应的特征向量，那么求解之后的x约为这个特征向量乘以lambda的倒数（非常大），同时误差也会被lambda的倒数放大得很厉害。如果右端向量与lambda对应的特征向量垂直的话，就不会有这样的问题。 你的例子里是什么情况我不清楚，但是现实问题中右端向量往往不是随意的，尤其在求解偏微分方程的时候，右端向量常常包含一定的特征空间的信息。不知这样说对解释你遇到的问题是否有所启示？ |
8楼2012-01-12 00:14:37
saladin983
铁杆木虫 (正式写手)
- 数学EPI: 4
- 应助: 6 (幼儿园)
- 贵宾: 0.2
- 金币: 9197.9
- 红花: 3
- 帖子: 687
- 在线: 115.7小时
- 虫号: 448168
- 注册: 2007-11-01
- 专业: 计算数学与科学工程计算
2楼2012-01-10 00:21:47
wangxn06
金虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 576.7
- 散金: 35
- 帖子: 509
- 在线: 250.5小时
- 虫号: 1156835
- 注册: 2010-11-27
- 专业: 流体力学
3楼2012-01-10 10:24:06
math2000
铁杆木虫 (职业作家)
- 数学EPI: 2
- 应助: 239 (大学生)
- 金币: 5846.2
- 红花: 18
- 帖子: 4810
- 在线: 458.7小时
- 虫号: 235375
- 注册: 2006-04-01
- 专业: 概率论与随机分析
4楼2012-01-10 16:16:01
