m=20000;???????????????%樣本數量
n=20;????????????????%迭代次數
x=?[normrnd（20.5[1m/2]）normrnd（40.5[1m/2]）];????%生成1*m的矩陣，分別為均值為3和5，標準差為1的正態分布
a=[0.30.7];u=[16];e=[0.850.9];????????????????%初始化，a為輸出的概率，u為均值，e為方差
w=zeros（m2）;??????????????%w為每個元素屬于哪個高斯分布的概率
p=zeros（12）;???????????????%p為兩個高斯分布的概率密度

for?k=1:n
????%E?step
????for?i=1:m
????????for?j=1:2
????????????p（1）=exp（-（（x（i）-u（1））^2/（2*e（1））））/（sqrt（2*pi）*sqrt（e（1）））;?????????????????