| 查看: 1242 | 回复: 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 |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 各位大神 小弟初拥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
- 专业: 计算数学与科学工程计算
5楼2015-01-11 18:57:32
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
