我不知道你想做什么。
尝试这个:
mu0 = @(a1,a2,t,tg,e)(E. * exp(a1。*(t-tg)./(a2 +(t-tg))))。*(t> tg)+ e。*(t <= tg);
mu = @(a1,a2,t,tg,e,gammaDot,tau,n)mu0(a1,a2,t,tg,e)./(1+(gammadot。* mu0(a1,a2,t,tg,e)/ tau)。^(1-n));
tau = 2.9e4;
n = 0.33;
E = 3.63e11;
A1 = 27.21;
A2 = 92.85;
tg = 373;
%%% t = b(1),γ= b(2)
[B,FVAL] = Fminunc(@(b)mu(a1,a2,b(1),tg,e,b(2),tau,n),randi(500,2,1));
FPRINTF(1,'t =%9.3f \ ngammadot =%9.3f \ n',b)
进行适当的更改以获得您想要的结果。