clc

[ij]=find（c2>0.5）;%相關系數矩陣
a3=c2（find（c2>0.5））;
A=[ija3];
A1=A（:3）;
A2=A1‘;
totalnum=length（A2）;
x=linspace（-11100）;??%將最大最小區間分成20個等分點（19等分）然后分別計算各個區間的個數
yy=hist（A2x）;??%計算各個區間的個數
yy=yy/totalnum;?%計算各個區間的個數
bar（xyy）?%畫出概率密度分布圖
axis?tight;
xlabel（‘股票相關系數‘）
ylabel（‘股票相關系數概率密度‘）
figure;
%%%%%------構造鄰接矩陣--------------------%%%%
matsize=88;
mat=zeros（matsizematsize）;
%%%
k?=?1;
%%?here?k?is?the?cout;for?example?geshu（1）?is?the?cout?of?>?0.001
?
??for?j?=?-1:0.1:1
?????for?i=1:totalnum
???????if?A2（i）>j
??????????nrow=A（i1）;
??????????ncol=A（i2）;
??????????mat（nrowncol）=1;
??????????mat（ncolnrow）=1;??
??????????
???????????end
?????end?????
?????[bm]=sort（sum（mat））;%%%%%%%%%%%%%%%按度排序
????dele=m（79:88）;%先刪除度排在前十的節點
????mat（:dele）=[];
????mat（dele:）=[];
?????[sci?sizes]?=?scomponents（mat）;??????
???????size=sort（sizes‘descend‘）;
???????ge1（k）=max（size）;%最大連通子圖的數目
???????size（11）=0;
??????geshu1（k）=max（size）;%第二大連通子圖的數目
??????k?=?k?+1;
????mat=zeros（matsizematsize）;

??end
??i=?-1:0.1:1;
plot（ige1‘-r<‘‘LineWidth‘0.5‘MarkerSize‘2）;
hold?on;
plot（igeshu1‘--bp‘‘LineWidth‘0.5‘MarkerSize‘2）;
xlabel（‘閾值‘）ylabel（‘節點個數‘）;
legend（‘最大連通子圖節點個數‘‘第二大連通子圖節點個數‘）;



%?%%%%首先給定一個鄰接矩陣mat
%?[bm]=sort（sum（mat））;%%%%%%%%%%%%%%%按度排序
%?????dele=m（81:88）;%先刪除度排在前八的節點，也就是
%?????mat（:dele）=[];
%?????mat（dele:）=[];

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2020-01-08?20:55??matlab模擬復雜網絡攻擊源代碼\
?????目錄???????????0??2020-01-08?20:49??matlab模擬復雜網絡攻擊源代碼\47651506xuyigongji\
?????文件????????1584??2016-01-08?19:46??matlab模擬復雜網絡攻擊源代碼\47651506xuyigongji\xuyigongji.m
?????文件???????70656??2018-06-21?21:21??matlab模擬復雜網絡攻擊源代碼\complex.doc
?????文件????????4099??2018-06-21?21:16??matlab模擬復雜網絡攻擊源代碼\新建文本文檔.txt