Освоил еще один способ, через преобразование Фурье. Самый быстрый и вероятно самой точный. Мы просто берем и исследуем период времени 2Пи. Подаем 2-х тоновый сигнал, с искажениями или без, на вход преобразования и по сути все. Частота отсчетов по теореме Котельникова должна быть минимум вдвое выше максимальной частоты сигнала (в т.ч. составляющей IMD). Далее, если 2 тона имеют частоты 999900 Гц и 1 Мгц, т.е. разница 100 Гц, просто вычисляем гармонику под номером 1000000, это будет один из исходных сигналов, потом 1000100, это будет составляющая IMD 3-го порядка, 1000200, это это будет составляющая IMD 5-го порядка и так далее. При этом совершенно не обязательно вычислять все на свете гармоники от нуля, достаточно вычислить только нужные, т.е. всего несколько штук. Это очень просто, быстро и качественно. Ну а далее по вышеприведенным формулам через логарифм (в формулах дробь можно перевернуть, что бы с минусом получалось) переводим в дб.