数学建模社区-数学中国

标题: 为什么一直出错呀。。。。那位高手帮帮忙。。 [打印本页]

作者: gh.sha 时间: 2010-12-17 11:11
标题: 为什么一直出错呀。。。。那位高手帮帮忙。。
function [A,nx,ny]=aaaa(n)
ny=n;
nx=3*(ny-1)+1;
A=zeros(nx-1,ny-1);
for i=1:(nx-1)*(ny-1)
p=mod(i,nx-1);
q=fix(i./(nx-1))+1;
if 0<p<ny && p<q<ny
      k=q+(p-1)*(nx-1);
   A(i,i)=1;
   A(i,k)=-1;
end
end

作者: gh.sha 时间: 2010-12-17 11:12
高手在哪里呀？？？

作者: gh.sha 时间: 2010-12-17 11:12

作者: gh.sha 时间: 2010-12-17 11:13
那位高手帮帮忙喽。。。。

作者: madio 时间: 2010-12-17 11:15
提示哪里有错？

作者: gh.sha 时间: 2010-12-17 11:16
??? Attempted to access A(60,-58); index must be a positive integer or logical.

Error in ==> aaaa at 10
A(i,q+(p-1)*(nx-1))=-1;

作者: gh.sha 时间: 2010-12-17 11:16
运行就出问题。。。。

作者: gh.sha 时间: 2010-12-17 11:21
高手在哪里呀？高手在哪里？

作者: gh.sha 时间: 2010-12-17 11:22
高手在哪里呀？高手在哪里？

作者: gh.sha 时间: 2010-12-17 11:23
高手在哪里呀？高手在哪里？

作者: gh.sha 时间: 2010-12-17 11:27
回复 madio 的帖子
??? Attempted to access A(60,-58); index must be a positive integer or logical.
Error in ==> aaaa at 10
A(i,q+(p-1)*(nx-1))=-1;
运行就出错。。。很无语。。

作者: 红速手 时间: 2010-12-17 11:44
我也不是很懂，提供点小意见供参考，A为一个25*26的矩阵，在这矩阵中是不可能找到A(60，-58)这个数的，所以说你下面的运算关系存在错误，矩阵的标号都为正。

作者: gh.sha 时间: 2010-12-17 11:48
回复 红速手 的帖子
改了，还是不对。。
function [A,nx,ny]=aaaa(n)
ny=n;
nx=3*(ny-1)+1;
A=zeros((nx-1)*(ny-1),(nx-1)*(ny-1));
for i=1

nx-1)*(ny-1)
p=mod(i,nx-1);
q=fix(i./(nx-1))+1;
if 0<p<ny && p<q<ny
   k=q+(p-1)*(nx-1);
   A(i,i)=1;
   A(i,k)=-1;
end
end

作者: 小草远在天涯 时间: 2010-12-17 12:06
??

作者: 小草远在天涯 时间: 2010-12-17 12:08
>> aaaa
??? Input argument 'n' is undefined.

Error in ==> F:\mat\work\aaaa.m
On line 2 ==> ny=n;

我的运行怎么错误是这样？

作者: gh.sha 时间: 2010-12-17 12:28
回复 小草远在天涯
>>aaaa(21)
应该是这样

作者: gh.sha 时间: 2010-12-17 12:44
走过路过的高手帮帮忙喽。。。。

作者: linmatsas 时间: 2010-12-17 15:10
0<p<ny那个程序可以用这样的符号的………………肯定分开写呀………………

作者: _math 时间: 2010-12-17 15:16
function [A,nx,ny]=aaaa(n)
ny=n;
nx=3*(ny-1)+1;
A=zeros((nx-1)*(ny-1),(nx-1)*(ny-1));
for i=1 : (nx-1)*(ny-1)
p=mod(i,nx-1);
q=fix(i/(nx-1))+1;
if 0<p&p<ny & p<q & q<ny
   k=q+(p-1)*(nx-1);
   A(i,i)=1;
   A(i,k)=-1;
end
end
楼主貌似if语句那里写错了。。。
我的可以运行的。。
aaaa(21)
运行后矩阵中有190个1和-1

作者: 伤感的幸福 时间: 2010-12-17 16:54
应该有问题啊
实现什么功能啊

作者: yujianyuanhaha 时间: 2010-12-17 22:26
真正的高手是不上网的···

作者: 柯雷 时间: 2010-12-18 23:39
dingyixiaxiaxiaxia

作者: zsy312 时间: 2010-12-21 21:11
function [A,nx,ny]=f(n)
ny=n;
nx=3*(ny-1)+1;
A=zeros(nx-1,ny-1);
for i=1: (nx-1)*(ny-1)
p=mod(i,nx-1);
q=fix(i./(nx-1))+1;
if 0<p&&p<ny && p<q&&q<ny
      k=q+(p-1)*(nx-1);
      A(i,i)=1;
      A(i,k)=-1;
end
end

作者: gh.sha 时间: 2010-12-26 22:27
回复 _math 的帖子
高手呀！！！

作者: 水木年华zzu 时间: 2010-12-27 12:19
查查这个值q+(p-1)*(nx-1)，出现了负数 p=mod(i,nx-1);
q=fix(i./(nx-1))+1;这两句前面是不是要加上ABS

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)