T=1000;
J=2;
mu=(ones(J,1)*0.09)'; %1x15 (1xd)
sigma=eye(J)*0.2;
X=mvnrnd(mu,sigma,T); %156x15 matrix (nxd)
options=optimset('Display','iter','TolFun',10^(-25),'TolX',10^-8,'MaxFunEvals',1000,'MaxIter',1000);
vp0=[rand(4,1)*0.15];
vpML=fminsearch(@(vp) tester(vp,X),vp0,options)

请大神看看
先谢谢啦！！