Tmin=500; Tmax=2500; Xs=0.15; Tm=1500;
for i=1:30
options.TolX = 1;
[Topt,yy]=fminbnd(@fun_simuhitran,Tmin,Tmax,options);
if((abs(Xs-Xopt)<5.0E-3)&&(abs(Topt-Tm)<2))
count=count+1; %A = count(str,pattern) 返回 pattern 在 str 中的出现次数
else
count=0; %非趋于收敛状态,重置count
end
Xs=Xopt;
Tm=Topt;
disp([i,count,Tm,Xs]);
if(count>3) % 连续3次以上迭代变化较小,退出迭代
break;
end
end