作者: 净心、精心 时间: 2013-8-25 18:32
if n<=0
disp( ‘error n必须为正数 ’);
end
细心一点就好了 作者: 每根头发都失眠 时间: 2013-8-25 19:51
净心、精心 发表于 2013-8-25 18:32
if n
还是不行哦 还是报错 这是整个程序 帮我啊 GOD
function RI=ri(n)
n=ceil(n);
if n<=0,error('n必须为正数');
end
if n==0 || n==1,RI=0;
return;
end
times=1000;
scaler=[9 8 7 6 5 4 3 2 1 1/2 1/3 1/4 1/5 1/6 1/7 1/8 1/9];
A=zeros(n);
lamda=zeros(times,1);
for num=1:times
rank=ceil(17*rand(0));
for i=1:n
for j=i:n
A(i,j)=scaler(rank(i,j));
A(j,i)=1/A(i,j);
A(i,i)=1;
end
end
rigenvector=eig(A);
lamda(num)=max(rigenvector);
end