%蔬菜油剔除2，7，8數(shù)據(jù)標(biāo)準(zhǔn)化后主成分分析

clc;
clear;
a =[1.0000    9.7000    5.2000   31.0000   52.7000    0.4000    0.4000    0.1000
    1.0000   11.1000    5.0000   32.9000   49.8000    0.3000    0.4000    0.1000
    1.0000   11.5000    5.2000   35.0000   47.2000    0.2000    0.4000    0.1000
    1.0000   10.0000    4.8000   30.4000   53.5000    0.3000    0.4000    0.1000
    1.0000   12.2000    5.0000   31.1000   50.5000    0.3000    0.4000    0.1000
    1.0000    9.8000    4.2000   43.0000   39.2000    2.4000    0.4000    0.5000
    1.0000   10.5000    5.0000   31.8000   51.3000    0.4000    0.4000    0.1000
    1.0000   10.5000    5.0000   31.8000   51.3000    0.4000    0.4000    0.1000
    1.0000   11.5000    5.2000   35.0000   47.2000    0.2000    0.4000    0.1000
    1.0000   10.0000    4.8000   30.4000   53.5000    0.3000    0.4000    0.1000
    1.0000   11.1000    5.0000   32.9000   49.8000    0.3000    0.4000    0.1000
    1.0000    9.3000    4.4000   43.3000   39.2000    2.3000    0.4000    0.5000
    1.0000   12.2000    5.0000   31.1000   50.5000    0.3000    0.4000    0.1000
    1.0000    9.7000    5.6000   35.2000   47.8000    0.5000    0.4000    0.2000
    1.0000   11.1000    4.6000   28.3000   50.6000    4.2000    0.4000    0.2000
    1.0000   11.0000    4.7000   29.8000   49.6000    3.4000    0.4000    0.2000
    1.0000   11.6000    6.0000   35.4000   45.9000    0.2000    0.4000    0.1000
    1.0000    9.8000    5.3000   31.7000   51.3000    0.8000    0.4000    0.2000
    2.0000    6.1000    4.1000   24.0000   64.3000    0.1000    0.3000    0.1000
    2.0000    6.2000    3.9000   27.1000   59.9000    1.3000    0.3000    0.3000
    2.0000    6.0000    4.9000   22.8000   64.4000    0.3000    0.3000    0.2000
         0    6.7000    3.4000   35.2000   49.3000    3.6000    0.4000    0.5000
         0    3.8000    3.9000   77.8000   12.4000    0.2000    0.3000    0.2000
         0   10.1000    4.7000   26.1000   52.6000    5.2000    0.4000    0.2000
         0    4.9000    2.2000   58.2000   25.2000    7.1000    0.6000    1.2000
         0    6.1000    2.7000   45.7000   37.4000    5.9000    0.5000    0.8000
    2.0000    7.2000    4.5000   25.6000   61.1000    0.2000    0.3000    0.2000
         0    5.9000    2.4000   55.5000   27.7000    6.2000    0.6000    1.1000
         0   10.9000    3.8000   25.3000   53.0000    5.6000    0.4000    0.2000
    2.0000    6.2000    4.1000   26.8000   61.4000    0.1000    0.3000    0.2000
    2.0000    6.1000    4.0000   25.2000   63.2000    0.1000    0.2000    0.2000
    5.0000   10.9000    3.6000   26.0000   52.6000    5.5000    0.4000    0.2000
    5.0000    9.6000    3.5000   30.3000   49.2000    5.9000    0.4000    0.3000
    7.0000   10.0000    2.3000   36.9000   47.1000    2.2000    0.5000    0.5000
    6.0000    5.1000    2.3000   55.9000   27.4000    6.8000    0.5000    0.5000
    3.0000    9.7000    3.4000   59.3000   20.5000    0.1000    1.5000    1.2000
         0    5.2000    2.0000   59.3000   22.7000    8.2000    0.6000    1.3000
         0    8.3000    3.5000   33.6000   45.6000    6.3000    0.4000    0.6000
    2.0000    6.1000    4.1000   26.7000   61.0000    0.6000    0.3000    0.2000
    2.0000    6.2000    4.0000   25.8000   62.2000    0.4000    0.3000    0.2000
    2.0000    6.0000    3.8000   29.6000   57.7000    1.2000    0.3000    0.3000
    2.0000    6.5000    2.8000   23.2000   66.1000    0.1000    0.2000    0.2000
    2.0000    7.0000    3.4000   23.2000   64.9000    0.1000    0.2000    0.2000
    2.0000    6.0000    4.0000   28.3000   60.1000    0.1000    0.3000    0.2000
    1.0000   11.2000    6.2000   37.7000   43.6000    0.2000    0.5000    0.2000
    1.0000    7.6000    4.7000   29.1000   56.1000    0.7000    0.4000    0.2000
    1.0000    7.7000    4.7000   31.1000   54.3000    0.8000    0.4000    0.2000
    1.0000    9.9000    5.2000   37.4000   44.2000    2.1000    0.4000    0.3000
    1.0000   11.5000    5.1000   27.8000   54.5000    0.2000    0.4000    0.1000
    1.0000   11.3000    5.8000   35.2000   46.7000    0.2000    0.4000    0.1000
    4.0000   14.9000    2.6000   68.2000   12.8000    0.6000    0.4000    0.3000
    4.0000    9.3000    2.8000   65.0000   17.0000    3.9000    0.5000    0.7000
    5.0000   10.5000    4.2000   25.5000   52.0000    7.8000    0.0500    0.0500
    5.0000   10.0000    4.2000   24.9000   53.2000    3.9000    0.4000    0.0500
    5.0000   10.4000    4.2000   25.9000   50.8000    7.5000    0.4000    0.4000
    5.0000   10.5000    4.2000   24.4000   52.1000    7.5000    0.4000    0.0500
    5.0000   10.5000    4.3000   24.6000   53.1000    7.6000    0.0500    0.0500
    5.0000   10.2000    4.0000   23.1000   55.1000    7.1000    0.5000    0.0500
    5.0000   10.9000    3.8000   27.2000   49.5000    6.4000    0.7000    0.9000
    5.0000   11.9000    3.8000   25.7000   52.7000    5.8000    0.0500    0.0500
         0    5.1000    2.0000   57.2000   26.1000    6.9000    0.7000    1.4000
         0    4.8000    1.9000   56.6000   25.7000    8.1000    0.8000    1.5000
         0   10.9000    3.2000   25.3000   53.2000    5.8000    0.4000    0.3000
         0   11.1000    3.8000   24.4000   52.8000    6.0000    0.4000    0.3000
         0   10.4000    4.8000   24.5000   51.8000    5.7000    0.5000    0.5000
         0   11.0000    4.3000   22.9000   54.1000    7.6000    0.0500    0.0500
         0   10.8000    4.3000   22.7000   54.6000    7.7000    0.0500    0.0500
         0   10.8000    4.3000   23.1000   54.1000    7.6000    0.0500    0.0500
         0    8.0000    4.3000   32.1000   44.8000    6.4000    1.5000    1.8000
         0    5.9000    5.2000   24.7000   62.7000    0.5000    1.1000    0.0500
         0    4.1000    4.3000   26.1000   62.8000    0.5000    0.0500    0.0500
         0    3.9000    3.3000   80.6000   11.3000    0.1000    0.0500    0.0500
         0    5.1000    2.0000   58.0000   25.0000    6.9000    0.7000    1.4000
         0    5.0000    1.8000   59.1000   23.2000    8.0000    0.7000    1.4000
    4.0000   10.9000    2.7000   76.7000    7.9000    0.8000    0.0500    0.0500
    4.0000   10.5000    2.8000   75.8000    8.0000    0.7000    0.0500    0.0500
    4.0000   12.0000    2.7000   75.1000    8.5000    0.8000    0.0500    0.0500
    4.0000   11.7000    2.9000   74.6000   10.1000    0.6000    0.0500    0.0500
    4.0000   11.4000    3.0000   73.0000   10.6000    0.7000    0.0500    0.0500
    6.0000    4.8000    1.8000   62.6000   20.0000    9.5000    0.0500    1.4000
    6.0000    5.5000    1.7000   59.0000   21.3000    9.3000    0.6000    1.5000
    6.0000    5.1000    1.9000   59.2000   22.3000    9.3000    0.7000    1.6000
    6.0000    4.8000    1.9000   61.6000   20.9000    8.0000    0.8000    1.5000
         0    6.2000    4.1000   31.1000   56.7000    0.8000    0.3000    0.0500
         0    6.2000    4.1000   31.1000   57.3000    0.5000    0.8000    0.0500
         0    6.7000    5.5000   24.1000   60.8000    0.9000    0.7000    0.7000
         0    8.3000    4.0000   34.0000   52.4000    0.3000    0.4000    0.0500
    3.0000   10.0000    3.3000   60.0000   21.3000    0.2000    1.5000    1.3000
    3.0000    9.6000    3.3000   57.7000   20.7000    0.2000    1.5000    1.8000
    1.0000   12.2000    5.4000   29.4000   53.0000    1.0000    0.0500    0.0500
    1.0000   11.0000    5.3000   35.0000   45.2000    1.3000    1.3000    0.7000
    1.0000   11.9000    5.6000   33.6000   48.9000    1.0000    0.0500    0.0500
    1.0000   11.4000    5.8000   34.5000   48.3000    1.0000    0.0500    0.0500
    1.0000   10.7000    5.4000   39.3000   43.2000    1.4000    0.0500    0.0500
    1.0000   11.2000    6.2000   35.8000   44.1000    0.7000    1.1000    0.0500
    1.0000   11.4000    5.8000   33.9000   48.1000    0.8000    0.8000    0.0500
    1.0000   11.5000    6.2000   39.7000   42.6000    0.8000    0.0500    0.0500
    1.0000   13.0000    6.2000   25.8000   55.0000    0.8000    0.0500    0.0500
    1.0000   13.0000    6.7000   29.6000   50.8000    0.5000    0.0500    0.0500
    1.0000   13.1000    6.3000   26.5000   53.6000    0.5000    0.5000    0.0500
    1.0000   11.6000    6.5000   38.4000   42.8000    0.5000    0.7000    0.0500
    1.0000   13.1000    5.7000   31.7000   49.5000    0.6000    0.0500    0.0500
    2.0000    6.1000    4.1000   25.1000   63.5000    0.5000    0.0500    0.0500
         0    5.4000    4.2000   41.4000   48.2000    0.5000    0.0500    0.0500
         0    4.6000    1.7000   61.8000   20.3000    9.3000    0.7000    1.5000
    2.0000    6.3000    4.2000   27.4000   61.4000    0.8000    0.0500    0.0500
    2.0000    6.2000    4.2000   27.1000   61.8000    0.8000    0.0500    0.0500
    2.0000    6.2000    4.2000   27.0000   60.9000    0.5000    0.3000    0.0500
    2.0000    6.2000    4.0000   28.3000   59.7000    0.9000    0.0500    0.0500
    2.0000    5.6000    4.2000   25.7000   58.9000    1.7000    2.8000    0.9000
    2.0000    6.4000    3.9000   26.0000   63.7000    0.5000    0.0500    0.0500
    2.0000    6.8000    4.3000   26.4000   60.5000    1.3000    0.0500    0.0500
    2.0000    6.0000    4.4000   27.1000   60.9000    0.9000    0.0500    0.0500
    2.0000    6.4000    4.8000   25.3000   61.8000    1.0000    0.0500    0.0500
    2.0000    5.9000    4.5000   24.1000   61.7000    0.9000    0.6000    0.6000
    2.0000    6.2000    4.1000   29.9000   57.8000    1.3000    0.0500    0.0500
    2.0000    6.6000    4.7000   24.5000   62.8000    0.3000    0.4000    0.0500
    2.0000    6.4000    4.4000   24.4000   63.7000    0.4000    0.4000    0.0500
    6.0000    5.4000    2.0000   53.2000   28.9000    7.3000    0.6000    1.3000
    6.0000    5.1000    1.9000   59.2000   22.4000    9.3000    0.6000    1.5000
    7.0000   10.7000    1.8000   30.2000   55.5000    0.9000    0.5000    0.3000
         0    7.5000    5.1000   29.5000   56.3000    0.1000    0.0500    0.0500
         0    8.0000    3.1000   40.1000   39.6000    6.9000    0.6000    1.0000
         0    8.2000    3.4000   34.4000   45.8000    5.7000    0.6000    0.8000
         0   11.0000    3.8000   24.6000   52.5000    6.1000    0.5000    0.4000
         0   10.5000    4.1000   22.3000   53.6000    8.3000    0.4000    0.0500
         0    9.4000    4.7000   34.9000   49.1000    0.2000    0.0500    0.0500
         0    9.8000    4.9000   34.8000   48.4000    1.7000    0.5000    0.0500
    6.0000    4.5000    1.7000   64.9000   18.6000    8.3000    0.0500    0.0500
    6.0000    5.7000    2.1000   54.6000   26.8000    8.0000    0.0500    0.0500
    6.0000    6.2000    2.2000   52.2000   29.0000    8.0000    0.0500    0.0500
    5.0000    9.7000    3.9000   25.1000   54.2000    5.9000    0.0500    0.0500]
b=a(:,3:6)
mean(b)
u=ans
var(b)
w=ans
i=1;
while i<=4
    v(i)=w(i)^(1/4)
    i=i+1
end
i=1;
while i<=132
   b(i,1:4)-u
   bb(i,:)=ans
   i=i+1
end
i=1;
while i<=132
    vv(i,:)=v
    i=i+1
end;
z=bb./vv
[pc,score,latent,tsquare]=princomp(z)
s(:,2:5)=score
s(:,1)=a(:,1)
i=1;
while i<=132
if s(i,1)<1
plot(s(i,2),s(i,3),'g+');
elseif s(i,1)<2
plot(s(i,2),s(i,3),'c+');
elseif s(i,1)<3
plot(s(i,2),s(i,3),'r+');
elseif s(i,1)<4
plot(s(i,2),s(i,3),'b+');
elseif s(i,1)<5
plot(s(i,2),s(i,3),'m+');
elseif s(i,1)<6
plot(s(i,2),s(i,3),'y+');
elseif s(i,1)<7
plot(s(i,2),s(i,3),'k+');
elseif s(i,1)<8
plot(s(i,2),s(i,3),'k*');
end
i=i+1;
hold on;
end;
grid on