Page 23 - 《中国药房》2021年第6期

EUS1＝Vs－Cs1 … … … … … … … … … … … … … … ⑦ 稳定策略，可以得到如图5C所示的复制动态相位图。
EUS2＝Vs－Cs2+RT s－μ（ksRT s+DS ）－（1－μ）[α（ksRT s+ 由 图 5 可 见 ，当 生 产 企 业 合 规 生 产 的 概 率 φ ＞
Ds ）+Is] … … … … … … … … … … … … … … … … … … … ⑧ （1－μ）kjIs－Cg
时，药政机关会逐渐减小监管力度α至
EUS＝φEUS1+（1－φ）EUS2＝Vs－φCs1－（1－φ）{Cs2－ （1－μ）kjIs
RT s+μ（ksRT s+Ds ）+（1－μ）[α（ksRT s+Ds ）+Is]} … … … … ⑨ 0；反之，药政机关会逐渐增大监管力度α至1。
2.3 行为选择演化稳定分析 dμ
（2）对于药检机构而言，令 F（μ）＝ ＝0，代入公
运用复制动态演化方式 ，得到三方主体的复制动 dt
[32]
态方程分别为： Is+Cj
式计算可得：μ＝0、μ＝1、α＝ 等3个均衡点。
dα kjIs+Dj
＝α（EUG1－EUG ）＝α（1－α）[（1－φ）（1－μ）kjIs－
dt 对F（μ）求导：F（μ）′＝（1－2μ）（1－φ）[α（kjIs+Dj ）－
Cg]… … … … … … … … … … … … … … … … … … … … … ⑩ Is－Cj]
dμ
＝μ（EUJ1－EUJ ）＝μ（1－μ）（1－φ）[α（kjIs+Dj ）－ 当α＝ Is+Cj 时，令 F（μ）′＝0，此时μ轴上的所有
dt kjIs+Dj
Is－Cj] … … … … … … … … … … … … … … … … … … … 点都是稳定状态，可以得到如图 6A 所示的复制动态相
dφ 位图。
＝φ（EUS1 －EUS ）＝φ（1－φ）{Cs2 －RT s －Cs1 +
dt dμ
dμ dμ
μ（ksRT s+Ds ）+（1－μ）[α（ksRT s+Ds ）+Is]} … … … … … … dt dt dt
dα
（1）对于药政机关而言，令F（α）＝ ＝0，代入公式
dt
0 μ 0 μ 0 μ
1 1 1
（1－μ）kjIs－Cg
⑩计算可得α＝0、α＝1、φ＝ 等 3 个均
Is+Cj Is+Cj Is+Cj
（1－μ）kjIs A. α ＝ 时 B. α ＞ 时 C. α ＜ 时
kjIs+Dj kjIs+Dj kjIs+Dj
衡点。
图6 药检机构的复制动态相位图
对F（α）求导：F（α）′＝（1－2α）[（1－φ）（1－μ）kjIs－
Fig 6 Replicator dynamic phase diagrams of drug
Cg]
inspection
（1－μ）kjIs－Cg
当φ＝ 时，令 F（α）′＝0，此时α轴上
Is+Cj 1
当α＞ 时，令 F（μ）′＝0，得μ＝ 。当 0≤
（1－μ）kjIs
的所有点都是稳定状态，可以得到如图 5A 所示的复制 kjIs+Dj 2
动态相位图。 1 1
μ＜ 时，F（μ）′＞0；当 ＜μ≤1时，F（μ）′＜0，此时μ＝1
2 2
dα dα dα
dt dt dt 是稳定策略，可以得到如图6B所示的复制动态相位图。
1
Is+Cj
当α＜ 时，令 F（μ）′＝0。当 0≤μ＜ 时，
0 α 0 α 0 α kjIs+Dj 2
1 1 1
1
F（μ）′＜0；当 ＜μ≤1时，F（μ）′＞0，此时μ＝0是稳定策
（1－μ）kjIs－Cg （1－μ）kjIs－Cg （1－μ）kjIs－Cg 2
A. φ＝ 时 B. φ＞ 时 C. φ＜ 时
略，可以得到如图6C所示的复制动态相位图。
（1－μ）kjIs （1－μ）kjIs （1－μ）kjIs
图5 药政机关的复制动态相位图 由 图 6 可 见 ，当 药 政 机 关 实 施 监 管 的 概 率 α ＞
Fig 5 Replicator dynamic phase diagrams of drug Is+Cj
时，药检机构会逐渐增加严格检测的概率μ至1；
administration kjIs+Dj
1 反之，药检机构会逐渐减小严格检测的概率μ至0。
（1－μ）kjIs－Cg
当φ＞ 时，令 F（α）′＝0，得α＝ 。
（1－μ）kjIs 2 dφ
（3）对于药品生产企业而言，令 F（φ）＝ ＝0，
1 1 dt
当 0≤α＜ 时，F（α）′＜0；当 ＜α≤1 时，F（α）′＞0，此
2 2 代 入 公 式 计 算 可 得 φ ＝0、φ ＝1、μ ＝
时α＝0 是稳定策略，可以得到如图 5B 所示的复制动态 －Cs2+RT s－Cs1－α（ksRT s+Ds ）－Is
等3个均衡点。
相位图。 ksRT s+Ds－α（ksRT s+Ds ）－Is
对 F（φ）求导：F（φ）′＝（1－2φ）{Cs2－RT s－Cs1+
（1－μ）kjIs－Cg
当φ＜ 时，令 F（α）′＝0。当 0≤α＜
（1－μ）kjIs μ（ksRT s+Ds ）+（1－μ）[α（ksRT s+Ds ）+Is]}
1 1 －Cs2+RT s－Cs1－α（ksRT s+Ds ）－Is
时，F（α）′＞0；当 ＜α≤1时，F（α）′＜0，此时α＝1是 当 μ ＝ 时 ，令
2 2 ksRT s+Ds－α（ksRT s+Ds ）－Is
中国药房 2021年第32卷第6期 China Pharmacy 2021 Vol. 32 No. 6 ·657 ·