Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
|
![]() |
#1 |
初级会员
注册日期: 2012-04-09
年龄: 35
帖子: 1
声望力: 0 ![]() |
![]()
1,关于GT函数(广义t分布)参数的极大似然估计。
函数形式:function y= GT( x,u,w,p,q) y=p./((2.*w.*q.^(1./p).*beta(1./p,q)).*(1+abs(x-u).^p./(q.*w.^p).^(-q-1./p))); end 问题的描述:其中x是一个很大的数据向量(1761*1),u,w,p,q是参数,我们要根据已知的x来估计它们。 尝试过的解决办法: Mle:(失败:在迭代过程中密度函数出现负值。可能是matlab中自带的mle与本问题不相容。Ps:对于其他密度函数改方法可以很好解决。) 乘积后取对数再对各个变量求偏导:似乎计算量太大了?我把变量降到两维都不行。 Fmincon:不太会用。。 2,求解laplace函数0.05分位数。 function y = Laplace( x,u,s ) y=1/(sqrt(2)*s)*exp(-sqrt(2)*abs(x-u)/s); end syms t double(solve(int(Laplace(x,0.1109,1.8574),-inf,t)==0.05)) 问题:matlab报错:Warning: Explicit solution could not be found. 3,t-location 密度函数调用在matlab中是什么命令? syms t solve(int(tpdf((x-0.118931)/1.28819,3.54706),-inf,t)==0.05) |
![]() |
![]() |