EM算法在混合高斯分布中的应用MATLAB代码

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）））;