Por ejemplo, tengo densidad de espectro de potencia calculada por pwelch
f = EEG.data;
f = f(2, :);
N = EEG.pnts;
Fs = EEG.srate;
NFFT = 2^nextpow2(N);
NOVERLAP = 0;
WINDOW = 512;
[spectra, freqs] = pwelch(x, WINDOW, NOVERLAP, NFFT,Fs);
deltaIdx = find(freqs>1 & freqs<4);
deltaPower = mean(spectra(deltaIdx));
El código proviene de EEGLab y se puede encontrar aquí: https://sccn.ucsd.edu/pipermail/eeglablist/2014/008043.html
¿Cómo debo calcular el poder absoluto de él? ¿Cuál es la relación entre la potencia absoluta y la densidad del espectro de potencia?
Tuve que buscar un poco para poder responder a esta pregunta. De hecho, una vez se hicieron preguntas similares sobre el procesamiento de señales .
Para hacer la explicación más sencilla, creo que la respuesta significa lo siguiente:
Con pwelch o un análisis FFT, puede calcular la amplitud de las sinusoides con frecuencias particulares (consulte esta respuesta ). Estas amplitudes al cuadrado dan como resultado la potencia absoluta dentro de estas frecuencias específicas. La potencia resultante por frecuencia es la densidad espectral de potencia (PSD). En neurociencia, las personas no suelen trabajar con frecuencias individuales, sino con bandas de frecuencia, como la banda alfa (8-14 Hz). Esto se puede calcular tomando el promedio o la integral.
La potencia absoluta (W), a la que se refirió, es la potencia de toda la señal. Por lo tanto, no hace una distinción entre las diferentes frecuencias. Esto se puede calcular sumando la potencia de cada frecuencia (es decir, tomando la integral de la señal). Al sumar, tiene la cantidad total de potencia dentro de la señal. La potencia absoluta se puede utilizar para normalizar la PSD, dividiendo la PSD por la potencia absoluta (como se describe en la respuesta sobre procesamiento de señales).
Promedio vs integral
Parece que es posible calcular el promedio o la integral de una banda de frecuencia. Pero, ¿cuándo tomarías la integral y cuándo tomarías el promedio? El promedio es una manera fácil de calcular qué tan alta es la potencia en una banda de frecuencia específica. De esta manera, no importa cuán grande sea el rango de una banda de frecuencia. El promedio de la banda alfa no estará más inflado que la banda delta (0-4 Hz) debido a un rango mayor.
Cuando tomas la integral, existirán diferencias que se ven afectadas por el tamaño de la banda de frecuencia, porque no divides la suma por la cantidad de muestras (es decir, frecuencias). Sin embargo, le permite normalizar el PSD con el poder absoluto. Digamos que la integral de la banda alfa da como resultado 2W
, la integral de la banda delta da como resultado 1W
y la potencia absoluta de todo el espectro es igual a 10W
. Entonces sabemos que la banda alfa juega una parte del 20% y la banda delta una parte del 10% en la señal total. Con un mero promedio, esto no sería posible porque 'despreciamos' el tamaño del rango.
Creo que depende de lo que busques exactamente. No he trabajado con el poder absoluto y no estoy familiarizado con su usuario. Quizás tenga algo que ver con las comparaciones entre exp. condiciones (dentro de bandas de frecuencia, es decir, alfa frente a alfa) o dentro de exp. condiciones (entre bandas de frecuencia). Si un experto real pudiera explicar por qué un método es más beneficioso que el otro, sería fantástico.
honi
Robin Kramer-diez tienen
federica