Compute the term structure of interest rates

Assignment Help MATLAB Programming
Reference no: EM131003897

Execution of the home task on Term Structure (Pricing, Vasicek Model, Monte Carlo Simulation, Bond Option,) with all the results produced in Matlab program including the doc report on the project Vasicek model

1) Vasicek model. Compute the price of a zero-coupon bonds with the Vasicek model. Use the function vbond.m with the following initial values:
• ro - initial short rate = 0.05
• k - speed of mean reversion = 0.025
• theta - long run risk-neutral mean short rate = 0.15
• sigma - annual volatility of short rate = 0.012
• t - time to evaluate price at = 0
• T - maturity of bond = 1, 2, 3„ 30 years

function P=vbond(ro,k,theta,sigma,t,T)
B=(1-exp(-k.*(T-t)))/k; A=exp((theta-sigma*sigma/2/k/k).*(B-(T-t))-sigma*sigma.*B.*B/4/k); V]=vrate(ro, k, theta, sigma, t);
P=A.*exp(-B.*M +.5.*B.*B.*V);

Exercise 1 Plot the term structure of interest rates - remember that the y(t, T) = -1/T-t log P(t, T). Describe the term structure.

Exercise 2 Compute the term structure of interest rates with the value of ro= 0.17. What is the difference between the two term structures? Why?

Exercise 3 Compare the 10-year bond price obtained with the exact formula and that with Monte Carlo. Do the same for the 5-year bond.

Exercise 4 Price with 200 Monte Carlo simulations the price of a 10-year bond with the CIR model, i.e.

drt = k(0 - rt)dt + σ√rtdWt

Remember that for the CIR we have

E[rt] = roe-kt + θ(1 - e-kt)

Var[rt] = rooσ2/k (e-kt e-2kt) + θσ2/2k(1 e-kt)2

and

rt˜c.Χ2 (4kθ/σ2, rte-kt/c)

where c = (1-ekT)σ2/4k, and Χ2 (.) is a non-central chi-Squared distribution with 4kθ/σ2 degrees of freedom and non-centrality parameter rte-kt/c. Use the Matlab function ncx2rnd(4kθ/σ2, rte-kt/c).

Exercise 5 ADDITIONAL POINTS. Simulate the distribution. Compare the exact price

P(t, T) = A(t, T) exp(- B (t,T) . rt)

where

A(t,T) = ((2h exp(k + h)(T - t)/2)/(2h + (k + h)(exp(T - t)h) - 1)2kθ/σ2

B(t,T) = (2 (exp((T - t)h) - 1))/ (2h+ (k + h)(exp(T - t)h) - 1)

h = √k2 +2σ2

you can also use the Matlab built-in functions, with the Monte Carlo price ob¬tained in the previous exercise.

Exercise 6 Consider a call option on a 10-year bond, with strike equal to the 5-year bond price (see exercise above). The option expires in 5 years. Use a Monte Carlo simulation to price this call option and compare it with the exact pricing (use the function vbondop . m) .

Exercise 7 Do the same exercise with strikes = 0.5, 0.7, 0.9 (i.e. in-the-money and out-of-the money). See what happens with different strikes.

Reference no: EM131003897

Questions Cloud

Explain the relationship among the variables : Do the following based on your scenario, - List and label the variables as independent, dependent, intervening, or moderating,  Explain the relationship among the variables, - Develop a problem statement for the situation.
Dataline case study : James Dougan, chief executive of Dataline Oil Services (DOS), faced a tough decision. A bid, which he and his fellow managers had put together to buy DOS from its parent company, had just been topped unexpectedly by another company. He and his team n..
Determine the price you are willing to pay : Write a brief paper (two pages maximum) supporting your rationale for the valuation you have chosen, and provide appropriate financial appendices (in addition to the two pages of text) showing your analyses.
From a prospective investor point of view : From a prospective investor's (or buyer's) point of view, determine the price you are willing to pay, if you are to invest. Write a brief paper (two pages maximum) supporting your rationale for the valuation you have chosen, and provide appropriate f..
Compute the term structure of interest rates : Plot the term structure of interest rates - remember that the y(t, T) = -1/T-t log P(t, T). Describe the term structure and Compute the term structure of interest rates with the value of ro= 0.17. What is the difference between the two term structur..
What are the x-intercepts of the graph : What are the x-intercepts of the graph? At what values of x does the graph have a vertical asymptote? At what values of y does the graph have a horizontal asymptote?
Autumn 1999 exam : A survey of a magazine's subscribers indicates that 60% own a home and 75% own a car. Ninety percent of the home owners who subscribe to the magazine, also own a car. What proportion of subscribers i. own both a car and a house? ii. own a car or a..
What is the index of refraction of this material : A laser beam shines along the surface of a block of transparent material (see figure below). Half of the beam goes straight to a detector, while the other half travels through the block and then hits the detector. The time delay between the arriva..
Calculate the second derivative of given function : Calculate the second derivative of f find where f is concave up concave down and has inflection points. Find any horizontal and vertical asymptotes of f.

Reviews

Write a Review

 

MATLAB Programming Questions & Answers

  Division by zero is not possible

You must also ensure that division by zero is not possible. If the user tries to divide by zero, the calculation should not be attempted, and an error message should appear in the Result text box.

  Find the answer through binomial distribution

The probability to get 3 times of "5" in 8 throws of a fair die is ? Using Binomial Distribution, and Using Matlab to simulate it.

  Finds the largest integer in the list recursively

Write a program that, given a list of 20 integers, finds the largest integer in the list recursively

  Computing project - beam analysis

Computing Project - Beam Analysis- The computing project Beam Analysis deals with load placement at various points along a beam. The concept is to model the relationship between location of load placement and the reaction forces at the end of the ..

  Is the interest rate expressed as a fractional amount

Create a MATLAB® function called future_value with three inputs: the investment (present value), the interest rate expressed as a fraction, and the number of compounding periods.

  Problems related to matlab

solve the system for the least-squares estimates of Δρi.

  Open a named pipe and to read data from the pipe

Open a named pipe and to read data from the pipe in matlab

  Impact of zigbee on telecommunication systems

Report on Impact of zigbee on telecommunication systems

  Generate and plot a continuous time signal

1.Generate and plot a continuous-time signal and a train of unit impulses using MATLAB 2. Generate samples of a given signal using a train of unit impulses using  MATLAB 3. Investigate the effect of sampling analogue signals

  Random processes

Theme: Random Processes, Compute and plot separate correlation functions of white noise data and of sinusoidal data. Do this using both the biased and unbiased autocorrelation functions. Are the results what you expected?

  Beginning with the row number passed

Write a MATLAB function that takes a matrix, a row number and a column number. Beginning with the row number passed to the function, scan down the column passed to the function and return the row number that contains the largest absolute value in ..

  Determine the level of absolute approximate relative error

2006ENG Numerical and Data Analysis. Convert the equation into f(T) = 0 when cp of 1.1kJ/kgK. In our MATLAB program, develop a plot of f(T) versus a range of T = 0 to 1200 K. MATLAB. On the axes, put proper labelling.

Free Assignment Quote

Assured A++ Grade

Get guaranteed satisfaction & time on delivery in every assignment order you paid with us! We ensure premium quality solution document along with free turntin report!

All rights reserved! Copyrights ©2019-2020 ExpertsMind IT Educational Pvt Ltd