MATLAB自带Weibull参数估计的函数:weibfit()
还有,你的似然函数写得不对