Есть посуточный ряд мощности некого сигнала в db
Нужно оценить средний ход за месяц
Такой способ не кажется мне верным
select month(tm), avg(db) from alias() group by 1
Дает очень мелкие вариации, в то время как
график исходного ряда по суткам довольно
волосатый.
Попробовал усреднить исходные данные,
ничего не выходит
Дб там такие ~-175 - ~165
В фоксе ни при каком set decimal
не считаются
Пример верха такого ряда
Левая колонка временной лаг, далее в колонках мощности по частотам
Пытался прицепить картинку, размер не пускает
1996-03-22, -73, -73, -74, -87, -92, -105, -114, -27
1996-03-23, -174, -172, -74, -130, -99, -115, -114, -27
1996-03-24, -175, -174, -171, -124, -101, -115, -127, -32
1996-03-25, -178, -177, -178, -133, -104, -116, -133, -33
1996-03-26, -176, -174, -165, -135, -105, -114, -125, -30
1996-03-27, -176, -177, -169, -137, -105, -116, -126, -32
1996-03-28, -172, -169, -161, -132, -105, -114, -120, -28
1996-03-29, -172, -174, -168, -136, -101, -113, -122, -31
1996-03-30, -177, -178, -179, -126, -103, -113, -124, -28
1996-03-31, -174, -170, -175, -113, -101, -112, -123, -32
1996-04-01, -175, -174, -175, -126, -108, -116, -130, -32
1996-04-02, -175, -173, -168, -132, -109, -116, -126, -32
1996-04-03, -177, -177, -179, -121, -108, -117, -128, -33
1996-04-04, -176, -177, -176, -116, -105, -116, -127, -33
1996-04-05, -175, -176, -176, -134, -109, -117, -130, -34
1996-04-06, -177, -177, -177, -132, -111, -123, -125, -35
1996-04-07, -178, -178, -179, -119, -111, -124, -139, -34
1996-04-08, -178, -179, -174, -131, -114, -122, -128, -34
1996-04-09, -176, -178, -179, -132, -115, -121, -137, -35
1996-04-10, -176, -175, -168, -137, -108, -118, -128, -33
1996-04-11, -177, -175, -166, -135, -107, -117, -127, -33
1996-04-12, -176, -175, -176, -137, -108, -117, -135, -33



Исправлено 2 раз(а). Последнее : boba, 10.10.19 17:01

Между логарифмическими величинами считать среднее арифметическое несколько неправильно
Тут кой-какие методики расчёта приведены.
transform.ru

Может быть надо сравнивать среднеквадратичное отклонение? Или вообще все три показателя, среднее и дисперсия +/-.

А стабильное изменение мощности покажет, скорее, медиана за одинаковый период. Но это только, если продолжительность процессов сопоставима по времени. То-есть, если стабильное нарастание и падение мощности продолжается несколько месяцев, а мы хотим сравнить разные месяцы, то медиана что-то должна показать.

Игорь, спасибо. Я ведь так и пытался вычислить, это в посте написано
Именно вторую формулу и нельзя вычислить по тем значениям,
что в данных.
См вторую формулу из приведенной вами ссылки.
Но клипбордом она суда попадает искаженной--

Lср.= 10 1g (100.1L1+100.1L2+100.1L3+...+100.1Ln) – 10 1gn, дБА (2)
Они выпадают за предел точности , нечего суммировать.
Наверное придется двинуть эти дб до чисел около нуля,
провести все расчеты, а потом вернуть назад порядок величины.
Программирование совершенно атрофирует мозги,
забываешь не только институт , но и школьную математику.
Сейчас не соображу, верно это или нет.

Немного лирики невпопад.
Я работал в НИИ где занимались ядерными исследованиями.
И там статистика - это все)
Потому, среднеквадратичные, и пр. отклонения, просто как пирожки на языке.

Но всегда крайне важен, физический смысл. Именно из него исходя, интерпретируются результаты.
Результат нужен для чего то.
Именно исходя из этого, строится интерпретация.

Можно говорить о теории, можно об эксперименте.
Но результат подгоняется всегда под что то. Потом теория пересматривается, и результат подгоняется оп новой.

В общем, надо понимать, что хочешь получить. Потом искать способы получения. Тут думаю, в инете, полно разных подходов, на любой вкус.

Владимир, там же вроде как алгоритм предложен где не надо вычислять все эти дикие 10^... А просто разности дБ по табличке брать, складывать (по идее числа будут очень скромные, ну разве что измерений 100500) потом вычитать опять же по таблице взятый логарифм от "числа измерений" - ну т.е. этих самых N

Игорь
Да, надо так
10 в степени среднего алгоритма выносится за сумму
10 только в разностных дб
Посмотрю, хватает ли точности фокса на разностные
степени.

boba
Программирование совершенно атрофирует мозги,

Правда? Бросить что ли программировать, пока совсем дауном не стал...

Я тоже не понял, в чем "секрет"... -174 дБ
и так тоже 0:
В хелпе пишут про Numreric диапазон: - .9999999999E+19 to .9999999999E+20
Отрицательные степени (числа ок. нуля) входят в диапазон, но их нельзя указать? Ну и засада...

Доб. Double:
Числа двойной точности с плавающей запятой обеспечивают точность в 15—17 десятичных цифр и масштабы в диапазоне примерно от 10−308 до 10+308. Под разрядность 10 разрядов (1024), под мантиссу 52 разряда, делим на пи (~LOG2(10)) = 28 десятичных разрядов! (почему-то у фокса 15). Double должно работать и в фоксе!
ru.wikipedia.org


Доб. Интересно, если каждый сигнал (отсчет) увеличить на произвольное кол дБ (прибавить громкость), потом вычислить по обычной формуле (развернув из логарифма в обычный уровень) среднее, или ср.квадратичное, или любое другое среднее, и потом опять уменьшить громкость результата (перевести в логарифм и вычесть это же кол дБ), то результат будет правильный?
...Скорее да, т.к. это равносильно все составляющие уровни в сумме (10^дб) умножить на одно число, просуммировать, потом на него же поделить, среднее должно остаться правильным. Про ср.квадратичное так может не сработать - надо подумать... (ср.кв., и среднее тоже - это все интеграл от уровня, или уровня^2, причем уровень задан в виде логарифма от истиного уровня...)

Для себя, кепстр, встретил давно чуть не в институтском возрасте, операции с логарифмами, энергетика излучения:
ru.wikipedia.org



Исправлено 5 раз(а). Последнее : of63, 12.10.19 21:53

Пивнул?

Вливайся в формулу, трезвеник х.

of63
Под разрядность 10 разрядов (1024), под мантиссу 52 разряда, делим на пи (~LOG2(10)) = 28 десятичных разрядов! (почему-то у фокса 15).

Мантисса на самом деле представляет 53 разряда точности см. про "неявную единицу" для нормализованных чисел.
И это позволяет хранить чуть меньше 16 десятичных чисел "точности". Хреново ты логарифмы считаешь Понятия не имею зачем там Пи приплёл...

> Понятия не имею зачем там Пи приплёл...
Это юмор, ЛОГ2(10 есть примерно пи, 2 в в 3й степени - это 10...

> Мантисса на самом деле представляет 53 разряда точности см. про "неявную единицу" для нормализованных чисел.
Не понял, 52 разряда - это сколько десятичных разрядов точности?

53 разряда это 15.95 десятичных знака. Т.е. 15 всегда, 16-й иногда сбивается - ну т.е. два разных десятичных числа различающиеся на единицу в 16 знаке представляются одним и тем же бинарным double.
Естественно, при проведении математических вычислений точность будет падать - в зависимости от разницы в "масштабе" чисел участвующих в операции.
К примеру 1e10 и 1е-10 сами по себе хранятся довольно таки "точно", но если их сложить или отнять, то меньшее потеряется полностью.

Вопрос простой, с какого ху фокс не видит величину (10^(-17))? это так заложено в фоксе?

Если че, психологическое вдрук, то вопрос не к тебе, Игорь, просто спросил, сам себя )



Исправлено 1 раз(а). Последнее : of63, 13.10.19 00:19

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

Аспид
Немного лирики невпопад.
Я работал в НИИ где занимались ядерными исследованиями.
...

С этого места по подробнее, как НИИ назывался, п/я Р-6575?

"Помедленнее, пожалуйста. Я записываю..."