作者quietfly (黑白照片)
看板MATLAB
标题[问题] fminsearch函数问题
时间Thu Nov 3 23:39:00 2011
我想再请问板上的高手
以下这个模型
因为目前已知的数据就是r 跟 d
其中12个待解参数 跟矩阵u,h都是在优化过程中解出来的
如果我想再加个条件 h(:)>0 应该怎麽写??
function [llf,LogLikelihood]=LLF(para,r,d)
T=length(r);
h=zeros(T,1);
c=para(1);
h(1)=mean((r-ones(T,1)*mean(r)).^2);
u(1:3)=r(1:3)-ones(3,1)*c;
for t=4:T
u(t)=r(t)-c-r(t-1:-1:t-3)'*para(2:4)-para(R+2)*r(t-3)*d(t);
h(t)=para(5)+para(6)*u(t-1)^2+para(7)*h(t-1)+para(8)*d(t)+
para(9)*d(t)*u(t-1)^2+para(10)*d(t)*h(t-1)+
para(11)*d(t)*u(t-1)^2*h(t-1);:
end
llf = -sum(log(h(4:T))+((u(4:T)).^2)'./h(4:T));
end
objectionfun = @(para) -LLF(para,r,d);
para0=[-0.001688;0.010819;-0.129009;-0.057876;0;0.15;0.75];
options = optimset('MaxFunEvals',10000000,'MaxIter',100000,'TolX',1e-20,'Hessian','bfgs','Algorithm','interior-point');
[parae,llfs,exitflag,output] = fminsearch(objectionfun,para0,options);
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 222.29.34.110
※ 编辑: quietfly 来自: 222.29.34.110 (11/03 23:43)