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

40、kjIs＝550、M＝800、Cj＝15、N＝50、Cs2＝120、RT s＝15、
1.0
α Cs1＝25、Is＝30，即药品生产企业合规生产的成本 Cs1小
μ
0.8 φ
于不合规生产的成本 Cs2加上惩罚和名誉损失之和 N 减
概率 0.6 去所获得的额外收益RT s。绘制药政机关的监管概率α、
药检机构严格检测的概率μ以及药品生产企业合规生产
0.4
的概率φ随着时间变化的动态演化仿真图，详见图11。
0.2
1.0
0 α
μ
0.8
φ
0 0.1 0.2 0.3 0.4 0.5 0.6
时间 0.6
图 9 药政机关、药检机构与药品生产企业的行为策略 概率
0.4
演化仿真图[E2 （1，0，0）]
Fig 9 Simulation diagram of behavior strategy evolu- 0.2
tion of drug administration，drug inspection 0
and manufacturing enterprises[E2 （1，0，0）] 0 0.025 0.050 0.075 0.100 0.125 0.150 0.175 0.200
时间
（3）在均衡点E3 （0，0，1）的状态下，试验参数的初始
图11 药政机关、药检机构与药品生产企业的行为策略
化赋值需满足－Cs2+RT s+Cs1－Is＜0，具体赋值如下：Cg＝
演化仿真图[E7 （0，1，1）]
60、kjIs＝20、M＝30、Cj＝40、N＝50、Cs2＝20、RT s＝25、
Fig 11 Simulation diagram of behavior strategy evolu-
Cs1＝25、Is＝35，即药品生产企业合规生产的成本 Cs1小
tion of drug administration，drug inspection
于不合规生产的成本 Cs2与贿赂花费 Is之和减去所获得
and manufacturing enterprises[E7 （0，1，1）]
的额外收益RT s。绘制药政机关的监管概率α、药检机构
从图 11 可以看出，随着时间的推移，药品生产企业
严格检测的概率μ以及药品生产企业合规生产的概率φ
合规生产的概率φ和药检机构严格检测的概率μ逐渐升
随着时间变化的动态演化仿真图，详见图10。
高，最终无限趋近于1；药政机关的监管概率α逐渐降低，
1.0 最终无限趋近于0。也就是说，当满足上述条件时，药政
α
μ
0.8 机关、药检机构以及药品生产企业的行为策略最终会稳
φ
定在{不监管，严格检测，合规}的状态。这与上述演算
0.6
概率 结果相一致，充分验证了模型结论的科学可信度。
0.4
2.6 模型结论
0.2
综合上述三方演化博弈模型的建立、求解与仿真模
0
拟试验，可以得到如下模型结论：
0 0.2 0.4 0.6 0.8 1.0 1.2
时间 Is+Cj
（1）当药政机关实施监管的概率α＞ 时，药
图10 药政机关、药检机构与药品生产企业的行为策略 kjIs+Dj
演化仿真图[E3 （0，0，1）] 检机构会逐渐增加严格检测的概率μ至 1。当药检机构
Fig 10 Simulation diagram of behavior strategy evolu- －Cs2+RT s－Cs1－α（ksRT s+Ds ）－Is
严格检测的概率μ＞
tion of drug administration，drug inspection ksRT s+Ds－α（ksRT s+Ds ）－Is
and manufacturing enterprises[E3 （0，0，1）] 时，药品生产企业合规生产的概率φ会逐渐增加至 1；当
从图 10 可以看出，随着时间的推移，药品生产企业 （1－μ）kjIs－Cg
药品生产企业合规生产的概率φ＞ 时，
合规生产的概率φ逐渐升高，最终无限趋近于 1；药政机 （1－μ）kjIs
关的监管概率α和药检机构严格检测的概率μ逐渐降低， 药政机关会逐渐减小监管的概率α至0；反之亦然。
最终无限趋近于0。也就是说，当满足上述条件时，药政 （2）药政机关实施监管的概率α与监管成本 Cg成反
机关、药检机构以及药品生产企业的行为策略最终会稳 比。药检机构严格检测的概率μ与日常合理收益 Vj无
定在{不监管，不严格检测，合规}的状态，这与上述演算 关，与监管强度kj、处罚罚款kjIs和名誉损失Dj成正比，与
结果相一致。 检测成本Cj成反比。药品生产企业合规生产的概率φ与
（4）在均衡点E7 （0，1，1）的状态下，试验参数的初始 日常合理收益 Vs无关，与不合规生产成本 Cs2、贿赂花费
化赋值需满足－Cs2+RT s+Cs1－N＜0，具体赋值如下：Cg＝ Is、监管强度ks以及处罚罚款ksRT s和名誉损失Ds成正比，

·660 · China Pharmacy 2021 Vol. 32 No. 6 中国药房 2021年第32卷第6期