function S=myfun8(nst,Snul,r,q,T,n,vg)
dt=T/n;
tt=[0:dt:T];
omega=log(1-nst(2)^2*nst(1)/2-nst(3)*nst(1))/nst(1);
for s = 1:n
g= gamrnd(dt/nst(1),nst(1));
vg(s+1) = vg(s) + nst(3)*g+ nst(2)*normrnd(0,sqrt(nst(3)));
%S(s+1) = Snul*exp(((r-q)+log(1-nst(2)^2*nst(1)/2-nst(3)*nst(1))/nst(1)).*tt(s+1)+vg(s+1));
end;

S1 = Snul*exp((r-q+omega).*tt+vg);

S=mean(S1);
这个是我写的VG model代码，但是怎么用数据，例如（ S&P index data),  去得到优化的三个系数呢？

[此贴子已经被作者于2008-4-5 18:01:23编辑过]

 谢谢楼上的，因为初学matlab，但如何实现呢？但是怎么用数据，例如（ S&P index data),  去得到优化的三个系数呢？
用fminsearch命令怎么有效得到，并且和数据结合呢？