在matlab中如何用MLE估计混合正态分布的参数？

手头有一组数据，大约190个，具备混合正态分布的特征（φ( x) = ρφ1 ( x) + (1 - ρ)φ2 ( x)），如何在matlab7.0中用MLE估计出参数（μ1 ,σ1 ,μ2 ,σ2 ,ρ）？谢谢！
原始数据如下：
-0.8428-0.7903-0.786-0.7825-0.7825-0.7778-0.7578-0.7547-0.7422-0.7264-0.7263-0.7245-0.7222-0.7189-0.7043-0.7039-0.6994-0.6971-0.6864-0.6828-0.6744-0.6729-0.669-0.6648-0.6624-0.6434-0.6431-0.6417-0.6417-0.6398-0.63608-0.6255-0.6232-0.6229-0.6224-0.6219-0.6203-0.6174-0.6028-0.6014-0.5956-0.5942-0.5939-0.58018-0.5612-0.5562-0.551-0.5268-0.4995-0.4917-0.4804-0.4509-0.4151-0.4036-0.3771-0.3657-0.2723-0.2492-0.2403-0.1981-0.1622-0.1271-0.0972-0.0928-0.086-0.0826-0.0741-0.0685-0.0671-0.0616-0.0517-0.0503-0.0493-0.0476-0.0456-0.0417-0.0415-0.0267-0.0222-0.0199-0.0191-0.0164-0.0151-0.0118-0.0073-0.0035-0.0031-0.0026-0.0023-0.0012-0.00040.0080.0083480.00920.00990.0120.01230.01240.01450.01460.01480.01870.02070.02170.02280.0230.02390.02440.02460.02530.02680.02730.02930.030.0310.03190.03370.03490.03550.03580.0370.03880.03880.04190.0420.04220.0440.04460.0450.0450.04920.05020.05060.05150.05160.05290.0550.05640.05840.05910.05940.05950.05970.06370.06460.06560.06850.06860.06860.07370.07560.07590.07780.07980.08210.0850.08590.08790.08850.09060.09190.09380.09840.0990.09940.10080.10220.1030.10620.10650.11260.11340.11750.12320.13460.13620.13730.17290.17360.17790.17820.21390.21680.23330.31260.31480.31780.3178

1# yuyi3860

详见：

http://blog.macro2.org/2011-07/用mle函数估计随机数分布的参数.html

写似然函数然后最优化 只能这么处理哈

Xaero 的答案准确、全面、及时，雪中送炭，谢谢！

很有用

Xaero  给力

liuxin9023 发表于 2011-7-20 07:28
写似然函数然后最优化 只能这么处理哈

为什么我的MATLAB里总是有错啊？

??? Error using ==> mlecustom>llf_pdfcdf at 443
The PDF function returned negative or zero values.

Error in ==> fminsearch at 310
    x(:) = xr; fxr = funfcn(x,varargin{:});

Error in ==> mlecustom at 186
        [phat,nll,err,output] = ...

Error in ==> mle at 227
            [phat, pci] = mlecustom(data,varargin{:});

+.軽描淡写 发表于 2013-3-27 10:15
为什么我的MATLAB里总是有错啊？

??? Error using ==> mlecustom>llf_pdfcdf at 443

您好,我也现在也遇到了和你一样的错误,请问你的错误解决了吗?能不能帮我解决一下啊

Xaero 发表于 2011-7-19 13:25
1# yuyi3860

详见：

谢谢你，真的很受用！！

dongtao1010 发表于 2013-5-6 14:21
您好,我也现在也遇到了和你一样的错误,请问你的错误解决了吗?能不能帮我解决一下啊

写似然函数然后最优化,有没有写好的程序借鉴一下啊？刚刚开始接触matlab，现在正要用到这个，但是一窍不通那~谢谢呐