| 查看: 1216 | 回复: 7 | ||||
duinichixin金虫 (著名写手)
|
[求助]
小弟求教几个程序-关于迭代法的 已有2人参与
|
|
小弟编写了Table4.1SOR中的r=0.01及n=8的程序，现有如下几个疑问：1、我运行出it为8，cpu=0.042，w(exp)=1.02，为什么与作者不一样？2、作者终止条件是F范数，我的是2范数，应该怎么修改？3.当n=64时出现‘Out of memory. Type HELP MEMORY for your options.’请问是我程序太复杂了吗？ 注：test为确定w(exp)程序，SOR为调用程序，test2为根据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; %迭代次数 %迭代过程 while norm(x-x0)>=eps x0=x; x=B*x0+f; n1=n1+1; if(n1>=M) disp('Warining:迭代次数太多，可能不收敛！'); 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 |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
fileCode有新解读？
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 各位大神 小弟初拥origin 做循环伏安曲线 不懂求教
已经有5人回复
- 数学的大神们，求解一个数学问题，大神们帮帮忙，小女子感激万分
已经有2人回复
- 求大神用matlab给我处理一组数据，有人没有啊？
已经有15人回复
- SVM 利用matlab仿真 求大神
已经有10人回复
- 求Matlab解非线性方程代码
已经有12人回复
- 关于四川大学考博专业课参考书的求助
已经有4人回复
- 怎么才能做出三个9的线性方程啊？？！！愁死我了
已经有11人回复
- ArcGIS文件迭代器的使用问题
已经有5人回复
duinichixin
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 892
- 散金: 596
- 红花: 4
- 帖子: 1009
- 在线: 338.5小时
- 虫号: 2475163
- 注册: 2013-05-18
- 性别: GG
- 专业: 计算数学与科学工程计算
2楼2015-01-11 18:35:51
duinichixin
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 892
- 散金: 596
- 红花: 4
- 帖子: 1009
- 在线: 338.5小时
- 虫号: 2475163
- 注册: 2013-05-18
- 性别: GG
- 专业: 计算数学与科学工程计算
3楼2015-01-11 18:37:21
mathstudy
金虫 (正式写手)
- 数学EPI: 2
- 应助: 161 (高中生)
- 金币: 4320.1
- 散金: 38
- 红花: 16
- 帖子: 446
- 在线: 141.1小时
- 虫号: 2515489
- 注册: 2013-06-20
- 性别: GG
- 专业: 计算数学与科学工程计算
4楼2015-01-11 18:42:26
duinichixin
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 892
- 散金: 596
- 红花: 4
- 帖子: 1009
- 在线: 338.5小时
- 虫号: 2475163
- 注册: 2013-05-18
- 性别: GG
- 专业: 计算数学与科学工程计算
5楼2015-01-11 18:57:32
zaq123321
专家顾问 (著名写手)
-
专家经验: +343 - 数学EPI: 6
- 应助: 298 (大学生)
- 贵宾: 0.247
- 金币: 11171.3
- 散金: 295
- 红花: 29
- 帖子: 1223
- 在线: 539.9小时
- 虫号: 405284
- 注册: 2007-06-17
- 性别: MM
- 专业: 生物大分子结构与功能
- 管辖: 数学
【答案】应助回帖
感谢参与，应助指数 +1
|
Chang your 2 norm to frobenius norm by replacing 2 with fro. Test it to make sure it works. U don't need to care about cpu time in the testing since your computer is different with his. Your computer memory may be not big enough to run n=64 case. Otherwise check whether there is some place to adjust your memory allocation in your code. Good luck [ 发自手机版 http://muchong.com/3g ] |
6楼2015-01-11 21:30:40
duinichixin
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 892
- 散金: 596
- 红花: 4
- 帖子: 1009
- 在线: 338.5小时
- 虫号: 2475163
- 注册: 2013-05-18
- 性别: GG
- 专业: 计算数学与科学工程计算
7楼2015-01-12 08:25:56
mathstudy
金虫 (正式写手)
- 数学EPI: 2
- 应助: 161 (高中生)
- 金币: 4320.1
- 散金: 38
- 红花: 16
- 帖子: 446
- 在线: 141.1小时
- 虫号: 2515489
- 注册: 2013-06-20
- 性别: GG
- 专业: 计算数学与科学工程计算
8楼2015-01-12 09:04:35
100