To close the series of articles on how to compute fractals using NumPy array operations, I want to show you one more thing: how to create animated GIFs from those NumPy arrays. The first point to make is that…

As a variation on the Mandelbrot fractal that I’ve computed before using NumPy array operations, I’d like to show you a few Julia fractals with some color added. I’ll concentrate on the Julia sets that are generated using the same quadratic polynomial that is used to…

If you use NumPy for numerical computations, which I recommend if you program in Python, one of the key things to do is to process entire arrays in one go as much as possible (this is also true for MATLAB). Using these so-called vectorized operations makes sure that…

Say that you have implemented a low-pass FIR filter with the correct cut-off frequency and transition band, but that you are not quite happy with the suppression of frequencies in the stop band. A straightforward trick to put in your toolbox is that you can double the suppression by…

Do you know what a quine is? Here is an example in Python…

The article Why use Symmetrical FIR Filters with an Odd Length? mentiones that symmetrical FIR filters are linear phase. In this follow-up article, I simply want to show you what the phase response of these filters actually looks like…

In my articles on filter design I mostly focus on a rather specific subset of all possible filters, namely symmetrical FIR filters with an odd number of coefficients. Each of these properties (1. FIR, 2. symmetrical, and 3. an odd number of coefficients) was chosen for…

This article is a wrap-up of what happened on TomRoelandts.com in 2017.

ʻOumuamua is the first ever interstellar body that was observed while it passed through our solar system, which makes it a really cool object. Hundreds of thousands of asteroids and comets have been observed, but all of them originated from within our own solar system. It was known that…

A Hadamard matrix is a matrix with all elements equal to \(+1\) or \(-1\), and for which the rows are mutually orthogonal. If you pick two rows from the matrix and write it as vectors \(\bf x\) and \(\bf y\), then these are orthogonal if…