Linear Algebra for Quantitative Finance › Matrices in Finance

A matrix is a rectangular array of numbers. While vectors represent single entities (one portfolio, one day of returns), matrices represent collections and relationships. These include return history across time, covariances between assets, and transformations between coordinate systems.

The jump from vectors to matrices is conceptually significant. A vector captures one snapshot, such as today's returns or a single portfolio's weights. A matrix can store an entire dataset or a specified set of pairwise summaries. When we write down a covariance matrix, we're encoding each asset's variance and every pairwise covariance in the investment universe. When we write down a return matrix, we're capturing the complete history of how multiple assets performed over multiple time periods. This compression of information into a structured rectangular array is what makes quantitative finance computationally tractable.

An $m \times n$ matrix has $m$ rows and $n$ columns. In finance, we commonly organize data with rows as time periods and columns as assets:

Here $r_{t,i}$ is the return of asset $i$ on day $t$. This $4 \times 3$ matrix contains 4 days of returns for 3 assets.

This chapter uses time as rows and assets as columns. Under that convention, each row represents the cross-section of the represented asset universe at one moment, while each column represents the complete time series of one asset. Extracting a row gives you all represented assets' returns on a specific day; extracting a column gives you one asset's return history. Other sources and libraries may use the transpose, so check the documented orientation before applying a matrix operation.

Matrix multiplication is the workhorse operation of linear algebra. For matrices $\mathbf{A}$ (size $m \times n$) and $\mathbf{B}$ (size $n \times p$), the product $\mathbf{C} = \mathbf{A}\mathbf{B}$ is an $m \times p$ matrix where:

where:

• $c_{ij}$: element in row $i$, column $j$ of the result matrix $\mathbf{C}$
• $a_{ik}$: element in row $i$, column $k$ of matrix $\mathbf{A}$
• $b_{kj}$: element in row $k$, column $j$ of matrix $\mathbf{B}$

Each element of the result is a dot product of a row from $\mathbf{A}$ with a column from $\mathbf{B}$.

Each element of the result shows how much the $i$-th row of $\mathbf{A}$ aligns with the $j$-th column of $\mathbf{B}$. In financial terms, when multiplying a return matrix by a weight vector, each resulting element aggregates the weighted contributions of all assets for that time period. You can think of matrix multiplication as performing many dot products simultaneously. The $(i,j)$ entry of the product answers the question "how does row $i$ of the first matrix relate to column $j$ of the second?"

This interpretation illuminates why matrix multiplication has the dimensional requirements it does. The dot product requires vectors of equal length, so for each row-column pair to produce a dot product, the row length (number of columns in $\mathbf{A}$) must equal the column length (number of rows in $\mathbf{B}$). The result matrix takes its row count from $\mathbf{A}$ and its column count from $\mathbf{B}$ because we compute one number for each possible row-column pairing.

A financial application is computing portfolio returns across multiple days. If $\mathbf{R}$ is a $T \times N$ matrix of simple returns ($T$ days, $N$ assets) and $\mathbf{w}$ is an $N \times 1$ weight vector, then $\mathbf{R}\mathbf{w}$ performs one dot product per day. These are portfolio simple returns when the stated weights apply at the beginning of every day, as in a portfolio rebalanced between days without modeled trading costs. A buy-and-hold allocation generally has different beginning weights after prices change, so its weights must be updated for each period. The result has shape $T\times1$ in this column-vector notation:

For real-valued asset returns with finite second moments, the covariance matrix captures how assets move together. For $N$ assets, it's an $N \times N$ symmetric matrix where element $(i,j)$ is the covariance between assets $i$ and $j$:

where:

• $\Sigma_{ij}$: the $(i,j)$ element of the covariance matrix $\Sigma$
• $r_i, r_j$: returns of assets $i$ and $j$
• $\mu_i, \mu_j$: expected (mean) returns of assets $i$ and $j$
• $\mathbb{E}[\cdot]$: the expectation operator

The formula shows what covariance measures: we're looking at the product of deviations from the mean. When asset $i$ is above its average and asset $j$ is also above its average, the product is positive. When both are below average, the product is again positive. But when one is above and the other below, the product is negative. By averaging these products across many observations, covariance measures positive or negative linear co-movement. Its absolute magnitude depends on the return scale, so a small covariance alone does not establish weak linear association; compare it with the product of the standard deviations, or use correlation. Zero covariance does not imply independence without additional assumptions; for example, a variable and its square can be dependent while having zero covariance under a symmetric distribution.

The diagonal elements $\Sigma_{ii}$ are variances of individual assets. The off-diagonal elements measure co-movement. Positive covariance means assets tend to move together. Negative covariance means they move oppositely.

Portfolio variance demonstrates the power of matrix notation. For portfolio weights $\mathbf{w}$ fixed over the modeled return distribution and an asset covariance matrix $\Sigma$ for that same distribution, the portfolio variance is:

where:

• $\sigma_p^2$: the variance of the portfolio's returns
• $\mathbf{w}$: the $N \times 1$ vector of portfolio weights
• $\mathbf{w}^T$: the transpose of $\mathbf{w}$ (a $1 \times N$ row vector)
• $\Sigma$: the $N \times N$ covariance matrix of asset returns

This compact formula packs a lot of computation: it accounts for each asset's variance and all pairwise covariances, weighted by the portfolio allocations. The expression $\mathbf{w}^T \Sigma \mathbf{w}$ is called a quadratic form because if you expand it, you get a polynomial where each term involves products of two weights. It's quadratic in the portfolio allocations.

Expanding for two assets:

where:

• $\sigma_1^2, \sigma_2^2$: the variances of assets 1 and 2
• $\sigma_{12}$: the covariance between assets 1 and 2

The first two terms are variance contributions from each asset. The third term is $2w_1w_2\sigma_{12}$, so its sign depends on both covariance and weights. Negative covariance lowers variance relative to zero covariance when the nonzero weights have the same sign; for opposite-sign positions, that effect reverses.

For nonnegative weights summing to one, portfolio variance is no greater than the weighted average of the individual variances. With both weights positive and covariance below the geometric mean of positive individual variances, this inequality is strict. The code below uses a different benchmark: the weighted average of individual volatilities. That is the portfolio volatility when all positively weighted, nonzero-volatility assets are perfectly positively correlated. Actual portfolio volatility is no greater than this benchmark, and is strictly lower if at least one such pair has correlation below one. The variance and volatility comparisons should not be interchanged.

The matrix formula has the same form for any number of assets. For 500 stocks, it remains $\mathbf{w}^T \Sigma \mathbf{w}$ while the ordered double sum contains 250,000 terms: 500 variance terms and 249,500 off-diagonal terms. Because the covariance matrix is symmetric, only 125,250 entries are unique. Compact notation keeps the calculation manageable as a portfolio grows from a toy example to a production system.

The portfolio volatility is lower than the weighted average of individual volatilities because the assets are imperfectly correlated. This is the mathematical basis of diversification.

The transpose of matrix $\mathbf{A}$, written $\mathbf{A}^T$, flips rows and columns: $(A^T)_{ij} = A_{ji}$.

Several special matrix types appear frequently in finance:

• Symmetric matrices: $\mathbf{A} = \mathbf{A}^T$. Covariance matrices are always symmetric. This property reflects a fundamental reality: the covariance between assets A and B must equal the covariance between B and A, since we're measuring the same relationship from both directions.
• Diagonal matrices: Non-zero elements only on the diagonal. Used to represent variance contributions when assets are uncorrelated, or to scale different variables by different amounts.
• Identity matrix: Diagonal matrix with 1s on the diagonal. Acts as the "1" of matrix multiplication: $\mathbf{I}\mathbf{A} = \mathbf{A}$. It's the matrix equivalent of multiplying by one, leaving any matrix unchanged.
• Positive definite matrices: A real symmetric matrix $\mathbf{A}$ is positive definite when $\mathbf{x}^T \mathbf{A} \mathbf{x} > 0$ for every non-zero real vector $\mathbf{x}$. Valid covariance matrices must be positive semi-definite (allowing zero). This ensures that every portfolio variance $\mathbf{w}^T\Sigma\mathbf{w}$ is non-negative, as variance must be by definition.

The eigenvalue calculation below previews Eigenvalue Decomposition, where we explain these values and their connection to variance. For now, use the nonnegative quadratic-form condition above as the definition of positive semi-definiteness.