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 ;