©copyright Raimo Olavi Toivonen 19852019. All rights reserved. Last updated 1 Jun 2019.

Home  Signal displays  Multi-envelope displays  FFT displays  LPC displays  Auditory displays 1983−  Loudness displays  Loudness curves  Timbre spectrum 1983−  LTAS displays  Harmonic FFT displays  Waterfall displays 1983−  F0 displays 1983−  F0 histograms with semi scale 1989−  Jitter/shimmer 1990−  Computer Voice Fields with Hz scale 1983−  FFT spectrogram displays  LPC spectrogram displays  Auditory spectrogram displays  Bark, ERB, mel, st ja Hz scales formant charts 1988−  Bark scale 1983−  ERB scale 1996−  Semi scale 1989−  Sone scale 1995− Works 1972−  Other links  Maailman kielten Bark-asteikkoiset formanttikartat

ISA software

Kepstrit lin/Hz-asteikolla
Intelligent Analyser™ (ISA) -ohjelmassa

|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. Signaali x(n) = xa(nT), n=0...N-1, T on näytejakso, xa= analogiasignaali. Jos xa on mikrofonin jännitesignaali, niin xa(nT) = kPa(nT), jossa Pa on paine ja k on vakiokerroin.

Hamming-ikkuna w(n) = 0.54-0.46cos(2πn/(N-1)).

Blackman-Harris-ikkuna w(n) = 0.35875-0.48829cos(2πn/(N-1))+0.14128cos(4πn/(N-1))-0.01168cos(6πn/(N-1)).

Laskenta käyttää decimation-in-time (DIT) algoritmia.
Olen vuosien varrella DSP-miehenä koodannut kepstrin (kaksoisDFT) täysin nollasta
(1)
Texasin 16-bittiselle TMS320-signaaliprosessoriperheelle konekielellä,
(2)
Motorolan 16- ja 32-bittisille M68000-mikroprosessoriperheille konekielellä ja C-kielellä,
(3) IBM:n 600-sarjan 32-bittiselle PowerPC-mikroprosessoriperheelle konekielellä ja C-kielellä,
(4)
Intelin 32- ja 64-bittisille mikroprosessoriperheille C++-kielellä.

Signaalin ikkunointi tapahtuu Hamming-ikkunalla, suorakaide-ikkunalla tai yli 90 dB:n Blackman-Harris-ikkunalla.


Analyysikuvat olen alusta saakka koodannut älykkäiksi Neon olio-ohjelmointikielellä.

Kepstri lin/Hz-asteikolla.

YIN lin/Hz-asteikolla.



FFT dB/Hz/Bark-asteikolla
.