24小时热门版块排行榜    

查看: 1063  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lijunjiexm

铁虫 (小有名气)

[求助] 求助一简单matlab程序的出错问题,谢谢

如下m文件,为什么在命令窗口运行 出现错误提示:
>> A=mesh2d2(10)
??? Index exceeds matrix dimensions.
Error in ==> mesh2d2 at 8
A = sparse (ii(keep), jj(keep), xx(keep));
>>
请大家指点,谢谢!

%%%%%%%%%%%%%%%%%%%%%%%
function A = mesh2d2 (n)
% create an n-by-n 2D mesh for the 2nd difference operator
nn = 1:n^2 ;
ii = [nn-n ; nn-1 ; nn ; nn+1 ; nn+n] ;
jj = repmat (nn, 5, 1) ;
xx = repmat ([-1 -14-1-1]',1, n^2) ;
keep = find (ii >= 1 & ii <= n^2 & jj >= 1 & jj <= n^2) ;
A = sparse (ii(keep), jj(keep), xx(keep));
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖

★ ★
余泽成(金币+2): 谢谢参与应助! 2011-11-26 23:43:22
??? Index exceeds matrix dimensions.
下标越界错误

在mesh2d2中的第八行之前,第七行之后的分号去掉
运行一次,把keep这个下标打印出来看,是不是没有find到,keep是0啊
导致在做稀疏矩阵A的时候ii(keep)访问矩阵的0号元素(不存在),因为matlab的矩阵下标从1开始的
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
2楼2011-11-26 17:53:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lijunjiexm 的主题更新
信息提示
请填处理意见