| 查看: 428 | 回复: 6 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 最爱花诗雨 的 4 个金币 | |||
最爱花诗雨金虫 (正式写手)
穷鬼
|
[交流]
【求助】MATLAB的求解问题
|
||
求助:matlab 中的jacobi程序和松弛迭代的程序 [ Last edited by wuguocheng on 2009-10-21 at 20:20 ] |
回复此楼» 猜你喜欢
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有6人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
最爱花诗雨
金虫 (正式写手)
穷鬼
- 应助: 4 (幼儿园)
- 贵宾: 0.804
- 金币: 1088.2
- 散金: 250
- 红花: 3
- 帖子: 473
- 在线: 65.7小时
- 虫号: 643243
- 注册: 2008-11-01
- 专业: 生物化学
2楼2009-10-22 18:34:57
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 仿真EPI: 1
- 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555140
- 散金: 27828
- 红花: 2147
- 沙发: 66666
- 帖子: 1602243
- 在线: 65200.8小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
3楼2009-10-22 19:31:57
★ ★ ★
nono2009(金币+2,VIP+0):谢谢应助。欢迎常来交流。 10-22 20:07
最爱花诗雨(金币+1,VIP+0):谢谢支持 10-27 19:15
nono2009(金币+2,VIP+0):谢谢应助。欢迎常来交流。 10-22 20:07
最爱花诗雨(金币+1,VIP+0):谢谢支持 10-27 19:15
|
function [x,n]=jacobi(A,b,x0,eps,varargin) if nargin==3 eps= 1.0e-6; M = 200; elseif nargin<3 error return elseif nargin ==5 M = varargin{1}; end D=diag(diag(A)); %求A的对角矩阵 L=-tril(A,-1); %求A的下三角阵 U=-triu(A,1); %求A的上三角阵 B=D\(L+U); f=D\b; x=B*x0+f; n=1; %迭代次数 while norm(x-x0)>=eps x0=x; x=B*x0+f; n=n+1; if(n>=M) disp('Warning: 迭代次数太多，可能不收敛！'); return; end end |
4楼2009-10-22 20:04:16
★ ★ ★
nono2009(金币+1,VIP+0):再奖1个！ 10-22 20:08
最爱花诗雨(金币+2,VIP+0):正是我要找的sor迭代，不知道能不能调用？ 10-23 18:02
nono2009(金币+1,VIP+0):再奖1个！ 10-22 20:08
最爱花诗雨(金币+2,VIP+0):正是我要找的sor迭代，不知道能不能调用？ 10-23 18:02
|
function [x,n]=SOR(A,b,x0,w,eps,M) if nargin==4 eps= 1.0e-6; M = 200; elseif nargin<4 error return elseif nargin ==5 M = 200; end if(w<=0 || w>=2) error; return; end D=diag(diag(A)); %求A的对角矩阵 L=-tril(A,-1); %求A的下三角阵 U=-triu(A,1); %求A的上三角阵 B=inv(D-L*w)*((1-w)*D+w*U); f=w*inv((D-L*w))*b; x=B*x0+f; n=1; %迭代次数 while norm(x-x0)>=eps x0=x; x =B*x0+f; n=n+1; if(n>=M) disp('Warning: 迭代次数太多，可能不收敛！'); return; end end |
5楼2009-10-22 20:04:57
6楼2009-10-22 20:06:30
★ ★ ★ ★
nono2009(金币+2,VIP+0):谢谢分享！ 10-22 20:27
最爱花诗雨(金币+1,VIP+0):谢谢~~~ 10-23 18:02
最爱花诗雨(金币+1,VIP+0):3Q 10-27 19:15
nono2009(金币+2,VIP+0):谢谢分享！ 10-22 20:27
最爱花诗雨(金币+1,VIP+0):谢谢~~~ 10-23 18:02
最爱花诗雨(金币+1,VIP+0):3Q 10-27 19:15
|
function [x,n]=SSOR(A,b,x0,w,eps,M) if nargin==4 eps= 1.0e-6; M = 200; elseif nargin<4 error return elseif nargin ==5 M = 200; end if(w<=0 || w>=2) error; return; end D=diag(diag(A)); %求A的对角矩阵 L=-tril(A,-1); %求A的下三角阵 U=-triu(A,1); %求A的上三角阵 B1=inv(D-L*w)*((1-w)*D+w*U); B2=inv(D-U*w)*((1-w)*D+w*L); f1=w*inv((D-L*w))*b; f2=w*inv((D-U*w))*b; x12=B1*x0+f1; x =B2*x12+f2; n=1; %迭代次数 while norm(x-x0)>=eps x0=x; x12=B1*x0+f1; x =B2*x12+f2; n=n+1; if(n>=M) disp('Warning: 迭代次数太多，可能不收敛！'); return; end end |
7楼2009-10-22 20:07:12