Вычисление мощности шума, Средний логарифм/логарифм среднего :) Опции

[Шаманъ]

Приветствую всех!

Вычисляю среднеквадратичную мощность шума и логарифмирую:

p = log(sum(x^2)/n),

тут все понятно

Если при вычислении мощности шума производить усреднение после логарифмирования, вот так:

p = sum(log(x^2)/n),

то результат будет смещенным, и потребуется коррекция.

Теперь вопрос. Если я вычисляю мощность шума блоками по 512отсчетов, потом логарифмирую и еще раз усредняю. Будет ли последнее усреднение вносить смещение и как оценить ошибку (если она будет)?

Заранее спасибо!

[jorikdima]

p = sum(log(x^2)/n)
Это что-то странное. Тоогда уж
p = sum(log(x^2))/n
чтоб осреднять. Но сумма логарифмов не равна логарифму суммы. Цель в чем?

[Шаманъ]

p = sum(log(x^2)/n)
Это что-то странное. Тоогда уж
p = sum(log(x^2))/n
чтоб осреднять.

Да, описка вышла.

Но сумма логарифмов не равна логарифму суммы. Цель в чем?

Я в курсе, что не равна, оттого и вопрос. Цель собственно в теме топика. Если посмотреть на вопрос немного шире, то в древних анализаторах спектра логарифмирование было до усреднения (видеофильтр) и потом выводили коррекцию для точного измерения мощности шума. Мне очень неудобно усреднять блоками более 512 отсчетов, вот и возникла такая идея.

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

Что, сильно плохая мысль ?

[mcheb]

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

Что, сильно плохая мысль ?

Допустим, один отсчёт равен 0. Тогда log(0^2) == ну в общем, идея понятна

[jorikdima]

идея понятна

Нет. Логарифм среднего не может быть нулем в оригинальном алгоритме?
Это обычная процедура по анализу числа на около-ноль перед логарифмированием.

Автор, Цель мне все-равно не ясна. В чем принципиальная разница?

[Шаманъ]

Автор, Цель мне все-равно не ясна. В чем принципиальная разница?

В смысле почему просто не усреднить бОльший блок, а потом взять логарифм? Это ограничения железа - мне придется переходить на вычисления с разрядностью 72 или 96 бит, получается очень непроизводительно. Есть правда одна идея - попробую сегодня.