Page 25 - 《中国药房》2021年第6期
P. 25
由于-Cg<0,因此当满足条件-Cs2+RT s+Cs1-Is<0 化赋值需满足kjIs-Cg<0、Cs2-RT s-Cs1+Is<0,具体赋值
时,系统可以到达稳定状态。也就是说,当药品生产企 如下:Cg=60、kjIs=20、M=30、Cj=45、N=50、Cs2=20、
业合规生产的成本小于不合规生产的成本与贿赂花费 RT s=25、Cs1=35、Is=30,即药政机关监管成本Cg大于药
之和减去所获得的额外收益时,E2 (0,0,1)是稳定点,此 检机构受贿的处罚 kjIs,药品生产企业合规生产的成本
时药政机关、药检机构与药品生产企业的行为策略为 Cs1大于不合规生产的成本 Cs2与贿赂花费 Is之和减去所
{不监管,不严格检测,合规}。 获得的额外收益RT s。绘制药政机关的监管概率α、药检
(4)将E4 (1,0,1)代入矩阵J,可以得到矩阵: 机构严格检测的概率μ以及生产企业合规生产的概率φ
Cg 0 0 随着时间变化的动态演化仿真图,详见图8。
0 0 0 α
1.0 μ
0 0 -Cs2+RT s+Cs1-N-Is φ
0.8
由于Cg>0,因此E4 (1,0,1)是不稳定点。
0.6
(5)将E5 (0,1,0)代入矩阵J,可以得到矩阵: 概率
0 0 0.4
-Cg
0 Is+Cj 0 0.2
0 0 Cs2-RT s-Cs1+N 0
由于Is+Cj>0,因此E5 (0,1,0)是不稳定点。 0 0.1 0.2 0.3 0.4 0.5
时间
(6)将E6 (1,1,0)代入矩阵J,可以得到矩阵:
图 8 药政机关、药检机构与药品生产企业的行为策略
0 0
演化仿真图[E1 (0,0,0)]
Cg
0 -M+Is+Cj 0 Fig 8 Simulation diagram of behavior strategy evolu-
0 0 Cs2-RT s-Cs1+N tion of drug administration,drug inspection
由于Cg>0,因此E6 (1,1,0)是不稳定点。 and manufacturing enterprises[E1 (0,0,0)]
(7)将E7 (0,1,1)代入矩阵J,可以得到矩阵:
从图8可以看出,随着时间的推移,药品生产企业合
0 0
-Cg 规生产的概率φ、药政机关的监管概率α和药检机构严格
0 0 0 检测的概率μ均逐渐降低,最终无限趋近于 0。也就是
0 0 -Cs2+RT s+Cs1-N
说,当满足上述条件时,药政机关、药检机构以及药品生
由于-Cg<0,因此当满足条件-Cs2+RT s+Cs1-N<0
产企业的行为策略最终会稳定在{不监管,不严格检测,
时,系统可以到达稳定状态。也就是说,当生产企业合
不合规}的状态,这与上述演算结果相一致。
规生产的成本小于不合规生产的成本与所受处罚、名誉 (2)在均衡点E2 (1,0,0)的状态下,试验参数的初始
损失之和减去所获得的额外收益时,E7 (0,1,1)是稳定 化赋值需满足Cg-kjIs<0、M-Is-Cj<0、Cs2-RT s-Cs1+
点,此时药政机关、药检机构与生产企业的行为策略为 N+Is<0,具体赋值如下:Cg=20、kjIs=40、M=55、Cj=
{不监管,严格检测,合规}。 45、N=45、Cs2=15、RT s=40、Cs1=45、Is=20,即药政机关
(8)将E8 (1,1,1)代入矩阵J,可以得到矩阵: 监管成本 Cg小于药检机构受贿的处罚 kjIs,药检机构受
0 0
Cg 贿的惩罚和名誉损失M小于检测成本Cj和贿赂Is之和,
0 0 0 药品生产企业合规生产的成本Cs1大于不合规生产的成
0 0 -Cs2+RT s+Cs1-N 本Cs2、贿赂花费Is以及惩罚与名誉损失之和N减去所获
由于Cg>0,因此E8 (1,1,1)是不稳定点。 得的额外收益RT s。绘制药政机关的监管概率α、药检机
由此,可以得到 4 个趋于稳定的博弈均衡点:E1 (0, 构严格检测的概率μ以及药品生产企业合规生产的概率
0,0)、E2 (1,0,0)、E3 (0,0,1)、E7 (0,1,1)。 φ随着时间变化的动态演化仿真图,详见图9。
2.5 模拟仿真与分析 从图9可以看出,随着时间的推移,药政机关的监管
上述结论是在行为参数的基础上计算所得,为了提 概率α逐渐升高,最终无限趋近于 1;药品生产企业合规
高模型结论的科学可信度,需要运用 MATLAB R2017a 生产的概率φ和药检机构严格检测的概率μ逐渐降低,最
软件进行仿真模拟试验,以更加清晰地展现出药政机 终无限趋近于0。也就是说,当满足上述条件时,药政机
关、药检机构以及生产企业三方主体的行为策略演化 关、药检机构以及药品生产企业的行为策略最终会稳定
趋势。 在{监管,不严格检测,不合规}的状态,这与上述演算结
(1)在均衡点E1 (0,0,0)的状态下,试验参数的初始 果相一致。
中国药房 2021年第32卷第6期 China Pharmacy 2021 Vol. 32 No. 6 ·659 ·