# パターン認識課題(ソースコード)
### オプション課題2(うまくいかなかった)

clear all

load digit.mat X T
[d,n,nc]=size(X)
%dはマス番号の総数
%nは訓練用の文字データがいくつあるか
%ncは1~10
nc = 10;
C=zeros(nc,nc);
covi=cov(X(:,:,1)')
for s=1:200
    for ct=1:10
        flag = false;num=ct;
        for i=1:10
            if i==10
                num=i;
                break
            end
            mui=mean(X(:,:,i),2);covi=cov(X(:,:,i)');
            for j=i+1:10
                muj=mean(X(:,:,j),2);
                covj=cov(X(:,:,j)');
                S=(covi+covj)/2; invS=inv(S);
                t=T(:,s,ct);
                p1=t'*invS*mui-mui'*invS*mui/2;
                p2=t'*invS*muj-muj'*invS*muj/2;
                if p2>p1
                    break
                end
                if j==10
                    flag=true;
                end
            end
            if flag
                num=i;
                break
            end
        end
        C(ct,num)=C(ct,num)+1;
    end
end