ex2_4.m

来自「华东理工大学自动化系《控制系统分析、设计和应用》教材配套程序实例,有关说明： 」· M 代码 · 共 17 行

% ex2_4
% 二阶自衡振荡过程的参数估计
load oscdata;
st=osc;
t=st(1,:);
h=st(2,:);
tm=max(t);
tt=0:0.1:tm;
hh=spline(t,h,tt);
x0=[2,15,3];
opti=optimset('fminbnd');
val=optimset(opti,'MaxFunEvals',3000);
x=lsqcurvefit('fun3',x0,tt,hh,0,100,val)
y=fun3(x,tt);
err=sum((hh-y).^2)
plot(t,h,'b*',tt,y,'k');
grid;