sos1.m

文件是编辑好的m文件

%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;st=step_data';t=st(1,:);h=st(2,:);
load invdata1;y=yy;
[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,1,24];
x=curvefit('onedt1',x0,tt,hh)
y=onedt1(x,tt);err=sum((hh-y).^2)
save fitdata1 t h tt y x;
plot(t,h,'r*',tt,y,'b','linew',2);grid;%pause;
%plot(tt,hh-y);grid;maxerr=max(abs(hh-y));