iris_classification.asv

libsvm支持向量机

clear
close all
%Iris classification
% s=rgb2gray(imread('D:\MATLAB7\work\处理后1\44_1.bmp'));
% Iris_44(1,:)=feature_extract(s);
% s=rgb2gray(imread('D:\MATLAB7\work\处理后1\44_2.bmp'));
% Iris_44(2,:)=feature_extract(s);
% s=rgb2gray(imread('D:\MATLAB7\work\处理后1\44_3.bmp'));
% Iris_44(3,:)=feature_extract(s);
% s=rgb2gray(imread('D:\MATLAB7\work\处理后1\44_4.bmp'));
% Iris_44(4,:)=feature_extract(s);
% %***********************************************
% s=rgb2gray(imread('D:\MATLAB7\work\处理后1\51_1.bmp'));
% Iris_51(1,:)=feature_extract(s);
% s=rgb2gray(imread('D:\MATLAB7\work\处理后1\51_2.bmp'));
% Iris_51(2,:)=feature_extract(s);
% s=rgb2gray(imread('D:\MATLAB7\work\处理后1\51_3.bmp'));
% Iris_51(3,:)=feature_extract(s);
% s=rgb2gray(imread('D:\MATLAB7\work\处理后1\51_4.bmp'));
% Iris_51(4,:)=feature_extract(s);
% %***************************************************
% s=rgb2gray(imread('D:\MATLAB7\work\处理后1\83_1.bmp'));
% Iris_83(1,:)=feature_extract(s);
% s=rgb2gray(imread('D:\MATLAB7\work\处理后1\83_2.bmp'));
% Iris_83(2,:)=feature_extract(s);
% s=rgb2gray(imread('D:\MATLAB7\work\处理后1\83_3.bmp'));
% Iris_83(3,:)=feature_extract(s);
% s=rgb2gray(imread('D:\MATLAB7\work\处理后1\83_4.bmp'));
% Iris_83(4,:)=feature_extract(s);
% %*****************************************************
% s=rgb2gray(imread('D:\MATLAB7\work\处理后1\88_1.bmp'));
% Iris_88(1,:)=feature_extract(s);
% s=rgb2gray(imread('D:\MATLAB7\work\处理后1\88_2.bmp'));
% Iris_88(2,:)=feature_extract(s);
% s=rgb2gray(imread('D:\MATLAB7\work\处理后1\88_3.bmp'));
% Iris_88(3,:)=feature_extract(s);
% s=rgb2gray(imread('D:\MATLAB7\work\处理后1\88_4.bmp'));
% Iris_88(4,:)=feature_extract(s);
% %*******************************************
load D:\APPLICATION1\Matlab7\work\package_emd\code\Data1_3\1_1.mat
N1=3;
% s=Imodes(25,:,N1);
for i=21:30
    if i==21
        s=Imodes(21,:,N1);
    else
        s=[s Imodes(i,:,N1)];
    end
end
Iris_44(1,:)=s;
clear Imodes
load D:\APPLICATION1\Matlab7\work\package_emd\code\Data1_3\1_2.mat
for i=21:30
    if i==21
        s=Imodes(21,:,N1);
    else
        s=[s Imodes(i,:,N1)];
    end
end
% s=Imodes(25,:,N1);
Iris_44(2,:)=s;
clear Imodes
load D:\APPLICATION1\Matlab7\work\package_emd\code\Data1_3\1_3.mat
for i=21:30
    if i==21
        s=Imodes(21,:,N1);
    else
        s=[s Imodes(i,:,N1)];
    end
end
% s=Imodes(25,:,N1);
Iris_44(3,:)=s;
clear Imodes
load D:\APPLICATION1\Matlab7\work\package_emd\code\Data1_3\1_4.mat
for i=21:30
    if i==21
        s=Imodes(21,:,N1);
    else
        s=[s Imodes(i,:,N1)];
    end
end
% s=Imodes(25,:,N1);
Iris_44(4,:)=s;
clear Imodes
%***********************************************
load D:\APPLICATION1\Matlab7\work\package_emd\code\Data1_3\2_1.mat
for i=21:30
    if i==21
        s=Imodes(21,:,N1);
    else
        s=[s Imodes(i,:,N1)];
    end
end
% s=Imodes(25,:,N1);
Iris_51(1,:)=s;
clear Imodes
load D:\APPLICATION1\Matlab7\work\package_emd\code\Data1_3\2_2.mat
for i=21:30
    if i==21
        s=Imodes(21,:,N1);
    else
        s=[s Imodes(i,:,N1)];
    end
end
% s=Imodes(25,:,N1);
Iris_51(2,:)=s;
clear Imodes
load D:\APPLICATION1\Matlab7\work\package_emd\code\Data1_3\2_3.mat
for i=21:30
    if i==21
        s=Imodes(21,:,N1);
    else
        s=[s Imodes(i,:,N1)];
    end
end
% s=Imodes(25,:,N1);
Iris_51(3,:)=s;
clear Imodes
load D:\APPLICATION1\Matlab7\work\package_emd\code\Data1_3\2_4.mat
for i=21:30
    if i==21
        s=Imodes(21,:,N1);
    else
        s=[s Imodes(i,:,N1)];
    end
end
% s=Imodes(25,:,N1);
Iris_51(4,:)=s;
clear Imodes
%***************************************************
load D:\APPLICATION1\Matlab7\work\package_emd\code\Data1_3\3_1.mat
for i=21:30
    if i==21
        s=Imodes(21,:,N1);
    else
        s=[s Imodes(i,:,N1)];
    end
end
% s=Imodes(25,:,N1);
Iris_83(1,:)=s;
clear Imodes
load D:\APPLICATION1\Matlab7\work\package_emd\code\Data1_3\3_2.mat
for i=21:30
    if i==21
        s=Imodes(21,:,N1);
    else
        s=[s Imodes(i,:,N1)];
    end
end
% s=Imodes(25,:,N1);
Iris_83(2,:)=s;
clear Imodes
load D:\APPLICATION1\Matlab7\work\package_emd\code\Data1_3\3_3.mat
for i=21:30
    if i==21
        s=Imodes(21,:,N1);
    else
        s=[s Imodes(i,:,N1)];
    end
end
% s=Imodes(25,:,N1);
Iris_83(3,:)=s;
clear Imodes
load D:\APPLICATION1\Matlab7\work\package_emd\code\Data1_3\3_4.mat
for i=21:30
    if i==21
        s=Imodes(21,:,N1);
    else
        s=[s Imodes(i,:,N1)];
    end
end
% s=Imodes(25,:,N1);
Iris_83(4,:)=s;
clear Imodes
%*****************************************************
load D:\APPLICATION1\Matlab7\work\package_emd\code\Data1_3\4_1.mat
for i=21:30
    if i==21
        s=Imodes(21,:,N1);
    else
        s=[s Imodes(i,:,N1)];
    end
end
% s=Imodes(25,:,N1);
Iris_88(1,:)=s;
clear Imodes
load D:\APPLICATION1\Matlab7\work\package_emd\code\Data1_3\4_2.mat
for i=21:30
    if i==21
        s=Imodes(21,:,N1);
    else
        s=[s Imodes(i,:,N1)];
    end
end
% s=Imodes(25,:,N1);
Iris_88(2,:)=s;
clear Imodes
load D:\APPLICATION1\Matlab7\work\package_emd\code\Data1_3\t1.mat
for i=21:30
    if i==21
        s=Imodes(21,:,N1);
    else
        s=[s Imodes(i,:,N1)];
    end
end
% s=Imodes(25,:,N1);
Iris_88(3,:)=s;
clear Imodes
%*******************************************

% train=[Iris_44(1:3,:);Iris_51(1:3,:);Iris_83(1:3,:);Iris_88(1:3,:)];
% test=[Iris_44(4,:);Iris_51(4,:);Iris_83(4,:);Iris_88(4,:)];
% target=[ones(3,1);2*ones(3,1);3*ones(3,1);4*ones(3,1)];
% test_target=[1;2;3;4];

train=[Iris_44(1:2,:);Iris_51(1:2,:)];
test=[Iris_44(3,:);Iris_51(3,:)];
target=[ones(2,1);2*ones(2,1)];
test_target=[1;2];
% train=[Iris_44(1:2,:);Iris_51(1:2,:);Iris_83(1:2,:);Iris_88(1:2,:)];
% test=[Iris_44(3,:);Iris_51(3,:);Iris_83(3,:);Iris_88(3,:)];
% target=[ones(2,1);2*ones(2,1);3*ones(2,1);4*ones(2,1)];
% test_target=[1;2;3;4];
 %train10=(train-min(min(train)))/(max(max(train))-min(min(train)));%规一化（全归一，按列归一）
 %test10=(test-min(min(train)))/(max(max(train))-min(min(train)));

 for ii=1:1000
  scale=ii/10000;
  model = svmtrain(target, train/scale,[,'-c 1000 -g 4  '])
 [predicted_label, accuracy] = svmpredict(test_target,test/scale, model);
 result(ii)=accuracy(1);
 end