The University of Sydney
School of Mathematics and Statistics
Computer Project
MATH2070/2970: Optimisation and Financial Mathematics Semester 2, 2019
Web Page: http://www.maths.usyd.edu.au/u/IM/MATH2070/
Lecturer: Anna Aksamit and Georg Gottwald
Due on Sunday 10th November at 11:59pm on TurnitIn
Submit exactly two files: a pdf with your report and m file with your Matlab code.
Report should be pleasant to read and include project formulations, descriptions and outputs
(tables, plots, histograms etc), all answers and discussion should be there.
Marking will be based on: accuracy, programming and presentation.
Please do not write your name on any sheet.
The deadline is a hard deadline in the sense that in case of a late submission (maximum
up to 10 days), you will be deducted
5% of the total marks for each day of delay. This
is non-negotiable, so make sure you submit in time; a submission on Monday the
11th at
12:01am will be an automatic deduction of
5%. It is your responsibility to check that your
submission was successful.
MATH2070: Do all questions except Question 6.
MATH2970: Do all questions.
In this computer project you will be analysing real stock market data downloaded from
Yahoo!Finance.
The file Data_2013_2019.csv which you can download from Ed, contains the daily closing prices of
the 30 stocks which make up the Dow Jones Industrial Average Index and closing prices of two indexes,
Dow Jones Industrial Average Index and S&P 500 Index. Prices are recorded on a (business)-daily
basis between 2/01/2013 and 30/09/2019.
There is one particularity with this time series: On 31st of Ausgust 2017 Dow and DuPont merged and
were traded as a new entity DowDuPont, then in 2019 Dow spun off of DowDuPont and was added to
the Dow Jones Industrial Average. Therefore only consider the 29 stocks without Dow (due to a short
trading history).
All prices are in US dollars.
Correlations and the covariance matrix
1. Export the data into Matlab using csvread and/or readtable. This question investigates the
correlations of the return rates of the 29 stocks. When analysing return rate data one has
several choices. A commonly used variable is the logarithmic change of price or the so called
log
return rate
: Let Skt be the price of k-th stock at time t, then consider Ykt = log Skt log Sk(t1)
(wrt the natural base).
(i) Calculate the maximal correlation between the
Yk, name and plot the two stock prices associated with the highest correlation as a function of time. On the graph present normalised
time series so that they start from the same value 100 on 2/01/2013.
(ii) Calculate the minimal correlation between the
Yk, name and plot the two stock prices associated with the smallest correlation as a function of time. On the graph present normalised
time series so that they start from the same value 100 on 2/01/2013.
Copyright c 2019 The University of Sydney 1
(iii) Visualise the correlation matrices for two subperiods: 1/12/2014–1/09/2016 and 1/09/2016–
1/02/2018 (you may use Matlab’s command
imagesc). Can you spot differences?
Plot the price of Dow Jones Industrial Average in the whole period. Can you relate it to
your observations about the correlation matrices?
(iv) Plot the histogram of the correlation coefficients
ρij for the two periods from the previous
point. Comment on your result.
Portfolio Theory
2. In this section we consider simple return rates, that is Rkt = SktSk(Stk(1) t1) , where Skt is the price
of the
k-th stock at time t. Carry out the following computational tasks for an unrestricted
optimal portfolio
P consisting of the 29 stocks included in the Dow Jones for an agent who
wants to invest $200,000 and has a risk aversion parameter
t = 0:2.
(a) Compute the dollar investment in each of the stocks and the corresponding expected return
and risk of
P .
(b) Illustrate the problem graphically and plot on the same graph in the
µσ-plane :
(i) The 29 stocks of the Dow Jones.
(ii) The minimum variance and efficient frontiers. Use a
t-range jtj ≤ 0:35 for your display.
(iii) A plot of
1000 random feasible portfolios satisfying jxij ≤ 20 (for each of the 29 stocks)
and
σi 0:05 for i = 1; : : : ; 1000.
You might notice that the random points occupy some region well-separated from the
minimum variance frontier (MVF) – comment on this and explain why (This is a/the
major part of the question).
(iv) The indifference curve of an investor with
t = 0:2 and their optimal portfolio P .
3. Determine which investors shortsell in the market consisting of the 29 stocks, and which stocks
they shortsell. Are there any stocks which no-one will shortsell or which everyone will shortsell?
4. Three funds with different risk profiles: In this question you will divide 29 stocks with
respect to their risk profile into 3 funds. Sort stocks from highest to lowest risk (expressed via
variance or standard deviation). Assuming that each stock has the same contribution to a given
fund, form high-risk fund from the 9 most risky stocks, low-risk fund from the 10 least risky
stocks and mid-risk fund from the rest.
(a) Compute expected returns and covariance matrix of the 3 funds.
(b) Let
P^ be an unrestricted optimal portfolio consisting of the 3 funds for an agent who wants
to invest $200,000 and has a risk aversion parameter
t = 0:2.
(i) Compute the dollar investment in each of the stocks and the corresponding expected
return and risk of
P^.
(ii) Plot on the second
µσ-plane graph :
The 3 funds.
The minimum variance and efficient frontiers. Use a t-range jtj ≤ 0:35 for your
display.
The indifference curve of an investor with t = 0:2 and their optimal portfolio P^.
The minimum variance frontier and optimal portfolio P from Question 2. Compare
solutions
P and P^ to the two problems based on computations and graphs.
2

Capital Asset Pricing Model
5. Assume that the daily risk free rate in the studied period was 0.002906%. Suppose that Standard
& Poor’s 500 Index is the market portfolio (S&P 500 Index prices are included in the data file).
Make a new
µσ-plane graph showing the risk free asset, market portfolio, and the Security
Market Line. Compute the
β’s of all relevant assets in this project (29 stocks, 3 funds, and
two optimal portfolio
Pand P^ from Questions 2 and 4). Plot these assets on the same graph.
Identify assets with
β’s greater than 1 and lower than 1. Comment on the result and describe
what Portfolio Theory would recommend an investor to do.
Log Optimal Portfolio
In this section we consider logarithmic utility maximisation problem. As in Markowitz Portfolio Theory,
there is a random vector of returns on
n stocks R = [R1; R2; :::; Rn]. The objective is to maximise the
expected logarithmic utility of the final wealth
W1, i.e.,
E log(W1) = E log(W0xT (1 + R)) = E log(xT (1 + R)) + log(W0)
where W0 is initial wealth. To ensure well-posedness of the problem there is no-shortselling constraint
imposed on feasible portfolio vector
x, that is x = [x1; x2; :::; xn] 2 Rn satisfies xi 0 for i = 1; :::; n
and Pn i=1 xi = 1.
Let us introduce
R¯ := 1 + R = [1 + R1; 1 + R2; :::; 1 + Rn], and denote the cumulative distribution
function of
R¯ by F, i.e., F(y) := P(R¯ y). We assume that R¯ is a discrete random variable, and
therefore
F is of the form F (y) = Pk:yky pk with pk = P(R¯ = yk) for each yk in a countable set
fy1; y2; :::g.
The problem of logarithmic utility maximisation can be written as
maximise
v(x; F ) := Eh log xT R¯ i
subject to x 0;
n X i
=1
xi = 1
Note that the above expectation is computed with respect to the distribution of R¯ which is uniquely
determined by
F.
For a fixed
F, a feasible portfolio xthat achieves the maximum of v(x; F ) is called log optimal
portfolio
, i.e.,

v(x; F ) = v(F) := max
x0; Pn i=1 xi=1
v(x; F ):

6. (a) Show that
(i)
v(x; F ) is concave in x and linear in F,
(ii)
v(F) is convex,
(iii) the set of log optimal portfolios with respect to a fixed
F is convex.
(b) Prove the following theorem:
Theorem: The log optimal portfolio xfor a fixed distribution F satisfies the following
necessary and sufficient conditions:
E (xR¯)Ti R¯  = 1 if xi > 0,
E (xR¯)Ti R¯  1 if xi = 0.
Hint: Using 6(a)(i) argue that log optimal portfolio can be characterised locally by an
appropriate condition on the directional derivative of
v in the direction from xto any
other feasible portfolio
x.
3