Prerequisites: Linear Algebra, Probability & the Z-Transform
The mathematics MIDA1 speaks from its first lecture and never stops to teach — summation, set and quantifier notation; vectors, matrices, the inner product and the least-squares normal equations; expectation, variance and correlation; and the discrete-time signals and Z-transform every model in the course is written in — each tied to the chapter that spends it and the exam problem that grades it.
01 · Motivation
What MIDA1 assumes on day one
The MIDA1 exam is, on every structured paper, four problems: a process-analysis problem, a prediction problem, an identification problem, and open theory. Not one of them can be started without the mathematics on this page — and not one lecture stops to teach it. It is inherited from the Bachelor’s calculus, linear-algebra and probability courses, and from the signals-and-systems background. This chapter collects that inheritance and points each piece at the exact MIDA1 chapter that will spend it.
This material is used from the first lecture and taught in none of them. It carries no exam weight of its own — no past paper tests “take an expectation” directly, and every past paper depends on it.
Every row below is evidence, not a generic syllabus: each prerequisite is here because a later chapter spends it, and the audit cites the line where that first happens.
| you need to be able to… | first spent in | graded by |
|---|---|---|
| read and over an indexed set | ch.01 foundations | AR/MA variance, least squares (ch.03, ch.06) |
| state a property and chase an implication | ch.02 stochastic-processes | Stationarity & representations (ch.02) |
| use , , and a cardinality | ch.02 stochastic-processes | Power spectral density (ch.04) |
| take an inner product and a transpose | ch.01 foundations | Least squares estimation (ch.06) |
| form a linear combination of past values | ch.03 ma-ar-arma-models | AR(1) / MA(2) variance (ch.03) |
| take an expectation and a variance | ch.02 stochastic-processes | Mean, covariance & spectrum (ch.02) |
| read a Pearson correlation coefficient | ch.02 stochastic-processes | Mean, covariance & spectrum (ch.02) |
| recognise a random walk / integrated signal | ch.06 identification | Experiment design (ch.09) |
Two further primers follow that the audit did not surface but the course quietly assumes: the Z-transform and the delay operator (section 05) — the language every model from chapter 03 on is written in — and a one-page least-squares refresher (section 06), the minimisation chapter 06 then applies to data. They are authored background, flagged again where they begin.
How to read this chapter
Skim the ledger. Anything you can already do, skip — each section is self-contained. Anything you cannot, read in full and do its exercise; the chapter that spends it will not slow down for you.
The running example: a room-temperature sensor
One example runs through every section, chosen because it is the object the course studies. A thermostat samples the deviation of a room’s temperature from its setpoint, once a minute. The room has thermal inertia — it forgets slowly — and it is buffeted by random disturbances : a door opening, sun through a window. The simplest honest model is first order:
Memory (the term) is what makes it a dynamical system; the white is the part no model can predict.
By the end of the chapter you will have written this signal’s mean and variance, its correlation across time, its transfer function and pole, and estimated its coefficient from a handful of readings by least squares — every prerequisite exercised on one object. It is deliberately the AR(1) process that chapter 02 analyses, chapter 03 solves with Yule–Walker, and chapter 06 identifies from data.
The prerequisite errors that cost the most marks
- Treating the delay operator as a number. It is the shift , not “one over ” to cancel at will. Mishandle it and every transfer function, spectral factorisation and predictor in chapters 04–07 comes out wrong.
- Confusing the expectation with the sample average . Expectation is the theoretical (ensemble) operator; the average is its estimator. Stationarity, Yule–Walker and the asymptotic analysis of chapter 07 are all stated in — conflate the two and the derivations dissolve.
- Dropping a transpose in the regressor. Chapter 06 stacks data into a matrix and a vector ; one misplaced transpose makes non-conformable and the normal equations meaningless.
02 · Notation
The notation you read in
Three pieces of notation appear on nearly every slide, and none is ever defined: the summation, the set, and the quantifier. Fluency here is not optional — a stationarity definition is one short sentence in this language and a paragraph without it.
- Summation and product notation
adds ; multiplies them. The index and its range are part of the meaning — over “all ” is an infinite sum and only makes sense when it converges.
- Set notation
(membership), (subset), and (union, intersection), and the cardinality (how many elements). Set-builder form reads “the such that…”.
- First-order logic
(“for all”), (“there exists”), and (“implies”). A property that must hold at every time is a statement; a definition that trades one condition for another is an implication.
Watch these three combine. The sample mean of the sensor series over readings is a summation,
An estimator — a number computed from data. Section 04 contrasts it with the theoretical mean .
and weak stationarity, the property chapter 02 opens with, is a single quantified sentence: the mean is constant and the covariance depends only on the lag,
If you cannot read that line — the , the index set , the fact that is a function of alone — chapter 02 will feel like a wall. Read it as: no matter which instant you pick, the average is the same number , and the covariance between two samples depends on how far apart they are, not on where they sit.
Where this is graded
Every “Stationarity & representations” question (ch.02) turns on parsing exactly this quantified statement, and every “Power spectral density” question (ch.04) sums a series over the integers. The notation is half the answer.
A signal has E[y(t)] = 3 for every t, but its covariance Cov(y(t), y(t−τ)) depends on both t and τ. Is it weakly stationary?
03 · Core idea
Linear algebra — the machinery of every estimator
Identification turns data into a vector of numbers and a model into another vector, then measures the gap between them. That is linear algebra, and three operations carry all of it: the linear combination, the inner product, and the transpose.
- Vector
An ordered list , written as a column. In this course a vector is usually a stack of past measurements or a stack of model parameters.
- Linear combination
— scale each vector and add. An AR model is a linear combination: built from past outputs with coefficients .
- Inner product
, a single number. It measures alignment; it is zero exactly when the two vectors are orthogonal.
- Transpose
flips rows and columns; . A matrix transpose is what turns a stack of row-observations into the normal-equations product .
- Euclidean norm
— the length of a vector, and therefore the size of an error when is a residual.
The one identity to keep at your fingertips is that the squared length of a vector is an inner product with itself,
A sum of squares is a squared norm. Every least-squares cost in the course is one of these in disguise.
because the entire prediction-error framework minimises a sum of squared residuals, which is a squared norm, which is an inner product. When chapter 06 writes the least-squares cost as and chapter 05 defines the optimal predictor as the one with minimum mean-square error, they are both measuring a length with this formula.
The sensor, as vectors
Take the readings . To ask “how well does the previous sample explain the next one?”, stack the inputs (each ) against the targets (each ). Their inner product and the self-product are the only two numbers the estimate of will need — section 06 finishes the job.
Where this is graded
“Least squares estimation” and “ARX identification by least squares” (ch.06) are, mechanically, the assembly of and from data. A dropped transpose or a mis-shaped stack is the most common way to lose the whole problem before any statistics begin.
Φ is an N×p matrix of stacked observations (N rows, p parameters). Which product is well-formed and gives the p×p matrix the normal equations need?
04 · Mechanics
Probability — the language of a random signal
A stochastic process is a signal whose value at each instant is a random variable. To describe one you need four numbers-from-randomness: the expectation, the variance, the autocovariance, and the correlation coefficient that normalises it.
- Expectation
, the probability-weighted average of a random variable — its theoretical mean, the number the sample average estimates. Linear: .
- Variance
, the expected squared spread about the mean. Its square root is the standard deviation.
- Autocovariance
— how the signal co-varies with a delayed copy of itself. is the variance.
- Pearson correlation coefficient
— the autocovariance stripped of units, so is perfect alignment and is none.
- Gaussian distribution
, the bell curve fixed entirely by its mean and variance — which is why a stationary Gaussian process is fixed entirely by and .
- White noise
A zero-mean process with for every : no structure across time, the part of a signal no predictor can improve on.
Now watch expectation do real work on the sensor. Assuming it has settled to a stationary regime with zero mean, take the variance of both sides of . Because is white and independent of the past, the cross term vanishes and
This is a one-line Yule–Walker calculation (ch.03) — and it only works because E[·] is linear and E[y(t−1)e(t)] = 0.
Push the same idea one lag further and , so the correlation coefficient of the sensor is . At lag one, : for an AR(1) process the lag-1 correlation coefficient is the model coefficient itself. That single fact links the Pearson correlation you are revising here to the autocorrelation of chapter 02 and the Yule–Walker equations of chapter 03.
Random walk — the edge case that breaks stationarity
Set and the sensor becomes , a random walk: each value is the running sum of every disturbance so far. The variance formula above divides by and blows up — the process is not stationary, it wanders without bound. Chapter 09 removes the wandering by differencing (, which is white and stationary), the “I” in ARIMA. Recognising an integrated signal is a prerequisite you meet again in experiment design.
Where this is graded
“Mean, covariance & spectrum” (ch.02) and every “AR(1) / MA(2) variance” problem (ch.03) are the calculation above, run on the specific model in the paper. The marks are in the linearity of and the vanishing of the noise cross term — not in arithmetic.
For the sensor y(t) = 0.5·y(t−1) + e(t) in steady state, what is the correlation coefficient ρ(1) between consecutive samples?
05 · Worked example
Discrete-time signals & the Z-transform
Authored primer — not an audit finding
The prereq-audit could not surface this section: the delay operator is pure notation, invisible to a concept probe, and “transfer function” registers as taught later (chapter 03). But every model from chapter 03 onward is written in this language, so it is genuine assumed background. Treat it as an appendix the exam still leans on.
MIDA1 never manipulates a signal one sample at a time for long. It packs the whole history into an operator and does algebra instead. The tool is the unit delay operator , defined by
shifts back one step; shifts forward. Some texts write it . It is an operator, never a number.
With it, the sensor’s difference equation collapses. Move the past output to the left:
Now divide — formally, as an operator — to write the output as the noise passed through a transfer function :
A transfer function is a ratio of polynomials in z. Its roots carry all the dynamics.
- Pole
A root of the denominator of — here . Poles govern how the system responds and whether it settles.
- Zero
A root of the numerator. The sensor’s only zero is at .
- Asymptotic stability
A discrete-time system is asymptotically stable exactly when every pole lies strictly inside the unit circle, . The sensor’s pole at qualifies, so its response to a disturbance decays.
The unit circle is the whole stability story in discrete time, and it is not decoration: a pole at is the random walk from the previous section — marginally stable, non-stationary — and a pole outside it is a signal that explodes. Chapter 04 chooses the canonical spectral factor by demanding poles and zeros inside the circle; chapter 05 builds the whitening filter , which is only usable when ‘s zeros are inside the circle so the inverse’s poles are too. Every one of those moves is a statement about where roots sit relative to .
Where the course picks this up
Chapter 03 writes MA, AR and ARMA models as — exactly this ratio, with higher-degree polynomials. Chapter 04 factors spectra into . Chapter 05 predicts by long-division of these polynomials. If is comfortable now, all three are bookkeeping; if it is not, they are fog.
A process has transfer function W(z) = 1 / (1 − 1.2·z⁻¹). Is it asymptotically stable, and what does that imply?
06 · Mechanics
Least squares in one page
Authored primer — the math, not chapter 06
Least squares is taught in chapter 06, so the audit correctly treats it as course content, not a prerequisite. What you need before chapter 06 is the underlying minimisation — set a gradient to zero, get a linear system. That refresher is here; the identification it powers (ARX, the meaning of the regressors) belongs to chapter 06 and is left there.
Before any algebra, the shape of the answer. Least squares is a projection: the residual sticks out perpendicular to everything the model can reach, and that perpendicularity is the normal equations. Every estimator in this course is a version of this picture.
Stack the data. Write each measurement as a target and a row of knowns (the regressor); collect of them into a vector and a matrix . A parameter vector predicts , and least squares picks the that makes the residual shortest — smallest squared norm, the identity from section 03:
Minimise a sum of squared errors — a squared Euclidean norm — over the parameters.
The cost is a quadratic in , so it has a single minimum where its gradient vanishes. Differentiating and setting the result to zero gives the normal equations:
Solvable exactly when ΦᵀΦ is invertible — chapter 06 calls that condition persistent excitation.
Geometrically, is the projection of onto the space its columns span, and the residual is orthogonal to every column — that orthogonality is the normal equations, .
1 · Set up the one-parameter regression
The sensor has a single unknown , so is a scalar and each regressor row is just the previous sample: . The normal equations collapse to one scalar equation.
2 · Assemble the two inner products
Reuse the readings from section 03. The four regression pairs are , so
3 · Solve
With scalars the inverse is a division:
Close to the true from five noisy points — and it would tighten as grows. That convergence-with-data is exactly the consistency question chapter 07 makes precise.
Why can the least-squares solution θ̂ = (ΦᵀΦ)⁻¹Φᵀy fail to exist?
07 · Practice
Work these before chapter 1
Do each on paper before opening its solution. They are ordered to mirror the chapter, and every one is shaped like the MIDA1 problem that will spend the prerequisite.
Try it Exercise 1 — Read a stationarity statement
Task. A process satisfies and . Using the quantified definition, decide whether it is weakly stationary and say which clause decides it.
Solution. Weak stationarity needs both a constant mean and a lag-only covariance . The covariance clause holds — depends on alone. But the mean varies with , so the first clause fails and the process is not weakly stationary. The mean clause is the decider; a lag-only covariance is not sufficient on its own.
Try it Exercise 2 — An inner product and a norm
Task. For and , compute , decide whether they are orthogonal, and give .
Solution. . It is non-zero, so the vectors are not orthogonal (orthogonality means inner product exactly zero). The norm is .
Try it Exercise 3 — Variance of an AR(1)
Task. A stationary zero-mean process obeys with . Find , and the correlation .
Solution. From with : . Then . The correlation coefficient is , so . Note needs no knowledge of — the variance cancels.
Try it Exercise 4 — Difference equation to transfer function
Task. Write in the delay operator, give its transfer function , and state its pole and zero. Is it asymptotically stable?
Solution. Group each side with : . Hence
The pole is and the zero is . The single pole satisfies , so the system is asymptotically stable. (Both roots are inside the unit circle, so this is also a canonical form in the sense of chapter 04 — worth noticing now.)
Try it Exercise 5 — A scalar least-squares fit
Task. You believe (no intercept) and observe the pairs . Use the normal equations to find .
Solution. With one parameter the normal equation is . Numerator: . Denominator: . So . The fitted line is the projection of the data onto the single regressor direction , and the residual is orthogonal to it.
Try it Exercise 6 — Spot the prerequisite in a MIDA1 problem
Task. A prediction problem says: “Given the ARMA process , verify it is in canonical form and compute the optimal one-step predictor.” Which prerequisites from this chapter are you tested on before writing a single predictor term?
Solution. Three. (1) The delay operator (section 05) — reading as a signal-shaping filter rather than a fraction of numbers. (2) Pole/zero location versus the unit circle (section 05) — “canonical form” means the pole and the zero both lie inside , which you check here as and . (3) Expectation and whiteness (section 04) — the predictor is optimal when its error is white, an statement about the residual. The long division that produces the predictor is the easy half; these three decide whether you set it up correctly at all.
Chapter 01 models the sensor as y(t) = 0.5·y(t−1) + e(t) and calls a predictor 'optimal' when its residual is white noise. Which two prerequisites from this chapter does that single sentence rest on?
08 · Recap
Chapter 00 — what you must carry forward
- Weak stationarity is one quantified sentence: a constant mean and a covariance that depends on the lag alone. Read it fluently and chapter 02 opens; stumble on the and it does not.
- A sum of squares is a squared norm is an inner product. is the identity under every least-squares cost and every mean-square error in the course.
- Expectation is linear and the noise cross term vanishes — that pair of facts turns into in one line, and gives . For an AR(1), the lag-1 correlation is the coefficient.
- The delay operator is an operator, not a number. It packs a difference equation into a transfer function , and the dynamics live in the roots.
- Stability is a compass reading against the unit circle. Poles inside decay; a pole on it is a random walk; a pole outside explodes. Canonical forms, whitening filters and spectral factors are all statements about roots versus that circle.
- Least squares is a projection. Minimising gives the normal equations ; they are solvable exactly when is invertible — the condition chapter 06 names persistent excitation.
- Exam radar: this chapter is never asked directly and is assumed in every problem. When a process-analysis or prediction problem feels impossible, the missing piece is usually on this page — most often a mishandled or a confusion between and the sample average.
Next: chapter 01 picks up exactly here — it writes the sensor’s difference equation, frames prediction as an optimisation, and calls a predictor optimal precisely when its residual is the white noise you just defined.