<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JMF</journal-id><journal-title-group><journal-title>Journal of Mathematical Finance</journal-title></journal-title-group><issn pub-type="epub">2162-2434</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jmf.2019.94034</article-id><article-id pub-id-type="publisher-id">JMF-96081</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Business&amp;Economics</subject><subject> Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Portfolio Mathematics with General Linear and Quadratic Constraints
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>David</surname><given-names>L. Stowe</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Department of Finance, Ohio University, Athens, Ohio, USA</addr-line></aff><pub-date pub-type="epub"><day>17</day><month>10</month><year>2019</year></pub-date><volume>09</volume><issue>04</issue><fpage>675</fpage><lpage>690</lpage><history><date date-type="received"><day>12,</day>	<month>July</month>	<year>2019</year></date><date date-type="rev-recd"><day>27,</day>	<month>October</month>	<year>2019</year>	</date><date date-type="accepted"><day>30,</day>	<month>October</month>	<year>2019</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  This paper explores the mathematics behind optimal portfolio construction when relative utility and risk are considered together in a general sense. I derive the portfolio optimization problems when subject to both a general liner constraint and a constraint to tracking error (a quadratic constraint), the most pervasive constraint placed on delegated portfolio managers. This unifies three very influential papers from the evolution of optimal portfolio theory. In addition, I also analyze the general linear constraint when applied to Sharpe Ratio maximization. When applied together, these formulations can allow principals and agents to better analyze alternatives and negotiate contracting in order to ensure that the constraints generate proper utility maximization.
 
</p></abstract><kwd-group><kwd>Portfolio Mathematics</kwd><kwd> Optimization</kwd><kwd> Delegated Portfolio Management</kwd><kwd> Tracking Error</kwd><kwd> Utility Maximization</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The purpose of this essay is to generalize the Markowitz portfolio choice problem and derive solutions to the variance minimization and utility maximization problems considering both general linear constraints and a single quadratic constraint. The Markowitz portfolio choice problem (introduced in [<xref ref-type="bibr" rid="scirp.96081-ref1">1</xref>] ) has been studied from many different perspectives. The problem is essential to minimize the variance of a linearly combined set of correlated random variables (returns), conditioned on the mean of these variables, with a constraint on the scalers or weights. The set of all such points is called the efficient frontier, i.e. the minimum variance envelope. The modern expression of the optimization problem outlined by Markowitz is due to Merton [<xref ref-type="bibr" rid="scirp.96081-ref2">2</xref>] who derives the solutions analytically using linear algebra. His solution considers the problem as a quadratic programming problem with two predefined linear constraints, the weight constraint (scalar constraint) and the expected return constraint (the mean). The statement of the problem is as follows. Matrices and vectors are in boldface.</p><p>w - n &#215; 1 weight vector of the portfolio.</p><p>Ω - n &#215; n covariance matrix (symmetric and positive definite).</p><p>r - n &#215; 1 vector of expected returns.</p><p>1 - n &#215; 1 vector of 1’s.</p><p>min w w ′ Ω w     s .t .     w ′ 1 = 1     and     w ′ r = r . (1.1)</p><p>However, the problem can also be solved with a general linear constraint as demonstrated by Best and Grauer [<xref ref-type="bibr" rid="scirp.96081-ref3">3</xref>] . The general linear constraint opens up possibilities to make considerations for and constraints on portfolio beta, sector or industry weights, multi-stage optimization, etc. I derive the solution to the general linear constraint problem in part 2 of this essay to familiarize the reader with my notation and so I can use it in my derivation of the maximum Sharpe Ratio problem.</p><p>However, since the original introduction of the model, modern thinking on portfolio construction (and financial theory in general) has evolved to consider the differences in incentives between the individual making the portfolio choices, the agent, and the individual hiring the delegated manager, the principal. This spawned a theory of active portfolio management, or delegated portfolio management, which centers on a time series of excess returns, essentially the returns of a security or portfolio in excess of a benchmark. This is a space transformation that causes agents to optimize in excess return space rather than raw return space (or even another excess return space with a different benchmark). Principals are then potentially motivated to constrain agent behavior in this space by putting limits on the variance of the delegated manager’s excess return, i.e. their tracking error variance.</p><p>There are many modern examples of papers analyzing delegated or relative problems. Bisburgen, Brandt, and Koijen [<xref ref-type="bibr" rid="scirp.96081-ref4">4</xref>] revisit the problem of multiple external agents relaxing the assumption about the certainty of the agents’ risk appetites and conclude that serious inefficiency exists without this assumption. Blake et al. [<xref ref-type="bibr" rid="scirp.96081-ref5">5</xref>] use the BBK framework and apply it to the delegated relationships in the pension industry and, among many other things, show that the delegated incentive is pervasive in professional investment management. Cuoco and Kaniel [<xref ref-type="bibr" rid="scirp.96081-ref6">6</xref>] analyze the linkage between compensation contracts and excess performance in delegated portfolio management. However, the delegated contracting problem can probably best be summed up by analyzing the tracking error constraint.</p><p>The first to analyze a constraint to tracking error in the context of a Markowitz portfolio choice model was Roll [<xref ref-type="bibr" rid="scirp.96081-ref7">7</xref>] . In this paper, Roll demonstrates the differences in portfolio choice for a decision-maker optimizing over variance versus one optimizing over tracking error. In a further refinement, Jorion [<xref ref-type="bibr" rid="scirp.96081-ref8">8</xref>] , applies a tracking error constraint to a mean/variance optimization problem, thus introducing a quadratic constraint into the problem. Whereas the Markowitz problem, under linear constraints, forms a minimum variance set shaped like a parabola in mean and variance, the tracking error constraint causes the efficient set to be shaped like an ellipse, with the feasible portfolios inside. See <xref ref-type="fig" rid="fig1">Figure 1</xref> for a visual example of these curves.</p><p>The closest application of a utility problem directly to Roll’s TEV frontier by Bertrand [<xref ref-type="bibr" rid="scirp.96081-ref9">9</xref>] . In this paper, he considers the problem of a fixed risk aversion constraint in mean/variance and its ability to generate preferred portfolios.</p><p>The variance minimization problem is one approach to analyzing portfolio choice. However, it is also reasonable to consider the problem of maximizing utility. The popular utility function considered in this problem is a “quadratic-style” utility where the investor’s utility increases in the expected return, w ′ r , and decreases in the expected risk or variance, w ′ Ω w . And, the parameter, a, measures the sensitivity or the tradeoff: the level of risk aversion. The relationship is expressed as follows:</p><p>max w w ′ r − a w ′ Ω w (1.2)</p><p><xref ref-type="fig" rid="fig2">Figure 2</xref> shows a few of the curves depicted in <xref ref-type="fig" rid="fig1">Figure 1</xref> with iso-utility curves corresponding to the same level of risk aversion, a.</p><p>Modern extensions of the Markowitz portfolio choice problem tend to focus more on altering something fundamental about the model or analyzing a specific constraint. For example, Chen et al. [<xref ref-type="bibr" rid="scirp.96081-ref10">10</xref>] optimize over semivariance but the model is adjusted to account for uncertainty in variable estimation based on historical data. Or authors may use modern techniques to analyze traditional problems like Kar et al. [<xref ref-type="bibr" rid="scirp.96081-ref11">11</xref>] that optimizes over the Sharpe Ratio and the Value at Risk simultaneously. These papers are very interesting and practical considerations for real world applications of this theory.</p><p>The primary contribution of this essay is in solving the problem proposed by Jorion [<xref ref-type="bibr" rid="scirp.96081-ref8">8</xref>] , optimizing with a tracking error constraint, under the conditions of the general linear constraint of Best and Grauer [<xref ref-type="bibr" rid="scirp.96081-ref3">3</xref>] . This expands and unifies two important developments in the mathematical analysis of the original Markowitz portfolio optimization model. I solve this problem with both a variance optimization approach and a utility maximization approach. The remainder of this paper precedes as follows. In section 2, I revisit the Best and Grauer [<xref ref-type="bibr" rid="scirp.96081-ref3">3</xref>] analysis with only a general linear constraint. In section 3, I derive the solution to the problem with the addition of a quadratic constraint (the tracking error constraint). This section contains this paper’s primary contribution and is a significant development in the ability to apply the general linear constraint to analyze modern delegated contracting with constraints and considerations to relative risk. In section 4, I consider an alternative utility function, the maximization of the Sharpe Ratio, Sharpe [<xref ref-type="bibr" rid="scirp.96081-ref12">12</xref>] , given these constraints. Lastly, section 5 concludes and suggests potential paths of future research.</p></sec><sec id="s2"><title>2. Mean/Variance Optimization with a General Linear Constraint</title><p>The general linear constraint problem analyzed in Best and Grauer [<xref ref-type="bibr" rid="scirp.96081-ref3">3</xref>] comes in two forms, variance minimization given linear constraints and utility maximization given linear constraints. Although these problems are already solved, I repeat this calculation to familiarize the reader with my notation and so I can refer to my solutions throughout the paper. Additionally, I leave the expressions in the purest form possible, without making extraneous substitution. As before, the notation used throughout this paper is as follows:</p><p>w - n &#215; 1 weight vector of the portfolio.</p><p>Ω - n &#215; n covariance matrix (symmetric and positive definite).</p><p>E - n &#215; k matrix of constraint conditions.</p><p>c - 1 &#215; k vector of constraint constants.</p><sec id="s2_1"><title>2.1. Minimize Variance with a General Linear Constraint</title><p>The variance function with a weight vector and a covariance matrix is w ′ Ω w . We want to minimize this function given a general linear constraint on the weights as follows:</p><p>min w w ′ Ω w     s .t .     w ′ E = c .</p><p>We introduce a new vector, Λ , to use as Lagrange multipliers:</p><p>Λ - k &#215; 1 vector of Lagrange multipliers.</p><p>The Lagrangian is:</p><p>L = w ′ Ω w − ( w ′ E − c ) Λ .</p><p>We differentiate the Lagrangian and set it equal to zero to find the critical values:</p><p>2 Ω w − E Λ = 0 .</p><p>Solving for the weight vector yields:</p><p>w = 1 2 Ω − 1 E Λ .</p><p>Next, we plug this weight vector into the condition to solve for the Lagrange multipliers:</p><p>w ′ E = 1 2 Λ ′ E ′ Ω − 1 E = c .</p><p>Λ = 2 ( E ′ Ω − 1 E ) − 1 c ′ .</p><p>Plugging this value into the expression for the weight vector above reveals the optimal weight vector given the general linear constraint:</p><p>w = Ω − 1 E ( E ′ Ω − 1 E ) − 1 c ′ . (2.1.1)</p><p>This is a beautiful expression without substitution. We can apply this to the variance function to calculate the variance given the constraints:</p><p>w ′ Ω w = c ( E ′ Ω − 1 E ) − 1 c ′ .</p><p>This is an equally beautiful expression. To solve this with the Markowitz constraints from Equation (1.1), we would make:</p><p>E = [ 1 r ] ,</p><p>and</p><p>c = [ 1 r ] .</p><p>E is essentially an augmented matrix with a column of 1s and a column of expected returns for each security. c is a 2 &#215; 1 matrix where the constraint on the top represents the sum of all the weights equaling 1 and the bottom r representing the sum of all the weights times the expected returns equaling a target weight, the weighted average expected return. The expression of the Markowitz minimum variance envelope with these constraints is:</p><p>w ′ Ω w = [ 1 r ] ′ [ [ 1 r ] ′ Ω − 1 [ 1 r ] ] − 1 [ 1 r ] .</p><p>And, we can vary the constraint, r, to trace out the entire frontier as depicted in <xref ref-type="fig" rid="fig1">Figure 1</xref>. <xref ref-type="table" rid="table1">Table 1</xref> includes the parameters used to make the calculations in the figures.</p></sec><sec id="s2_2"><title>2.2. Maximize Utility with a General Linear Constraint</title><p>The utility function in Equation (1.2) can also be computed with a general linear constraint. This form of utility is particularly convenient because it is expressed as a quadratic programming problem. We introduce a new constant into this utility, the risk aversion level, which is essentially the slope of the utility function (in variance) and reintroduce the return vector:</p><p>a - level of risk aversion.</p><p>r - n &#215; 1 vector of expected returns.</p><p>The problem is as follows:</p><p>max w w ′ r − a w ′ Ω w     s .t .     w ′ E = c .</p><p>We set up the Lagrangian:</p><p>L = w ′ r − a w ′ Ω w − ( w ′ E − c ) Λ .</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Numerical values used to build the figures</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >r</th><th align="center" valign="middle" >b</th><th align="center" valign="middle"  colspan="4"  >Ω</th></tr></thead><tr><td align="center" valign="middle" >0.06 0.10 0.13 0.04</td><td align="center" valign="middle" >0 1 0 0</td><td align="center" valign="middle" >0.0111 0.0064 0.0052 0.0013</td><td align="center" valign="middle" >0.0064 0.0184 0.0095 0.0084</td><td align="center" valign="middle" >0.0052 0.0095 0.0178 0.0068</td><td align="center" valign="middle" >0.0013 0.0084 0.0068 0.0125</td></tr></tbody></table></table-wrap><p>Every figure and calculation done in this paper were built with these three matrices. These values were chosen somewhat arbitrarily for ease of exposition. The final examples are mildly representative of a realistic scenario even though the values are not calculated from real world data.</p><p>Differentiate the function and set the conditions to zero:</p><p>r − 2 a Ω w − E Λ = 0 .</p><p>Solve for the weight vector:</p><p>w = 1 2 a Ω − 1 ( r − E Λ ) .</p><p>Reapply the constraint condition to calculate Λ :</p><p>w ′ E = 1 2 a ( r − E Λ ) ′ Ω − 1 E = c .</p><p>Λ = ( E ′ Ω − 1 E ) − 1 ( E ′ Ω − 1 r − 2 a c ′ ) .</p><p>And lastly, plug Λ into the weight vector, do some rearranging, and this yields the optimal weights:</p><p>w = Ω − 1 E ( E ′ Ω − 1 E ) − 1 c ′ + 1 2 a Ω − 1 ( r − E ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 r ) . (2.2.1)</p><p>This expression is messier than above but after the rearranging, one can see that the final solution in Equation (2.2.1) is simply the minimum variance solution from Equation (2.1.1) shifted by an additional term.</p></sec><sec id="s2_3"><title>2.3. Optimize in Excess Return/Tracking Error</title><p>The problem proposed by Roll [<xref ref-type="bibr" rid="scirp.96081-ref7">7</xref>] suggests that delegated managers optimize in excess return/tracking error volatility (TEV) rather than mean/variance. The solutions to these problems are identical to the solutions derived in sections 2.1 and 2.2 albeit with a different weight vector and different constraints. We need to introduce another weight vector into the process:</p><p>b - n &#215; 1 constant vector of weights (for a benchmark).</p><p>Tracking error variance and excess return are then defined as follows:</p><p>( w − b ) ′ Ω ( w − b ) - Tracking error variance (TEV).</p><p>( w − b ) ′ r - Excess return.</p><p>Since b is a constant vector, we could substitute it into any of the expressions in 2.1 and 2.2 and still have a valid derivation. Thus, similar to Equation (2.1.1), the weight vector that solves the TEV minimization problem:</p><p>min w ( w − b ) ′ Ω ( w − b )     s .t .     ( w − b ) ′ E = c .</p><p>is:</p><p>( w − b ) = Ω − 1 E ( E ′ Ω − 1 E ) − 1 c ′ . (2.3.1)</p><p>And, the weight vector that solves a utility maximization problem:</p><p>max w ( w − b ) ′ r − a ( w − b ) ′ Ω ( w − b )     s .t .     ( w − b ) ′ E = c .</p><p>is:</p><p>( w − b ) = Ω − 1 E ( E ′ Ω − 1 E ) − 1 c ′ + 1 2 a Ω − 1 ( r − E ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 r ) . (2.3.2)</p><p>Because we are in excess return space, the constraints on this problem need to change slightly. To replicate the solutions in Roll (1992), the constraint vector needs to be set as follows:</p><p>c = [ 0 r w − r b ] .</p><p>Taken the two vectors, w and b, which should each normalize to 1, their difference is 0. Also the second constraint is now the excess return constraint, r w − r b , in order to draw the curves appropriately. See <xref ref-type="fig" rid="fig3">Figure 3</xref> for a representation of these curves drawn in tracking error space.</p></sec></sec><sec id="s3"><title>3. Optimize with an Additional Tracking Error Constraint</title><p>The variance optimization problem given a tracking error constraint was originally solved in the context of portfolio choice in Jorion [<xref ref-type="bibr" rid="scirp.96081-ref8">8</xref>] . In this framework, he used the standard Markowitz constraints, the length of the weight vector and the weighted average expected return, to derive an ellipse in mean/variance space inside of which contains all portfolios with tracking error less than the constraint. In this section, which is the primary contribution of this essay, I derive “Jorion’s ellipse” with the general linear constraint of Best and Grauer [<xref ref-type="bibr" rid="scirp.96081-ref3">3</xref>] . This contribution unifies the analysis of linear constraints on the portfolio optimization model and the significance of relative portfolio optimization to modern delegated contracting. The closed-form nature of these solutions gives analysts a conceptual, theoretical framework to perform comparative analysis without having to rely on numerical simulation.</p><sec id="s3_1"><title>3.1. Optimize Var with General Linear and Tracking Error Constraints</title><p>As with before, we will need a couple new variable definitions:</p><p>T 2 - Tracking error variance constraint.</p><p>λ - Lagrange multiplier for the tracking error constraint.</p><p>We, once again, begin by stating the problem:</p><p>min w w ′ Ω w     s .t .     w ′ E = c ,   ( w − b ) ′ Ω ( w − b ) = T 2 .</p><p>The Lagrangian then follows:</p><p>L = w ′ Ω w − ( w ′ E − c ) Λ − λ ( ( w − b ) ′ Ω ( w − b ) − T 2 ) .</p><p>Next we differentiate by w as set the resulting equations equal to zero:</p><p>2 Ω w − E Λ − λ ( 2 Ω w − 2 Ω b ) = 0 .</p><p>Solving for w yields the following vector:</p><p>w = 1 2 ( 1 − λ ) Ω − 1 ( E Λ − 2 λ Ω b ) .</p><p>Applying the first constraint allows us to calculate one of the Lagrangian terms, Λ :</p><p>w ′ E = 1 2 ( 1 − λ ) ( Λ ′ E ′ − 2 λ b ′ Ω ) Ω − 1 E = c .</p><p>Λ = ( E ′ Ω − 1 E ) − 1 ( 2 ( 1 − λ ) c ′ + 2 λ E ′ b ) .</p><p>Substituting Λ back into the vector for w yields:</p><p>w = Ω − 1 E ( E ′ Ω − 1 E ) − 1 c ′ + λ 1 − λ ( Ω − 1 E ( E ′ Ω − 1 E ) − 1 E ′ b − b ) .</p><p>Because the second condition is easier to calculate if we consider the vector (w – b), we can subtract b from the vector above:</p><p>( w − b ) = ( Ω − 1 E ( E ′ Ω − 1 E ) − 1 c ′ − b ) + λ 1 − λ ( Ω − 1 E ( E ′ Ω − 1 E ) − 1 E ′ b − b ) .</p><p>Applying the condition yields the following equation, after some minor simplification:</p><p>( w − b ) ′ Ω ( w − b ) = ( λ 1 − λ ) 2 ( σ b 2 − b ′ E ( E ′ Ω − 1 E ) − 1 E ′ b )     + 2 ( λ 1 − λ ) ( σ b 2 − b ′ E ( E ′ Ω − 1 E ) − 1 E ′ b )     + c ( E ′ Ω − 1 E ) − 1 c ′ − 2 b ′ E ( E ′ Ω − 1 E ) − 1 c ′ + σ b 2 = T 2 .</p><p>where σ b 2 = b ′ Ω b , the variance of the benchmark.</p><p>Lastly, after some profuse and tedious elementary algebra, and an application of the quadratic equation, the constant can be solved for:</p><p>( λ 1 − λ ) = − 1 &#177; ( 1 − ( c + 2 b ′ E ) ( E ′ Ω − 1 E ) − 1 c ′ − σ b 2 + T 2 σ b 2 − b ′ E ( E ′ Ω − 1 E ) − 1 E ′ b ) 1 / 2 .</p><p>And the optimal solution for w is obtained by plugging in the constant:</p><p>w = Ω − 1 E ( E ′ Ω − 1 E ) − 1 c ′     − ( 1 &#177; ( 1 − ( c + 2 b ′ E ) ( E ′ Ω − 1 E ) − 1 c ′ − σ b 2 + T 2 σ b 2 − b ′ E ( E ′ Ω − 1 E ) − 1 E ′ b ) 1 / 2 )     ⋅ ( Ω − 1 E ( E ′ Ω − 1 E ) − 1 E ′ b − b ) . (3.1.1)</p><p>This weight vector is the solution to the problem of optimizing variance given a general linear and a tracking error constraint. Additionally, as before, this vector is a transformation of the general minimum variance solution shifted by a term in two directions. The &#177; draws out half of the curve in one direction and the other half in the other direction. The linear constraints in combination with the tracking error constraint can cause there to be an infinite number, two, one or no real solutions to the problem. It is certainly notable that the two constraints can work in conflict and generate imaginary vectors. Essentially, if the linear constraints generate a curve that doesn’t pass through the tracking error constrained ellipse, then there will be no real solutions to this problem.</p></sec><sec id="s3_2"><title>3.2. Optimize Utility Given General Linear &amp; Tracking Error Constraints</title><p>The utility optimization problem can also be important to principals. Knowing the best possible outcome given a tracking error constraint could be a useful metric when deciding how strict to make this constraint. No new variables are needed. The statement of the problem is as follows:</p><p>max w w ′ r − a w ′ Ω w     s .t .     w ′ E = c ,   ( w − b ) ′ Ω ( w − b ) = T 2 .</p><p>The Lagrangian is then set up:</p><p>L = w ′ r − a w ′ Ω w − ( w ′ E − c ) Λ − λ ( ( w − b ) ′ Ω ( w − b ) − T 2 ) .</p><p>Then it is differentiated and set to zero:</p><p>r − 2 a Ω w − E Λ − 2 λ ( Ω w − Ω b ) = 0 .</p><p>Before moving to the next step, it is convenient to add 2 a Ω b to both sides:</p><p>r − 2 a ( Ω w − Ω b ) − E Λ − 2 λ ( Ω w − Ω b ) = 2 a Ω b .</p><p>Then we should solve for w:</p><p>w = 1 2 ( a + λ ) Ω − 1 ( r − 2 a Ω b − E Λ ) + b .</p><p>Applying the first constraint yields:</p><p>w ′ E = 1 2 ( a + λ ) ( r ′ − 2 a Ω b − E Λ ) ′ Ω − 1 E + b ′ E = c .</p><p>And solving for Λ gives the following:</p><p>Λ = ( E ′ Ω − 1 E ) − 1 ( E ′ Ω − 1 r − 2 a E ′ b − 2 ( a + λ ) ( c ′ − E ′ b ) ) .</p><p>Next we plug Λ back into the equation for w:</p><p>( w − b ) = ( 1 2 ( a + λ ) ) Ω − 1 ( ( r − 2 a Ω b ) − E ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 ( r − 2 a Ω b ) )     + Ω − 1 E ( E ′ Ω − 1 E ) − 1 ( c ′ − E ′ b ) .</p><p>Applying the final constraint yields the following after some algebra and cancellation:</p><p>( w − b ) ′ Ω ( w − b ) = ( 1 2 ( a + λ ) ) 2 ( r − 2 a Ω b ) ′ ( Ω − 1 − Ω − 1 E ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 ) ( r − 2 a Ω b )         + ( c − b ′ E ) ( E ′ Ω − 1 E ) − 1 ( c ′ − E ′ b ) = T 2 .</p><p>Solving for the constant is a simple process:</p><p>( 1 2 ( a + λ ) ) = &#177; ( T 2 − ( c − b ′ E ) ( E ′ Ω − 1 E ) − 1 ( c ′ − E ′ b ) ( r − 2 a Ω b ) ′ ( Ω − 1 − Ω − 1 E ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 ) ( r − 2 a Ω b ) ) 1 / 2 .</p><p>And plugging this back into the equation for w yields the optimal weight vectors for utility when it is constrained by a general linear constraint and a tracking error constraint:</p><p>w = Ω − 1 E ( E ′ Ω − 1 E ) − 1 c ′ − Ω − 1 E ( E ′ Ω − 1 E ) − 1 E ′ b + b     &#177; ( T 2 − ( c − b ′ E ) ( E ′ Ω − 1 E ) − 1 ( c ′ − E ′ b ) ( r − 2 a Ω b ) ′ ( Ω − 1 − Ω − 1 E ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 ) ( r − 2 a Ω b ) ) 1 / 2     ⋅ Ω − 1 ( ( r − 2 a Ω b ) − E ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 ( r − 2 a Ω b ) ) . (3.2.1)</p><p>One can see that this solution is also a shift from the minimum variance portfolio. Also, because the utility function is a parabola and the feasible set is an ellipse, the &#177; gives two solutions, one for the utility minimization problem and one for the utility maximization problem. This problem suffers from the same concerns regarding conflicting constraints. It is possible for this problem to have no real solutions if the constraints are chosen carelessly.</p></sec></sec><sec id="s4"><title>4. Maximize the Sharpe Ratio with General Linear Constraints</title><p>With the introduction of the Capital Asset Pricing Model, finance theory went the direction of assuming a risk free asset when analyzing linear asset pricing problems. The Sharpe Ratio, introduced in Sharpe [<xref ref-type="bibr" rid="scirp.96081-ref12">12</xref>] and defined below:</p><p>Sharpe   Ratio = w ′ r − r f ( w ′ Ω w ) 1 / 2 ,</p><p>is popular even today and assumes a risk free asset, which has expected return, r f . This ratio is simply the slope of a line connecting the risk free asset to any risky point in the Markowitz feasible set. See <xref ref-type="fig" rid="fig4">Figure 4</xref> for representations of this line.</p><p>Maximizing the slope of this line in mean/variance space generates what is known as the Capital Market Line (CML) and its intersection with the minimum variance envelope is called the market portfolio.</p><p>The problem with the risk free asset in the context of a general solution to the variance minimization problem is that the covariance matrix, Ω , is no longer invertible because the risk free asset creates a row and column of zeros. Thus, Equation (2.1.1) is not a valid solution in this scenario. Finding the CML under the standard linear constraints of portfolio optimization is rather straight forward but with a general linear constraint, it becomes slightly more complicated.</p>Maximizing the Sharpe Ratio with a General Linear Constraint<p>Below is the solution to the Sharpe Ratio maximization problem. It is important to note that the solution to the following problem gives the weight vector for what is called the Market Portfolio if the constraints are the standard weight and return constraints. The statement of the problem is below:</p><p>max w S = w ′ r − r f ( w ′ Ω w ) 1 / 2     s .t .     w ′ E = c .</p><p>This formulation doesn’t lend itself well to quadratic programming, so in order to solve it we will use a well-known result called the mutual fund theorem, from Merton [<xref ref-type="bibr" rid="scirp.96081-ref2">2</xref>] , or the two fund theorem. Essentially, a linear combination of any two portfolios on the efficient envelope also resides on the efficient envelope. Thus, the solution to the maximization problem above is a linear combination of any other two portfolios on the efficient set. Therefore, I will use two of the results I already derived, w ∗ the variance minimization solution, Equation (2.1.1), and w U the utility maximization solution, Equation (2.2.1).</p><p>w ∗ = Ω − 1 E ( E ′ Ω − 1 E ) − 1 c ′ .</p><p>w U = Ω − 1 E ( E ′ Ω − 1 E ) − 1 c ′ + Ω − 1 ( r − E ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 r ) .</p><p>And, their linear combination is below with weight x:</p><p>w = ( 1 − x ) w ∗ + x w U .</p><p>For simplicity, I define a new weight vector w H ,</p><p>w H = w U − w ∗ = Ω − 1 ( r − E ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 r ) .</p><p>And, the weight vector, w, becomes:</p><p>w = w ∗ + x w H .</p><p>Plugging this vector into our function for S yields:</p><p>S = w ′ ∗ r + x w ′ H r − r f ( w ′ ∗ Ω w ∗ − 2 x w ′ ∗ Ω w H + x 2 w ′ H Ω w H ) 1 / 2 .</p><p>w H is particularly convenient because it can be shown that</p><p>w ′ ∗ Ω w H = c ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 r − c ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 r = 0.</p><p>Thus, S becomes:</p><p>S = w ′ ∗ r + x w ′ H r − r f ( w ′ ∗ Ω w ∗ + x 2 w ′ H Ω w H ) 1 / 2 .</p><p>Now basic calculus can be applied to find the maximum with respect to x.</p><p>d S d x = − x w ′ H Ω w H ⋅ ( w ′ ∗ r − r f ) − w ′ H r ⋅ w ′ ∗ Ω w ∗ ( w ′ ∗ Ω w ∗ + x 2 w ′ H Ω w H ) 3 / 2 = 0.</p><p>And the solution follows easily:</p><p>x = w ′ H r ⋅ w ′ ∗ Ω w ∗ w ′ H Ω w H ⋅ ( w ′ ∗ r − r f ) .</p><p>Next, we need to calculate the expected returns and variances to reparametrize the expression into the proper variables and find the constant:</p><p>x = r ′ Ω r − r ′ Ω E ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 r r ′ Ω r − r ′ Ω E ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 r ⋅ c ( E ′ Ω − 1 E ) − 1 c ′ c ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 r − r f = c ( E ′ Ω − 1 E ) − 1 c ′ c ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 r − r f .</p><p>And lastly, plug this constant back into the function for w = w ∗ + x w H :</p><p>w = Ω − 1 E ( E ′ Ω − 1 E ) − 1 c ′ + ( c ( E ′ Ω − 1 E ) − 1 c ′ c ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 r − r f )     ⋅ Ω − 1 ( r − E ( E ′ Ω − 1 E ) − 1 E ′ Ω − 1 r ) . (4.1.1)</p><p>Thus, this is the weight vector that maximizes the Sharpe Ratio given a general linear constraint, or equivalently, if the constraints are the standard Markowitz constraints, this is the weight vector of the market portfolio. It is expressed in the familiar and convenient form of a shift to the minimum variance weight vector.</p></sec><sec id="s5"><title>5. Conclusions</title><p>In this paper, I derive the solutions to portfolio choice problems with a general linear constraint when they are also subject to a tracking error constraint. The primary method used to derive the solutions is quadratic programming. This essay connects the research of some of the most fundamental papers in the field of portfolio choice and generalizes a solution to the Tracking Error Volatility constraint problem of Jorion [<xref ref-type="bibr" rid="scirp.96081-ref8">8</xref>] . I set up the objective function in both variance minimization and utility maximization. Although Jorion solves the variance minimization problem for specific constraints, my method solves the problem with a general constraint. Moreover, Jorion does not derive the utility maximization problem with a tracking error constraint, which I have done in the essay with a general linear constraint.</p><p>Additionally, I have revisited the problems of Best and Grauer [<xref ref-type="bibr" rid="scirp.96081-ref3">3</xref>] and reformulated them so the similarity is more apparent between their solutions and mine. I have discussed how the general constraint matrices need to be set in order to replicate the problems of Markowitz [<xref ref-type="bibr" rid="scirp.96081-ref1">1</xref>] and Roll [<xref ref-type="bibr" rid="scirp.96081-ref7">7</xref>] . I use my derivations of the general linear constrained problems to tackle the problem of maximizing the Sharpe Ratio and derive the optimal solution to this problem. Maximizing the Sharpe ratio can be considered an alternative form of utility and is a very important concept in the foundation of linear asset pricing, particularly when a risk free rate is assumed to exist.</p><p>There are limitations to empirical research given these derivations. Real world data is oftentimes too volatile to apply the methods of modern portfolio theory directly. These derivations are best used as a conceptual framework to make theoretical assertions about the behavior of investment managers. The implications from these theories can, however, be creatively tested and existing empirical data can be used to discover nuances of the decision-making or real world investment managers. One example in particular that this model exposes is that managers that are incentivized by relative risk and return metrics will usually build portfolios that are riskier than a principal would desire. Although these equations likely couldn’t test that. This result could be tested in numerous other ways.</p><p>There are many other examples of constrained portfolio choice problems that fall into the framework outlined in this paper. Although my final expressions are cumbersome, they are also left in their purist form so the reader can see the form unobstructed by substitution. The general linear constraint is very useful when considering real world constraints to asset management, and the tracking error constraint is a pervasive, explicit and implicit constraint to which almost all asset managers abide. This paper will hopefully allow researchers to compare the effects of differing constraints on portfolios in a more approachable manner, with the solutions readily available.</p></sec><sec id="s6"><title>Conflicts of Interest</title><p>The author declares no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s7"><title>Cite this paper</title><p>Stowe, D.L. (2019) Portfolio Mathematics with General Linear and Quadratic Constraints. Journal of Mathematical Finance, 9, 675-690. https://doi.org/10.4236/jmf.2019.94034</p></sec></body><back><ref-list><title>References</title><ref id="scirp.96081-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Markowitz, H. (1952) Portfolio Selection. The Journal of Finance, 7, 77-91. https://doi.org/10.1111/j.1540-6261.1952.tb01525.x</mixed-citation></ref><ref id="scirp.96081-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Merton, R.C. (1972) An Analytic Derivation of the Efficient Portfolio Frontier. Journal of Financial and Quantitative Analysis, 7, 1851-1872. https://doi.org/10.2307/2329621</mixed-citation></ref><ref id="scirp.96081-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Best, M.J. and Grauer, R.R. (1990) The Efficient Set Mathematics When Mean-Variance Problems Are Subject to General Linear Constraints. Journal of Economics and Business, 42, 105-120. https://doi.org/10.1016/0148-6195(90)90027-A</mixed-citation></ref><ref id="scirp.96081-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Binsbergen, J., Brandt, M. and Koijen, R. (2008) Optimal Decentralized Investment Management. Journal of Finance, 63, 1849-1895. https://doi.org/10.1111/j.1540-6261.2008.01376.x</mixed-citation></ref><ref id="scirp.96081-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Blake, D.A., Rossi, G., Timmermann, A., Tonks, I. and Wermers, R. (2013) Decentralized Investment Management: Evidence from the Pension Fund Industry. The Journal of Finance, 68, 1133-1178. https://doi.org/10.1111/jofi.12024</mixed-citation></ref><ref id="scirp.96081-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Cuoco, D. and Kaniel, R. (2011) Equilibrium Prices in the Presence of Delegated Portfolio Management. Journal of Financial Economics, 101, 264-296. https://doi.org/10.1016/j.jfineco.2011.02.012</mixed-citation></ref><ref id="scirp.96081-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Roll, R. (1992) A Mean/Variance Analysis of Tracking Error. The Journal of Portfolio Management, 18, 13-22. https://doi.org/10.3905/jpm.1992.701922</mixed-citation></ref><ref id="scirp.96081-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Jorion, P. (2003) Portfolio Optimization with Tracking-Error Constraints. Financial Analysts Journal, 59, 70-82. https://doi.org/10.2469/faj.v59.n5.2565</mixed-citation></ref><ref id="scirp.96081-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Bertrand, P. (2010) Another Look at Portfolio Optimization under Tracking Error Constraints. Financial Analysts Journal, 66, 78-90. https://doi.org/10.2469/faj.v66.n3.2</mixed-citation></ref><ref id="scirp.96081-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Chen, L., Peng, J., Zhang, B. and Rosyida, I. (2017) Diversified Models for Portfolio Selection Based on Uncertain Semivariance. International Journal of Systems Science, 48, 637-648. https://doi.org/10.1080/00207721.2016.1206985</mixed-citation></ref><ref id="scirp.96081-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Kar, M.B., Kar, S., Guo, S., Li, X. and Majumder, S. (2019) A New Bi-Objective Fuzzy Portfolio Selection Model and Its Solution through Evolutionary Algorithms. Soft Computing, 23, 4367-4381. https://doi.org/10.1007/s00500-018-3094-0</mixed-citation></ref><ref id="scirp.96081-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Sharpe, W.F. (1966) Mutual Fund Performance. The Journal of Business, 39, 119-138. https://doi.org/10.1086/294846</mixed-citation></ref></ref-list></back></article>