%目標(biāo)函數(shù)
function F = tracklsq(input)
a=input(1);
b=input(2);
opt = simset('solver','ode5','SrcWorkspace','Current');
[tout,xout,yout] = sim('trackoptim',[0 5],opt);
%計(jì)算誤差信號(hào)
F = yout-1;       
%調(diào)用優(yōu)化函數(shù)
%調(diào)入仿真模型
trackoptim   
%初始化變量                   
input=[8.163 1];  
%設(shè)置優(yōu)化參數(shù)
options = optimset('LargeScale','off','Display','iter', 'TolX',0.001,'TolFun',0.001);
%調(diào)用優(yōu)化函數(shù)
[output,fval]=lsqnonlin('tracklsq',input,[],[],options)
% 輸出參數(shù)
a=output(1)
b=output(2)