re.m

复杂网络建随机网模拟建网代码

clear all;
st=cputime;
clc;
disp('请不要进行其它操作,程序正在运行....');
head1='oldBA';
head2='newBA';
jindu='已完成';
jin='%....';
m=['1','2','3','4','5','6','7','8','9','0'];
mm='0';
h=waitbar(0,'程序正在运行，请等待。。。');
for i=1:10
    filename=[head1,m(i),'.txt'];
    fp=fopen(filename,'r');
    n=fscanf(fp,'%6d',1);
    [A,COUNT]=fscanf(fp,'%4d',[n,n]);
    fclose(fp);
    [M,N]=eig(A);
    alpha1(i)=N(1,1)/N(n-1,n-1);
    filename=[head2,m(i),'.txt'];
    fp=fopen(filename,'r');
    n=fscanf(fp,'%6d',1);
    [A,COUNT]=fscanf(fp,'%4d',[n,n]);
    fclose(fp);
    [M,N]=eig(A);
    alpha2(i)=N(1,1)/N(n-1,n-1);
    B=[jindu,m(i),mm,jin];
    disp(B);
    waitbar(i/10);
    
    
end
old=sum(alpha1)/i;
new=sum(alpha2)/i;
et=cputime-st;
fp=fopen('equal.txt','w');
fprintf(fp,'运行时间%f秒',et);
fprintf(fp,'未加权矩阵结果:');
fprintf(fp,'%12.8f',alpha1);
fprintf(fp,'平均值:');
fprintf(fp,'%12.8f',old);
fprintf(fp,'加权矩阵结果:');
fprintf(fp,'%12.8f',alpha2);
fprintf(fp,'平均值:');
fprintf(fp,'%12.8f',new);
fclose(fp);
alpha1
alpha2
old
new