| 查看: 1569 | 回复: 31 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
大神看一下这个程序有什么问题已有3人参与
|
||
|
format long A=[10 -1 -2;-1 10 -2;-1 -1 5]; b=[7.2;8.3;4.2]; omega=1.15;X0=[0;0;0;0];eps=1e-5;N=300; x=SOR(A,b,omega,X0,eps,N); function x=SOR(A,b,omega,x0,eps,N) n=length(x0);x=ones(n,1);k=0; while k<=N for i=1:n x(i)=omega*(- A(i,1:i-1)*x(1:i-1)... - A(i,i+1:n)*x0(i+1:n)+b(i))/A(i,i)+(1-omega)*x0(i); end k=k+1; if norm(x-x0,inf)<eps break; end x0=x; end if k>N warning('算法超出最大迭代次数！'); else disp(['迭代次数= ',num2str(k)]) x end 运行结果是： Index exceeds matrix dimensions. Error in SOR (line 5) x(i)=omega*(- A(i,1:i-1)*x(1:i-1)... 不知道出现什么问题请指教！ |
回复此楼» 猜你喜欢
导师想让我从独立一作变成了共一第一
已经有8人回复
-
博士读完未来一定会好吗
已经有23人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
myyounger
至尊木虫 (文坛精英)
- 应助: 8 (幼儿园)
- 贵宾: 0.105
- 金币: 23288.5
- 散金: 500
- 红花: 144
- 帖子: 16619
- 在线: 484.6小时
- 虫号: 638546
- 注册: 2008-10-28
- 专业: 脑血管结构、功能异常及相
11楼2014-11-05 18:32:58
hytao2012
铁杆木虫 (正式写手)
木头虫子
- 应助: 53 (初中生)
- 金币: 6303.3
- 散金: 115
- 红花: 11
- 帖子: 479
- 在线: 205.4小时
- 虫号: 2050091
- 注册: 2012-10-08
- 性别: GG
- 专业: 固体力学
12楼2014-11-05 18:36:58
 |
13楼2014-11-05 21:32:05