无标度网络 matlab,Matlab构建无标度网络的程序,运行出错,高手帮忙看怎么改!!

N=10000;m0=3;m=3;

adjacent_matrix=sparse(m0,m0);

for i=1:m0

for j=1:m0

if j==i

adjacent_matrix(i,j)=1;

end

end

end

adjacent_matrix=sparse(adjacent_matrix);

node_degree=zeros(1,m0+1);

node_degree(2:m0+1)=sum(adjacent_matrix);

for iter=4:N

iter

total_degree=2*m*(iter-4)+6;

cum_degree=cumsum(node_degree);

choose=zeros(1,m)

%%%

r1=rand(1)*total_degree;

for i=1:iter-1

if(r1>=cum_degree(i))&(r1

choose(1)=i;

break

end

end

choose(1)%%

r2=rand(1)*total_degree;

for i=1:iter-1

if(r2>=cum_degree(i))&(r2

choose(2)=i;

break

end

end

while choose(2)==choose(1)

r2=rand(1)*total_degree;

for i=1:iter-1

if(r2>=cum_degree(i))&(r2

choose(2)=i;

break

end

end

end

choose(2) %%%

r3=rand(1)*total_degree;

for i=1:iter-1

if(r3>=cum_degree(i))&(r3

choose(3)=i;

break

end

end

while(choose(3)==choose(1))|(choose(3)==choose(2))

r3=rand(1)*total_degree;

for i=1:iter-1

if(r3>=cum_degree(i))&(r3

choose(3)=i;

break

end

end

end

choose(3)%%%

for k=1:m

z=choose(k)

adjacent_matrix(z,iter)=1;

adjacent_matrix(iter,z)=1;

end

node_degree=zeros(1,iter+1);

node_degree(2:iter+1)=sum(adjacent_matrix);

end

save data adjacent_matrix;

??? Subscript indices must either be real positive integers or logicals.

Error in ==> d:\MATLAB6p5\work\Untitled3.m

On line 63  ==>         adjacent_matrix(z,iter)=1;

[本帖最后由 mooni 于 2009-4-22 10:50 编辑]