function [bzh,  remarix, dwhpcs, var, per_ex] = kk1(X)
[m,n]=size(X);
% 繪制第一個圖形
subplot(1,2,1);
boxplot(X,0,'+',0); 
xlabel('變量值的范圍');
ylabel('五個指標（變量）') ;

%計算樣本的標準差
stdr=std(X,0,1);
%計算樣本各列的均值
meanr=mean(X);
%標準化原始數據 
sr=(X-meanr(ones(m,1),:))./stdr(ones(m,1),:);
%樣本相關矩陣
R=sr'*sr/(m-1);
 
%進行主成分分析
[pcs,newdata,variances,t2]=princomp(sr) ;
%方差貢獻率
percent_explained=100*variances/sum(variances); 
 
%繪制第二個圖形
subplot(1,2,2); 
pareto(percent_explained);
xlabel('主成分');
ylabel('方差占的比重(%)');
bzh = sr;
remarix = R;
dwhpcs = pcs;
var = variances;
per_ex = percent_explained;