
输入:t=[0.25 0.5 1 1.5 2 3 4 6 8]
t =
0.2500 0.5000 1.0000 1.5000 2.0000 3.0000 4.0000 6.0000 8.0000
输入:c=[19.21 18.15 15.36 14.10 12.89 9.32 7.45 5.24 3.01]
c =
19.2100 18.1500 15.3600 14.1000 12.8900 9.3200 7.4500 5.2400 3.0100
输入:y=log(c)
y =
2.9554 2.8987 2.7318 2.6462 2.5565 2.2322 2.0082 1.6563 1.1019
输入:a=polyfit(t,y,1)
a =
-0.2347 2.9943
/a(2)为a返回值中的常数lnc0输入:c0=exp(a(2))
c0 =
19.9709
/a(1)为a返回值中的x的系数a(1)=-k输入:k=-a(1)
k =
0.2347
输入:t1=0.25:0.01:8
输入:c1=c0exp(-kt1)
输入:plot(t,c,’k+’,’LineWidth’, 2)
输入:hold on
输入:plot(t1,c1,’LineWidth’, 2)

方法2:
输入:t=[0.25 0.5 1 1.5 2 3 4 6 8]
t =
0.2500 0.5000 1.0000 1.5000 2.0000 3.0000 4.0000 6.0000 8.0000
输入:c=[19.21 18.15 15.36 14.10 12.89 9.32 7.45 5.24 3.01]
c =
19.2100 18.1500 15.3600 14.1000 12.8900 9.3200 7.4500 5.2400 3.0100
输入:y=log(c)
y =
2.9554 2.8987 2.7318 2.6462 2.5565 2.2322 2.0082 1.6563 1.1019
输入:a=polyfit(t,y,1)
a =
-0.2347 2.9943
输入:k=-a(1)
k =
0.2347
输入:d=300
输入:v=d/exp(a(2))
v =
15.0219
发表回复