Расчет мощности сигнала по спектру

[alex_m]

Пытаюсь прикрутить к SDR приемнику S-метр.
Сделал измерение пиковой мощности сигнала путем поиска бина с максимальным сигналом.
Магнитуду бинов перевожу в децибеллы по формуле db = 20 * Log10(magnitude);
К результату добавляю константу, чтобы получить уровень в dBm.
Из dBm уже считаю S уровень. Вроде работает.

Но возникло два вопроса:

1) Правильно ли я считаю пиковую мощность, выбирая уровень одного бина с максимальным уровнем?
Можно ли это считать PEP мощностью?

2) Как правильно посчитать RMS мощность в заданной полосе?

[EU1SW]

1. нет
2. Мощность в полосе вычисляется суммой всех бинов FFT, попадающих в заданную полосу. Мощность - это площадь фигуры, тут все просто.
Поэтому для вычисления показаний S-метра нужно брать 20 логарифмов от суммы магнитуд FFT после ФОС

[alex_m]

суммой всех бинов FFT

пробовал суммировать квадраты мощностей всех бинов, но получилась ерунда...

[EU1SW]

не, сначала просуммировать, а потом "взять", я там подумал и дописал в пред. пост, но опоздал )

[Genadi Zawidowski]

суммировать квадраты мощностей

Не первыполняйте план, не надо квадраты мощностей, достаточно квадраты амплитуд.

[EU1SW]

если суммировать квадраты амплитуд то нужно брать 10 логарифмов, если просто суммировать амплитуды - 20 логарифмов.
результат одинаковый, но отпадает необходимость брать дофига квадратов )

[alex_m]

попробовал через амплитуду, вроде работает
Уровень как пиковый, так и RMS меняется соответственно изменению аттенюатора на выходе NWT

[EU1SW]

если вы имеете ввиду под "пиковым" магнитуду одной палки FFT - то это не правильно

[alex_m]

результат одинаковый, но отпадает необходимость брать дофига квадратов )

убрать квадраты не выходит, а вот если складывать квадраты и потом 10*Log10(sum), результат соответствует.
Проверял от 0 до 50 дБ с шагом 10 дБ - все вроде ок, если не считать небольших флуктуаций

если вы имеете ввиду под "пиковым" магнитуду одной палки FFT - то это не правильно

вот и у меня сомнения в правильности такого подхода, а как правильно?

[EU1SW]

если формально, то PEP вычислить можно выполнив обратное FFT и найдя максимум амплитуды во временном домене...
но что то мне подсказывает, что просуммировав все бины FFT в частотном домене мы в результате получим то же самое...
надо вспоминать... )))

UPD: для RMS явно где то должно быть усреднение по времени