Acoustic Time-frequency analysis

This function generates a waterfall plot in dB of an audiofile

function [time_freq_mat,frequency,sound_pressure_level,time_vector]=acoustic_Wf(pressure_signal,Fs,window_name,window_length,window_overlap,plot_curve,varargin)
% This function perform an acoustic time-frequency analysis
%
% SYNTHAXIS :
% [time_freq_mat,frequncy_sound_level,sound_level]=acoustic_Wf(signal,Fs,window_name,window_length,window_overlap,plot_curve)
% [time_freq_mat,frequncy_sound_level,sound_level]=acoustic_Wf(signal,Fs,window_name,window_length,window_overlap,plot_curve,weigthing)
% [time_freq_mat,frequncy_sound_level,sound_level]=acoustic_Wf(signal,Fs,window_name,window_length,window_overlap,plot_curve,weigthing,octave)
%
% INPUT :
% pressure_signal : N-length vector : vector time domain signal (Best practice
% supposes that the N-length vector is a power of 2)
% Fs : double : sampling frequency
% window_name : string : should be one of the above list
% rectangular : rectangular window
% triangular : triangular with L parameter set to N
% bartlett : Bartlett window
% hanning : hanning window
% hamming : hamming window
% blackmann : blackmann window
% flattop : flattop window
% window_length : double : length of the window in seconds
% window_overlap : double : overlap between twoe windows in seconds
% plot_curve : logical : 1 to plot the checking figure; 0 don't
% weighting : string : apply the standard acoustic weighting 'dBa', 'dBb'
%
%
% OUTPUT :
% time_freq_mat :N x ceil(0.5*N-length) matrix : time-frequency sound
% pressure level matrix
% sound_pressure_level : ceil(0.5*N-length) vector : vector of sound
% pressure level
% frequency : ceil(0.5*N-length) vector : window shape for checking purposes
%
% Autor : fganny
% Version : 1.0