capbyblk

Price caps using Black option pricing model

Syntax

[CapPrice,Caplets] = capbyblk(RateSpec,Strike,Settle,Maturity,Volatility)

[CapPrice,Caplets] = capbyblk(___,Name,Value)

Description

[CapPrice,Caplets] = capbyblk(RateSpec,Strike,Settle,Maturity,Volatility) price caps using the Black option pricing model. capbyblk computes prices of vanilla caps and amortizing caps.

[CapPrice,Caplets] = capbyblk(___,Name,Value) adds optional name-value pair arguments.

Examples

Price a Cap Using the Black Option Pricing Model

Consider an investor who gets into a contract that caps the interest rate on a $100,000 loan at 8% quarterly compounded for 3 months, starting on January 1, 2009. Assuming that on January 1, 2008 the zero rate is 6.9394% continuously compounded and the volatility is 20%, use this data to compute the cap price. First, calculate the RateSpec:

ValuationDate = 'Jan-01-2008';
EndDates ='April-01-2010';
Rates = 0.069394;
Compounding = -1; 
Basis = 1;

RateSpec = intenvset('ValuationDate', ValuationDate, ...
'StartDates', ValuationDate,'EndDates', EndDates, ...
'Rates', Rates,'Compounding', Compounding,'Basis', Basis)

RateSpec = struct with fields:
           FinObj: 'RateSpec'
      Compounding: -1
             Disc: 0.8554
            Rates: 0.0694
         EndTimes: 2.2500
       StartTimes: 0
         EndDates: 734229
       StartDates: 733408
    ValuationDate: 733408
            Basis: 1
     EndMonthRule: 1

Compute the price of the cap.

Settle = 'Jan-01-2009'; % cap starts in a year
Maturity = 'April-01-2009';
Volatility = 0.20;
CapRate = 0.08;
CapReset = 4;
Principal=100000;

CapPrice = capbyblk(RateSpec,  CapRate, Settle, Maturity, Volatility,...
'Reset',CapReset,'ValuationDate',ValuationDate,'Principal', Principal,...
'Basis', Basis)

CapPrice = 51.6125

Price a Cap Using a Different Curve to Generate the Future Forward Rates

Define the OIS and Libor rates.

Settle = datenum('15-Mar-2013');
CurveDates = daysadd(Settle,360*[1/12 2/12 3/12 6/12 1 2 3 4 5 7 10],1);
OISRates = [.0018 .0019 .0021 .0023 .0031 .006  .011 .017 .021 .026 .03]';
LiborRates = [.0045 .0047 .005 .0055 .0075 .0109  .0162 .0216 .0262 .0309 .0348]';

Create an associated RateSpec for the OIS and Libor curves.

OISCurve = intenvset('Rates',OISRates,'StartDate',Settle,'EndDates',CurveDates,'Compounding',2,'Basis',1);
LiborCurve = intenvset('Rates',LiborRates,'StartDate',Settle,'EndDates',CurveDates,'Compounding',2,'Basis',1);

Define the Cap instruments.

Maturity = {'15-Mar-2018';'15-Mar-2020'};
Strike = [0.04;0.05];
BlackVol = 0.2;

Price the cap instruments using the term structure OISCurve both for discounting the cash flows and generating future forward rates.

[Price, Caplets] = capbyblk(OISCurve, Strike, Settle, Maturity, BlackVol)

Price = 2×1

    0.7472
    0.9890

Caplets = 2×7

         0    0.0000    0.0033    0.2996    0.4443       NaN       NaN
         0    0.0000    0.0003    0.1134    0.2112    0.2292    0.4349

Price the cap instruments using the term structure LiborCurve to generate future forward rates. The term structure OISCurve is used for discounting the cash flows.

[PriceLC, CapletsLC] = capbyblk(OISCurve, Strike, Settle, Maturity, BlackVol,'ProjectionCurve',LiborCurve)

PriceLC = 2×1

    1.3293
    1.6329

CapletsLC = 2×7

         0    0.0000    0.0337    0.4250    0.8706       NaN       NaN
         0    0.0000    0.0052    0.1767    0.4849    0.3663    0.5998

Compute the Price of Two Amortizing Caps Using the Black Model

Define the RateSpec.

Rates = [0.0358; 0.0421; 0.0473; 0.0527; 0.0543];
ValuationDate = '15-Nov-2011';
StartDates = ValuationDate;
EndDates = {'15-Nov-2012';'15-Nov-2013';'15-Nov-2014' ;'15-Nov-2015';'15-Nov-2016'};
Compounding = 1;
RateSpec = intenvset('ValuationDate', ValuationDate,'StartDates', StartDates,...
             'EndDates', EndDates,'Rates', Rates, 'Compounding', Compounding)

RateSpec = struct with fields:
           FinObj: 'RateSpec'
      Compounding: 1
             Disc: [5x1 double]
            Rates: [5x1 double]
         EndTimes: [5x1 double]
       StartTimes: [5x1 double]
         EndDates: [5x1 double]
       StartDates: 734822
    ValuationDate: 734822
            Basis: 0
     EndMonthRule: 1

Define the cap instruments.

Settle ='15-Nov-2011';
Maturity = '15-Nov-2015';
Strike = [0.03;0.035];
Reset = 1;
Principal ={{'15-Nov-2012' 100;'15-Nov-2013' 70;'15-Nov-2014' 40;'15-Nov-2015' 10}};

Price the amortizing caps.

Volatility = 0.10;  
Price = capbyblk (RateSpec, Strike, Settle, Maturity, Volatility,...
'Reset', Reset,'Principal', Principal)

Price = 2×1

    3.0339
    2.0141

Price a Cap Using the Shifted Black Model

Create the RateSpec.

ValuationDate = 'Mar-01-2016';
EndDates = {'Mar-01-2017';'Mar-01-2018';'Mar-01-2019';'Mar-01-2020';'Mar-01-2021'};
Rates = [-0.21; -0.12; 0.01; 0.10; 0.20]/100;
Compounding = 1;
Basis = 1;

RateSpec = intenvset('ValuationDate',ValuationDate,'StartDates',ValuationDate, ...
'EndDates',EndDates,'Rates',Rates,'Compounding',Compounding,'Basis',Basis)

RateSpec = struct with fields:
           FinObj: 'RateSpec'
      Compounding: 1
             Disc: [5x1 double]
            Rates: [5x1 double]
         EndTimes: [5x1 double]
       StartTimes: [5x1 double]
         EndDates: [5x1 double]
       StartDates: 736390
    ValuationDate: 736390
            Basis: 1
     EndMonthRule: 1

Price the cap with a negative strike using the Shifted Black model.

Settle = 'Jun-01-2016'; % Cap starts in 3 months.
Maturity = 'Sep-01-2016';
ShiftedBlackVolatility = 0.31;
CapRate = -0.003;  % -0.3 percent strike.
CapReset = 4;
Principal = 100000;
Shift = 0.01; % 1 percent shift.

CapPrice = capbyblk(RateSpec,CapRate,Settle,Maturity,ShiftedBlackVolatility,...
'Reset',CapReset,'ValuationDate',ValuationDate,'Principal',Principal,...
'Basis',Basis,'Shift',Shift)

CapPrice = 26.0733

Input Arguments

RateSpec — Interest-rate term structure

structure

Interest-rate term structure (annualized and continuously compounded), specified by the RateSpec obtained from intenvset. For information on the interest-rate specification, see intenvset.

Data Types: struct

Strike — Rate at which cap is exercised

decimal

Rate at which cap is exercised, specified as a NINST-by-1 vector of decimal values.

Data Types: double

Settle — Settlement date for cap

serial date number | date character vector

Settlement date for the cap, specified as a serial date number or a date character vector.

Data Types: double | char

Maturity — Maturity date for cap

serial date number | date character vector | cell array of date character vectors

Maturity date for the cap, specified as a serial date number or date character vector.

Data Types: double | char

Volatility — Volatilities values

numeric

Volatilities values, specified as a NINST-by-1 vector of numeric values.

The Volatility input is not intended for volatility surfaces or cubes. If you specify a matrix for the Volatility input, capbyblk internally converts it into a vector. capbyblk assumes that the volatilities specified in the Volatility input are flat volatilities, which are applied equally to each of the caplets.

Data Types: double

Name-Value Pair Arguments

Specify optional comma-separated pairs of Name,Value arguments. Name is the argument name and Value is the corresponding value. Name must appear inside quotes. You can specify several name and value pair arguments in any order as Name1,Value1,...,NameN,ValueN.

Example: [CapPrice,Caplets] = capbyblk(RateSpec,Strike,Settle,Maturity,Volatility,'Reset',CapReset,'Principal',100000,'Basis',7)

'Reset' — Reset frequency payment per year

1 (default) | numeric

Reset frequency payment per year, specified as a NINST-by-1 vector.

Data Types: double

'Principal' — Notional principal amount

100 (default) | numeric

Notional principal amount, specified as a NINST-by-1 vector or a NINST-by-1 cell array. When Principal is a NINST-by-1 cell array, each element is a NumDates-by-2 cell array, where the first column is dates and the second column is associated principal amount. The date indicates the last day that the principal value is valid.

Use Principal to pass a schedule to compute the price for an amortizing cap.

Data Types: double | cell

'Basis' — Day-count basis of instrument

0 (actual/actual) (default) | integer from 0 to 13

Day-count basis representing the basis used when annualizing the input forward rate, specified as a NINST-by-1 vector of integers.

• 0 = actual/actual

• 1 = 30/360 (SIA)

• 2 = actual/360

• 3 = actual/365

• 4 = 30/360 (PSA)

• 5 = 30/360 (ISDA)

• 6 = 30/360 (European)

• 7 = actual/365 (Japanese)

• 8 = actual/actual (ICMA)

• 9 = actual/360 (ICMA)

• 10 = actual/365 (ICMA)

• 11 = 30/360E (ICMA)

• 12 = actual/365 (ISDA)

• 13 = BUS/252

For more information, see basis.

Data Types: double

'ProjectionCurve' — Rate curve used in generating future forward rates

if ProjectionCurve is not specified, then RateSpec is used both for discounting cash flows and projecting future forward rates (default) | structure

The rate curve to be used in generating the future forward rates. This structure must be created using intenvset. Use this optional input if the forward curve is different from the discount curve.

Data Types: struct

'Shift' — Shift in decimals for shifted Black model

0 (no shift) (default) | positive decimal

Shift in decimals for the shifted Black model, specified using a scalar or NINST-by-1 vector of rate shifts in positive decimals. Set this parameter to a positive rate shift in decimals to add a positive shift to the forward rate and strike, which effectively sets a negative lower bound for the forward rate. For example, a Shift of 0.01 is equal to a 1% shift.

Data Types: double

Output Arguments

CapPrice — Expected price of cap

vector

Expected price of the cap, returned as a NINST-by-1 vector.

Caplets — Caplets

array

Caplets, returned as a NINST-by-NCF array of caplets, padded with NaNs.

More About

Shifted Black

The Shifted Black model is essentially the same as the Black’s model, except that it models the movements of (F + Shift) as the underlying asset, instead of F (which is the forward rate in the case of caplets).

This model allows negative rates, with a fixed negative lower bound defined by the amount of shift; that is, the zero lower bound of Black’s model has been shifted.

Algorithms

Black Model

$$\begin{array}{l}dF={\sigma }_{Black}Fdw\\ call=e^{-\gamma T}\left[FN(d_1)-KN(d_2)\right]\\ put=e^{-\gamma T}\left[KN(-d_2)-FN(-d_1)\right]\\ d_1=\frac{\ln \left(\frac{F}{K}\right)+\left(\frac{{\sigma }_B{}^2}{2}\right)T}{{\sigma }_B\sqrt{T}},d_2=d_1-{\sigma }_B\sqrt{T}\\ {\sigma }_B={\sigma }_{Black}\end{array}$$

Where F is the forward value and K is the strike.

Shifted Black Model

$$\begin{array}{l}dF={\sigma }_{Shifted\_Black}\left(F+Shift\right)dw\\ call=e^{-\gamma T}\left[\left(F+Shift\right)N(d_{s1})-\left(K+Shift\right)N(d_{s2})\right]\\ put=e^{-\gamma T}\left[\left(K+Shift\right)N(-d_{s2})-\left(F+Shift\right)N(-d_{s1})\right]\\ d_{s1}=\frac{\ln \left(\frac{F+Shift}{K+Shift}\right)+\left(\frac{{\sigma }_{sB}{}^2}{2}\right)T}{{\sigma }_{sB}\sqrt{T}},d_{s2}=d_{s1}-{\sigma }_{sB}\sqrt{T}\\ {\sigma }_{sB}={\sigma }_{Shifted\_Black}\end{array}$$

Where F+Shift is the forward value and K+Shift is the strike for the shifted version.