www.pudn.com > 实验7-线形拟合.rar > fun076.m

function f=fun076(x,t1) 
t=[0.083 0.167 0.25 0.5 0.75 1.0 1.5 2.25 3 4 6 8 10 12];  
c=[10.9 21.1 27.3 36.4 35.5 38.4 34.8 24.2 23.6 15.7 8.2 8.3 2.2 1.8]; 
if nargin==1 
   f=c-x(1)*x(3)*(exp(-x(2)*t)-exp(-x(1)*t))/(x(1)-x(2)); 
else 
   plot(t,c,'r+') 
   f=x(1)*x(3)*(exp(-x(2)*t1)-exp(-x(1)*t1))/(x(1)-x(2)); 
end