Hoppa till huvudinnehållet
Till KTH:s startsida Till KTH:s startsida

EQ2801 Optimal Filtering 7,5 hp

Course memo Autumn 2022-50845

Version 1 – 08/11/2022, 4:31:43 PM

Course offering

Autumn 2022-1 (Start date 29/08/2022, English)

Language Of Instruction

English

Offered By

EECS/Intelligent Systems

Course memo Autumn 2022

Course presentation

This course gives thorough knowledge of linear estimation theory. The main theme of the course is optimal linear estimation, Kalman and Weiner filtering, which are systematic methods to solve estimation problems with applications in several technical disciplines, for example in telecommunications, automatic control and signal processing but also in other disciplines, such as econometrics and statistics.The course also provides an introduction to optimal filtering for non-linear systems. The course assumes familiarity with basic concepts from matrix theory, stochastic processes, and linear systems theory. The course is directed towards the students who intend to work with development and research within these fields.

The following topics are covered; Basic estimation theory, time discrete and time continuous Wiener filters, time discrete Kalman filters, properties of Wiener and Kalman filters, smoothing, Extended Kalman filters, sigma-point filters and particle filters.

Headings denoted with an asterisk ( * ) is retrieved from the course syllabus version Spring 2019

Content and learning outcomes

Course contents

This course gives thorough knowledge of linear estimation theory. The main theme of the course is optimal linear estimation, Kalman and Weiner filtering, which are systematic methods to solve estimation problems with applications in several technical disciplines, for example in telecommunications, automatic control and signal processing but also in other disciplines, such as econometrics and statistics.The course also provides an introduction to optimal filtering for non-linear systems. The course assumes familiarity with basic concepts from matrix theory, stochastic processes, and linear systems theory. The course is directed towards the students who intend to work with development and research within these fields.

The following topics are covered; Basic estimation theory, time discrete and time continuous Wiener filters, time discrete Kalman filters, properties of Wiener and Kalman filters, smoothing, Extended Kalman filters, sigma-point filters and particle filters.

Intended learning outcomes

After successfully completing the course, the student should be able to

•                     Understand to which type of estimation problems linear estimation can be applied.

•                     Understand the relationship between computational complexity, filter structure, and performance.

•                     Understand the relationship between optimal filtering, linear estimation, and Wiener/Kalman filtering.

•                     Approach estimation problems in a systematic way.

•                     Compute, analyze, and modify state space models.

•                     Derive and manipulate the time discrete and time continuous Wiener filter equations and compute the Wiener filter for a given estimation problem.

•                     Derive and manipulate the time discrete Kalman filter equations and compute the Kalman filter for a given estimation problem.

•                     Analyze properties of optimal filters.

•                     Implement Wiener and Kalman filters (time discrete) and state space models using Matlab.

•                     Simulate state space models and optimal filters, analyze the results, optimize the filter performance, and provide a written report on the findings.

•                     Know about common methods for optimal filtering in the case of non-Gaussian noise or non-linear models, such as Extended Kalman filter, sigma point filtering and particle filtering.

•                    Formulate logical arguments, orally and in writing, in a way that is considered valid in scientific publications and presentations within the topic area.

Learning activities

The course is organized jointly with the PhD level course FEM3200.

Lectures

During the 8 weekly lectures, we go through the theory and together discuss important concepts, using clicker questions and group discussions. Therefore, the lectures are primarily offered in real life at campus, but (somewhat primitively) recorded lectures will also be available as a complement.

Weekly homework

As a student, you are expected to spend a large part of the effort in the course, on solving weekly homework problems (in total 6 sets of homework problems). You may discuss the problems with other students in the course, but the handed in solutions should be formulated individually.

Projects

Two home projects should be solved during the course. The tasks involve a combination of pen/paper derivations and Matlab implementations. The project assignments may be solved and reported by groups of no more than 2 students. For Project 1, a short written report should be handed in, whereas Project 2 is reported orally.

Detailed plan

Learning activities Content Preparations
Lec. 1 Least squares estimation, conditional mean Read Anderson&Moore: Sect. 2.3 and 5.2
Lec. 2 Spectral factorization, causal Wiener filtering (continuous+discrete time) Read Simon: Sect. 3.4
Lec. 3 Recursive Bayesian filtering, Kalman Read Theodoridos: Sect. 17.3, Theodoridos: App. 12, Simon: Sect. 15.1
Lec. 4 Innovations and Kalman, Extended Kalman Read Anderson&Moore: Sect. 5.3–5.4, 8.1–8.2, “An innovations approach to least-squares
estimation–Part I: Linear filtering in additive white
noise,” T.Kailath, IEEE Transactions on Automatic
Control, Volume: 13, Issue: 6, Dec 1968., Simon: Sect. 13.2.3
Lec. 5 Innovations form Kalman, time invariant Kalman, frequency domain expressions. Read Anderson&Moore: Sect. 6.3, Chapt. 4
Lec. 6 Smoothing, square root algorithms Read Anderson&Moore: Chapt. 7, Sect. 6.5,
Simon: Chapt. 9, Sect. 6.3, I Bierman, Gerald J., and Catherine L. Thornton. "Numerical compari on of Kalman filter algorithms: Orbit determination case study." Automatica 13.1
(1977): 23-35.
Lec. 7 Importance sampling, particle filtering Read Theodoridis, Chapt. 14.1–14.5 and 17, Simon, Chapt. 13-15
Lec. 8 Particle filtering cont. Unscented Kalman Read Simon, Chapt. 14-15,  Theodoridis, Chapt. 17


Time and location, see Schema HT-2022-276

Preparations before course start

Recommended prerequisites

EQ1220/EQ1270 Signal Theory/Stochastic Signals and Systems, or equivalent

EQ2300 Digital Signal Processing, EQ2401 Adaptive Signal Processing

Literature

The material in the course will be taken from several text books, all of them
being available as E-books through KTH Library:

  • D. Simon, Optimal State Estimation, John Wiley & Sons, 2006
  • B. D. O. Anderson, J. B. Moore, Optimal Filtering, Dover Publications, 2012.
  • S. Theodoridis, Machine Learning; A Bayesian and Optimization Perspective, 2nd ed., Academic press, 2020.

In addition, a few research articles will be used (also available through KTH Library).

Software

The projects and homework require some implementation work, preferably in Matlab.

Examination and completion

Grading scale

A, B, C, D, E, FX, F

Examination

  • INL1 - Homework assignments, 4.5 credits, Grading scale: A, B, C, D, E, FX, F
  • PRO1 - Project assignment, 1.5 credits, Grading scale: A, B, C, D, E, FX, F
  • PRO2 - Project assignment, 1.5 credits, Grading scale: A, B, C, D, E, FX, F

Based on recommendation from KTH’s coordinator for disabilities, the examiner will decide how to adapt an examination for students with documented disability.

The examiner may apply another examination format when re-examining individual students.

The section below is not retrieved from the course syllabus:

Homework assignments ( INL1 )

There are 6 sets of weekly homework. 

Homework number Available from Due
1 Lecture 1 Lecture 2
2 Lecture 2 Lecture 4
3 Lecture 4 Lecture 5
4 Lecture 5 Lecture 6
5 Lecture 6 Lecture 7
6 Lecture 7 Lecture 8

You may discuss the problems with other students in the course, but the handed in homework solutions should be formulated individually.

Project assignments 1 and 2 (PRO1 and PRO2)

The task of Project 1, will be based on the material from Lectures 1-4 and will therefore be made available at Lecture 4. The task of Project 2 will be based on material from Lecture 7 (and 8) and will be made available at Lecture 7. 

The project assignments may be solved and reported by groups of 2 students (or individually, though it's recommended to work two and two). The projects involves a combination and paper/pen calculations and numerical implementation and evaluations.

For Project 1, a short written report documenting the solution and results should be handed in, before Lecture 7. For Project 2, the implementation should be handed in and the results are reported orally, at latest during the first week of study period 2.

Other requirements for final grade

•                     PRO1 – Project assignment, 1.5, grade scale: A, B, C, D, E, FX, F

•                     PRO2 – Project assignment, 1.5, grade scale: A, B, C, D, E, FX, F

•                     INL1 – Homework assignments, 4.5, grade scale: A, B, C, D, E, FX, F

Final grade based on 70% from INL1 and 15% each from PRO1 and PRO2, respectively.

The course requires significant individual effort. Solving the homework problems requires good familiarity with the theory but also an ability to formulate a practical problem using suitable mathematical models and applying the theory to these. The written presentation of solutions and project also provide training in the ability to formulate logical arguments in a way that is considered valid in scientific publications. One of the project assignments is presented in a technical report, the other one in an oral presentation.

Grading criteria/assessment criteria

The different homework sets and projects cover different parts of the intended learning outcomes. The basic idea of the grading is therefore that an acceptable level has to be reached for all the homework sets and both projects, in order to obtain a passing grade. Also, since the ability to formulate valid mathematical arguments (orally and written) is an important part of the intended learning outcomes, this is judged separately and contributes to the final grade. The details are as follows.

For each homework set and for each of the projects, you obtain points according to the following criteria.

Mathematical/technical content, 1-3 points:

#points Criterion
1 point More than 50% of the difficulties of the problem set/project have been solved, at least by providing a valid approach and working out some details. 
2 points More than 75% of the difficulties of the problem set/project have been solved, at least by providing a valid approach and working out some details. In addition, at least two of the problems have been solved fully. 
3 points Full solutions are provided for all problems and only minor details are flawed/missing. 

Presentation, 1-2 points:

(only judged if at least 1 point was obtained for the technical content, and only judged for the problems or parts of the project that have been ``solved'' according to above).

#points Criterion
1 point The presentation manages to convey the main ideas of the solved problems, but not fully convincing in all details.
2 points Clear and convincing presentation of the solved problems.

These points are summed up and result in a grade according to the table below. When calculating the number of points for the final grade, each of the project points are weighted by a factor two, so that the maximum is 6·5+2(5+5)=50 points. Please note that at least 1+1 points is needed for each of the 6 homework sets and for each of the projects, to obtain a passing grade. Note also that only the grading scale A-D is used for the projects, so that a passing grade for the projects give a D (or higher). 

Grade Homeworks (INL1) Projects (PRO1&2) Final grade
E 12-15 - 20-25
D 16-19 2 26-32
C 20-23 3 33-39
B 24-27 4 40-45
A 28-30 5 46-50

 

Opportunity to complete the requirements via supplementary examination

Students who have passed all but one of the homeworks+projects get the chance to do an extra assignment to pass the course. The corresponding homework/project will then get the minimum passing number of points, whereafter the final grade is calculated as described above.

Ethical approach

  • All members of a group are responsible for the group's work.
  • In any assessment, every student shall honestly disclose any help received and sources used.
  • In an oral assessment, every student shall be able to present and answer questions about the entire assignment and solution.

Further information

No information inserted

Round Facts

Start date

29 Aug 2022

Course offering

  • Autumn 2022-50845

Language Of Instruction

English

Offered By

EECS/Intelligent Systems

Contacts

Course Coordinator

Teachers

Examiner