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 ·
   20   21   22   23   24   25   26   27   28   29   30