Этот пример показывает явным образом, как вычислить общее гармоническое искажение в дБн для сигнала, состоящего из основного принципа и двух гармоник. Явное вычисление проверяется по результату, возвращенному thd.

Создайте сигнал, выбранный на уровне 1 кГц. Сигнал состоит из основного принципа на 100 Гц с амплитудой 2 и две гармоники на уровне 200 и 300 Гц с амплитудами 0.01 и 0.005. Получите общее гармоническое искажение явным образом и использование thd.

t = 0:0.001:1-0.001;
x = 2*cos(2*pi*100*t)+0.01*cos(2*pi*200*t)+0.005*cos(2*pi*300*t);
tharmdist = 10*log10((0.01^2+0.005^2)/2^2)

tharmdist = -45.0515

r = thd(x)

r = -45.0515

Создайте сигнал, выбранный на уровне 1 кГц. Сигнал состоит из основного принципа на 100 Гц с амплитудой 2 и три гармоники в 200, 300, и 400 Гц с амплитудами 0.01, 0.005, и 0.0025.

Определите номер гармоник к 3. Это включает основной принцип. Соответственно, степень в 100, 200, и 300 Гц используется в вычислении THD.

t = 0:0.001:1-0.001;
x = 2*cos(2*pi*100*t)+0.01*cos(2*pi*200*t)+ ...
    0.005*cos(2*pi*300*t)+0.0025*sin(2*pi*400*t);
r = thd(x,1000,3)

r = -45.0515

Определение количества гармоник, равных 3, игнорирует степень на уровне 400 Гц в вычислении THD.

Создайте сигнал, выбранный на уровне 1 кГц. Сигнал состоит из основного принципа на 100 Гц с амплитудой 2 и три гармоники в 200, 300, и 400 Гц с амплитудами 0.01, 0.005, и 0.0025.

Получите периодограмму оценка PSD сигнала и используйте оценку PSD в качестве входа к thd. Определите номер гармоник к 3. Это включает основной принцип. Соответственно, степень в 100, 200, и 300 Гц используется в вычислении THD.

t = 0:0.001:1-0.001;
fs = 1000;
x = 2*cos(2*pi*100*t)+0.01*cos(2*pi*200*t)+ ...
    0.005*cos(2*pi*300*t)+0.0025*sin(2*pi*400*t);
[pxx,f] = periodogram(x,rectwin(length(x)),length(x),fs);
r = thd(pxx,f,3,'psd')

r = -45.0515

Определите THD путем введения спектра мощности, полученного с Окном Хэмминга и пропускной способностью разрешения окна.

Создайте сигнал, выбранный на уровне 10 кГц. Сигнал состоит из основного принципа на 100 Гц с амплитудой 2 и три нечетных гармоники в 300, 500, и 700 Гц с амплитудами 0.01, 0.005, и 0.0025. Задайте количество гармоник к 7. Определите THD.

fs = 10000;
t = 0:1/fs:1-1/fs;
x = 2*cos(2*pi*100*t)+0.01*cos(2*pi*300*t)+ ...
    0.005*cos(2*pi*500*t)+0.0025*sin(2*pi*700*t);
[sxx,f] = periodogram(x,hamming(length(x)),length(x),fs,'power');
rbw = enbw(hamming(length(x)),fs);
r = thd(sxx,f,rbw,7,'power')

r = -44.8396

Сгенерируйте сигнал, который напоминает вывод слабо нелинейного усилителя с тоном на 2,1 кГц, как введено. Сигнал выбирается в течение 1 секунды на уровне 10 кГц. Вычислите и постройте спектр мощности сигнала. Используйте окно Kaiser с β = 38 для вычисления.

Fs = 10000;
f = 2100;

t = 0:1/Fs:1; 
x = tanh(sin(2*pi*f*t)+0.1) + 0.001*randn(1,length(t));

periodogram(x,kaiser(length(x),38),[],Fs,'power')

Гармоники терпят от шума на частотах 4,2 кГц, 6,3 кГц, 8,4 кГц, 10,5 кГц, 12,6 кГц и 14,7 кГц. Все частоты за исключением первой больше, чем частота Найквиста. Гармоники искажаются соответственно в 3,7 кГц, 1,6 кГц, 0,5 кГц, 2,6 кГц и 4,7 кГц.

Вычислите общее гармоническое искажение сигнала. По умолчанию thd обрабатывает искаженные гармоники как часть шума.

thd(x,Fs,7);

Повторите вычисление, но теперь обработайте искаженные гармоники как часть сигнала.

thd(x,Fs,7,'aliased');

Создайте сигнал, выбранный на уровне 10 кГц. Сигнал состоит из основного принципа на 100 Гц с амплитудой 2 и три нечетных гармоники в 300, 500, и 700 Гц с амплитудами 0.01, 0.005, и 0.0025. Задайте количество гармоник к 7. Определите THD, степень в гармониках и соответствующие частоты.

fs = 10000;
t = 0:1/fs:1-1/fs;
x = 2*cos(2*pi*100*t)+0.01*cos(2*pi*300*t)+ ...
    0.005*cos(2*pi*500*t)+0.0025*sin(2*pi*700*t);
[r,harmpow,harmfreq] = thd(x,10000,7);
[harmfreq harmpow]

ans = 7×2

  100.0000    3.0103
  201.0000 -321.0983
  300.0000  -43.0103
  399.0000 -281.9259
  500.0000  -49.0309
  599.0000 -282.1066
  700.0000  -55.0515

Степени в четных гармониках находятся на порядке $$-\mathrm{}\mathrm{}300$$ дБ, который соответствует амплитуде $$10^{-\mathrm{}15}$$.

Сгенерируйте синусоиду частоты 2,5 кГц, выбранные на уровне 50 кГц. Добавьте Гауссов белый шум со стандартным отклонением 0.00005 к сигналу. Передайте результат через слабо нелинейный усилитель. Постройте THD.

fs = 5e4;
f0 = 2.5e3;
N = 1024;
t = (0:N-1)/fs;

ct = cos(2*pi*f0*t);
cd = ct + 0.00005*randn(size(ct));

amp = [1e-5 5e-6 -1e-3 6e-5 1 25e-3];
sgn = polyval(amp,cd);
thd(sgn,fs);

График показывает, что спектр раньше вычислял отношение и область, обработанную как шум. Уровень DC исключен из вычисления. Основной принцип и гармоники маркированы.