реклама на сайте
подробности

 
 
> Операция обратная к "скользящему среднему"
_pv
сообщение Jan 6 2014, 12:36
Сообщение #1


Гуру
******

Группа: Свой
Сообщений: 2 563
Регистрация: 8-04-05
Из: Nsk
Пользователь №: 3 954



есть измеренные данные про которые известно, что по ним прошлись скользяшим средним с неким известным размером D.
каким образом можно восстановить исходные данные которые были до усреднения?
из простого на ум приходит только пройтись по этим уже усреднённым данным еще раз таким же скользящим средним и получить оценку ошибки вызванную усреднением, которую потом добавить обратно к первоначальным данным. Оно вроде как работает, но понятно что это хоть и очень простой, но не совсем честный способ.
Как это делать математически правильно?
понятно что можно сделать Фурье, поделить спектр на характеристику скользящего среднего (по сути КИХ фильтра с одинаковыми единичными коэффициентами) и преобразовать обратно. но при делении на нули в том месте где характеристика имеет провалы пожалуй будут проблемы, оно и понятно так как частоты кратные размеру фильтра давятся в ноль и обратному восстановлению не подлежат.
А если без преобразования в частотную область, с какой функцией надо сделать свёртку чтобы получить фильтр обратный скользящему среднему? или даже в общем случае, каким образом преобразовать коэффициенты КИХ фильтра, чтобы произведение исходного фильтра и пробразованного давало 1.
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
bookd
сообщение Jan 11 2014, 20:32
Сообщение #2


Частый гость
**

Группа: Свой
Сообщений: 151
Регистрация: 4-02-09
Из: Поволжье
Пользователь №: 44 403



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


--------------------
Всеобщая дебилизация не повод наносить ущерб своему здоровью.
Go to the top of the page
 
+Quote Post
Tarbal
сообщение Jan 13 2014, 12:11
Сообщение #3


Профессионал
*****

Группа: Свой
Сообщений: 1 351
Регистрация: 21-05-10
Пользователь №: 57 439



Цитата(bookd @ Jan 12 2014, 00:32) *
Ну не знаю, читаю уже четвертую страницу, реальных советов ноль. ТС решил проблему еще на первой странице, а все остальное просто хлам на мой взгляд. А так хочется увидеть качественное решение, с исходными данными, конкретным результатом и списком налагаемых ограничений...

Да ну. Я показал, что его решение можно использовать вопреки высказанному скептицизму и как произвести коррекцию ошибки измерения начального отсчета.

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


Сообщение отредактировал Tarbal - Jan 13 2014, 12:18
Go to the top of the page
 
+Quote Post
Stanislav
сообщение Jan 13 2014, 18:39
Сообщение #4


Гуру
******

Группа: Свой
Сообщений: 4 363
Регистрация: 13-05-05
Из: Москва
Пользователь №: 4 987



Понеслась опять нелёгкая, похоже...

Цитата(Tarbal @ Jan 13 2014, 16:11) *
Да ну. Я показал, что его решение можно использовать вопреки высказанному скептицизму и как произвести коррекцию ошибки измерения начального отсчета.

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

,

верно? biggrin.gif

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

Цитата(Tanya @ Jan 13 2014, 19:50) *
Удовлетворю Ваше любопытство - вот каждый день вижу хроматограф газовый.
..........
Вот вспомнилось. Первое (почти) электронное устройство, вышедшее из моих рук, было для хроматографа.
А вот и хымыки пожаловали.
Жалобу на вас, что ли, накатать?
Или предложить поискать площадку, где ваши таланты раскроются в полную силу?
Пока попробую второе.
Не соблаговолите ли сходить на ххымию, господа флудеры?
http://www.ximuk.ru/


--------------------
Самонадеянность слепа. Сомнения - спутник разума. (с)
Go to the top of the page
 
+Quote Post
AndrewN
сообщение Jan 13 2014, 20:25
Сообщение #5


Местный
***

Группа: Участник
Сообщений: 336
Регистрация: 7-03-07
Из: Петербург
Пользователь №: 25 961



QUOTE (Stanislav @ Jan 13 2014, 22:39) *
А вот и хымыки пожаловали
Аппаратная функция (как пример интегрального преобразования) есть у всех измерительных приборов. Чем плох хроматограф в качестве примера?
Go to the top of the page
 
+Quote Post
Stanislav
сообщение Jan 13 2014, 20:43
Сообщение #6


Гуру
******

Группа: Свой
Сообщений: 4 363
Регистрация: 13-05-05
Из: Москва
Пользователь №: 4 987



Цитата(AndrewN @ Jan 14 2014, 00:25) *
Аппаратная функция (как пример интегрального преобразования) есть у всех измерительных приборов. Чем плох хроматограф в качестве примера?
Примера чего именно?
И, самое интересное, какое этот пример имеет отношение к вопросу темы?

-------------------------

Уважаемые модераторы. По-моему, тема не стоит и выеденного яйца, поскольку, по сути, посвящена "доказательству" (ну, или "опровержению") совершенно банального тождества (см. выше). Ответ получен топикстартером на первой же странице.
Посему считаю, что её лучше закрыть, во избежание многостраничного флуда. Ну, или перенести в "курилку", на худой конец всё, что здесь написано после ответа на вопрос темы.


--------------------
Самонадеянность слепа. Сомнения - спутник разума. (с)
Go to the top of the page
 
+Quote Post
AndrewN
сообщение Jan 13 2014, 20:57
Сообщение #7


Местный
***

Группа: Участник
Сообщений: 336
Регистрация: 7-03-07
Из: Петербург
Пользователь №: 25 961



QUOTE (Stanislav @ Jan 14 2014, 00:43) *
Уважаемые модераторы. По-моему, тема не стоит и выеденного яйца, поскольку, по сути, посвящена "доказательству" (ну, или "опровержению") совершенно банального тождества (см. выше). Ответ получен топикстартером на первой же странице.
Посему считаю, что её лучше закрыть, во избежание многостраничного флуда. Ну, или перенести в "курилку", на худой конец всё, что здесь написано после ответа на вопрос темы.

Оки-доки. Вы высказались. А теперь не мешайте другим разговаривать. Люди сами разберутся, как, где и о чём им разговаривать.

Если вам не понятны теория, проблемы и области применений интегральных преобоазований - почитайте учебники.
Go to the top of the page
 
+Quote Post
Stanislav
сообщение Jan 13 2014, 21:23
Сообщение #8


Гуру
******

Группа: Свой
Сообщений: 4 363
Регистрация: 13-05-05
Из: Москва
Пользователь №: 4 987



AndrewN, на мои вопросы ответа нет, я правильно понял?

Хорошо, писать здесь больше не буду, если не попросят.
И опять же, против поговорки: "мели, Емеля, - твоя неделя" ничего особо не имею.
Только после этого называть Электроникс техническим форумом как-то язык не поворачивается...

А о теории, уважаемый AndreyN, у нас настолько разные представления, что могу вам посоветовать держать свои советы при себе.
Называть "интегрированием" фильтрацию скользящим средним можно только в понятийном пространстве, ортогональном к общепринятому.
Желаю удачи на сей ниве, тем не менее. sm.gif
Правила форума ещё бы не мешало почитать...


У всех посетителей темы прошу прощения за оффтоп.


--------------------
Самонадеянность слепа. Сомнения - спутник разума. (с)
Go to the top of the page
 
+Quote Post
rudy_b
сообщение Jan 14 2014, 14:26
Сообщение #9


Знающий
****

Группа: Свой
Сообщений: 888
Регистрация: 25-09-08
Из: Питер
Пользователь №: 40 458



Цитата(Stanislav @ Jan 14 2014, 01:23) *
...Называть "интегрированием" фильтрацию скользящим средним можно только в понятийном пространстве, ортогональном к общепринятому.

Скользящее среднее это именно интегрирование на заданном временном интервале. В самом общепринятом смысле.

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

Попробую немного по другому выразить уже сказанное раньше. Разность текущей и предыдущей сумм S[i] - S[i-1] равна разности текущего отсчета Х[i] и отсчета N тактов ранее X[i-N], где N - длина суммирования.
Т.е. X[i]=S[i]-S[i-1]+X[i-N].
Если мы теперь откатим на N точек назад и сделаем то же самое, то увидим, что разность S[i-N] - S[i-N-1]=X[i-N]-X[i-2N], т.е. X[i-N]=S[i-N]-S[i-N-1]+X[i-2N]. Подставив в предыдущее выражение получим
X[i]=S[i]-S[i-1]+X[i-N]=S[i]-S[i-1] + S[i-N]-S[i-N-1] + X[i-2N].
Т.е. если мы будем суммировать разности сумм (S[i]-S[i-1]) с шагом N (такую сумму разностей назовем SS), то сможем вычислить текущее значение X[i] зная значение в точке k*N (k-целое) тактов назад
X[i]=SS + X[i-k*N].
Если провести подобную операцию с N последовательными точками, т.е. вести N сумм SS (SS[m], где m=0...N-1), то мы в любой момент времени можем вычислить текущее значение зная значение в точке, с которой мы начали вести суммы SS (SS[0],m=0).
X[i]=SS[m] + X[m], где m - остаток от деления i/N, а X[m]-константа (истинный отсчет, значение буфера сумматора...)
Это полностью эквивалентно знанию всех значение буфера сумматора в какой-то момент времени с дальнейшим отслеживанием его содержимого.

Если вдуматься, то уже можно с высокой степенью достоверности утверждать, что задача восстановления исходных значений сводится к задаче восстановления содержимого буфера сумматора. Любые иные методы дадут дополнительные ошибки.

А сделать это можно или получая пропись с самого начала (нулевой буфер) или, как предложил ViKo, подав калибровочный сигнал. Если форма и время начала сигнала известны, то буфер можно восстановить практически точно (до погрешности АЦП). И вот как именно это сделать и стоит подумать. Как пример, в процессе работы вместо сигнала с датчика подать некоторый медленный и достаточно длинный сигнал (период заведомо много больше N). Если время его включения неизвестно, можно попытаться поиграться с корреляцией (по самим суммам) и максимально точно определив его время подогнать значения в буфере (уже по отсчетам). Можно сделать несколько итераций.

Да, давайте считать все вычисления в целых, без деления и т.д. Понятно, что любая плавучка даст погрешности и они приведут к ошибкам, но это следующий уровень приближения. Если калибровочный сигнал велик, то отосительные погрешности восстановления буфера будут порядка погрешностей АЦП, а это уже неистребимые никакими цифровыми методами ошибки.
Go to the top of the page
 
+Quote Post
GetSmart
сообщение Jan 14 2014, 21:49
Сообщение #10


.
******

Группа: Участник
Сообщений: 4 005
Регистрация: 3-05-06
Из: Россия
Пользователь №: 16 753



Цитата(rudy_b @ Jan 14 2014, 20:26) *
Если вдуматься, то уже можно с высокой степенью достоверности утверждать, что задача восстановления исходных значений сводится к задаче восстановления содержимого буфера сумматора. Любые иные методы дадут дополнительные ошибки.

А сделать это можно или получая пропись с самого начала (нулевой буфер) или, как предложил ViKo, подав калибровочный сигнал. Если форма и время начала сигнала известны, то буфер можно восстановить практически точно (до погрешности АЦП). И вот как именно это сделать и стоит подумать.

Слишком абстрактная задача. Если какая-то информация о входном сигнале доступна "обходным путём" или наоборот (воздействием), то это сразу меняет дело. Калибровочный сигнал во-первых прервёт основной сигнал, а во-вторых достаточно на вход скользящего подать 0 и через N тактов буфер скользящего самоочистится. Или любую постоянку/НЧ.

Или другой вариант. Входной сигнал сглаживается через два скользящих с разным N, некратным друг другу (простые числа). На выходе имеем два отдельных сигнала, функционально связанных с исходным. Тогда входной сигнал наблюдать или воздействовать на него не нужно.

ЗЫ
Ранее, восстанавливая сигнал через ДПФ, наверное надо было слепые зоны умножать на 1. Это для практического восстановления сигнала. Для математического конечно же не хватает информации о состоянии фильтра. Практическое восстановление будет вполне правдоподобным.

Сообщение отредактировал GetSmart - Jan 14 2014, 22:31


--------------------
Заблуждаться - Ваше законное право :-)
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- _pv   Операция обратная к "скользящему среднему"   Jan 6 2014, 12:36
- - Fat Robot   Для решения задачи нужно прежде всего почитать про...   Jan 6 2014, 13:11
|- - iiv   Цитата(Fat Robot @ Jan 6 2014, 18:11) Для...   Jan 6 2014, 13:47
|- - Tarbal   Цитата(iiv @ Jan 6 2014, 17:47) и конечно...   Jan 6 2014, 14:14
|- - iiv   Цитата(Tarbal @ Jan 6 2014, 19:14) Умноже...   Jan 6 2014, 15:09
- - delaver   Цитата(_pv @ Jan 6 2014, 16:36) есть изме...   Jan 6 2014, 13:30
|- - _pv   Цитата(delaver @ Jan 6 2014, 20:30) Не мо...   Jan 6 2014, 14:17
|- - GetSmart   Цитата(_pv @ Jan 6 2014, 20:17) да, всё о...   Jan 8 2014, 09:04
|- - _pv   Цитата(GetSmart @ Jan 8 2014, 16:04) Отку...   Jan 8 2014, 09:33
- - Tarbal   Возможно прямоугольная форма импульсной характерис...   Jan 6 2014, 14:20
- - GetSmart   Цитата(_pv @ Jan 6 2014, 18:36) каким обр...   Jan 6 2014, 23:58
- - AndrewN   QUOTE (_pv @ Jan 6 2014, 16:36) есть изме...   Jan 7 2014, 14:24
|- - _pv   Цитата(AndrewN @ Jan 7 2014, 21:24) Никак...   Jan 7 2014, 14:42
|- - AndrewN   QUOTE (_pv @ Jan 7 2014, 17:42) QUOTE X[i...   Jan 7 2014, 15:36
- - TSerg   Цитата(_pv @ Jan 6 2014, 16:36) ..есть из...   Jan 7 2014, 17:27
- - Fat Robot   Вот и синтезируйте фильр, обратный к вашему скольз...   Jan 7 2014, 20:19
- - RHnd   Я может что-то не понимаю, но в чем загвоздка? Взя...   Jan 7 2014, 21:10
|- - Fat Robot   Посмотрите, где лежат полюсы вашего "обратног...   Jan 7 2014, 21:31
|- - RHnd   Цитата(Fat Robot @ Jan 8 2014, 01:31) Пос...   Jan 8 2014, 16:10
- - thermit   Топикстартер задал тривиальный вопрос, на который ...   Jan 7 2014, 21:57
|- - AndrewN   Верующие могут отправлять свои естественные культо...   Jan 7 2014, 23:22
|- - Fat Robot   Особо неверующим будет вдвойне интересно, если в к...   Jan 8 2014, 07:25
|- - ViKo   Цитата(thermit @ Jan 8 2014, 00:57) Вот с...   Jan 8 2014, 09:10
- - thermit   ЦитатаFat Robot: Особо неверующим будет вдвойне ин...   Jan 8 2014, 07:52
- - thermit   ЦитатаViKo: о достаточно потерять хотя бы одно вхо...   Jan 8 2014, 09:24
|- - ViKo   Цитата(thermit @ Jan 8 2014, 12:24) Добро...   Jan 8 2014, 09:31
|- - ViKo   Цитата(ViKo @ Jan 8 2014, 12:31) Но, так ...   Jan 8 2014, 11:58
- - GetSmart   По существу же сказано, если имеются только данные...   Jan 8 2014, 09:29
- - GetSmart   С рекурсивным алгоритмом в принципе "нечестно...   Jan 8 2014, 10:30
- - rudy_b   Задачка, неожиданно, оказалась интересной. Смотрит...   Jan 8 2014, 11:34
|- - GetSmart   Цитата(rudy_b @ Jan 8 2014, 17:34) Как пр...   Jan 8 2014, 12:15
|- - rudy_b   Цитата(GetSmart @ Jan 8 2014, 16:15) Из ч...   Jan 9 2014, 13:18
|- - GetSmart   Цитата(rudy_b @ Jan 9 2014, 19:18) Да, со...   Jan 9 2014, 15:16
- - thermit   Дык, если задаться вопросом "что можно сделат...   Jan 8 2014, 12:25
|- - fontp   QUOTE (thermit @ Jan 8 2014, 16:25) Дык, ...   Jan 8 2014, 17:21
- - GetSmart   Скользящее среднее нелинейно искажает высокие част...   Jan 8 2014, 12:50
|- - _pv   Цитата(GetSmart @ Jan 8 2014, 19:50) Скол...   Jan 8 2014, 16:17
||- - GetSmart   Цитата(_pv @ Jan 8 2014, 22:17) но ничего...   Jan 8 2014, 16:19
|- - Tarbal   Цитата(GetSmart @ Jan 8 2014, 16:50) Скол...   Jan 10 2014, 14:51
|- - GetSmart   Цитата(Tarbal @ Jan 10 2014, 20:51) Разве...   Jan 10 2014, 15:33
||- - Tarbal   Цитата(GetSmart @ Jan 10 2014, 19:33) Свё...   Jan 10 2014, 16:16
||- - GetSmart   Цитата(Tarbal @ Jan 10 2014, 22:16) Сверт...   Jan 10 2014, 20:16
||- - Tarbal   Цитата(GetSmart @ Jan 11 2014, 00:16) Она...   Jan 10 2014, 23:31
|- - GetSmart   Цитата(rudy_b @ Jan 10 2014, 22:55) Отсюд...   Jan 10 2014, 17:26
- - thermit   Цитатаfontp: В присутствии шума ничем нельзя эти н...   Jan 8 2014, 17:29
|- - fontp   QUOTE (thermit @ Jan 8 2014, 21:29) С эти...   Jan 8 2014, 17:56
- - GetSmart   Можно предложить сперва поиск стартовых значений д...   Jan 9 2014, 00:59
- - rudy_b   Т.е. получается как с алиасом. Среднее по N точкам...   Jan 9 2014, 20:29
|- - GetSmart   Цитата(rudy_b @ Jan 10 2014, 02:29) Т.е. ...   Jan 10 2014, 08:16
|- - Tanya   Интересный велосипед тут изобретают... Автор прого...   Jan 10 2014, 09:03
|- - AndrewN   QUOTE (Tanya @ Jan 10 2014, 12:03) спектр...   Jan 10 2014, 20:03
- - TSerg   ТС проговорился, что "есть некие измерения пр...   Jan 10 2014, 10:35
|- - Tanya   Цитата(TSerg @ Jan 10 2014, 14:35) Собств...   Jan 10 2014, 16:55
|- - TSerg   Цитата(Tanya @ Jan 10 2014, 20:55) Вот до...   Jan 10 2014, 18:05
|- - Tanya   Цитата(TSerg @ Jan 10 2014, 22:05) Ничего...   Jan 10 2014, 18:39
- - rudy_b   При нормальной работе, если мы получаем пропись с ...   Jan 10 2014, 16:55
|- - AndrewN   QUOTE (rudy_b @ Jan 10 2014, 19:55) Т.е. ...   Jan 10 2014, 23:36
- - Tarbal   Предположим, что начальное значение было измерено ...   Jan 10 2014, 19:32
- - GetSmart   А был ли мальчик? Если рассуждать логически, то и...   Jan 11 2014, 13:55
|- - GetSmart   Цитата(bookd @ Jan 12 2014, 02:32) Ну не ...   Jan 12 2014, 00:38
|- - Tanya   Цитата(Tarbal @ Jan 13 2014, 16:11) и еще...   Jan 13 2014, 12:49
||- - iiv   Цитата(Tanya @ Jan 13 2014, 18:49) Вот по...   Jan 13 2014, 13:16
||- - Tanya   Цитата(iiv @ Jan 13 2014, 17:16) в хромот...   Jan 13 2014, 13:39
||- - iiv   Цитата(Tanya @ Jan 13 2014, 19:39) В коло...   Jan 13 2014, 15:20
||- - Tanya   Цитата(iiv @ Jan 13 2014, 19:20) интересн...   Jan 13 2014, 15:50
|- - ViKo   Цитата(Stanislav @ Jan 13 2014, 21:39) Сп...   Jan 13 2014, 19:39
||- - Stanislav   Цитата(ViKo @ Jan 13 2014, 23:39) Пусть о...   Jan 13 2014, 19:43
||- - ViKo   Цитата(Stanislav @ Jan 13 2014, 22:43) Ды...   Jan 13 2014, 19:44
||- - AndrewN   QUOTE (Stanislav @ Jan 14 2014, 00:23) На...   Jan 13 2014, 23:14
|- - Tarbal   Цитата(Stanislav @ Jan 13 2014, 21:39) По...   Jan 13 2014, 22:27
- - thermit   Цитатаbookd: А так хочется увидеть качественное р...   Jan 11 2014, 21:00
|- - Stanislav   Цитата(thermit @ Jan 12 2014, 01:00) ...в...   Jan 12 2014, 06:45
- - Stanislav   Как это что? Тождество. Вместо m подставьте выраже...   Jan 13 2014, 20:04
- - Rst7   Moderator: Я смотрю, тут началось отнюдь не конст...   Jan 13 2014, 20:10
- - Tarbal   Станислав, куда вы пропали? Я с нетерпением жду во...   Jan 14 2014, 14:07
- - thermit   ЦитатаTarbal: Можете написать импульсную характери...   Jan 14 2014, 19:44
- - Tarbal   Чем вам дискуссия не нравится? Как из этой импуль...   Jan 14 2014, 19:52
- - thermit   Цитатаrudy_b: Скользящее среднее это именно интегр...   Jan 14 2014, 20:16
|- - AndrewN   QUOTE (thermit @ Jan 15 2014, 00:16) Нуда...   Jan 14 2014, 20:54
|- - Stanislav   Цитата(AndrewN @ Jan 15 2014, 00:54) Help...   Jan 14 2014, 21:02
- - Stanislav   Коль попросили, отвечу. Цитата(Tarbal @ Jan ...   Jan 14 2014, 20:26
- - thermit   ЦитатаTarbal: Как из этой импульсной характеристик...   Jan 14 2014, 20:40
|- - Tarbal   Цитата(thermit @ Jan 14 2014, 23:40) Если...   Jan 15 2014, 04:16
- - thermit   ЦитатаAndrewN: А чего вы, собственно, доказываете?...   Jan 14 2014, 20:56
|- - AndrewN   QUOTE (thermit @ Jan 14 2014, 23:56) Я ни...   Jan 14 2014, 21:04
|- - rudy_b   Цитата(thermit @ Jan 15 2014, 00:56) Я ни...   Jan 14 2014, 21:08
|- - Stanislav   Цитата(thermit @ Jan 15 2014, 00:56) А чи...   Jan 14 2014, 21:17
|- - rudy_b   Цитата(Stanislav @ Jan 15 2014, 01:13) Та...   Jan 14 2014, 21:20
|- - Stanislav   Цитата(rudy_b @ Jan 15 2014, 01:20) Дык п...   Jan 14 2014, 21:24
- - thermit   Цитатаrudi_b: Красиво врут теоретики, просто прият...   Jan 14 2014, 21:47
- - thermit   ЦитатаTarbal: Я не это спросил. Ваши попытки Я рас...   Jan 15 2014, 09:55
2 страниц V   1 2 >


Reply to this topicStart new topic
2 чел. читают эту тему (гостей: 2, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 23rd July 2025 - 20:30
Рейтинг@Mail.ru


Страница сгенерированна за 0.03053 секунд с 7
ELECTRONIX ©2004-2016