平均相对权益损失由以下公式得出:Eloss=PiE*我- Ei(0)PjEprej(5),其中损失是相对于初始预冲击值Eprei计算的。或者,对于动力学的给定步骤t,我们考虑权益的相对损失hi(t)和每家银行的银行间杠杆矩阵∧ij(t):hi(t)=Ei(0)- Ei(t)Ei(0)(6)∧ij(t)=wi公司→j(0)Ei(0),如果银行j在时间(t)之前未违约- 1) 0否则(7),其中Ei(0)是i的初始权益。假设违约损失为100%,相对权益损失hi(t)的动态方程为:hi(t+1)=min1,hi(t)+NXj=1∧ij(t)[pDj(t+1)- pDj(t)](8) 式中,pDj(t)=hj(t)eα[hj(t)-1] 是步骤t中银行j的违约概率。控制参数α∈ (0, ∞) 允许从线性债务等级(α=0,意味着违约概率与股权损失成正比)(Bardocia等人,2015)持续切换到最终算法(α→ ∞, i、 例如,只有当股本耗尽时,才会发生违约,否则银行不会传染)(Fur fine,2003)。冲击传播动力学t结束时的平均相对权益损失*is:Eloss=Xi[高(t*) - hi(1)]Ei(0)PjEj(0)(9),其中hi(1)是对i的初始冲击。因此,Eloss不考虑对系统的初始冲击,而只考虑网络对系统损失的影响。我们使用两种停止条件进行模拟:当差值k[h(t)时-h(t- 1) ]E(0)kbecomes小于公差tol=10-5,或当相互作用的数量等于maxiter=10.1031021010.10.20.30.40.50.50.60.60.80.40.520130.00.20.40.60.81.0ElossFigure 2:Eloss(使用线性DebtRank计算)作为链路密度和均匀冲击幅度的函数,对于使用公式(3)构建的2008年(左)和2013年(右)数据,φ=1。
|