求助，MATLAB程序运行时出错，怎么改？

langia · 发表于 2008-6-2 17:32:09

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

您需要登录才可以下载或查看，没有账号？立即加入

×

各位前辈们，有谁特别熟悉MATLAB软件？帮我看下这个程序，为什么运行的时候，总是提示我这个错误：
Warning: Divide by zero.
> In lllllllll at 17
Warning: Divide by zero.
> In lllllllll at 18
??? Attempted to access G(0,0); index must be a positive integer or logical.
Error in ==> lllllllll at 29
G(i,i)=G(i,i)+gij/(k^2);

程序如下：
%牛顿——拉夫逊法进行潮流计算
n=input('请输入节点数：n=')
nl=input('请输入支路数：nl=')
isb=input('请输入平衡母线节点号：isb=')
A=input('请输入由支路参数形成的矩阵： A=')
A=zeros(nl,5);
G=zeros(n)
B=zeros(n)
for a=1:nl
i=A(a,1);
j=A(a,2);
r=A(a,3);
x=A(a,4);
k=A(a,5);
z=r*r+x*x;
gij=r/z;
bij=-x/z;
   if(k<0)
      k=abs(k);
      G(i,i)=G(i,i)+gij;
      B(i,i)=B(i,i)+bij+k;
      G(j,j)=G(j,j)+gij;
      B(j,j)=B(j,j)+bij+k;
      G(i,j)=-gij;
      B(i,j)=-bij;
elseif(k==0)
k=1;
G(i,i)=G(i,i)+gij/(k^2);
B(i,i)=B(i,i)+bij/(k^2);
G(j,j)=G(j,j)+gij;
B(j,j)=B(j,j)+bij;
G(i,j)=-gij/k;
B(i,j)=-bij/k;
elseif(k>0)
G(i,i)=G(i,i)+gij/(k*k);
B(i,i)=B(i,i)+bij/(k*k);
      G(j,j)=G(j,j)+gij;
      B(j,j)=B(j,j)+bij;
      G(i,j)=-gij/k;
B(i,j)=-bij/k;
end
if(i>j)
         G(j,i)=G(i,j);
         B(j,i)=B(i,j);
end
end
for a=1:n
   for b=1:n
      if(a<b)
         G(b,a)=G(a,b);
         B(b,a)=B(a,b);
      end
   end
end
n=input('请输入节点数：n=')
nl=input('请输入支路数：nl=')
isb=input('请输入平衡母线节点号：isb=')
g=0%迭代次数
e=input('请输入各节点参数形成的矩阵： e=')%%%%%%%%%%原始数据输入
f=input('请输入各节点参数形成的矩阵： f=')
Pis=input('请输入各节点参数形成的矩阵： Pis=')
Qis=input('请输入各节点参数形成的矩阵： Qis=')
e=zeros(n,1)
f=zeros(n,1)
Pis=zeros((n-1),1)
Qis=zeros((n-1),1)
max=1
while(max>1.0e-5)
A=input('请输入由支路参数形成的矩阵： A=')
G=zeros(n,n)
B=zeros(n,n)
for a=1:nl
i=A(a,1)
j=A(a,2)
r=A(a,3)
x=A(a,4)
k=A(a,5)
z=r*r+x*x
gij=r/z
   bij=-x/z
   if(k<0)
      k=abs(k);
      G(i,i)=G(i,i)+gij;
      B(i,i)=B(i,i)+bij+k;
      G(j,j)=G(j,j)+gij;
      B(j,j)=B(j,j)+bij+k;
      G(i,j)=-gij;
         B(i,j)=-bij;
elseif(k==0)
k=1;
G(i,i)=G(i,i)+gij/(k*k);
B(i,i)=B(i,i)+bij/(k*k);
      G(j,j)=G(j,j)+gij;
      B(j,j)=B(j,j)+bij;
      G(i,j)=-gij/k;
B(i,j)=-bij/k;
elseif(k>0)
G(i,i)=G(i,i)+gij/(k*k);
B(i,i)=B(i,i)+bij/(k*k);
      G(j,j)=G(j,j)+gij;
      B(j,j)=B(j,j)+bij;
      G(i,j)=-gij/k;
B(i,j)=-bij/k;
end
if(i>j)
         G(j,i)=G(i,j);
         B(j,i)=B(i,j);
end
end
for a=1:n
   for b=1:n
      if(a<b)
         G(b,a)=G(a,b);
         B(b,a)=B(a,b);
      end
   end
end
这是求节点导纳矩阵。请问，为什么我输进去参数之后，老是提示那个错误？？
还有，定义矩阵的时候，定义完之后要不要加;号？？到底怎么区分？？急死了。
大家快来帮帮我。。。。。。。。。。。。。。。。。。
先谢谢啦。。。。。。。

qiuwei_002 · 发表于 2008-6-2 17:54:57

能否提供一下参数？

sadever · 发表于 2008-6-2 19:14:17

A=input('请输入由支路参数形成的矩阵： A=')" i9 a, L# ^- E) y1 j5 e E/ Y
A=zeros(nl,5);2 m- H+ g% m- c2 Y
这两句是不是有问题，既然赋值了，为什么又给零矩阵。
这两个错误是很简单的错误，一个是除零，一个是矩阵下标为零，你只要给个断点就可以了。

hhusunguoqiang · 发表于 2008-6-18 10:43:55

楼上的说的很对，做除法的时候有0元素出现！

账号		自动登录	找回密码
密码			立即加入

求助，MATLAB程序运行时出错，怎么改？

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

相关帖子

浏览过的版块