代码如下：根据最小二乘原理计算。t=[10,15,20,30];y=[130,250,967,1823];res=lsqcurvefit(@(kar,t) kar(1)/1+kar(2)*exp(-kar(3)*t),[1 1 1],t,y)res = 792.5000 616.0802 73.0729意思是K=792.5000a= 616.0802r=73.0729
