作者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)