xueyuanxo
2012-04-09, 22:07
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)
函数形式：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)