%設置全局變量 meas species bad numobs
global meas species bad numobs quadclass

%散點圖
gscatter(meas(:,1), meas(:,2), species,'mgb','vsp');
xlabel('Sepal 長度');
ylabel('Sepal 寬度');

%二次分類
quadclass= classify(meas(:,1:2),meas(:,1:2),species);

%計算錯誤分類      strcmp函數比較字符串
bad = ~strcmp(quadclass,species);
%計算meas第一列大小
numobs = size(meas,1);
%sum函數求元素的和
sum(bad) / numobs;

%在散點圖上顯示錯誤分類 用x表示
hold on;
plot(meas(bad,1), meas(bad,2), 'kx');
hold off;