clear
clc
load('M:\HNAU\实验\酸牛奶荧光\处理数据\2018-04-13\特征谱\data.mat');
%标准化数据（Z标准化方法）
Std=std(s);M=mean(s);s=(s-repmat(M,size(s,1),1))./repmat(Std,size(s,1),1);
Std=std(t);M=mean(t);t=(t-repmat(M,size(t,1),1))./repmat(Std,size(t,1),1);
[COEF,SCORE latent]=princomp(s); %主成分得分，即原始数据在新生成的主成分空间里的坐标值,COEF的每一列代表一个新生成的主成分的系数
X=COEF(:,1:6);
x1=X(1:20,:);
x2=X(21:40,:);
x3=X(41:60,:);
x4=X(61:80,:);
x5=X(81:100,:);
xmean = mean(X);
xstd = std(X);
X = (X - xmean(ones(100,1),:))./xstd(ones(100,1),:);
%分类器
%class=ones(50,1);class(11:20)=2*class(10,1);class(21:30)=3*class(10,1);class(31:40)=4*class(10,1);class(41:50)=5*class(10,1);
class=ones(75,1);class(16:30)=2*class(15,1);class(31:45)=3*class(15,1);class(46:60)=4*class(15,1);class(61:75)=5*class(15,1);
% 将数据分为训练集和测试集
idxTrain=1:15;
idxTest=11:20;
P1 = (x1(idxTrain,:) - xmean(ones(15,1),:))./xstd(ones(15,1),:);
P2 = (x2(idxTrain,:)- xmean(ones(15,1),:))./xstd(ones(15,1),:);
P3 = (x3(idxTrain,:)- xmean(ones(15,1),:))./xstd(ones(15,1),:);
P4 = (x4(idxTrain,:)- xmean(ones(15,1),:))./xstd(ones(15,1),:);
P5 = (x5(idxTrain,:)- xmean(ones(15,1),:))./xstd(ones(15,1),:);
input1=[P1; P2; P3; P4; P5];
%input1=X;
[model,k,ClassLabel]=LDATraining(input1,class)
Q1 = (x1(idxTest,:)- xmean(ones(10,1),:))./xstd(ones(10,1),:);
Q2 = (x2(idxTest,:)- xmean(ones(10,1),:))./xstd(ones(10,1),:);
Q3 = (x3(idxTest,:)- xmean(ones(10,1),:))./xstd(ones(10,1),:);
Q4 = (x4(idxTest,:)- xmean(ones(10,1),:))./xstd(ones(10,1),:);
Q5 = (x5(idxTest,:)- xmean(ones(10,1),:))./xstd(ones(10,1),:);
input2=[Q1; Q2; Q3; Q4; Q5];
Class=LDATesting(input2,k,model,ClassLabel);
figure
h1=bar(1:10,Class(1:10),'b');
hold on
h2=bar(11:20,Class(11:20),'g');
h3=bar(21:30,Class(21:30),'r');
h4=bar(31:40,Class(31:40),'c');
h5=bar(41:50,Class(41:50),'m');
hold off
legend([h1 h2 h3 h4 h5],'A','B','C','D','E')
set(gca,'linewidth',2);
set(gca,'xtick',[0:10:50],'ytick',[0:1:5]);
xlabel('Sample number','FontName','Times New Roman','FontSize',14)
ylabel('Predicted class','FontName','Times New Roman','FontSize',14)
set(gca,'linewidth',2,'FontName','Times New Roman','FontSize',12)
