标题:[求助]fsolve可不可以解我这个二元非线性方程组
只看楼主
linger1983
Rank: 1
等 级:新手上路
帖 子:2
专家分:0
注 册:2006-5-12
 问题点数:0 回复次数:1 
[求助]fsolve可不可以解我这个二元非线性方程组
function f=myfunction2(alpha,zta,y,n,r)
z=-y.^alpha;
z=1-exp(z);
h=r*log(alpha*zta)+(alpha-1).*sum(log(y(1:r)))-sum(-z(1:r))+(zta-1).*sum(log(z(1:r)))+(n-r).*log(1-z(r).^zta);
f(1)=diff(h,alpha);
f(2)=diff(h,zta);



n=30;r=20;
alpha0=0.3;zta0=2;
for j=1:n
y(j)=-log(1-rand(1).^(1/zta0));
y(j)=y(j).^(1/alpha0);
end
y=sort(y);
x0=[0.3,2];
syms alpha zta;
[x,fval]=fsolve(@myfunction2,x0,alpha,zta,y,n,r)

是一个参数的最大似然估计,不知道怎么解出alpha,zta.请高手指点,谢谢!!!

搜索更多相关主题的帖子: 方程 fsolve alpha 非线性 log 
2006-05-14 15:49
linger1983
Rank: 1
等 级:新手上路
帖 子:2
专家分:0
注 册:2006-5-12
得分:0 
没人告诉我啊
2006-05-15 12:55



参与讨论请移步原网站贴子:https://bbs.bccn.net/thread-64125-1-1.html




关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.434292 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved