%此文件是Bayes最小誤差法的實現函數，主文件直接對其進行調用，無參數，返回值是識別結果。

function?[number]=Bayes_minerror（）
?????X=[];
?????A=[];
????number=0;
?????currPath?=?fileparts（mfilename（‘fullpath‘））;
?????manu_picture=strcat（currPath‘\‘‘manu.bmp‘）;
?????manu_picture=pre_figure（manu_picture）;
?????for?m=1:10
?????????Y=[];
?????????folder_name=sprintf（‘%d‘m-1）;
?????????my_str1=strcat（currPath‘\手寫數字\‘folder_name）;??%這幾句主要進行模板數字的路徑處理
?????????for?n=1:10
?????????????Z=[];
?????????filefolder_name=sprintf（‘%d‘n）;?????
?????????????file_name=strcat（folder_name‘-‘filefolder_name）;
?????????????model_picture=strcat（my_str1‘\‘file_name‘.bmp‘）;?????????????
?????????????A=pre_figure（model_picture）;???%調用預處理函數，獲得特征
?????????????for?z=1:5??????????????????%將5*5的特征矩陣轉換成1*25的矩陣
?????????????????Z=[ZA（z:）];???%矩陣重組
?????????????end
?????????????Y（n:）=Z;???????????%矩陣重組
?????????end
?????????X（::m）=Y;??????????%X是3維矩陣，第一維分別代表每個樣本文件夾里的10個數字。
?????end??????????????????????%第二維是每個樣本的25個特征，第三維分別代表10個文件夾。
?????
?????mymanu_picture=[];
?????for?z=1:5
?????????mymanu_picture=[mymanu_picturemanu_picture（z:）];?%把待識別的手寫圖片5*5轉換成1*25
?????end

?%%%%%%%%???下面是Bayes最小誤差法，基本原理見圖像模式識別—VC++技術實現.pdf???%%%%%%%%%%
?????
???hx=zeros（110）;
???for?n=1:10
???????Xmeans=zeros（125）;
???????%求樣本平均值
?????for?k=1:10
?????????for?i=1:25?????????????
?????????????Xmeans（i）=Xmeans（i）+X（kin）;
?????????end
?????end
?????for?i=1:25
?????Xmeans（i）=Xmeans（i）/10;
?????end
?????
?????%求協方差矩陣
?????mode=zeros（1025）;
?????for?i=1:10
?????????for?j=1:25
?????????????mode（ij）=X（ijn）;
?????????end?????????
?????end
?????for?i=1:25
?????????for?j=1:25
?????????????s=0.0;
?????????????for?k=1:10
?????????????????s=s+（mode（ki）-Xmeans（i））*（mode（kj）-Xmeans（j））;
?????????????end
?????????????s=s/9+0.0000000001*rand（）;??%這里加了個隨機數，因為會有很多項s為0，導致S中出現全0列，使下面對S求逆進行不了
?????????????%s=s/9;
?????????????S（ij）=s;
?????????end????
?????end
?????
?????%求先驗概率
?????Pw=0.1;???????%總樣本數為100，每類樣本數都是10，故此為0.1
?????invS=inv（S）;??%求S的逆矩陣
?????detS=det（S）;??%求S的行列式
?????
?????%?求判別函數
?????Xx=mymanu_picture;
?????for?i=1:25
?????????Xx（i）=Xx（i）-Xmeans（i）;?????
?????end
?????t1=Xx*invS*Xx‘;
?????t2=log（Pw）;
?????t3=log（detS+1）;
?????hx（n）=-t1/2+t2-t3/2;
???end
?????????????
????%求判別函數的最大值
????maxval=hx（1）;
????for?i=1:10
????????if（hx（i）>maxval）
????????????maxval=hx（i）;
????????????number=i-1;
????????end
????end?????????????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2750??2012-04-20?15:19??手寫數字識別\Bayes_MinError.m

?????文件???????4118??2012-04-20?15:19??手寫數字識別\handwritten_numbers.m

?????文件????????540??2012-04-20?15:19??手寫數字識別\myoutput.m

?????文件???????2067??2012-04-20?15:38??手寫數字識別\pre_figure.m

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\0\0-1.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\0\0-10.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\0\0-2.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\0\0-3.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\0\0-4.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\0\0-5.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\0\0-6.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\0\0-7.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\0\0-8.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\0\0-9.bmp

????..A..H.?????36352??2012-04-19?14:54??手寫數字識別\手寫數字\0\Thumbs.db

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\1\1-1.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\1\1-10.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\1\1-2.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\1\1-3.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\1\1-4.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\1\1-5.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\1\1-6.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\1\1-7.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\1\1-8.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\1\1-9.bmp

?????文件???????1810??2012-04-19?14:54??手寫數字識別\手寫數字\1\quanzhi.txt

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\2\2-1.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\2\2-10.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\2\2-2.bmp

?????文件??????58678??2012-04-19?14:54??手寫數字識別\手寫數字\2\2-3.bmp

............此處省略218個文件信息