clear
%女生的身高體重的均值和協方差
fid=fopen（‘FEMALE.txt‘‘r‘）;
Cf?=?textscan（fid?‘%f?%f‘）;
mheightf=mean（Cf{1}）;
mweightf=mean（Cf{2}）;
sigme2f=cov（Cf{1}Cf{2}）;
%男生的身高體重的均值和協方差
fid2=fopen（‘MALE.txt‘‘r‘）;
Cm?=?textscan（fid2?‘%f?%f‘）;
mheightm=mean（Cm{1}）;
mweightm=mean（Cm{2}）;
sigme2m=cov（Cm{1}Cm{2}）;

%身高和體重相關時，判別函數
m_female=[mheightfmweightf]‘;
m_male=[mheightmmweightm]‘;
disp（‘請輸入身高和體重矩陣：‘）
R=input（‘‘）;
R=R‘?
%?%?%?%?syms?h?w
%?%?%?%?R=[h?w]‘;
g1=-1/2*（（R?-?m_female）‘*inv（sigme2f）*（R?-?m_female）-（R?-?m_male）‘*inv（sigme2m）*（R?-?m_male））;
g2=-1/2*log（det（sigme2f）/det（sigme2m））;
g3=log（0.9/0.1）;
gi=g1+g2+g3
%決策結果
disp（‘身高和體重相關時，經過Bayes最小錯誤率的判別結果為：‘）
if?gi>0
????disp（‘It?is?a?female!‘）
else?
????disp（‘It?is?a?male!‘）
end

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

?????文件???????1043??2010-04-26?13:19??test_error.m

?????文件???????1015??2010-03-29?22:22??xunlian_error.m

?????文件????????833??2010-04-06?23:23??bayes_error.m

?????文件????????864??2010-03-28?22:13??bayes_error_notc.m

?????文件????????426??2004-10-08?09:03??FEMALE.TXT

?????文件????????424??2004-10-08?09:04??MALE.TXT

?????文件????????348??2010-03-29?10:48??test1.txt

?????文件???????3018??2010-03-29?11:58??test2.txt

-----------?---------??----------?-----??----

?????????????????7971????????????????????8