极大似然估计的MATLAB程序，哪位大神帮忙看看，多谢了

你就把alpha，beta和lambda看成x1，x2，x3，x看成已知的数。然后fminunc就行了。（我觉得需要设置alpha和beta值的范围。）

magicsun 发表于 2015-6-25 10:26
你就把alpha，beta和lambda看成x1，x2，x3，x看成已知的数。然后fminunc就行了。（我觉得需要设置alpha和be ...

你的意思是不需要我做实验弄的那些观测值？

需要，这些就是x，已知的数

magicsun 发表于 2015-6-25 11:59
需要，这些就是x，已知的数

>> clear
>> fx=@(x)-log(x(1)/x(2)*((1.1101e+03-x(3))/x(2))^(x(1)-1)*exp(-(1.1101e+03-x(3))/x(2))^x(1))-log(x(1)/x(2)*((1.0035e+03-x(3))/x(2))^(x(1)-1)*exp(-(1.0035e+03-x(3))/x(2))^x(1))-log(x(1)/x(2)*((941.2414-x(3))/x(2))^(x(1)-1)*exp(-(941.2414-x(3))/x(2))^x(1))-log(x(1)/x(2)*((971.2918-x(3))/x(2))^(x(1)-1)*exp(-(971.2918-x(3))/x(2))^x(1))-log(x(1)/x(2)*((1.1419e+03-x(3))/x(2))^(x(1)-1)*exp(-(1.1419e+03-x(3))/x(2))^x(1));
[x,fval,exitflag,output]=fminsearch(fx,[1 200 800])
是这样吗？

求对数之后alpha和beta就不用限定了，他们是ln（alpha），自动保证大于0了。

应该就是这意思吧！但我不知道你的公式写的对不对！

magicsun 发表于 2015-6-25 19:23
应该就是这意思吧！但我不知道你的公式写的对不对！

x =

   1.0e+03 *

    0.0007    0.0651    1.4948


fval =

  -0.2028 + 4.4201i


exitflag =

     1


output =

    iterations: 41
     funcCount: 147
     algorithm: 'Nelder-Mead simplex direct search'
       message: [1x194 char]
这是我运行的结果，但是极小值是虚数？

还是错的，那可能还是需要设定alpha 和beta的范围。

magicsun 发表于 2015-6-25 20:46
还是错的，刚才犯错了。

？什么意思？

那可能还是需要设定alpha 和beta的范围。