1# clark1025
用fminsearch做
先在script file里面猜测一组数据(对要估计的参数猜测)eg:guess=[1 2 3](代表你要估计的三个参数,1 2 3是我自己随便写的)
然后用 estimate=fminsearch('myfct',guess);
调用函数
function[f]=myfct(x)
global y T; %把script里的sample设成global variable
T=length(y);
yita=x(1);
alpha1=x(2);
beta1=x(3);
然后设定你的log likelihood function
lnl=lnl=-(T-1)/2*log(2*pi)-slnh/2-sz/2;(随便举的例子)
然后
f=-lnl;(这个必须要,因为fminsearch是找最小值,加了负号才变成找最大值)
end;