asym_correlation_profile.m

可以用来产生复杂网络

Nstat=10;
bpd=2;

srand=dir_srand+dir_srand';

k2=full(sum(srand));
k2_max=max(k2);

b2_max=1+floor(bpd.*log10(k2_max));
b1_max=b2_max;

n_1_2_asym=zeros(b1_max+1,b2_max+1);
binned_srand;

n_1_2_asym_orig=n_1_2;
dir_generate_srand;
srand=dir_srand+dir_srand';

binned_srand;

aver_n_1_2_asym=n_1_2;
aver_sq_n_1_2_asym=n_1_2.^2;

for k=2:Nstat;
   k
   dir_generate_srand;
   srand=dir_srand+dir_srand';

   binned_srand;
   aver_n_1_2_asym=aver_n_1_2_asym+n_1_2;
   aver_sq_n_1_2_asym=aver_sq_n_1_2_asym+n_1_2.^2;
end;
aver_n_1_2_asym=aver_n_1_2_asym./Nstat;
aver_sq_n_1_2_asym=aver_sq_n_1_2_asym./Nstat;
err_n_1_2_asym=sqrt(aver_sq_n_1_2_asym-aver_n_1_2_asym.^2);


asym_ratio_1_2_asym=n_1_2_asym_orig./(aver_n_1_2_asym+0.0001.*(aver_n_1_2_asym==0));
dev_n_1_2_asym_orig=(n_1_2_asym_orig-aver_n_1_2_asym)./(err_n_1_2_asym+0.0001.*(aver_n_1_2_asym==0));


asym_ratio_1_2_asym(end,:)=asym_ratio_1_2_asym(end-1,:);
asym_ratio_1_2_asym(:,end)=asym_ratio_1_2_asym(:,end-1);
dev_n_1_2_asym_orig(end,:)=dev_n_1_2_asym_orig(end-1,:);
dev_n_1_2_asym_orig(:,end)=dev_n_1_2_asym_orig(:,end-1);

figure;  pcolor(asym_ratio_1_2_asym);  colorbar


shading interp; 
colorbar;

figure;  pcolor(dev_n_1_2_asym_orig);  



shading interp; 
colorbar;

num2str(asym_ratio_1_2_asym(1:end-1,1:end-1))
num2str(dev_n_1_2_asym_orig(1:end-1,1:end-1))