xs24
2008-12-27, 15:34
%--myfun()
clear;clc
X=[16.68 27.75 31.84 18.4 53 55 29 175;
19.7 27.75 32.94 19.2 59 55 33 287;
15.2 23.4 32.98 16.24 46 65 42 153;
7.25 8.97 21.3 4.76 34 62 39 163;
29.45 56.49 40.74 43.68 75 69 27 214;
32.93 42.78 49.98 33.87 66 50 33 260;
25.39 37.85 36.76 27.56 68 63 36 243;
15.05 19.49 27.21 14.21 56 76 36 175;
19.82 28.78 33.41 20.17 59 71 39 183;
21.13 35.2 39.16 26.52 52 62 35 173;
16.75 28.72 29.62 19.23 55 58 30 152;
15.83 28.03 26.4 17.43 61 61 33 160];
%由于原始数据单位不一致,需要进行标准化处理
STX=std(X);
MX=mean(X);
[n,m]=size(X);
for ii=1:1:m
SX(:,ii)=(X(:,ii)-MX(ii))./STX(ii);
end
[PC,SCORE,LATENT]=princomp(SX);
[num2str(m),'个综合因子PC(按列放置)’]
PC
[num2str(n),'个样本的Z得分(按行放置)’]
SCORE
[num2str(m),'个综合因子的特征值(按列放置)’]
LATENT'
[num2str(m),'个综合因子的方差占总方差百分比%(按列放置)’]
explained=100*LATENT'/sum(LATENT')
%----将n个样本的分类进行图形表达
figure(1);plot3(SCORE(:,1),SCORE(:,2),SCORE(:,3),'*');grid on;
xlabel('第一主成分’);ylabel('第二主成分’);zlabel('第三主成分’);title('三维样本分类图);
for ii=1:1:n;text(SCORE(:,1),SCORE(:,2),SCORE(:,3),['工厂’,num2str(ii)]);end
figure(2);plot(SCORE(:,1),SCORE(:,2),'*');grid on;
xlabel('第一主成分’);ylabel('第二主成分’);title('二维样本分类图’);
for ii=1:1:n;text(SCORE(:,1),SCORE(:,2),['工厂’,num2str(ii)]);end
运行不了,请问高手们,怎么改?谢谢:sweat:
clear;clc
X=[16.68 27.75 31.84 18.4 53 55 29 175;
19.7 27.75 32.94 19.2 59 55 33 287;
15.2 23.4 32.98 16.24 46 65 42 153;
7.25 8.97 21.3 4.76 34 62 39 163;
29.45 56.49 40.74 43.68 75 69 27 214;
32.93 42.78 49.98 33.87 66 50 33 260;
25.39 37.85 36.76 27.56 68 63 36 243;
15.05 19.49 27.21 14.21 56 76 36 175;
19.82 28.78 33.41 20.17 59 71 39 183;
21.13 35.2 39.16 26.52 52 62 35 173;
16.75 28.72 29.62 19.23 55 58 30 152;
15.83 28.03 26.4 17.43 61 61 33 160];
%由于原始数据单位不一致,需要进行标准化处理
STX=std(X);
MX=mean(X);
[n,m]=size(X);
for ii=1:1:m
SX(:,ii)=(X(:,ii)-MX(ii))./STX(ii);
end
[PC,SCORE,LATENT]=princomp(SX);
[num2str(m),'个综合因子PC(按列放置)’]
PC
[num2str(n),'个样本的Z得分(按行放置)’]
SCORE
[num2str(m),'个综合因子的特征值(按列放置)’]
LATENT'
[num2str(m),'个综合因子的方差占总方差百分比%(按列放置)’]
explained=100*LATENT'/sum(LATENT')
%----将n个样本的分类进行图形表达
figure(1);plot3(SCORE(:,1),SCORE(:,2),SCORE(:,3),'*');grid on;
xlabel('第一主成分’);ylabel('第二主成分’);zlabel('第三主成分’);title('三维样本分类图);
for ii=1:1:n;text(SCORE(:,1),SCORE(:,2),SCORE(:,3),['工厂’,num2str(ii)]);end
figure(2);plot(SCORE(:,1),SCORE(:,2),'*');grid on;
xlabel('第一主成分’);ylabel('第二主成分’);title('二维样本分类图’);
for ii=1:1:n;text(SCORE(:,1),SCORE(:,2),['工厂’,num2str(ii)]);end
运行不了,请问高手们,怎么改?谢谢:sweat: