 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (3187)
- >考研 (580)
- >導(dǎo)師招生 (466)
- >蟲(chóng)友互識(shí) (125)
- >文獻(xiàn)求助 (83)
- >碩博家園 (51)
- >考博 (45)
- >論文投稿 (40)
- >博后之家 (39)
- >招聘信息布告欄 (29)
- >休閑灌水 (28)
- >基金申請(qǐng) (22)
- >找工作 (17)
- >教師之家 (15)
- >文獻(xiàn)求助完結(jié)子版 (14)
- >論文道賀祈福 (13)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 610 | 回復(fù): 3 | ||
duinichixin金蟲(chóng) (著名寫(xiě)手)
|
[求助]
小弟求教幾個(gè)程序-關(guān)于迭代法的（多謝） 已有1人參與
|
|
小弟編寫(xiě)了Table4.1SOR中的r=0.01及n=8的程序，現(xiàn)有如下幾個(gè)疑問(wè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時(shí)出現(xiàn)‘Out of memory. Type HELP MEMORY for your options.’請(qǐng)問(wèn)是我程序太復(fù)雜了嗎？ 注：test為確定w(exp)程序，SOR為調(diào)用程序，test2為根據(jù)w(exp)計(jì)算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的對(duì)角矩陣 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的對(duì)角矩陣 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ù) %迭代過(guò)程 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ù)此樓» 猜你喜歡
291求調(diào)劑
已經(jīng)有15人回復(fù)
-
材料學(xué)碩333求調(diào)劑
已經(jīng)有17人回復(fù)
-
070300化學(xué)求調(diào)劑
已經(jīng)有5人回復(fù)
-
070300一志愿211，312分求調(diào)劑院校
已經(jīng)有6人回復(fù)
-
物理學(xué)調(diào)劑
已經(jīng)有3人回復(fù)
-
08工科，295，接受跨專業(yè)調(diào)劑
已經(jīng)有5人回復(fù)
-
材料科學(xué)與工程調(diào)劑
已經(jīng)有12人回復(fù)
-
本2一志愿C9-333分，材料科學(xué)與工程，求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿西交288化工專碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
308求調(diào)劑
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 各位大神 小弟初擁origin 做循環(huán)伏安曲線 不懂求教
已經(jīng)有5人回復(fù)
- 數(shù)學(xué)的大神們，求解一個(gè)數(shù)學(xué)問(wèn)題，大神們幫幫忙，小女子感激萬(wàn)分
已經(jīng)有2人回復(fù)
- 求大神用matlab給我處理一組數(shù)據(jù)，有人沒(méi)有啊？
已經(jīng)有15人回復(fù)
- SVM 利用matlab仿真 求大神
已經(jīng)有10人回復(fù)
- 求Matlab解非線性方程代碼
已經(jīng)有12人回復(fù)
- 關(guān)于四川大學(xué)考博專業(yè)課參考書(shū)的求助
已經(jīng)有4人回復(fù)
- 怎么才能做出三個(gè)9的線性方程�。�？��！愁死我了
已經(jīng)有11人回復(fù)
- ArcGIS文件迭代器的使用問(wèn)題
已經(jīng)有5人回復(fù)
duinichixin
金蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 892
- 散金: 596
- 紅花: 4
- 帖子: 1009
- 在線: 338.5小時(shí)
- 蟲(chóng)號(hào): 2475163
- 注冊(cè): 2013-05-18
- 性別: GG
- 專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
2樓2015-01-11 18:39:45
匿名
用戶注銷(xiāo) (職業(yè)作家)
Matlab專家
- 應(yīng)助: 223 (大學(xué)生)
- 金幣: 2891.3
- 散金: 5302
- 紅花: 59
- 沙發(fā): 10
- 帖子: 3308
- 在線: 945小時(shí)
- 蟲(chóng)號(hào): 0
- 注冊(cè): 2012-05-01
- 性別: GG
- 專業(yè): 力學(xué)
3樓2015-01-11 19:20:45
duinichixin
金蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 892
- 散金: 596
- 紅花: 4
- 帖子: 1009
- 在線: 338.5小時(shí)
- 蟲(chóng)號(hào): 2475163
- 注冊(cè): 2013-05-18
- 性別: GG
- 專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
4樓2015-01-11 19:55:33
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 duinichixin 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293分求調(diào)劑，外語(yǔ)為俄語(yǔ) +4 | 加一一九 2026-03-31 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +11 | 多久上課 2026-03-27 | 12/600 |
|
|
|
[考研] 286分調(diào)劑 +10 | Faune 2026-03-30 | 11/550 |
|
|
|
[考研] 各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩 +10 | Koxui 2026-03-28 | 10/500 |
|
|
|
[考研] 求調(diào)劑323材料與化工 +10 | 1124361 2026-03-24 | 10/500 |
|
|
[考研] 284求調(diào)劑 +14 | junqihahaha 2026-03-26 | 15/750 |
|
|
|
[考研] 085701求調(diào)劑初試286分 +5 | secret0328 2026-03-28 | 5/250 |
|
|
|
[考研] 337求調(diào)劑 +6 | 《樹(shù)》 2026-03-29 | 6/300 |
|
|
|
[考研] 085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富 +6 | 來(lái)財(cái)aa 2026-03-25 | 6/300 |
|
|
|
[考研] 343求調(diào)劑085601 +3 | 要努力學(xué)習(xí)x 2026-03-29 | 3/150 |
|
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
|
[考研] 壓國(guó)家一區(qū)線，求導(dǎo)師收留，有恩必謝！ +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
|
[考研] 08開(kāi)頭275求調(diào)劑 +4 | 拉誰(shuí)不重要 2026-03-26 | 4/200 |
|
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|