%%?數據讀取?生成隨機矩陣
clear;
list?=?importdata（‘list.txt‘）;??%%call?by?list{n}
list?=?char（list）;
data?=?load（“data.txt“）;
[m?n]?=?size（data）;
flog?=?fopen（‘log.txt‘‘w‘）;????%日志文件
fprintf（flog‘%s\r\n‘‘20180726最大負荷日?支路數據‘）;

%%?數據預處理
zero_num?=?0;???%?除零
fid?=?fopen（‘all_zero.txt‘‘w‘）;
new_data?=?zeros（1n）;
for?i?=?1:m
????if?sum（data（i:））?==?0????%檢測存儲全零數據
????????zero_num?=?zero_num+1;
????????fprintf（fid‘%s\r\n‘list（i:））;
????else
????????new_data（i-zero_num:）?=?data（i:）;
????????new_list（i-zero_num:）?=?list（i:）;?%#ok
????end
end
%disp（[‘濾除全零數據共‘num2str（zero_num）‘條‘]）;
fprintf（flog‘%s\r\n‘[‘濾除全零數據共‘num2str（zero_num）‘條‘]）;
fclose（fid）;
clear?data?list?zero_num;

rng?default;????%?取樣
[bootstatsub]?=?bootstrp（1@meannew_data）;
[m?n]?=?size（new_data）;


%%?%曲線平滑??異常值處理
inv_num=0;
for?i?=?1:m
????new_data（i:）?=?smooth（new_data（i:）10‘rlowess‘）;??%曲線平滑
????%moving?rlowess?rloess
????
????%new_data（i:）?=?zscore（new_data（i:））;???%歸一化處理
????[new_data（i:）~]?=?mapminmax（new_data（i:）01）;
????
????if?abs（max（new_data（i:））?-?min（new_data（i:）））?????????new_data（i:）?=?ones（1?n）;
????????inv_num?=?inv_num+1;
????end
end
%disp（[‘處理恒不變數據共‘num2str（inv_num）‘條‘]）;
fprintf（flog‘%s\r\n‘[‘處理恒不變數據共‘num2str（inv_num）‘條‘]）;
clear?inv_num?;

%%?自帶函數進行mean聚類
cluster?=?zeros（m1）;
for?K?=?1:sqrt（m/2）
????[cluster（:K）]?=?kmeans（new_dataK）;???%保存每次聚類結果‘replicate‘5?迭代次數
????%?[idxcenter]?=?kmeans（new_data?K）;
end
eva?=?evalclusters（new_datacluster‘DaviesBouldin‘）;??%DB參數評估最佳聚類數
[~DB_max]?=?max（eva.CriterionValues）;

fprintf（flog‘%f\r\n‘eva.CriterionValues）;
%?eva2?=?evalclusters（new_datacluster‘CalinskiHarabasz‘）;
%?[~CH_max]?=?max（eva2.CriterionValues）;
%?eva3?=?evalclusters（new_datacluster‘Silhouette‘）;
%?[~SH_max]?=?max（eva3.CriterionValues）;

pattern=zeros（m?n+1）;???????%結果存儲
pattern（:1:n）?=?new_data（::）;
%pattern（:1+n）?=?cluster（:?