sos.m

文件是编辑好的m文件

figure('pos',[100,100,300,200],'color','w');
axes('pos',[.1,.1,.8,.8]);
%t=[0,10,20,40,60,80,100,150,200,300,400,500];
%h=[0,-0.5,9.5,18,33,45,55,63,78,86,95,98];
load d:\yzpta\stepdata;t=tt;h=yy;tm=max(t);
%load d:\yzpta\invdata;[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;
tt=0:0.5:tm;hh=spline(t,h,tt);
x0=[2,15,3];
FUN1= inline('x0(1)*(1-exp(-((tt-x0(3)).*(tt>x0(3)))/x0(2)))','x0','tt');
x=lsqcurvefit(FUN1,x0,tt,hh)
y=FUN1(x,tt);err=sum((hh-y).^2)
save d:\yzpta\fitdata1 t h tt y x;
plot(t,h,'r*',tt,y,'k');grid;%pause;
axis([0 50 0 1.5])
%plot(tt,hh-y);grid;maxerr=max(abs(hh-y));