Cómo calcular la THD en MATLAB

distorsión armónica total o THD, es una medida de la cantidad de distorsión está presente en una señal debido a las oscilaciones de los armónicos de la señal de entrada. THD y THD + N - una variante incluido el ruido, son mediciones comunes para evaluar la calidad de audio . MATLAB es un paquete de programación matemática que puede ser utilizado para calcular THD con sólo unas pocas líneas de código . Instrucciones
1

importar los datos de la señal, o generar una señal de prueba para fines de demostración

t = 0:1000 ; señal . = Cos ( 2 * pi * 30 * t) + 0.04 * sin ( 2 * pi * 60 * t) + 0.01 * cos ( 2 * pi * 90 * t);
2

Crear un espectro utilizando el dspdata.msspectrum () función < br . >

periodogram_obj = spectrum.periodogram ; mass_spec = dspdata.msspectrum ( periodogram_obj , señal , ' Fs ' , 1000, ' Nfft ' , longitud ( señal) ) ;
3

Extract la magnitud del espectro en las frecuencias armónicas

my_data = mass_spec.Data ( 31,61,91 ) .
4

calcular el THD, la suma de las magnitudes menos el señal de entrada dividida por la señal de entrada

THD = sqrt (sum ( my_data ( 2 : fin ))) /sqrt ( my_data ( 1 ));
5

Multiplicar por 100 . para expresar THD, o THD + N, como porcentaje

THD_pct = 100 * THD ; .