 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 611 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
duinichixin金蟲 (著名寫手)
|
[求助]
小弟求教幾個程序-關(guān)于迭代法的（多謝） 已有1人參與
|
|
|
小弟編寫了Table4.1SOR中的r=0.01及n=8的程序，現(xiàn)有如下幾個疑問：1、我運(yùn)行出it為8，cpu=0.042，w(exp)=1.02，為什么與作者不一樣？2、作者終止條件是F范數(shù)，我的是2范數(shù)，應(yīng)該怎么修改？3.當(dāng)n=64時出現(xiàn)‘Out of memory. Type HELP MEMORY for your options.’請問是我程序太復(fù)雜了嗎？ 注：test為確定w(exp)程序，SOR為調(diào)用程序，test2為根據(jù)w(exp)計算it程序，各位大神不吝賜教。 test.m: n=8; r=0.01; a=ones(n,1); b=ones(n-1,1); I=eye(n); MM=2*diag(a)-diag(b,-1)-diag(b,1); N=0.5*diag(b,-1)-0.5*diag(b,1); A=MM+2*r*N+100*I/(n+1); B=A; X=rand(n,n); F=A*X+X*B; E=kron(I,A)+kron(B',I); b=F(  ;%x=E\b; x0=b; ww=0.01:0.01:1.99; D=diag(diag(E)); %求E的對角矩陣 L=-tril(E,-1); %求E的下三角矩陣 U=-triu(E,1); %求E的上三角矩陣 for i=1:length(ww) w=ww(i); B=inv(D-L*w)*((1-w)*D+w*U); p(i)=max(max(abs(eig(B)))); %plot(ww,p); [x,n1]=SOR(E,b,x0,w); mm(i)=n1; end plot(ww,mm); Xe=reshape(x,n,n); errror=norm((Xe-X),2) SOR.m: function [x,n1]=SOR(E,b,x0,w,eps,M) if nargin==4 eps=1.0e-6; M=2000; elseif nargin<4 error return elseif nargin==5 M=200; end if(w<=0 || w>=2) %收斂條件要求 error; return; end D=diag(diag(E)); %求E的對角矩陣 L=-tril(E,-1); %求E的下三角矩陣 U=-triu(E,1); %求E的上三角矩陣 B=inv(D-L*w)*((1-w)*D+w*U); f=w*inv((D-L*w))*b; x=B*x0+f; n1=1; %迭代次數(shù) %迭代過程 while norm(x-x0)>=eps x0=x; x=B*x0+f; n1=n1+1; if(n1>=M) disp('Warining:迭代次數(shù)太多，可能不收斂！'); return; end end test2.m w=1.02 n=8; r=0.01; a=ones(n,1); b=ones(n-1,1); I=eye(n); M=2*diag(a)-diag(b,-1)-diag(b,1); N=0.5*diag(b,-1)-0.5*diag(b,1); A=M+2*r*N+100*I/(n+1); B=A; X=rand(n,n); F=A*X+X*B; E=kron(I,A)+kron(B',I); b=F( ;%x=E\b; x0=b; [x,n1]=SOR(E,b,x0,w) Xe=reshape(x,n,n); errror=norm((Xe-X),2) 1.png 2.png |
回復(fù)此樓» 猜你喜歡
334分 一志愿武理-080500 材料求調(diào)劑
已經(jīng)有10人回復(fù)
-
324分 085600材料與化工
已經(jīng)有7人回復(fù)
-
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩
已經(jīng)有11人回復(fù)
-
309求調(diào)劑
已經(jīng)有6人回復(fù)
-
調(diào)劑推薦
已經(jīng)有7人回復(fù)
-
一志愿西交大080500材料學(xué)碩349
已經(jīng)有7人回復(fù)
-
一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑
已經(jīng)有4人回復(fù)
-
352分-085602-一志愿985
已經(jīng)有6人回復(fù)
-
273求調(diào)劑
已經(jīng)有5人回復(fù)
-
349求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 各位大神 小弟初擁origin 做循環(huán)伏安曲線 不懂求教
已經(jīng)有5人回復(fù)
- 數(shù)學(xué)的大神們，求解一個數(shù)學(xué)問題，大神們幫幫忙，小女子感激萬分
已經(jīng)有2人回復(fù)
- 求大神用matlab給我處理一組數(shù)據(jù)，有人沒有��？
已經(jīng)有15人回復(fù)
- SVM 利用matlab仿真 求大神
已經(jīng)有10人回復(fù)
- 求Matlab解非線性方程代碼
已經(jīng)有12人回復(fù)
- 關(guān)于四川大學(xué)考博專業(yè)課參考書的求助
已經(jīng)有4人回復(fù)
- 怎么才能做出三個9的線性方程�。�？！！愁死我了
已經(jīng)有11人回復(fù)
- ArcGIS文件迭代器的使用問題
已經(jīng)有5人回復(fù)
匿名
用戶注銷 (職業(yè)作家)
Matlab專家
- 應(yīng)助: 223 (大學(xué)生)
- 金幣: 2891.3
- 散金: 5302
- 紅花: 59
- 沙發(fā): 10
- 帖子: 3308
- 在線: 945小時
- 蟲號: 0
- 注冊: 2012-05-01
- 性別: GG
- 專業(yè): 力學(xué)
3樓2015-01-11 19:20:45
duinichixin
金蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 892
- 散金: 596
- 紅花: 4
- 帖子: 1009
- 在線: 338.5小時
- 蟲號: 2475163
- 注冊: 2013-05-18
- 性別: GG
- 專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算
2樓2015-01-11 18:39:45
duinichixin
金蟲 (著名寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 892
- 散金: 596
- 紅花: 4
- 帖子: 1009
- 在線: 338.5小時
- 蟲號: 2475163
- 注冊: 2013-05-18
- 性別: GG
- 專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算
4樓2015-01-11 19:55:33
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 調(diào)劑 +4 | 好好讀書。 2026-03-28 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 復(fù)試調(diào)劑 +7 | 雙馬尾痞老板2 2026-03-31 | 7/350 |
|
|
[考研] 304求調(diào)劑 +8 | 素年祭語 2026-03-31 | 11/550 |
|
|
|
[考研] 318求調(diào)劑 +10 | 陳晨79 2026-03-30 | 10/500 |
|
|
|
[考研] 求調(diào)劑 生物學(xué) 377分 +6 | zzll03 2026-03-31 | 6/300 |
|
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +20 | Mr. Z 2026-03-25 | 21/1050 |
|
|
|
[考研] 330分求調(diào)劑 +6 | qzenlc 2026-03-29 | 6/300 |
|
|
|
[考研] 362求調(diào)劑 +9 | 西南交材料專碩3 2026-03-31 | 9/450 |
|
|
|
[考研] 282求調(diào)劑 不挑專業(yè) 求收留 +4 | Yam. 2026-03-30 | 5/250 |
|
|
|
[考研] 一志愿085600中科院寧波所276分求調(diào)劑 +12 | 材料學(xué)257求調(diào)劑 2026-03-28 | 13/650 |
|
|
|
[考研] 334求調(diào)劑 +7 | Trying] 2026-03-31 | 7/350 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-31 | 4/200 |
|
|
|
[考研] 085601一志愿西北工業(yè)大學(xué)初試346 +4 | 085601初試346 2026-03-30 | 4/200 |
|
|
|
[考研] 327求調(diào)劑 +5 | 小卡不卡. 2026-03-29 | 5/250 |
|
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 321求調(diào)劑 +7 | 璞玉~~ 2026-03-25 | 8/400 |
|
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|