 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2849 | 回復: 17 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
wangjunyahit銀蟲 (初入文壇)
|
[求助]
求Matlab高手解決線性方程組的迭代求解問題
|
|
|
遇到一個線性方程組Ax=b的求解問題，使用Matlab直接求解但不確定結(jié)果是否正確， 想使用迭代解法求解驗證結(jié)果，但苦于找不到合適的迭代解法和程序，希望各位大神 給予幫助，感激涕零�。�！ 方程組是這樣的： 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 ] |
回復此樓» 猜你喜歡
086502化學工程342求調(diào)劑
已經(jīng)有4人回復
-
材料求調(diào)劑一志愿哈工大324
已經(jīng)有7人回復
-
328求調(diào)劑
已經(jīng)有7人回復
-
340求調(diào)劑
已經(jīng)有5人回復
-
一志愿南昌大學324求調(diào)劑
已經(jīng)有4人回復
-
求調(diào)劑推薦 材料 304
已經(jīng)有15人回復
-
291求調(diào)劑
已經(jīng)有19人回復
-
331環(huán)境科學與工程求調(diào)劑
已經(jīng)有3人回復
-
289求調(diào)劑
已經(jīng)有8人回復
-
315分求調(diào)劑
已經(jīng)有7人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- matlab求解非線性方程組，錯誤提示怎么解決
已經(jīng)有5人回復
- 用不動點迭代求非線性方程組，求得的不是想要的解，怎么辦
已經(jīng)有12人回復
- matlab求解非線性方程組，求助！
已經(jīng)有6人回復
- 求助！matlab用fsolve函數(shù)求解非線性方程組的問題！
已經(jīng)有19人回復
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(shù)值比較
已經(jīng)有4人回復
- 如何使用matlab求解非線性方程組的所有整數(shù)解？
已經(jīng)有9人回復
- Matlab求解偏微分方程組
已經(jīng)有22人回復
- matlab求解非線性方程組
已經(jīng)有16人回復
- matlab數(shù)值求解邊界條件微分方程組
已經(jīng)有7人回復
- 求matlab求解一個方程組的問題
已經(jīng)有6人回復
- 求高人指點用matlab求解非線性方程組，解決了追加100金幣；
已經(jīng)有11人回復
- 請教Matlab高手解多元非線性方程組中的參數(shù)
已經(jīng)有6人回復
- matlab的fsove 命令求解非線性方程組
已經(jīng)有6人回復
- MATLAB求解一個比較復雜的方程
已經(jīng)有14人回復
- 都不知道怎么命名標題：迭代法解方程組？
已經(jīng)有16人回復
- 【求助】MATLAB 解方程組【已解決】
已經(jīng)有7人回復
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題
已經(jīng)有21人回復
- 【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實數(shù)解！
已經(jīng)有18人回復
- 【求助】求教 多元指數(shù)方程如何解 有軟件否？
已經(jīng)有6人回復
- 【求助】matlab怎么求解偏微分方程組啊，先謝謝了
已經(jīng)有13人回復
- 【求助】100金幣求助高人一個程序解方程
已經(jīng)有15人回復
- 非線性方程組的迭代法（數(shù)值計算高手請進）
已經(jīng)有7人回復
- 【求助】求教matlab解非線性方程組
已經(jīng)有9人回復
- 【求助】求一個用牛頓迭代法解非線性方程組的程序
已經(jīng)有10人回復
- 【求助】如何用Runge-Kutta迭代求解二階常微分方程組【已解決】
已經(jīng)有9人回復
9樓2012-11-25 17:38:11
【答案】應助回帖
★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
xiegangmai: 金幣+2, 謝謝應助 2012-11-24 20:54:11
wangjunyahit: 金幣+3, ★★★很有幫助 2012-11-25 11:12:49
感謝參與，應助指數(shù) +1
xiegangmai: 金幣+2, 謝謝應助 2012-11-24 20:54:11
wangjunyahit: 金幣+3, ★★★很有幫助 2012-11-25 11:12:49
|
矩陣的rcond顯然可以改好啊。例如將第一列乘以4e27。 >> Anew = A; Anew(:,1) = A(:,1) *4e27;rcond(Anew) 0.0103 >>x_1= Anew\b; x_1_real = x_1; x_1_real(1) = 4e27 * x_1(1) 不過結(jié)果解還是一樣的。 |
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ù)值分析沒有,這個可以用高斯迭代法和雅可比迭代法求解,不過雖然用迭代法,跑程序時還是有警告,可能求解還是不是非常可靠,我求解了一個,下面我把高斯迭代法那個我求解用的程序和運行結(jié)果貼在下面: 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('迭代結(jié)果:'); x=x0' disp('逆矩陣求解結(jié)果：'); 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 迭代結(jié)果: x = 1.029393177754839e+020 6.565808888911417e-008 2.259822851544943e-008 3.842185487752933e-009 4.113382036189012e-009 逆矩陣求解結(jié)果： 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è): 工程熱物理相關(guān)交叉領(lǐng)域
4樓2012-11-25 09:28:56
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +6 | 閆旭東 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 322求調(diào)劑 +5 | 宋明欣 2026-03-27 | 5/250 |
|
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
|
[考研] 286求調(diào)劑 +3 | 丟掉懶惰 2026-03-27 | 6/300 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +12 | minht 2026-03-26 | 12/600 |
|
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
|
[考研] 070300求調(diào)劑306分 +3 | 26要上岸 2026-03-27 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
|
[考研] 329求調(diào)劑 +7 | 鈕恩雪 2026-03-25 | 7/350 |
|
|
|
[考研] 294分080500材料科學與工程求調(diào)劑 +4 | 柳溪邊 2026-03-26 | 4/200 |
|
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
|
[考研] 321求調(diào)劑 +5 | 材料cailiao 2026-03-21 | 5/250 |
|
|
|
[考研] 329求調(diào)劑 +5 | 1() 2026-03-22 | 5/250 |
|
|
|
[考研]
 085602 289分求調(diào)劑
+8
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +6 | jiajunX 2026-03-22 | 6/300 |
|
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +10 | 紙魚ly 2026-03-21 | 11/550 |
|
|
|
[考研]
 085404電子信息284分求調(diào)劑
+4
|
13659058978 2026-03-24 | 4/200 |
|