Vitajte na [www.pocitac.win] Pripojiť k domovskej stránke Obľúbené stránky
Importovať dáta signálu , alebo generovať skúšobný signál pre demonštračné účely
t = 0:1000 ; . Signál = cos ( 2 * pi * 30 * t ) + 0,04 * sin ( 2 * pi * 60 * t ) + 0,01 * cos ( 2 * pi * 90 * t ) ; .
2
Vytvoriť spektrum pomocou dspdata.msspectrum funkcie ( )
periodogram_obj = spectrum.periodogram ; mass_spec = dspdata.msspectrum ( periodogram_obj , signál , " Fs " , 1000 , " NFFT " , dĺžka ( signál ) ) ;
3
Extract . veľkosť spektra na harmonických frekvenciách
my_data = mass_spec.Data ( 31,61,91 ) ;
4
Vypočítajte THD , súčet veľkostí mínus vstupný signál rozdelený podľa vstupného signálu
THD = sqrt ( sum ( my_data ( 2 : end ) ) ) /sqrt ( my_data ( 1 ) ) ;
5
Multiply 100. . . vyjadriť THD , THD + N alebo ako percento
THD_pct = 100 * THD ;