面板数据聚类分析MATLAB代码.TXT

clc
clear all


%% 导入电子表格中的数据
% 用于从以下电子表格导入数据的脚本:
%
%    工作簿: E:\Fw_ wslsldmm\面板聚类数据.xlsx 工作表: 工作表1
%
% 要扩展代码以供其他选定数据或其他电子表格使用，请生成函数来代替脚本。

% 由 MATLAB 自动生成于 2020/01/10 18:46:02

% %% 导入数据
% [~, ~, raw] = xlsread('E:\Fw_ wslsldmm\面板聚类数据.xlsx','工作表1','C2:C28');
% raw(cellfun(@(x) ~isempty(x) && isnumeric(x) && isnan(x),raw)) = {''};
% cellVectors = raw(:,1);
% 
% %% 创建表
% B = table;
% 
% %% 将导入的数组分配给列变量名称
% B.area = cellVectors(:,1);
% 
% %% 清除临时变量
% clearvars area raw cellVectors;
A=xlsread('E:\Fw_ wslsldmm\面板聚类数据.xlsx','D2:K325');
%27个地区 12个年份 8个变量,将其转化为截面数据
for g=1:12
    for m=1:27

        for l=1:8
            DD(m,l,g)=A(27*(g-1)+m,l);
         end
    end
end
%% 计算每个时间节点的主成分
for i=1:12
     %得到的数据矩阵的行数和列数
     a=size(DD(:,:,1),1);
     b=size(DD(:,:,1),2);
    %计算系数矩阵:CM
    CM(:,:,i)=corrcoef(DD(:,:,i));
    %计算CM的特征值和特征向量
    [V(:,:,i),D(:,:,i)]=eig(CM(:,:,i));
    %将特征值按降序排列到DS中
     for j=1:b
         DS(j,1,i)=D(b+1-j,b+1-j,i);
     end
     %计算贡献率
     for t=1:b
         DS(t,2,i)=DS(t,1,i)/sum(DS(:,1,i));%单个贡献率
         DS(t,3,i)=sum(DS(1:t,1,i))/sum(DS(:,1,i));%累计贡献率
     end
     %假定主成分的信息保留率
            T=1;
            for k=1:b
                if DS(k,3,i) >= T
                com_num=k;
                 break;
                end