第2期 曲卫华，等：企业、政府与公众公共健康提升激励机制演化分析 .245 研究第4)种情形。 2.1 选择某种策略初始情况变化对演化结果影响 数值实验 数值实验结果如图5，其中x、y和z分别表示污染 2 企业选择实行能源转型的初始比例、政府选择监管的 由于所列式子的复杂性，针对企业实行能源转 比例以及公众选择参与环境管理的比例。基于山西省 型后，政府对其减免税收的比率α、政府给予实行能 某发电企业和专家咨询，满足上文系统博弈收敛于平 源转型企业成本的补贴B,以及公众公共健康损失 衡点(1,1,1)的第4)种情形下参数值如下（单位：万 赔偿E,、政府对不实行能源转型企业的罚金F等政 元)：a=0.2,B=0.3,A=150,F=100,C1=520,C2= 策参数采用Matlab软件（版本：R2012b)对第4）中 420,C3=5,C,=8,E,=6,T=0.3,πe=1000,T。=40 情形做数值模拟，以期得到有益的结论为政府监管 T2=8,T3=3,T1=2,R1=1.5,E2=1,01=0.8,02=0.75, 相关部门及公共卫生部门的决策提供参考。 R2=2。从图5中得知排污企业群体、政府监管机构与 公众群体之间策略交往行为演化的路径依赖性。 1.0 1.0 ◇0aa网o866560000 1.0 oooooooooooc 0g学4*** 0.8 ▣=0.8 ★x0.6 ◆x=0.4 ，幸安马 。。职 x=0.2 0.2doo000oo0o0oo0000 06★ ▣=0.8 会以华称草学学区学草量单华华华草单的的的的论务 ▣1-0.8 0.4 04 ★x=0.6 0.2d00000000000000000nm000000000m0 0.4 ★x=0.6 ◆1x-0.4 =0.4 0x0.2 0x-0.2 0 0.0050.0100.0150.0200.0250.030 00.0050.0100.0150.0200.0250.030 0 0.0050.0100.0150.0200.0250.030 (a)y=0.2,0.8 (b)x=0.6z=0.8 (c)x-0.23-0.8 ◇0oooo98888688060 1.0 1.0 48 1.0 0.9 0.8 口1-0.8 0.8D 0.7 1=0.6 104 4444华溜 01=0.8 0.6 01=0.2 ooooooooooooooooooooooocpoo00 0.6 t=0.6 0.5 ★=0.6 ◆J=0.4 0.4 0.4 ◆x0.4000 0.4 0J-0.2 0.3 0-0.2 0.2 050.0100.0150.0200.0250.300.26eo00o0000000000 0.005 00.0050.0100.0150.0200.0250.030 0.0050.0100.0150.0200.0250.030 (a7y0.2,z0.3 (b)x0.3,z0.8 (c)x-0.9J=0.8 1.0 1.0 Lba4a66663 1.0 2000080oooooooooooxoooo ▣2=0.8 ▣-0.8 020.8 ★ -0.6 ★x=0.6 0.6 ★2=0.6 ◆20.4 ◆x-0.4 ◆2=0.4 0 2=0.2 SSSSO000 0 01x-0.2 。。8 0.4 0z=0.2 0 0.0050.01000150.0200.0250.030 00.0050.0100.0150.0200.0250.030 00.0050.0100.0150.0200.0250.030 (a")=0.6,2-=0.8 (b"-0.6.-0.2 (c")-0.2,=0.2 图5选择某一策略的不同初始情况对演化结果的影响 Fig.5 The effect of various initial conditions choosing a policy on evolution results 从图5中的(a)、(b)和(c)看出，x的收敛方向与 在x固定的条件下，z值越大，y收敛于平衡状态的速 速度和y、z有关，在y固定的条件下，z越大，x收敛速 度也越快。从图5中的(a)、(b)和(c)看出，在y固定 度越快，可以很快收敛于理想的稳定状态：但在：值较 的条件下，x取值越大，z收敛于平衡状态的速度要稍 高且固定的条件下，如z=0.8,从(a)、(c)可知，y越大， 慢：在x固定的条件下，y取值越大，z收敛于平衡状态 x收敛于稳定状态的速度越慢。从图5中的(a)、(b) 的速度也稍慢。从图5中的9个子图可看出：污染企 和(c)看出，y的收敛方向和速度与x、z有关，在z固 业、政府与公众群体各自行为的演化结果和收敛时间 定的条件下，x越大，y很快收敛于理想的稳定状态。 不仅受本群体中选择“能源转型”、“实行监管”与“参与研究第 ４）种情形。 ２ 数值实验 由于所列式子的复杂性，针对企业实行能源转 型后，政府对其减免税收的比率 α 、政府给予实行能 源转型企业成本的补贴 β ，以及公众公共健康损失 赔偿 Ｅ１ 、政府对不实行能源转型企业的罚金 Ｆ 等政 策参数采用 Ｍａｔｌａｂ 软件（版本：Ｒ２０１２ｂ） 对第 ４） 中 情形做数值模拟，以期得到有益的结论为政府监管 相关部门及公共卫生部门的决策提供参考。 ２．１ 选择某种策略初始情况变化对演化结果影响 数值实验结果如图５，其中 ｘ 、 ｙ 和 ｚ 分别表示污染 企业选择实行能源转型的初始比例、政府选择监管的 比例以及公众选择参与环境管理的比例。 基于山西省 某发电企业和专家咨询，满足上文系统博弈收敛于平 衡点（１，１，１）的第 ４）种情形下参数值如下（单位：万 元）： α ＝０．２， β ＝ ０．３， Ａ ＝ １５０， Ｆ ＝ １００， Ｃ１ ＝ ５２０， Ｃ２ ＝ ４２０， Ｃ３ ＝５， Ｃ４ ＝８， Ｅ１ ＝ ６， τ ＝ ０．３， πＦ ＝ １０００， π０ ＝ ４０， π２ ＝ ８， π３ ＝３， π１ ＝２， Ｒ１ ＝１．５， Ｅ２ ＝１， θ１ ＝０．８， θ２ ＝０．７５， Ｒ２ ＝２。 从图 ５ 中得知排污企业群体、政府监管机构与 公众群体之间策略交往行为演化的路径依赖性。 图 ５ 选择某一策略的不同初始情况对演化结果的影响 Ｆｉｇ．５ Ｔｈｅ ｅｆｆｅｃｔ ｏｆ ｖａｒｉｏｕｓ ｉｎｉｔｉａｌ ｃｏｎｄｉｔｉｏｎｓ ｃｈｏｏｓｉｎｇ ａ ｐｏｌｉｃｙ ｏｎ ｅｖｏｌｕｔｉｏｎ ｒｅｓｕｌｔｓ 从图 ５ 中的（ａ）、（ｂ）和（ｃ）看出， ｘ 的收敛方向与 速度和 ｙ 、 ｚ 有关，在 ｙ 固定的条件下， ｚ 越大， ｘ 收敛速 度越快，可以很快收敛于理想的稳定状态；但在 ｚ 值较 高且固定的条件下，如 ｚ ＝０．８，从（ａ）、（ｃ）可知， ｙ 越大， ｘ 收敛于稳定状态的速度越慢。 从图 ５ 中的（ａ）、（ｂ） 和（ｃ）看出， ｙ 的收敛方向和速度与 ｘ 、 ｚ 有关，在 ｚ 固 定的条件下， ｘ 越大， ｙ 很快收敛于理想的稳定状态。 在 ｘ 固定的条件下， ｚ 值越大， ｙ 收敛于平衡状态的速 度也越快。 从图 ５ 中的（ａ）、（ｂ）和（ｃ）看出，在 ｙ 固定 的条件下， ｘ 取值越大， ｚ 收敛于平衡状态的速度要稍 慢；在 ｘ 固定的条件下， ｙ 取值越大， ｚ 收敛于平衡状态 的速度也稍慢。 从图 ５ 中的 ９ 个子图可看出：污染企 业、政府与公众群体各自行为的演化结果和收敛时间 不仅受本群体中选择“能源转型”、“实行监管”与“参与 第 ２ 期 曲卫华，等：企业、政府与公众公共健康提升激励机制演化分析 ·２４５·