Nothing
# This code is based on the Matlab implementations of PLP and Rasta
# feature calculations by Daniel P. W. Ellis of Columbia University /
# International Computer Science Institute. For more details, see:
# http://www.ee.columbia.edu/~dpwe/resources/matlab/rastamat/
hz2mel <- function(f, htk=FALSE){
if(!is.numeric(f) || f < 0)
stop("frequencies have to be non-negative")
if(htk){
z <- 2595 * log10(1 + f/700)
} else {
# Mel calculation like in Slaney's Auditory Toolbox
f_0 <- 0
f_sp <- 200/3
brkfrq <- 1000
# Starting Mel value for log region
brkpt <- (brkfrq - f_0)/f_sp
# The magic 1.0711703 which is the ratio needed to get from 1000 Hz to
# 6400 Hz in 27 steps.
logstep <- exp( log(6.4)/27 )
linpts <- (f < brkfrq)
z <- 0 * f
z[linpts] <- (f[linpts] - f_0)/f_sp
z[!linpts] <- brkpt + (log(f[!linpts]/brkfrq))/log(logstep)
}
return(z)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.