24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2494  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wangxn06

金虫 (正式写手)

[求助] 解线性代数方程组

问题描述:有一个线性方程组AX=b(A是方阵)(系数矩阵A的大小为200*200以上)
                  系数矩阵A 是对称矩阵,对角不占优
                  不同的自由项b所计算得到的剩余残差不同,即的到r=b-AX'的量级上有很大差别。
问题:1:同样的系数矩阵,不同的自由项,为什么得到的残差量级上会有很大的差别,这是不是由于系数矩阵的某些特点引起的——比如说矩阵A的条件数过大?
      2.是否有合适的方法可以使得残差达到比较合适的量级。
         我用过全主元高斯消去法、列主元高斯消去法、并且在全主元高斯消去法的基础上,应用了解的迭代改进法、解对称正定线性方程组的楚列斯基法,但是这些都没有提高解的精确度(这些方法对某些自由向量b,解的结果精度都比较高,但是对某些自由向量b解的精度都不足,而这些自由向量b之间到底有什么差别也不好说)。Sample Text
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangxn06

金虫 (正式写手)

引用回帖:
: Originally posted by saladin983 at 2012-01-12 00:14:37:
如果使用迭代法的话,可以利用preconditioning来改善条件数,具体的需要从系数矩阵的特性来考虑,常用的比如iLU。

关于后面一个问题,我想我们可以打个比方考虑这样的情况:系数矩阵有很小的特征值lambda(接 ...

我试着说说你的意思:
x1为系数矩阵的特征向量,α为x1对应的特征值。 即有:Ax1=αx1
同时有 Ax=b
如果b和x1接近,可以近似有Ab=αb
因此有Ab=αAx
固有b=αx
因此x-b/α
不知道 您是不是这个意思。
9楼2012-01-12 10:24:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

saladin983

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
计算残量的话还需要给定x,不然怎么比较不同b时的残量?你是怎么设定x的呢?
2楼2012-01-10 00:21:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangxn06

金虫 (正式写手)

引用回帖:
: Originally posted by saladin983 at 2012-01-10 00:21:47:
计算残量的话还需要给定x,不然怎么比较不同b时的残量?你是怎么设定x的呢?

我不太明白您的意思,但是我试着按我的理解回答一下。
AX=b,我先用某种解线性代数方程组的方法,比如高斯全主元消去法,得到X的解,比如说为X'. 然后计算AX', 最后残余向量r=b-AX'.
我的问题是,理论而言这个残余向量为0时,计算得到的X'才是这个方程的真解。但是,我在计算的过程中,自由项不同,残余向量的二范数的值差别很大。
对于某些自由向量,残余向量的二范数可达到10e-12的以上
而对有些自由向量,残余向量的二范数只能达到10e-2。
是否有好的方法使得残余向量的二范数尽可能的接近零呢?即,使得得到的解x‘接近真解。
3楼2012-01-10 10:24:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

math2000

铁杆木虫 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
小雨萌萌(金币+1): 3Q~ 2012-01-11 14:53:11
方程组AX=b是不相容的,即没有解,需要解最小二乘解!
因为A是对称矩阵,可以考虑正交分解的方法。
4楼2012-01-10 16:16:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211 0703化学 346分求调剂 +19 土豆er? 2026-04-09 19/950 2026-04-09 20:11 by maddjdld
[考研] 301求调剂 +3 静静想想 2026-04-05 3/150 2026-04-09 19:12 by 诗与自由
[考研] 312求调剂 +3 李鸿飞飞 2026-04-06 3/150 2026-04-09 17:32 by wp06
[考研] 337求调剂 +4 Gky09300550, 2026-04-09 4/200 2026-04-09 17:18 by 帕尔马拉特
[考研] 085600,321分求调剂 +19 大馋小子 2026-04-04 20/1000 2026-04-09 14:12 by Delta2012
[考研] 招收有机化学、化工,药学,食品灯专业学生 +3 yrfhjgdj 2026-04-08 3/150 2026-04-09 10:15 by QYQX_123
[考研] 293调剂 +21 yj1221 2026-04-08 22/1100 2026-04-09 08:34 by 5268321
[考研] 材料考研求调剂总分280 +30 mkjlz1 2026-04-06 35/1750 2026-04-08 21:25 by cyh—315
[考研] 324求调剂 +17 想上学求调 2026-04-03 17/850 2026-04-08 20:04 by 我减肥1
[考研] 304求调剂 +10 素年祭语 2026-04-06 17/850 2026-04-08 09:05 by 蓝云思雨
[考研] 288环境专硕,求调材料方向 +35 lllllos 2026-04-04 39/1950 2026-04-07 23:24 by 一只好果子?
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[论文投稿] Decision: Revise for Editor还会送审吗 100+3 CccccccccFD 2026-04-04 5/250 2026-04-07 10:58 by 北京莱茵润色
[考研] 0857大类环境工程B区求调剂 +3 龚禹铭 2026-04-05 3/150 2026-04-06 10:22 by 蓝云思雨
[考研] 277求调剂 +5 考研调剂lxh 2026-04-05 5/250 2026-04-05 19:03 by chy09050039
[考研] 341求调剂 +3 学无止境,冲 2026-04-05 3/150 2026-04-05 09:40 by lbsjt
[考研] 可跨专业调剂 +3 周的得地 2026-04-04 6/300 2026-04-04 22:21 by barlinike
[考研] 342求调剂 +3 Liang7111 2026-04-04 5/250 2026-04-04 19:47 by dongzh2009
[考研] 372分材料与化工(085600)一志愿湖南大学求调剂 +3 蓝笺片 2026-04-03 4/200 2026-04-03 17:58 by Jimmyandyou
[考研] 一志愿北京交通大学材料工程总分358 +4 cs0106 2026-04-03 4/200 2026-04-03 13:41 by 百灵童888
信息提示
请填处理意见