importare i dati di segnale , o generare un segnale di prova a scopo dimostrativo
t = 0:1000 ; segnale . = Cos ( 2 * pi * 30 * t) + 0.04 * sin ( 2 * pi * 60 * t ) + 0,01 * cos ( 2 * pi * 90 * t ) ;
2
Creare uno spettro utilizzando il dspdata.msspectrum () la funzione < br .
periodogram_obj > = spectrum.periodogram ; mass_spec = dspdata.msspectrum ( periodogram_obj , segnale , ' F ' , 1000, ' NFFT ' , lunghezza ( segnale) ) ;
3
estratto l'ampiezza dello spettro alle frequenze armoniche
my_data = mass_spec.Data ( 31,61,91 ) .
4
Calcola il THD , la somma delle grandezze meno il segnale di ingresso diviso per il segnale di ingresso
THD = sqrt ( sum ( my_data ( 2 : end ) ) ) /sqrt ( my_data ( 1) ) ;
5
moltiplicare per 100 . per esprimere THD o THD + N , in percentuale
THD_pct = 100 * THD ; .
Programmazione © www.354353.com