我想用3个泊松分布的总和拟合一个实验分布:
数据=Cal_H2;
Lambda_Zero=pd2.Lambda;
Lambda_One=ResultOne(2);
TwoDist=@(数据,Ampl_1,Ampl_2,Lambda_2)(1-Ampl_1-Ampl_2)。*poisspf(数据,Lambda_零)+Ampl_1.*poisspf(数据,Lambda_1)+Ampl_1.*poisspf(数据,Lambda_2);
StartParam=[0.02 0.9 20];
LowerBound=[0];
上限=[1 150];
oldopts=statset(“定制”);
newopts=statset(旧选项,“GradObj”,“开”);
[结果二,限定二]=mle(数据,“pdf”,两区,“开始”,星图,“Lowerbound”,LowerBound,“上限”,上界,“最佳乐趣”,“fmincon”,“选项”,newopts)
但由于结果必须归一化,我需要施加A+B+C=0,A,B,C在0和1之间,我不知道如何将这个方程发送到
铁铬镍铁合金
通过我们投入的选择
最大似然误差
. 此外,当我直接运行它时,我有一个问题:
未定义的函数或变量“Hstr”。
fmincon错误(第482行)
sfminbx(funfcn,X,l,u,详细,选项,defaultopt,computeLambda,f,GRAD,HESS,Hstr,varargin{:});
自定义中的错误(第221行)
fmincon(@llf_diff,start,[]、[]、[]、[]、[]、lb、ub、[]、opts、unsensdata、censtata、unsensrefq、censFreq、llf、fun1Args、fun2Args、checkFunVals、delta);
mle中的错误(第247行)
[phat,pci]=mlecustom(数据,varargin{:});
如果有人能帮我,我的MATLAB不是很流利,仅仅发送一个指向帮助/文档的链接可能不会有用,因为我很可能无法理解细节。非常感谢!