©copyright Raimo Olavi Toivonen 1983-2023. All rights reserved. Last updated on Jan 8, 2023. |
|X(k)|dB (k=0...N/2-1) = 20log10|SUM(w(n)x(n)e-jk(2π/N)n)|, n=0...N-1. |Y(k)|dB (k=0...N/4-1) = |SUM(|X(n)|dB e-jk(2π/N)n)|, n=0...N/2-1. Signal x(n) = xa(nT), n=0...N-1, T is the sample period, xa= analog signal. If xa is the voltage signal of the microphone, then xa(nT) = kPa(nT), where Pa is the pressure and k is the constant factor.
Hamming-window w(n) = 0.54-0.46cos(2πn/(N-1)).
Blackman-Harris-window w(n) = 0.35875-0.48829cos(2πn/(N-1))+0.14128cos(4πn/(N-1))-0.01168 cos(6πn/(N-1)).
The calculation uses a decimation-in-time (DIT) algorithm.
Signal
windowing is done with a
Hamming window, a rectangular window or a
Blackman-Harris window with more than 90 dB.
Analysis
pictures I have coded from the very beginning in
Neon object-oriented
programming language.
Cepstrum on the lin/Hz scale.