| 查看: 1665 | 回复: 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)... 不知道出现什么问题请指教！ |
回复此楼» 猜你喜欢
自荐读博
已经有9人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
-
自然科学基金委宣布启动申请书“瘦身提质”行动
已经有4人回复
-
求个博导看看
已经有18人回复
linccnu
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 2436.5
- 红花: 1
- 帖子: 191
- 在线: 52.2小时
- 虫号: 2568402
- 注册: 2013-07-27
- 性别: GG
- 专业: 信息处理方法与技术
25楼2014-11-09 15:17:32
myyounger
至尊木虫 (文坛精英)
- 应助: 8 (幼儿园)
- 贵宾: 0.105
- 金币: 23288.5
- 散金: 500
- 红花: 144
- 帖子: 16619
- 在线: 484.6小时
- 虫号: 638546
- 注册: 2008-10-28
- 专业: 脑血管结构、功能异常及相
11楼2014-11-05 18:32:58
hytao2012
铁杆木虫 (正式写手)
木头虫子
- 应助: 53 (初中生)
- 金币: 6313.3
- 散金: 115
- 红花: 11
- 帖子: 479
- 在线: 205.8小时
- 虫号: 2050091
- 注册: 2012-10-08
- 性别: GG
- 专业: 固体力学
12楼2014-11-05 18:36:58
 |
13楼2014-11-05 21:32:05