1/3 octave spectra analysis of sustained phonation
Nasality is described using the parameters Fn_mean, Fn_std and Fn_trend. These parameters are evaluated by the function 'nasality.m' which is designed for prolonged vowels /i/ and may be called as:
[x,Fs]=audioread('prolonged_vowel_i.wav'); [Fn_mean, Fn_std, Fn_trend]=nasality(x,Fs)
Signalprocessing Toolbox Statistics and Machine Learning Toolbox
This algorithm is based on 1/3 the octave filter bank designed by Christophe Couvreur. Avalilable at: https://www.mathworks.com/matlabcentral/fileexchange/69-octave/content/octave/oct3bank.m
If you find this algorithm useful and wish to reference it in your publications, you can make a reference directly to this article:
M. Novotny, J. Rusz, R. Cmejla, H. Ruzickova, J. Klempír, E. Ruzicka. (2016). "Hypernasality associated with basal ganglia dysfunction: evidence from Parkinson’s disease and Huntington’s disease, PeerJ, 4: 2530.