 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2855 | 回復: 17 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
wangjunyahit銀蟲 (初入文壇)
|
[求助]
求Matlab高手解決線性方程組的迭代求解問題
|
|
|
遇到一個線性方程組Ax=b的求解問題，使用Matlab直接求解但不確定結果是否正確， 想使用迭代解法求解驗證結果，但苦于找不到合適的迭代解法和程序，希望各位大神 給予幫助，感激涕零�。�！ 方程組是這樣的： A=[-1.2e-21,1.4e6,1.4e6,-1.4e6,1.3e6; 8.9e-22,-1.4e6,4.9e3,3.8e3,0; 1.8e-22,1.4e3,-1.4e6,4.3e3,0; 4.4e-23,3.2e2,1.2e3,-5.1e6,4.3e3; 5.3e-23,0,0,1.3e3,-2.3e6] b=[0;-1.8e-4;-1.3e-2;-1.5e-2;-4.0e-3] 直接求解如下： x=A\b; x=[ 1.029393177754837e+20 6.565808888911404e-08 2.259822851544940e-08 3.842185487752932e-09 4.113382036189007e-09] rcond(A)=1.694258513618419e-29;矩陣A可逆的條件數(shù)估值接近0，矩陣A條件不好，得到的解有誤差 希望各位大神 給予幫助，感激涕零！�。� ~EZ9ND[(A@F3IUUY9EJVQ4C.jpg [ Last edited by wangjunyahit on 2012-11-24 at 16:10 ] |
回復此樓» 猜你喜歡
學碩274求調劑
已經有11人回復
-
一志愿鄭大材料工程290求調劑
已經有12人回復
-
環(huán)境工程 085701，267求調劑
已經有8人回復
-
083000學碩274求調劑
已經有10人回復
-
一志愿085600中科院寧波所276分求調劑
已經有4人回復
-
材料專業(yè)申博
已經有3人回復
-
307求調劑
已經有20人回復
-
285求調劑
已經有8人回復
-
材料與化工（0856）304求B區(qū)調劑
已經有8人回復
-
370求調劑
已經有3人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- matlab求解非線性方程組，錯誤提示怎么解決
已經有5人回復
- 用不動點迭代求非線性方程組，求得的不是想要的解，怎么辦
已經有12人回復
- matlab求解非線性方程組，求助！
已經有6人回復
- 求助！matlab用fsolve函數(shù)求解非線性方程組的問題！
已經有19人回復
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(shù)值比較
已經有4人回復
- 如何使用matlab求解非線性方程組的所有整數(shù)解？
已經有9人回復
- Matlab求解偏微分方程組
已經有22人回復
- matlab求解非線性方程組
已經有16人回復
- matlab數(shù)值求解邊界條件微分方程組
已經有7人回復
- 求matlab求解一個方程組的問題
已經有6人回復
- 求高人指點用matlab求解非線性方程組，解決了追加100金幣；
已經有11人回復
- 請教Matlab高手解多元非線性方程組中的參數(shù)
已經有6人回復
- matlab的fsove 命令求解非線性方程組
已經有6人回復
- MATLAB求解一個比較復雜的方程
已經有14人回復
- 都不知道怎么命名標題：迭代法解方程組？
已經有16人回復
- 【求助】MATLAB 解方程組【已解決】
已經有7人回復
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題
已經有21人回復
- 【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實數(shù)解！
已經有18人回復
- 【求助】求教 多元指數(shù)方程如何解 有軟件否？
已經有6人回復
- 【求助】matlab怎么求解偏微分方程組啊，先謝謝了
已經有13人回復
- 【求助】100金幣求助高人一個程序解方程
已經有15人回復
- 非線性方程組的迭代法（數(shù)值計算高手請進）
已經有7人回復
- 【求助】求教matlab解非線性方程組
已經有9人回復
- 【求助】求一個用牛頓迭代法解非線性方程組的程序
已經有10人回復
- 【求助】如何用Runge-Kutta迭代求解二階常微分方程組【已解決】
已經有9人回復
wangjunyahit
銀蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 366.5
- 帖子: 19
- 在線: 11.2小時
- 蟲號: 1697137
- 注冊: 2012-03-17
- 性別: GG
- 專業(yè): 工程熱物理相關交叉領域
4樓2012-11-25 09:28:56
2樓2012-11-24 20:35:20
fabaoyi
新蟲 (初入文壇)
- 應助: 1 (幼兒園)
- 金幣: 28
- 帖子: 11
- 在線: 8.6小時
- 蟲號: 2066628
- 注冊: 2012-10-17
- 性別: GG
- 專業(yè): 應用數(shù)學方法
【答案】應助回帖
★ ★
感謝參與，應助指數(shù) +1
xiegangmai: 金幣+2, 謝謝參與 2012-11-25 19:58:30
感謝參與，應助指數(shù) +1
xiegangmai: 金幣+2, 謝謝參與 2012-11-25 19:58:30
|
你學過數(shù)值分析沒有,這個可以用高斯迭代法和雅可比迭代法求解,不過雖然用迭代法,跑程序時還是有警告,可能求解還是不是非�？煽�,我求解了一個,下面我把高斯迭代法那個我求解用的程序和運行結果貼在下面: function gsdd(A,b,x0,fs,wucha,ddcs) D=diag(diag(A)); U=-triu(A,1); L=-tril(A,-1); ljzx=inv(A)*b'; if det(D)==0 disp('對角矩陣D奇異，迭代不能實現(xiàn)'); return end B=(D-L)\U;f=(D-L)\b'; for k=1:ddcs x1=B*x0'+f; epsilon=norm(x1-x0',fs); wcx=norm(x1-x0',fs)/(norm(x1,fs)+eps); x0=x1'; if (epsilon end end if (epsilon>wucha)||(wcx>wucha)||k==ddcs disp(' 超過給定最大迭代次數(shù)，可能不收斂 ') return end disp('迭代次數(shù):'); k disp('迭代結果:'); x=x0' disp('逆矩陣求解結果：'); x=ljzx end >> A=[-1.2e-21,1.4e6,1.4e6,-1.4e6,1.3e6; 8.9e-22,-1.4e6,4.9e3,3.8e3,0; 1.8e-22,1.4e3,-1.4e6,4.3e3,0; 4.4e-23,3.2e2,1.2e3,-5.1e6,4.3e3; 5.3e-23,0,0,1.3e3,-2.3e6]; >> b=[0;-1.8e-4;-1.3e-2;-1.5e-2;-4.0e-3]'; >> x0=[1 1 1 1 1]; >> gsdd(A,b,x0,1,10^(-15),1000) Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 1.694259e-029. > In gsdd at 5 Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 2.352941e-028. > In gsdd at 10 Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 2.352941e-028. > In gsdd at 10 迭代次數(shù): k = 527 迭代結果: x = 1.029393177754839e+020 6.565808888911417e-008 2.259822851544943e-008 3.842185487752933e-009 4.113382036189012e-009 逆矩陣求解結果： x = 1.029393177754837e+020 6.565808888911401e-008 2.259822851544940e-008 3.842185487752931e-009 4.113382036189006e-009 |
3樓2012-11-25 01:28:13
wangjunyahit
銀蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 366.5
- 帖子: 19
- 在線: 11.2小時
- 蟲號: 1697137
- 注冊: 2012-03-17
- 性別: GG
- 專業(yè): 工程熱物理相關交叉領域
5樓2012-11-25 09:33:37
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 285求調劑 +5 | AZMK 2026-03-30 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 083000環(huán)境科學與工程調劑，總分281 +4 | 橙子（勝意） 2026-03-30 | 4/200 |
|
|
|
[考研] 285求調劑 +6 | AZMK 2026-03-29 | 9/450 |
|
|
|
[考研] 一志愿：西北大學，英一數(shù)一408-284分求調劑 +5 | 12.27 2026-03-27 | 5/250 |
|
|
|
[考研] 一志愿廈門大學材料工程專碩354找調劑�。�！ +5 | 貝唄鋇鋇 2026-03-30 | 5/250 |
|
|
|
[考研] 298求調劑 +3 | 什么是胖頭魚 2026-03-30 | 5/250 |
|
|
|
[考研] 本科211生物醫(yī)學工程085409求調劑339分 +3 | 里子木yy 2026-03-29 | 3/150 |
|
|
|
[考研] 化學308分調劑 +7 | 你好明天你好 2026-03-30 | 8/400 |
|
|
|
[考研] 337求調劑 +6 | 《樹》 2026-03-29 | 6/300 |
|
|
|
[考研] 317求調劑 +10 | 蛋黃咸肉粽 2026-03-26 | 10/500 |
|
|
|
[考研] 材料與化工304求B區(qū)調劑 +4 | 邱gl 2026-03-26 | 7/350 |
|
|
|
[考研] 295求調劑 +5 | wei-5 2026-03-26 | 5/250 |
|
|
|
[考研] 290求調劑 +3 | dfffsar 2026-03-29 | 3/150 |
|
|
|
[考研]
 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+11
|
nnnnnnn5 2026-03-25 | 11/550 |
|
|
|
[考研] 343求調劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
|
[考研] 286求調劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
|
[考研] 調劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
|
[考研] 求b區(qū)院校調劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
|
[考研] 300分，材料，求調劑，英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
|
[考研]
材料考研調劑生
+3
|
黃粱一夢千年 2026-03-24 | 3/150 |
|