1. Harmonik Bileşenlerin Ölçülmesi: Sinyalin harmonik bileşenleri ölçülür ve bu bileşenlerin RMS değerleri alınır 12.
2. Karelerin Toplamı: Ölçülen harmonik bileşenlerin kareleri toplanır 23.
3. Maksimum Değerle Bölme: Toplanan karelerin toplamı, sinyalin temel bileşeninin RMS değerine bölünür 12.
4. Yüzde Olarak İfade: Elde edilen değer karekök işlemine tabi tutulur ve yüzde olarak ifade edilir (genellikle % işareti kullanılır) 23.

import numpy as np
import scipy.fftpack

# Örnek veri
sampled_data = np.loadtxt('veri.txt', delimiter=',')

# FFT (Fast Fourier Transform) ile harmonik bileşenleri elde etme
abs_yf = np.abs(scipy.fftpack.fft(sampled_data))

# THD hesaplama
def thd(abs_data):
    sq_sum = 0
    for r in range(len(abs_data)):
        sq_sum += (abs_data[r]**2)
    sq_harmonics = sq_sum - (max(abs_data)**2)
    thd = 100 * np.sqrt(sq_harmonics) / max(abs_data)
    return thd

# Sonuç
print("Total Harmonic Distortion: %s" % thd(abs_yf))  # Yüzde olarak sonuç verir
``` [2](https://stackoverflow.com/questions/42077247/calculating-thd-in-python).
Kopyala