| 查看: 2330 | 回复: 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)... 不知道出现什么问题请指教！ |
回复此楼» 猜你喜欢
现代”学阀”该如何界定
已经有14人回复
-
求合成方法
已经有8人回复
-
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有5人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有5人回复
-
有机合成以后会不会被AI改变？做科研的虫友怎么看
已经有8人回复
-
我的奶奶
已经有3人回复
 |
28楼2014-11-11 07:20:49
myyounger
至尊木虫 (文坛精英)
- 应助: 8 (幼儿园)
- 贵宾: 0.105
- 金币: 23288.5
- 散金: 500
- 红花: 144
- 帖子: 16619
- 在线: 484.6小时
- 虫号: 638546
- 注册: 2008-10-28
- 专业: 脑血管结构、功能异常及相
11楼2014-11-05 18:32:58
hytao2012
铁杆木虫 (正式写手)
木头虫子
- 应助: 53 (初中生)
- 金币: 6771.3
- 散金: 115
- 红花: 11
- 帖子: 568
- 在线: 208.8小时
- 虫号: 2050091
- 注册: 2012-10-08
- 性别: GG
- 专业: 固体力学
12楼2014-11-05 18:36:58
|
|
13楼2014-11-05 21:32:05