利用聚类算法划分风电出力的典型日

clc
clear
close all
load('P_WindY5');
figure
plot(P_WindY5);
xlabel('时刻节点/5min');
ylabel('功率/MW');
title('风电出力年曲线图');
P_WindY365 = zeros(365,24*12);
for i=1:365
    P_WindY365(i,1:24*12) = P_WindY5(1,24*12*i-24*12+1:24*12*i);
end
%继而是调用K-MEANS聚类，得出8天的典型日                
%% 分区1的kmeans风电聚类
[IDW,~,~,~] = kmeans( P_WindY365 , 8 );         
E_Wind1=find(IDW==1);
E_Wind1=P_WindY365(E_Wind1,:);           
E_Wind2=find(IDW==2);
E_Wind2=P_WindY365(E_Wind2,:);       
E_Wind3=find(IDW==3);
E_Wind3=P_WindY365(E_Wind3,:);
E_Wind4=find(IDW==4);
E_Wind4=P_WindY365(E_Wind4,:);
E_Wind5=find(IDW==5);
E_Wind5=P_WindY365(E_Wind5,:);
E_Wind6=find(IDW==6);
E_Wind6=P_WindY365(E_Wind6,:);
E_Wind7=find(IDW==7);
E_Wind7=P_WindY365(E_Wind7,:);
E_Wind8=find(IDW==8);
E_Wind8=P_WindY365(E_Wind8,:);
figure(2)
subplot(2,4,1);
plot(E_Wind1(:,1:24*12)','-','LineWidth',1.5);xlabel('时间/h');ylabel('风电场群出力聚类W/m2');title(['典型日1频次',num2str(numel(E_Wind1(:,1)))]);axis([1 24*12 0 1600]);
subplot(2,4,2);
pl