function?[?pai?Yz?mu?E?]?=?GMM_function（?XKerrmax_time?）
%GMM函數?其中X表示n*d的樣本?K為高斯個數
%???err為閾值?max_time為最大迭代次數
tic
if?nargin?==?2
????err=1e-14;
????max_time=1000;
end
if?nargin?==?3
????max_time=1000;
end
%%%%%%主程序
X?=?X‘;??%X為D*N型數據，跟PRML對樣本數據描述相反
[dim?num_data]?=?size（X）;

%Initialize
%-------------------------------
%K?=?numel（unique（data.y））;
[inx?C]?=?kmeans（X‘K）;
mu?=?C‘;
for?i=1:dim%%%%%%%%%%%%%%%%%%%%求均值居然出現負值…………
????for?j=i:K
????????if（mu（ij）<0）
????????????mu（ij）=0;
????????end
????end
end
pai?=?zeros（1K）;
E?=?zeros（dimdimK）;
for?k=1:K
????pai（k）?=?sum（inx==k）;???
????%%%inition?E
????E（::k）=iniE（pai（k）inx?CX‘k）;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%畫圖測試
%?figure（1）;
%?for?i=1:K
%?[x?y]=meshgrid（linspace（0108

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????6231??2014-11-25?20:41??test.fig
?????文件???????10889??2014-11-25?20:41??test.m
?????文件????????2619??2014-11-14?13:02??GMM_function.m