交易对手信用风险与CVA - MATLAB与Simulink万博1manbetx - 万博1manbetx,s manbetx 845,万博尤文图斯

交易对手信用风险和CVA

打开生活的脚本

这个例子展示了如何计算单边信贷价值(估值)调整(CVA)的银行持有一个投资组合的香草利率掉期与几个交易对手。CVA是由于交易对手违约而导致的场外合约或合约组合的预期损失。特定交易对手的CVA定义为每一时刻的预期风险折现乘以该交易对手在该时刻违约的概率，所有这些乘以1减去回收率。CVA公式为:

$C V 一个＝（ 1 - R ） \int_{0}^{T} d 我年代 c E E （ t ） d P D （ t ）$

在哪里R经济复苏,discEE贴现的预期风险敞口t,PD违约概率分布。

预期风险敞口的计算首先是模拟给定合同或投资组合的许多未来风险因素情景。风险因素可以是利率，如本例所示，但会根据投资组合而异，包括外汇汇率、股票或商品价格，或任何会影响合约市场价值的因素。一旦模拟了足够多的情景，合约或投资组合就可以在每个情景的一系列未来日期进行定价。结果是一个由契约值组成的矩阵或“立方体”。

这些价格在考虑到银行可能存在的担保协议和净额协议(如在本例中，几份合同的价值可能相互抵消，从而降低它们的总风险)后转换为风险敞口。

对每个情景的合约值进行折现以计算折现的风险敞口。然后，可以通过每个模拟日期的贴现敞口的简单平均值来计算贴现预期敞口。

最后，交易对手违约概率一般由信用违约互换(CDS)市场报价得出，交易对手的CVA可根据上述公式计算。假设交易对手违约独立于其风险敞口(没有错误风险)。

这个例子演示了一个香草利率掉期的投资组合，目标是计算特定交易对手的CVA。

读交换组合

掉期投资组合的价值有时接近于零t＝0．每次掉期都与一个交易对手有关，可能包含在净额协议中，也可能不包含在净额协议中。

从电子表格读取交换swapFile =“cva-swap-portfolio.xls”；掉期= readtable (swapFile,“表”，“交换组合”）;互换。LegType =[互换。LegType~swaps.LegType]; swaps.LegRate = [swaps.LegRateReceiving swaps.LegRatePaying]; swaps.LegReset = ones(size(swaps,1),1); numSwaps = size(swaps,1);

的交换参数的更多信息CounterpartyID和NettingID,请参阅creditexposures．的交换参数的更多信息主要，成熟，LegType，LegRate，LatestFloatingRate，期,LegReset,请参阅swapbyzero．

从利率曲线创建利率

解决= datenum (的14 - 12月- 2007）;期限= [3 6 12 5*12 7*12 10*12 20*12 30*12]';ZeroRates = [0.033 0.034 0.035 0.040 0.042 0.044 0.048 0.0475]';男高音ZeroDates = datemnth(解决);复合= 2;基础= 0;RateSpec = intenvset (startdate可以的解决,“EndDates”ZeroDates,．..“利率”ZeroRates,“复合”复合,“基础”、基础);图;情节(ZeroDates ZeroRates,“啊——”）;包含(“日期”）;datetick (“keeplimits”）;ylabel (“零率”）;网格在；标题(“结算日的收益率曲线”）;

设置可变的仿真参数

您可以改变生成的模拟利率场景的数量。本示例首先将模拟日期设置为更频繁的日期，然后在将来变得更不频繁。

%蒙特卡罗模拟次数numScenarios = 1000;%计算每月的模拟日期，然后再计算季度的日期。0:12 simulationDates = datemnth(解决);simulationDates = [simulationDates datemnth(simulationDates(end)，3:3:74)]';numDates =元素个数(simulationDates);

计算浮动复位日期

对于每个模拟日期，计算每个交换之前的浮动重置日期。

floatDates = cfdates(解决- 360,swaps.Maturity swaps.Period);互换。FloatingResetDates = zeros(numSwaps,numDates);为i = numDates:-1:1 thisDate = simulationDates(i);floatDates(floatDates > thisDate) = 0;swaps.FloatingResetDates(:,我)= max (floatDates [], 2);结束

建立Hull-White单因素模型

模拟为合约估值的风险因素是零曲线。对于本例，您使用单因素赫尔-怀特模型对利率期限结构进行建模。这是一个短期利率模型，其定义为:

$d r ＝［ θ （ t ） - 一个 r ］ d t + σ d z$

在哪里

$d r$ :短期汇率在时间上发生微小变化后的变化， $d t$
$一个$ :平均回归率
$σ$ :短期利率的波动性
$d z$ :韦纳过程(标准的正常过程)
$θ （ t ）$ :漂移函数定义为:

$θ （ t ）＝ F_{t} （ 0 ， t ） + 一个 F （ 0 ， t ） + \frac{σ^{2}}{2 一个} （ 1 - e^{- 2 一个 t} ）$

$F （ 0 ， t ）$ :某一时刻的瞬时转发速率 $t$

$F_{t} （ 0 ， t ）$ :的偏导数 $F$ 关于时间

一旦你模拟了短期利率的路径，在每个模拟日期使用公式生成一个完整的收益率曲线:

$R （ t ， T ）＝ - \frac{1}{（ T - t ）} \ln 一个（ t ， T ） + \frac{1}{（ T - t ）} B （ t ， T ） r （ t ）$

$\ln 一个（ t ， T ）＝ \ln \frac{P （ 0 ， T ）}{P （ 0 ， t ）} + B （ t ， T ） F （ 0 ， t ） - \frac{1}{4 {一个}^{3.}} σ^{2} （ e^{- 一个 T} - e^{- 一个 t} ）^{2} （ e^{2 一个 t} - 1 ）$

$B （ t ， T ）＝ \frac{1 - e^{- 一个（ T - t ）}}{一个}$

$R （ t ， T ）$ :零利率在时间 $t$ 一段时间 $T - t$

$P （ t ， T ）$ :零息债券在同一时期的价格 $t$ 每次一美元 $T$

每个场景都包含了整个期限结构，并在我们选择的每个模拟日期进行建模。

指的是利用市场数据校准赫尔白模型有关Hull-White单因素模型校准的更多细节，请参见金融工具工具箱™用户指南中的示例。

α= 0.2;σ= 0.015;hw1 = HullWhite1F (RateSpec、ασ);

模拟的场景

对于每个场景，使用赫尔-怀特单因素利率模型模拟每个估值日期的未来利率曲线。

%使用可重复随机数生成器(改变种子以产生%不同的随机场景)。prevRNG = rng (0,“旋风”）;dt = diff (yearfrac(解决simulationDates 1));nPeriods =元素个数(dt);场景= hw1.simTermStructs (nPeriods,．..“nTrials”numScenarios,．..“deltaTime”, dt);%恢复随机数生成器状态rng (prevRNG);通过每个已实现的利率计算贴现率%的场景。dfactors = 1 (numDates numScenarios);为i = 2:numDates tenorDates = datemnth(simulationDates(i-1)，Tenor);rateAtNextSimDate = interp1 (tenorDates挤压(场景(张,,,:)),．..simulationDates(我),“线性”，“extrap”）;%计算D (t1, t2)dfactors(我)= zero2disc (rateAtNextSimDate,．..repmat (simulationDates (i), 1 numScenarios), simulationDates(张),1,3);结束dfactors = cumprod (dfactors, 1);

检查一个场景

创建一个特定场景的收益率曲线演变的表面图。

我= 20;图;surf(Tenor, simulationDates, scenario (:，:，i))轴紧datetick (“y”，“mmmyy”）;包含(的男高音(个月)）;ylabel (“观察日期”）;zlabel (“利率”）;甘氨胆酸ax =;斧子。6ew = [-49 32]; title(sprintf(“情景%d收益率曲线演变”\n”,我));

计算按市价计算掉期价格

对于每一种情形，掉期投资组合都在每个未来模拟日期定价。价格是用价格近似函数计算的，hswapapprox．由于蒙特卡罗模拟的性能要求，在CVA应用程序中使用简化的近似函数对合同进行定价是很常见的。

由于模拟日期不对应于互换现金流日期(其中浮动利率被重置)，估计最新的浮动利率与1年利率(所有互换的期限为1年)之间插值的最接近的模拟利率曲线。

然后将掉期价格聚合成一个值的“立方体”，其中包含每种情景的每个模拟日期的所有未来合约值。合同价格的立方体是一个三维矩阵，其中每一行代表一个模拟日期，每一列代表一个合同，每一“页”代表一个不同的模拟场景。

%计算此场景的所有按市价计算的值。使用一个%近似函数在这里改善性能。值= hcomputeMTMValues (simulationDates互换,场景,男高音);

检查情况下价格

创建一个特定情景下所有掉期价格演变的图表。

我= 32;图;情节(simulationDates值(:,:,i));datetick;ylabel (“盯市价格”）;标题(sprintf (“按照情景%d互换价格”,我));

可视化模拟投资组合价值

绘制模拟每个场景的总投资组合价值。随着每个情景的推移，合约的价值会根据模型利率期限结构的变化而上下波动。当掉期接近到期时，其价值将开始趋近于零，因为所有剩余现金流的总价值将在每个现金流日后减少。

%随时间查看投资组合价值图;totalPortValues = squeeze(sum(values, 2));情节(simulationDates totalPortValues);标题(“所有场景的MTM投资组合总价值”）;datetick (“x”，“mmmyy”) ylabel (“投资组合价值(美元)”)包含(“模拟日期”）

按交易对手计算风险敞口

个别合约的曝光(我)时间t为合约最高金额(6),0：

$E_{我} （ t ）＝马克斯｛ V_{我} （ t ）， 0 ｝$

而特定交易对手的风险敞口只是单个合同风险敞口的总和:

$E_{c p} （ t ）＝ \sum E_{我} （ t ）＝ \sum 马克斯｛ V_{我} （ t ）， 0 ｝$

但是，在有净额协定的情况下，合同是合计在一起的，可以相互抵消。因此，净额协议中所有合同的总风险敞口为:

$E_{n 一个} （ t ）＝马克斯｛ \sum V_{我} （ t ）， 0 ｝$

计算整个投资组合以及每个交易对手在每个模拟日期的风险敞口creditexposures函数。

未联网合同用a表示南在NettingID向量。未净合约的风险敞口，如果为正值，则等于该合约的市场价值，否则为零。

包括在净值协议中的合同的价值加在一起，可以相互抵消。有关按市价计算风险敞口的更多细节，请参阅参考资料。

[exposure, expcpty] =信用风险(价值，掉期。CounterpartyID,．..“NettingID”, swaps.NettingID);

在我们的模拟中绘制每个场景的总投资组合风险。与合约价值图相似，随着互换到期，每种情况下的风险敞口将接近于零。

%随着时间的推移查看投资组合的风险敞口图;totalPortExposure =挤压(sum(敞口,2));情节(simulationDates totalPortExposure);标题(“所有情况下的投资组合敞口”）;datetick (“x”，“mmmyy”) ylabel (敞口(美元))包含(“模拟日期”）

曝光资料

在分析银行对交易对手的潜在未来风险敞口时，有几种风险敞口是有用的。在这里，你可以计算每个交易对手的几个(非贴现)风险敞口，以及整个投资组合。

pfizer(未来潜在暴露量):在未来某一特定日期暴露量分布的高百分位数(95%)(也称为峰值暴露量(PE))
MPFE(最大潜在未来暴露):所有日期的最大PFE
EE:(预期曝光率):每个日期曝光率分布的平均值
EPE(预期正照射量):预期照射量随时间的加权平均值
EffEE(有效预期暴露量):任何时候的最大预期暴露量，t，或以前的时间
EffEPE(有效预期正向暴露):有效预期暴露的加权平均值

有关进一步的定义，请参见参考文献中的巴塞尔协议II文件。

计算全部投资组合风险portExposures =总和(曝光,2);%计算每个交易对手和整个投资组合的风险敞口cpProfiles = exposureprofiles (simulationDates,曝光);portProfiles = exposureprofiles (simulationDates portExposures);

首先是整个投资组合，然后是一个特定的交易对手。

%可视化投资组合暴露概况图;情节(simulationDates portProfiles。工业的．..simulationDates portProfiles。MPFE *的(numDates, 1),．..simulationDates portProfiles。情感表达,．..simulationDates portProfiles。EPE *的(numDates, 1),．..simulationDates portProfiles。EffEE,．..simulationDates portProfiles。EffEPE *的(numDates, 1));传奇({“PFE(95%)”，“马克斯pfizer”，“Exp曝光(EE)”，“Time-Avg EE (EPE)”，．..'Max past EE (EffEE)'，“Time-Avg EffEE (EffEPE) '}) datetick (“x”，“mmmyy”)标题(的投资组合风险敞口概要”）;ylabel (敞口(美元))包含(“模拟日期”）

可视化特定交易对手的暴露概况。

cpIdx = find(expcpty == 5);图;情节(simulationDates cpProfiles (cpIdx)。工业的．..simulationDates cpProfiles (cpIdx)。MPFE *的(numDates, 1),．..simulationDates cpProfiles (cpIdx)。情感表达,．..simulationDates cpProfiles (cpIdx)。EPE *的(numDates, 1),．..simulationDates cpProfiles (cpIdx)。EffEE,．..simulationDates cpProfiles (cpIdx)。EffEPE *的(numDates, 1));传奇({“PFE(95%)”，“马克斯pfizer”，“Exp曝光(EE)”，“Time-Avg EE (EPE)”，．..'Max past EE (EffEE)'，“Time-Avg EffEE (EffEPE) '}) datetick (“x”，“mmmyy”，“keeplimits”)标题(sprintf (“交易对手%d敞口概况”cpIdx));ylabel (敞口(美元))包含(“模拟日期”）

贴现风险敞口

使用每个模拟利率场景的贴现因子计算贴现预期风险敞口。给定场景中给定估值日期的贴现因子是从一个模拟日期到下一个模拟日期的增量贴现因子的乘积，以及该场景的利率路径。

%在每一种情况下，获得每个交易对手的贴现敞口discExp = 0(大小(曝光));为i = 1: numScenarios discExp(:,:我)= bsxfun (@times, dfactors(:,我),曝光(:,:我));结束预期风险折现率discProfiles = exposureprofiles (simulationDates discExp,．..“ProfileSpec”，“EE”）;

绘制总投资组合以及每个交易对手的贴现预期风险敞口。

将每个交易对手的折扣EE汇总成一个矩阵discEE = [discProfiles.EE];投资组合折现EE图;情节(simulationDates sum (discEE, 2)) datetick (“x”，“mmmyy”，“keeplimits”)标题(“投资组合预期风险折现”）;ylabel (的贴现风险敞口(美元))包含(“模拟日期”）

交易对手折扣EE图;情节(simulationDates discEE) datetick (“x”，“mmmyy”，“keeplimits”)标题(“每个交易对手的预期风险折现”）;ylabel (的贴现风险敞口(美元))包含(“模拟日期”）

每一方违约概率曲线的校准

给定交易对手的违约概率是由该交易对手CDS的当前市场价差所暗示的。使用的函数cdsbootstrap以生成每个模拟日期的累积违约概率。

为每一交易对手进口CDS市场信息cd = readtable (swapFile,“表”，“CDS利差”）;disp (CDS);

Date cp1 cp2 cp4 cp5 _____________ ____________ ___ {'3/20/2008'} 140 85 115 170 140 {'3/20/2009'} 185 120 150 205 175 {'3/20/2010'} 215 170 195 245 210 {'3/20/2011'} 275 215 240 285 265 {'3/20/2012'} 340 255 290 320 310

CDSDates = datenum (CDS.Date);CDSSpreads = table2array (CDS(:, 2:结束));(RateSpec ZeroData =。EndD一个te年代R一个teSpec.Rates];校准每个交易对手的违约概率defprobb = 0 (length(simulationDates)， size(CDSSpreads,2)); / /模拟日期为i = 1:size(defprobb,2) probData = cdsbootstrap(ZeroData， [CDSDates CDSSpreads(:，i)]，．..解决,“probDates”, simulationDates);DefProb(:,我)= probData (:, 2);结束%每个交易对手的累积违约概率图。图;情节(simulationDates DefProb)标题(“每个交易对手的违约概率曲线”）;包含(“日期”）;网格在；ylabel (“累积概率”) datetick (“x”，“mmmyy”) ylabel (违约概率的)包含(“模拟日期”）

CVA的计算

信贷价值(估值)调整(CVA)公式为:

$C V 一个＝（ 1 - R ） \int_{0}^{T} d 我年代 c E E （ t ） d P D （ t ）$

在哪里R经济复苏,discEE贴现的预期风险敞口t,PD违约概率分布。这假设敞口独立于违约(没有错误方向的风险)，并且还假设敞口是使用风险中性概率获得的。

在估值日期上用有限和近似积分如下:

$C V 一个（一个 p p r o x ）＝（ 1 - R ） \sum_{我＝ 2}^{n} d 我年代 c E E （ t_{我} ）（ P D （ t_{我} ） - P D （ t_{我 - 1} ））$

在哪里t_1是今天的日期,t_2、……t＿n未来估值日期。

假设CDS信息对应的是具有指数的交易对手cpIdx．计算出的CVA是我们对交易对手的信贷敞口的当前市场价值cpIdx．对于本例，将恢复速率设置为40%。

复苏= 0.4;CVA = (1-Recovery) * sum(discEE(2:end，:) .* diff(DefProb)); / /恢复为i = 1:numel(CVA) fprintf('交易对手的CVA %d = $%.2f\n',我,脑血管意外(i));结束

交易对手1的CVA = 2228.36美元，交易对手2的CVA = 2487.60美元，交易对手3的CVA = 920.39美元，交易对手4的CVA = 5478.50美元，交易对手5的CVA = 5859.30美元

图;酒吧(CVA);标题(“每个交易对手的CVA”）;包含(“对手”）;ylabel (“CVA的美元”）;网格在；

参考文献

Pykhtin, Michael和Steven Zhu，交易对手信用风险建模指南， GARP, 2007年7 / 8月，第37期，第16-22页。
Pykhtin, Michael和Dan Rosen，在交易水平和CVA上为交易对手风险定价, 2010年。
新巴塞尔协议:https://www.bis.org/publ/bcbs128.pdf256页。

另请参阅

cdsbootstrap|cdsprice|cdsrpv01|cdsspread

外部网站

信用违约互换的定价和估值(4分22秒)