diff --git a/libraries/dsp/doc/pfb.txt b/libraries/dsp/doc/pfb.txt
new file mode 100644
index 0000000000000000000000000000000000000000..cf131191fa887c726841e72082b86b48bc42b537
--- /dev/null
+++ b/libraries/dsp/doc/pfb.txt
@@ -0,0 +1,14 @@
+
+1) Quantization in LOFAR2 station: https://support.astron.nl/confluence/display/L2M/L4+SDPFW+Decision%3A+LOFAR2.0+SDP+Firmware+Quantization+Model
+2) Fixed point numbers: https://support.astron.nl/confluence/display/L2M/L3+SDP+Decision%3A+Definition+of+fixed+point+numbers
+3) Rounding: https://support.astron.nl/confluence/display/L2M/L4+SDPFW+Decision%3A+Number+representation%2C+resizing+and+rounding
+
+4) LOFAR station Polyphase Filterbank (PFB) model in Matlab: https://git.astron.nl/desp/apertif_matlab/-/blob/master/matlab/ The apertif_matlab_readme.txt gives a brief desciption of all files in this repository. The one_pfb.m runs the model.
+
+5) PFB implementation in Apertif and in LOFAR2 station firmware: https://git.astron.nl/desp/hdl/-/blob/master/libraries/dsp/fft/doc/ASTRON_SP_054_filterbank_spec_part2.pdf
+
+Het FIR filter als geheel is symmetrisch, maar per phase is het niet symmetrisch.
+
+Het FIR filter heet ook prototype filter en geeft de filtering aan per output bin (= subband). Alle bins ondergaan dezelfde filtering. In een rechtstreekse implementatie zou je steeds na elk input sample de FFT doen en dan voor elke FFT output bin het FIR filter. Dat is echter niet efficient, omdat je met een factor D kunt downsamplen, wat betekent dat je D-1 samples weggooid. Daarom worden de FIR filters aan de output bins als het ware door de FFT heen gedrukt naar de input. Dan wordt het een polyphase FIR filter en pas je de FFT niet toe na elk sample, maar na elk blok van D samples. Hiermee wordt de implementatie dan wel efficient, omdat alleen gedowsamplede outputs berekend worden. In ons geval is D = N_points = 1024 van de FFT, daarom heet dit een kritisch gesamplede filterbank. Als je D < N_points kiest krijg je een zogenaamde oversampled filterbank.
+
+De FFT kun je zien als een bank met mixers die het input signaal voor elke bin naar 0 Hz mixen. De bin waarde is een complex getal, omdat je amplitude en phase nodig hebt om de bin golf uniek te kunnen definieren. Meerdere bin waardes sequentieel vormen het gedownsamplede bin signaal in de tijd. De downsample factor D = N_points als de FFT per block van N_point.