wolffd@0: function varargout = mirkurtosis(orig,varargin)
wolffd@0: %   k = mirkurtosis(x) calculates the kurtosis of x, indicating whether
wolffd@0: %       the curve is peaked or flat relative to a normal distribution.
wolffd@0: %   x can be either:
wolffd@0: %       - a spectrum (spectral kurtosis),
wolffd@0: %       - an envelope (temporal kurtosis), or
wolffd@0: %       - any histogram.
wolffd@0: 
wolffd@0: 
wolffd@0: varargout = mirfunction(@mirkurtosis,orig,varargin,nargout,struct,@init,@main);
wolffd@0: 
wolffd@0: 
wolffd@0: function [x type] = init(x,option)
wolffd@0: if not(isamir(x,'mirdata')) || isamir(x,'miraudio')
wolffd@0:     x = mirspectrum(x);
wolffd@0: end
wolffd@0: type = 'mirscalar';
wolffd@0: 
wolffd@0: 
wolffd@0: function k = main(x,option,postoption)
wolffd@0: if iscell(x)
wolffd@0:     x = x{1};
wolffd@0: end
wolffd@0: y = peaksegments(@kurtosis,get(x,'Data'),...
wolffd@0:                            get(x,'Pos'),...
wolffd@0:                            get(mircentroid(x),'Data'),...
wolffd@0:                            get(mirspread(x),'Data'));
wolffd@0: if isa(x,'mirspectrum')
wolffd@0:     t = 'Spectral kurtosis';
wolffd@0: elseif isa(x,'mirenvelope')
wolffd@0:     t = 'Temporal kurtosis';
wolffd@0: else
wolffd@0:     t = ['Kurtosis of ',get(x,'Title')];
wolffd@0: end
wolffd@0: k = mirscalar(x,'Data',y,'Title',t,'Unit','');
wolffd@0: 
wolffd@0: 
wolffd@0: function k = kurtosis(d,p,c,s)
wolffd@0: k = sum((p-c).^4.*d) ./ sum(d) ./ s.^4;