Алгоритм расчета THD Опции

[stepit]

Здравствуйте, реализовал следующий алгоритм расчета THD:

1. Снимаем сигнал y со звуковой карты;
2. Накладываем окно Хенинга на сигнал: yw = hanning(y);
3. Делаем FFT и ищем частоту основной гармоники: S = fft(yw), f = max(abs(S));
4. Берем значения спектра в гармониках: h1 = S(f), h2=S(2*f).. hn = S(n*f);
5. Рассчитываем THD = 100 * norm([h2 .. hn]) / norm(h1).

В этом алгоритме меня беспокоит точность определения амплитуды гармоники, и вобще низнаю как оценить точность работы алгоритма. Может кто знает другие алгоритмы расчета THD, тогда пожалуйста поделитесь.

[V_G]

Про зависимость THD от частоты в общем не слышал (он же по определению total!), но от частоты дискретизации - вполне вероятно, как и отношения с/ш. Есть такая программа SpecLab (от звуковухи работает), посмотрите на ней THD и с/ш при разных частотах дискретизации. При увеличении частоты дискретизации THD будет возрастать за счет добавления амплитуд шума на частотах "гармоник". Так что полосу анализа (и частоту дискретизации соответственно) надо выбирать ОЧЕНЬ осмысленно, чтобы не копить лишние шумы