Consultas

Hola Maila,
Para hallar la f0 con fpeaks es necesario usar un objeto de tipo spec. No hay argumento "at" en esta función. Pero en la clase utilizamos el argumento nmax que nos da una cierta cantidad máxima de picos, que además son los de mayor amplitud.
Por otra parte, si miras en la ayuda de la función "fund", vas a ver que el único argumento necesario es el archivo .wav
Por lo que el primer argumento no puede ser fspec, que es en si un objeto creado por nosotros en R y si haces class(fspec) vas a ver que te dice "Matrix" "array".
Marcar la posición no es un argumento obligatorio, pero en este caso lo hacemos con at=center para seguir obteniendo los datos de la misma porción de audio en la que venimos trabajando.

Es muy importante que usen la ayuda de R, porque ahí pueden ver los argumentos necesarios, la clase a la que pertenecen, cuales son los argumentos opcionales etc.

Espero que esta respuesta te ayude en la tarea. Saludos,

lucia

Hola Maila y Franco. Me alegro que aprovechen el foro para preguntar e intercambiar entre ustedes. A continuación algunos comentarios.

a) Como dice Franco, usar localpeaks es una opción. Pero también podrían usar simplemente fpeaks, que es una función un poco más sencilla. El atributo "at" se usa al calcular el espectro, con la función "spec", e indica el punto en el tiempo en el que se calcula el espectro. Para controlar la cantidad de puntos detectados existe el atributo "nmax", de la función fpeaks. No se si te sirve la aclaración. Hay que explorar los parámetros que recibe cada función, con "help(nombre_de_función)". Y además les recomiendo ver los ejemplos de las transparencias y del libro.

En resumen, el procedimiento para calcular los picos del espectro es el siguiente:
1) Cargar un archivo de audio
Con la función "data" si es un archivo de la biblioteca, o con "readWave" si es un archivo .wav
2) Calcular el espectro
Con la función "spec". Pueden ajustar varios parámetros, como el largo de la ventana y el instante de tiempo donde se ubica esa ventana.
3) Detectar los picos
Con la función "fpeaks" o con "localpeaks", pero con "fpeaks" debería ser suficiente para el ejercicio.
También pueden ajustar varios parámetros para controlar la cantidad de picos detectados (vean las transparencias y el libro).
Tengan presente que esta función recibe como parámetro de entrada el espectro calculado en la parte anterior.

b) Como bien dice Franco hay un error en los argumentos que le estás pasando a la función "fund". Fijate que recibe un objeto del tipo wave (el que obtenemos en el punto 1 de mi anterior respuesta). Para ver lo que deberías pasar como parámetro te sugiero usar la ayuda de R, "help(nombre_de_función)".

En este caso dice lo siguiente (copio solo lo relevante para los parámetros que te generan dudas):
fund(wave, f, ..., fmax = f/2, ...)
wave - an R object.
f - sampling frequency of wave (in Hz). Does not need to be specified if embedded in wave.
fmax - the maximum frequency to detect (in Hz).

Ahí se puede ver que el primer argumento es un objeto R y si bien no lo especifica tiene que ser un "wave" como indica su nombre. Luego el argumento "f" corresponde a la frecuencia de muestreo. Verás que aclara que no es necesario especificarlo si está indicado en el archivo de audio, lo cual es lo habitual. Por último, "fmax" es la frecuencia máxima a considerar como candidato para la detección de f0. Es decir, la función va a devolver una frecuencia fundamental por debajo de este límite. Esto tiene sentido para que el cálculo de f0 sea más rápido y para evitar posibles errores, pero no es necesario especificarlo en todos los casos. Por defecto se usa el valor fmax=f/2, es decir la mitad de la frecuencia de muestreo, por lo que se considera todo el espectro. Recuerden que una señal digital muestreada a una frecuencia de muestreo f puede tener componentes de frecuencia hasta la frecuencia de Nyquist, es decir f/2.

Espero que haya aclarado.

Saludos