Cómo programar en matlab para calcular la densidad media de energía eólica
Matlab implementa la estimación clásica del espectro de potencia
El resultado de FFT es un espectro y el resultado de PSD es un espectro de potencia; El espectro de potencia de señales con diferentes espectros es. Pueden ser iguales; el espectro de potencia es el cuadrado del módulo de amplitud y el resultado es un número real.
En Matlab, la densidad espectral de autopotencia puede ser calculado directamente usando la función psd Según matlab, psd puede lograr la estimación del método Welch, es decir, es equivalente a utilizar el método de periodograma promedio mejorado para estimar la densidad espectral de potencia de señales aleatorias. El resultado obtenido por psd debería ser más suave.
1. Método directo:
El método directo también se llama método de periodograma y considera los N datos de observación de la secuencia aleatoria x (n) como una secuencia con energía limitada. y Calcule directamente la transformada discreta de Fourier de x (n) para obtener X (k), luego tome el cuadrado de su amplitud y divídalo por N como una estimación del verdadero espectro de potencia de la secuencia x (n).
Ejemplo de código Matlab:
clear;
Fs=1000; Frecuencia de muestreo
n=0: 1/Fs: 1;
Generar una secuencia que contenga ruido
xn=cos(2*pi*40*n) 3*cos(2*pi*100*n) randn(size(n)) ;
ventana=boxcar(longitud(xn)); ventana rectangular
nfft=1024;
[Pxx, f]=periodograma(xn, ventana, nfft, Fs); método directo
plot(f, 10*log10(Pxx));