k_l_bayesianclassifier.m

Karhunen-Loève特征变换 模式识别课程实验

% M-file function, K_L_BayesianClassifier.m 
% 用于K-L变换的两类正态分布模式的贝叶斯分类器
% P1 第1组先验概率
% P2 第2组先验概率
% c11   第1组判对总数
% c12   第1组判错总数
% c12_p 第1组判错率
% c21   第2组判对总数
% c22   第2组判错总数
% c22_p 第2组判错率

function [c11,c12,c12_p,c21,c22,c22_p] = K_L_BayesianClassifier

clear;
     load Y2_123.dat;
     load Y3_123.dat;
Y1(:,:) = Y2_123(:,:);
Y2(:,:) = Y3_123(:,:);
n = 3;


g1 = 1;
g2 = 2;
P1 = 0.5;
P2 = 0.5;
[m1,s1] = K_L_Parameter(Y1);
[m2,s2] = K_L_Parameter(Y2);
c11 = 0;
c12 = 0;
c21 = 0;
c22 = 0;
for k = 1:50
    c1 = K_L_Classifier(g1,g2,Y1,k,m1,s1,P1,m2,s2,P2,n);
    c2 = K_L_Classifier(g1,g2,Y2,k,m1,s1,P1,m2,s2,P2,n);
    if (c1 == g1)
        c11 = c11 + 1;
    else
        c12 = c12 + 1;
    end
    if (c2 == g2)
        c21 = c21 + 1;
    else
        c22 = c22 + 1;
    end    
end
c12_p = c12 / (c11 + c12);
c22_p = c22 / (c21 + c22);