| 5 | 1/1 | 返回列表 |
| 查看: 2478 | 回復(fù): 10 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
wangxn06金蟲 (正式寫手)
|
[求助]
解線性代數(shù)方程組
|
||
|
問題描述:有一個(gè)線性方程組AX=b(A是方陣)(系數(shù)矩陣A的大小為200*200以上) 系數(shù)矩陣A 是對(duì)稱矩陣,對(duì)角不占優(yōu) 不同的自由項(xiàng)b所計(jì)算得到的剩余殘差不同,即的到r=b-AX'的量級(jí)上有很大差別。 問題:1:同樣的系數(shù)矩陣,不同的自由項(xiàng),為什么得到的殘差量級(jí)上會(huì)有很大的差別,這是不是由于系數(shù)矩陣的某些特點(diǎn)引起的——比如說矩陣A的條件數(shù)過大? 2.是否有合適的方法可以使得殘差達(dá)到比較合適的量級(jí)。 我用過全主元高斯消去法、列主元高斯消去法、并且在全主元高斯消去法的基礎(chǔ)上,應(yīng)用了解的迭代改進(jìn)法、解對(duì)稱正定線性方程組的楚列斯基法,但是這些都沒有提高解的精確度(這些方法對(duì)某些自由向量b,解的結(jié)果精度都比較高,但是對(duì)某些自由向量b解的精度都不足,而這些自由向量b之間到底有什么差別也不好說)。Sample Text |
金蟲 (正式寫手)
|
我不太明白您的意思,但是我試著按我的理解回答一下。 AX=b,我先用某種解線性代數(shù)方程組的方法,比如高斯全主元消去法,得到X的解,比如說為X'. 然后計(jì)算AX', 最后殘余向量r=b-AX'. 我的問題是,理論而言這個(gè)殘余向量為0時(shí),計(jì)算得到的X'才是這個(gè)方程的真解。但是,我在計(jì)算的過程中,自由項(xiàng)不同,殘余向量的二范數(shù)的值差別很大。 對(duì)于某些自由向量,殘余向量的二范數(shù)可達(dá)到10e-12的以上 而對(duì)有些自由向量,殘余向量的二范數(shù)只能達(dá)到10e-2。 是否有好的方法使得殘余向量的二范數(shù)盡可能的接近零呢?即,使得得到的解x‘接近真解。 |
鐵桿木蟲 (正式寫手)
鐵桿木蟲 (職業(yè)作家)
金蟲 (正式寫手)
|
首先方程相容不相容,由系數(shù)矩陣陣A的秩以及擴(kuò)張矩陣的秩的大小有關(guān)系,當(dāng)系數(shù)矩陣的秩小于擴(kuò)展矩陣的秩的時(shí)候,方程組不相容。 首先假設(shè)是不相容的,這意味著,矩陣A肯定是不滿秩的。這個(gè)或許我可以試試。 但是問題又出來了,如果矩陣不滿秩,意味著矩陣A奇異,那么不管自由項(xiàng)是什么,方程組都得不到唯一的解。但是,我的問題中,在某些自由項(xiàng)下,還是可以有很好的結(jié)果的。所以,方程組不可能是不相容的。 同時(shí),您提到了正交分解法,我想問一下這個(gè)方法相比其他方法有什么優(yōu)點(diǎn),是否從理論上可以使得得到的解x'更加接近真解。 非常感謝! |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 江蘇蘇北高校誠(chéng)邀調(diào)劑同學(xué) +3 | zzll406 2026-03-31 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +6 | Faune 2026-03-30 | 6/300 |
|
|
[考研] 一志愿中海洋材料357 +4 | 麥恩莉. 2026-03-30 | 4/200 |
|
|
[考研] 272求調(diào)劑,接受跨專業(yè)調(diào)劑! +3 | 閑魚盧 2026-03-31 | 3/150 |
|
|
[考研] 土木304求調(diào)劑 +4 | 兔突突突, 2026-03-31 | 5/250 |
|
|
[考研] 085600材料與化工329分求調(diào)劑 +5 | 這是有太陽(yáng)哇 2026-03-25 | 5/250 |
|
|
[考研] 085600 286分 材料求調(diào)劑 +11 | 麻辣魷魚 2026-03-27 | 12/600 |
|
|
[考研] 295材料工程專碩求調(diào)劑 +10 | 1428151015 2026-03-27 | 10/500 |
|
|
[考研] 349求調(diào)劑 +6 | 李木子啊哈哈 2026-03-25 | 6/300 |
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-27 | 8/400 |
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|