matlab plv,脑电脑网络分析代码使用流程介绍

%% 计算网络属性

clear all

% close all

clc

SaveFile = 'F:\kangkang_data\properties\';

mkdir(SaveFile);

file = 'F:\kangkang_data\cohdata\S_cohdata\Re_Coh0\f_a\';%计算plv属性

Sub_File = dir([file,'*.mat']);

for Sub = 1:length(Sub_File)%循环被试

load([file,Sub_File(Sub).name]);

% eval(['load',' ',file,'Re_Coh0_',num2str(Sub),'.mat']);

[m,n]=size(Re_Coh0);

CorrMatrix1_1 = Re_Coh0;

CorrMatrix1_1 = Re_Coh0-eye(m,n); %去掉主对角线上的值

% Re_Coh0_1 = Re_Coh0; %不把主对角线置零，cluster和E_loc值有差别

%聚类系数

Cluster_1 = mean(clustering_coef_wu(CorrMatrix1_1));

%路径长度

Matrix11 = ones(size(CorrMatrix1_1))-CorrMatrix1_1;

Matrix11(1:n+1:end) = 0;

[D11 B11] = distance_wei(Matrix11);

CharpNet_1 = sum(sum(D11))/(m*(m-1));

%全局效率

E_glo_1 = efficiency_wei(CorrMatrix1_1);

%局部效率

E_loc_1 = mean(efficiency_wei(CorrMatrix1_1,1));

% %度

% degree = degrees_und(CorrMatrix1_1);

properties_1(Sub,:)=[Cluster_1 CharpNet_1 E_glo_1 E_loc_1];

disp(['将要计算第',num2str(Sub),'个被试','_第',num2str(1),'种刺激！！！！'])

end

% properties_1([4,5,7,14,15,19],:) = [];

% properties_1([5,14],:) = [];

save([SaveFile,'properties_S_FA','.mat'],'properties_1');

% %%

% file = 'F:\kangkang_data\cohdata\F_cohdata\Re_Coh0\';

% Sub_File = dir([file,'*.mat']);

% sch_data=[];

% %%

% for Sub = 1:length(Sub_File)%循环被试

% load([file,Sub_File(Sub).name]);

% %for i = 1:5

% % eval(['load',' ',file,'Re_Coh0_',num2str(Sub),'.mat']);

% [m,n]=size(Re_Coh0);

% CorrMatrix1_1 = Re_Coh0;

% % CorrMatrix1_1 = plvMatrix1-eye(m,n); %去掉主对角线上的值

% %聚类系数

% Cluster_2 = mean(clustering_coef_wu(CorrMatrix1_1));

% %路径长度

% Matrix11 = ones(size(CorrMatrix1_1))-CorrMatrix1_1;

% Matrix11(1:n+1:end) = 0;

% [D11 B11] = distance_wei(Matrix11);

% CharpNet_2 = sum(sum(D11))/(m*(m-1));

% %全局效率

% E_glo_2 = efficiency_wei(CorrMatrix1_1);

% %局部效率

% E_loc_2 = mean(efficiency_wei(CorrMatrix1_1,1));

% properties_2(Sub,:)=[Cluster_2 CharpNet_2 E_glo_2 E_loc_2];

% disp(['将要计算第',num2str(Sub),'个被试','_第',num2str