clear
syms?pw1?pw2;
pf=input（‘女生的先驗概率：‘）;
pm=1-pf;
X=input（‘請輸入身高體重矩陣：‘）;
X=X‘

%求出為女生的概率密度函數
fid1=fopen（‘FEMALE.txt‘‘r‘）;
rf=textscan（fid1?‘%f?%f‘）;
fheight=mean（rf{1}）;
fweight=mean（rf{2}）;
fm=[fheight;fweight];
fM1=cov（rf{1}rf{2}）;
?fM=det（fM1）;
fw1=1/（（2*pi）*sqrt（fM））*exp（-0.5*（X-fm）‘*inv（fM1）*（X-fm））;
%求出為男生的概率密度函數
fid2=fopen（‘MALE.txt‘‘r‘）;
rm=textscan（fid2?‘%f?%f‘）;
mheight=mean（rm{1}）;
fweight=mean（rm{2}）;
X2=[fheight;fweight];
mN1=cov（rm{1}rm{2}）;
mN=det（（mN1））;
fw2=1/（（2*pi）*sqrt（mN））*exp（-0.5*（X-X2）‘*inv（mN1）*（X-X2））;
%求后驗概率
hpf=pf*fw1/（pf*fw1+pm*fw2）;
hpm=pm*fw2/（pf*fw1+pm*fw2）;
?%設置風險矩陣
f=[0?7;3?0];
%計算風險及進行判決
Rf=f（11）*hpf+f（12）*hpm
Rm=f（21）*hpf+f（22）*hpm
disp（‘身高和體重相關時，最小風險Bayes判決結果：‘）
if?Rf>=Rm
?????disp（‘It?is?a?male!‘）;
else
?????disp（‘It?is?a?female!‘）;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2019-04-14?15:40??模式識別-貝葉斯matlab實例\
?????目錄???????????0??2019-04-14?15:39??模式識別-貝葉斯matlab實例\最小風險貝葉斯決策\
?????文件?????????915??2010-03-28?22:26??模式識別-貝葉斯matlab實例\最小風險貝葉斯決策\bayes_sike.m
?????文件????????1005??2010-03-28?22:25??模式識別-貝葉斯matlab實例\最小風險貝葉斯決策\bayes_sike_notc.m
?????文件?????????426??2004-10-08?09:03??模式識別-貝葉斯matlab實例\最小風險貝葉斯決策\FEMALE.TXT
?????文件?????????424??2004-10-08?09:04??模式識別-貝葉斯matlab實例\最小風險貝葉斯決策\MALE.TXT
?????文件?????????348??2010-03-29?10:48??模式識別-貝葉斯matlab實例\最小風險貝葉斯決策\test1.txt
?????文件????????3018??2010-03-29?11:58??模式識別-貝葉斯matlab實例\最小風險貝葉斯決策\test2.txt
?????文件????????1197??2010-03-29?12:40??模式識別-貝葉斯matlab實例\最小風險貝葉斯決策\test_sike.m
?????文件????????1171??2010-03-29?22:20??模式識別-貝葉斯matlab實例\最小風險貝葉斯決策\xunlian_sike.m
?????目錄???????????0??2017-11-12?10:17??模式識別-貝葉斯matlab實例\貝葉斯程序示例\
?????目錄???????????0??2017-11-12?10:17??模式識別-貝葉斯matlab實例\貝葉斯程序示例\最小錯誤率\
?????文件????????1596??2016-12-02?15:52??模式識別-貝葉斯matlab實例\貝葉斯程序示例\最小錯誤率\bayes_cell.m
?????文件???????14778??2016-12-02?15:37??模式識別-貝葉斯matlab實例\貝葉斯程序示例\最小錯誤率\參考程序.docx
?????目錄???????????0??2017-11-12?10:17??模式識別-貝葉斯matlab實例\貝葉斯程序示例\最小風險\
?????文件????????1736??2016-12-09?13:57??模式識別-貝葉斯matlab實例\貝葉斯程序示例\最小風險\bayes.asv
?????文件????????1736??2016-12-09?14:03??模式識別-貝葉斯matlab實例\貝葉斯程序示例\最小風險\bayes.m
?????文件?????????304??2016-12-02?15:53??模式識別-貝葉斯matlab實例\貝葉斯程序示例\最小風險\bayes_cell2.m
?????文件???????15456??2016-12-09?14:05??模式識別-貝葉斯matlab實例\貝葉斯程序示例\最小風險\參考程序-最小風險.docx
?????文件???????86016??2016-12-09?14:10??模式識別-貝葉斯matlab實例\貝葉斯程序示例\貝葉斯程序示例.doc