Matlab通過市場數(shù)據(jù)校準(zhǔn)Hull-White利率模型參數(shù)
原文鏈接:http://tecdat.cn/?p=17973
原文出處:拓端數(shù)據(jù)部落公眾號
?
利率衍生證券的定價(jià)依賴于描述基本過程的模型。這些利率模型取決于您必須通過將模型預(yù)測與市場上可用的現(xiàn)有數(shù)據(jù)進(jìn)行匹配來確定的一個(gè)或多個(gè)參數(shù)。在Hull-White模型中,有兩個(gè)與短期利率過程相關(guān)的參數(shù):均值回歸和波動(dòng)率。
對于Hull-White模型,關(guān)于均值回歸(α)和波動(dòng)率(σ)最小化是二維的。也就是說,校準(zhǔn)Hull-White模型可最大程度地減少模型的預(yù)測價(jià)格與觀察到的市場價(jià)格之間的差異。
Hull-White校準(zhǔn)案例
使用市場數(shù)據(jù)來識別為構(gòu)建工具定價(jià)的Hull-White樹所需的隱含波動(dòng)率(σ)和均值回歸(α)系數(shù),必須對市場數(shù)據(jù)進(jìn)行插值。
考慮具有以下參數(shù)的上限:
Strike = 0.0690;
Reset = 4;
Principal = 1000;
?
datestr(capletDates')
ans =
21-Mar-2008
21-Jun-2008
21-Sep-2008
21-Dec-2008
21-Mar-2009
21-Jun-2009
21-Sep-2009
21-Dec-2009
21-Mar-2010
21-Jun-2010
21-Sep-2010
21-Dec-2010
21-Mar-2011
在最佳情況下,查找?guī)в?code>Strike?=的 市場波動(dòng)率?0.0690
,以及列出的到期日,但找到這些確切數(shù)據(jù)的可能性很小。結(jié)果,使用市場上可用的數(shù)據(jù)并進(jìn)行插值。
根據(jù)市場數(shù)據(jù),您可以獲得不同日期和行使價(jià)的上限信息。
?
到期協(xié)議1 = 0.0590協(xié)議2 = 0.07902008年3月21日0.15330. 15262008年6月21日0.17310. 17302008年9月21日0. 17270. 17262008年12月21日0. 17520. 17472009年3月21日0. 18090. 18082009年6月21日0. 18090. 17922009年9月21日0. 18050. 17972009年12月21日0.18020. 17942010年3月21日0.18020. 17332010年6月21日0. 17570. 17512010年9月21日0. 17550. 17502010年12月21日0. 17550. 17452011年3月21日0. 17260. 1719
?要求日期,協(xié)議價(jià)和實(shí)際波動(dòng)分為三個(gè)變量:??MarketStrike
,??MarketMat
,和?MarketVol
。
MarketStrike = [0.0590; 0.0790];
MarketMat = {'21-Mar-2008';
'21-Jun-2008';
'21-Sep-2008';
MarketVol = [0.1533 0.1731 0.1727 0.1752 0.1809 0.1800 0.1805 0.1802 0.1735 0.1757 ...
使用此數(shù)據(jù)完成輸入?yún)?shù)?RateSpec
:
RateSpec =
FinObj: 'RateSpec'
Compounding: 4
Disc: [14x1 double]
Rates: [14x1 double]
EndTimes: [14x1 double]
StartTimes: [14x1 double]
EndDates: [14x1 double]
StartDates: 733428
ValuationDate: 733428
Basis: 0
EndMonthRule: 1
調(diào)用校準(zhǔn)例程以找到波動(dòng)率參數(shù)Alpha和Sigma的值
根據(jù)市場數(shù)據(jù)計(jì)算的值?Alpha
?,Sigma
?。例如,您可以將目標(biāo)函數(shù)的公差定義為?100*eps
?:
Local minimum possible.
> In hwcalbycapfloor at 93
In hwcalbycap at 75
Alpha =
1.0000e-06
Sigma =
0.0127
先前的警告表明轉(zhuǎn)換不是最佳的。使用的搜索算法未找到符合所有約束的解決方案。要辨別該解決方案是否可接受,請通過為以下項(xiàng)指定第三個(gè)輸出(OptimOut
)來查看?優(yōu)化結(jié)果:
hwcalb(RateSpec, MarketStrike, MarketMat,...
?結(jié)構(gòu)的?OptimOut.residual
?字段?OptimOut
是優(yōu)化殘差。此值包含Black caplet和優(yōu)化過程中計(jì)算出的值之間的差異。您可以使用該?OptimOut.residual
?值計(jì)算與Black Caplet價(jià)格相比的百分比差異(誤差),然后確定殘差是否可以接受,決定是否接受以單一數(shù)值參數(shù)化的市場?Alpha
?和?Sigma
。
使用市場數(shù)據(jù)和Black模型獲取價(jià)格
Black -76模型被認(rèn)為是世界范圍內(nèi)廣泛接受的互換期權(quán)定價(jià)模型,互換期權(quán)計(jì)算公式看起來和期權(quán)價(jià)格公式一 樣,這表示在使用Black模型對互換期權(quán)進(jìn)行定價(jià)時(shí),可以用和定價(jià)期權(quán)時(shí)同樣的方法進(jìn)行計(jì)算.為了使Black模型生效,只需輸入由市場決定的隱含波動(dòng)率即可。
為了確定優(yōu)化的有效性,請使用Black公式和市場數(shù)據(jù)計(jì)算參考值。請注意,您必須首先對市場數(shù)據(jù)進(jìn)行插值以獲得用于計(jì)算的數(shù)據(jù)量:
FlatVol = interp2(Mats, Strikes, MarketVol, datenum(Maturity),
使用Black模型計(jì)算價(jià)格:
capbybike, Settle, Maturity, FlatVol,...
Caplets =
0.3210
1.6355
比較優(yōu)化值和Black值并以圖形方式顯示
在計(jì)算了參考值之后,請通過分析和圖形方式比較這些值,以確定是否計(jì)算出的單個(gè)值?Alpha
?并?Sigma
?提供適當(dāng)?shù)慕浦担?/p>
plot(MarketMatNum(2:end), Caplets, 'or', MarketMatNum(2:end), O
ylabel('Caplet Price');
t
set(gcf, 'NumberT

使用Black,HW分析和HW樹模型比較上限價(jià)格

最受歡迎的見解
1.用機(jī)器學(xué)習(xí)識別不斷變化的股市狀況—隱馬爾科夫模型(HMM)的應(yīng)用
2.R語言GARCH-DCC模型和DCC(MVT)建模估計(jì)
3.R語言實(shí)現(xiàn) Copula 算法建模依賴性案例分析報(bào)告
4.R語言COPULAS和金融時(shí)間序列數(shù)據(jù)VaR分析
5.R語言多元COPULA GARCH 模型時(shí)間序列預(yù)測
6.用R語言實(shí)現(xiàn)神經(jīng)網(wǎng)絡(luò)預(yù)測股票實(shí)例
7.r語言預(yù)測波動(dòng)率的實(shí)現(xiàn):ARCH模型與HAR-RV模型
8.R語言如何做馬爾科夫轉(zhuǎn)換模型markov switching model
9.matlab使用Copula仿真優(yōu)化市場風(fēng)險(xiǎn)
使用Black模型,Hull-White分析模型和Hull-White樹模型比較相應(yīng)上限價(jià)格。