一类扩充Goodwin振子的动力学行为+源代码(4)
{█(G_2∙P_1 □(→┴(k_2^Q ) ̃ ) G_2+Q_2 " 促进" @G_2 □(→┴(k_2^Q ) ̃ ) G_2+Q_2 "阻抑" )┤ (7)
本文考虑P_1对G_2转录的阻抑作用和P_2对G_1转录的促进作用。在促进的情况下，调控蛋白与基因的复合物转录；在阻抑的情况下，复合物不活跃，自由基因转录。
同理，酶促反应的动力学机理可以表示为以下形式：
S∙P_1 □(→┴k_3^P ) P_3 (8)
本文仅考虑P_1对S的正向催化作用，即P_1与S结合后，将加快S的反应速率。
转录产物mRNA Q_1和Q_2、调控蛋白质P_1和P_2以及酶促反应的产物P_3，由于降解反应，都只有有限的寿命。对于mRNA的翻译，mRNA、蛋白质和酶促反应的产物的降解，有如下的反应式：
Q_i □(→┴k_i^P ) Q_i+P_i "，" i=1,2，"翻译"
〖 Q〗_i □(→┴d_i^Q ) 0，i=1,2 , mRNA"降解" (9)
P_i □(→┴d_i^P ) 0，i=1,2 ,3 蛋白质"降解"
翻译和降解反应为单步骤过程，并且简单的降解反应总是一阶的。现在，我们给出关于Q_1，Q_2，P_1，P_2，P_3的常微分动力方程：
(dq_i)/dt=q ̇_i=k_i^Q F_i (p_j )-d_i^Q q_i，i=1,2，j=2,1
(dp_i)/dt=p ̇_i=k_i^P q_i-d_i^P p_i，i=1,2 (10)
(dp_3)/dt=(p_3 ) ̇=k_3^P F_3 (p_1 )-d_3^P p_3
该动力系统包含十个动力学参数k_1^Q 、k_2^Q 、k_3^S 、k_1^P 、k_2^P 、d_1^Q 、d_2^Q 、d_1^P 、d_2^P 、d_3^P和三个结合函数F_1 (p_2 )、F_2 (p_1 )、F_3 (p_1 )，系统是线性的。这个特征对于分析Jacobi矩阵和确定稳定点的稳定性是很重要的。
2 定性分析