在求得了各个温度下氧气和二氧化碳浓度与时间的拟合方程后，就可以根据公式[5]~[8]，求得氧气的呼吸消耗速率ro2和二氧化碳的呼吸生成速率rco2，拟合曲线见图1；根据公式[7]、[8]通过多重线性回归分析，获得表3结果：
表3 香菇在273k、283k和293k下根据多重回归求得米氏方程各参数值
温度
k
km
%o2
vmax
mg/kghr
ki
%co2
r2
标准误差
（std.error）
293
0.27
981
1.33
0.985
0.0059
0.23
828
1.58
0.985
0.000034
283
0.4
743
2.86
0.989
0.060
0.14
314
8.20
0.989
0.000038
273
0.05
117
8.11
0.973
0.22
0.017
57.4
23.58
0.973
0.00015
0.017 57.4 23.58 0.973 0.00015
注：2、4、6行数据由公式[7]求得，3、5、7行由公式[8]求得。
表3中，每组方程各项的标准回归化系数（beta）在软件中显示高亮度（即p-level<0.05）说明线性显著。
图 1 o2、co2浓度对香菇呼吸速度的影响
fig.1 effects of o2,co2 concentration on respiratory rate of shiitake
由图1可以看出在293k条件下，香菇的由始至终呼吸熵在0.86~1.02之间变化，其中，当氧气浓度在1.0%和二氧化碳浓度在13.4%条件下呼吸熵约为1；在283k条件下，香菇的呼吸熵在0.4~1.6之间变化，当氧气浓度为5.2%和二氧化碳浓度为7.2%，rq约为1；273k条件下rq区间为0.4~1.03，当氧气浓度为0.7%和二氧化碳浓度为8.4%，rq约为1。由此表明本实验是在有氧呼吸的条件下进行的，计算有效。
2.3 温度对香菇呼吸作用的影响
三个温度下lnvm的值统计见表4：
表 4 不同1/t下，lnvm 的值
1/温度（k）
1/273
1/283
1/293
lnvm1
4.76
6.61
6.89
lnvm2
4.05
5.75
6.72
不同1/t下，lnvm线性回归预测值与观察值关系中氧气斜率=-0.93，复相关系数=0.927；二氧化碳斜率=0.99，复相关系数=0.990。
通过计算求得的在氧气最大消耗率活化能为7.73kj/mol；二氧化碳最大产生率的活化能为8.23kj/mol。根据公式[9]，在已知一定温度下氧气的最大消耗速率或二氧化碳的最大产出率，就可以求得其他温度条件下的最大变化率。例如：在有氧呼吸条件下，已知香菇在273k下氧气的最大消耗率为117mg/kghr，根据公式[9]代入相关参数值可求得在288k下氧气最大消耗率为139.7 mg/kghr。
3 结果与讨论
本次实验利用米氏方程机理，选取了具有一定特征的函数形式，拟合了一定温度下氧气、二氧化碳浓度随时间变化的方程，并且复相关系数都达到了0.95以上。
米氏方程的各个参数值为进一步气调包装的参数设定提供理论依据，并且方程原始模型仍为时间的函数，大大简化了由于呼吸作用而导致体系的复杂性。实验中同时也获得了对呼吸熵的动态变化。
温度因子对香菇呼吸作用的影响则通过arrhenius关系求得ea，同时求得在有氧条件下不同温度下的氧气的最大消耗率或二氧化碳的最大产出率，为以后在不同温度下进行气调包装提供数值参考。
作者：刘晓丹 谢晶