ex2_3.m

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

% ex2_3
% 具有反相特性过程的参数估计
figure('pos',[100,100,300,200],'color','w');
axes('pos',[.15,.1,.8,.8]);
load invdata2;
[r,tm]=size(y);
t=y(1,[1:10:200,200:100:tm]);
h=y(2,[1:10:200,200:100:tm]); 
tt=t;
hh=h;
x0=[1,12,2,24,3];
opti=optimset('fminbnd');
val=optimset(opti,'MaxFunEvals',3000);
x=lsqcurvefit('fun2',x0,tt,hh,0,100,val)
y=fun2(x,tt);
err=sum((hh-y).^2);
save c:\fitdata2 t h tt y x;
plot(t,h,'b*',tt,y,'k');
grid;
axis([0 100 -0.1 0.05]);
set(gca,'ytick',[-.1:0.05:.05],'ytickmode','manual');