用以下一组数据拟合上述模型中的参数k1、k2:
M=64000/490 =130.6122 (毫克/百毫升)
建立函数文件:
function f=curvefun1(k,t)
f=k(1)*64000/490*(exp(-k(2)*t)-exp(-k(1)*t))/(k(1)-k(2))
输入拟合数据:
t=[0 0.25 0.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5 6 7 8 9 10 11 12 13 14 15 16];
c=[0 30 68 75 82 82 77 68 68 58 51 50 41 38 35 28 25 18 15 12 10 7 7 4];
任取k1、k2的一组初始值:k0=[2,1];
时间(小时) 0.25 0.5 0.75 1 1.5 2 2.5 3 3.5 4 4.5 5
酒精含量 30 68 75 82 82 77 68 68 58 51 50 41
时间(小时) 6 7 8 9 10 11 12 13 14 15 16
酒精含量 38 35 28 25 18 15 12 10 7 7 4
1