24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yan0415

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 1439.4
散金: 25
帖子: 591
在线: 99.4小时
虫号: 1378778
注册: 2011-08-24
专业: 泛函分析

[求助] jacobi迭代法算法出错，请高手帮忙

jacobi迭代法算法：
myjacobi.m
function [x,relerr,niter]=myjacobi(A,b,x0,tol,maxiter)
relerr=inf;
niter=1;
S=diag( diag(A) );
T=A-S;
while relerr >=tol &niter < maxiter,
x=S ( b-T*x0 );
relerr=norm( x-x0,inf )/norm( x,inf );
x0=x;
niter=niter+1;
end
输入相关数据后，metlab提示第七行x=S ( b-T*x0 );有问题：
??? Subscript indices must either be real positive integers or logicals.

Error in ==> myjacobi at 7
x=S ( b-T*x0 );
这是怎么会是哪？
谢谢！

回复此楼

» 猜你喜欢

职称评审没过，求安慰已经有21人回复
垃圾破二本职称评审标准已经有15人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有20人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
聘U V热熔胶研究人员已经有10人回复
求助文献已经有3人回复
投稿返修后收到这样的回复，还有希望吗已经有8人回复
三无产品还有机会吗已经有6人回复

1楼2012-12-18 17:56:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yan0415

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 1439.4
散金: 25
帖子: 591
在线: 99.4小时
虫号: 1378778
注册: 2011-08-24
专业: 泛函分析

引用回帖:

2楼: Originally posted by gotokenya at 2012-12-18 20:23:49
是x=S\( b-T*x0 )吧

我输入的数据是：
h6=hilb(6);
>> x6=ones(6,1);
>> b6=h6*x6;
>> x0=ones(6,1)+1;
>> [y,relerr,inter]=myjacobi(h6,b6,x0,1e-3,1000)
但metlab提示：
??? Error using ==> mrdivide
Matrix dimensions must agree.

Error in ==> myjacobi at 7
x=S/( b-T*x0 );
麻烦您再给看一下问题出在哪啦，谢谢！