|
Операция обратная к "скользящему среднему" |
|
|
|
Jan 6 2014, 12:36
|
Гуру
     
Группа: Свой
Сообщений: 2 563
Регистрация: 8-04-05
Из: Nsk
Пользователь №: 3 954

|
есть измеренные данные про которые известно, что по ним прошлись скользяшим средним с неким известным размером D. каким образом можно восстановить исходные данные которые были до усреднения? из простого на ум приходит только пройтись по этим уже усреднённым данным еще раз таким же скользящим средним и получить оценку ошибки вызванную усреднением, которую потом добавить обратно к первоначальным данным. Оно вроде как работает, но понятно что это хоть и очень простой, но не совсем честный способ. Как это делать математически правильно? понятно что можно сделать Фурье, поделить спектр на характеристику скользящего среднего (по сути КИХ фильтра с одинаковыми единичными коэффициентами) и преобразовать обратно. но при делении на нули в том месте где характеристика имеет провалы пожалуй будут проблемы, оно и понятно так как частоты кратные размеру фильтра давятся в ноль и обратному восстановлению не подлежат. А если без преобразования в частотную область, с какой функцией надо сделать свёртку чтобы получить фильтр обратный скользящему среднему? или даже в общем случае, каким образом преобразовать коэффициенты КИХ фильтра, чтобы произведение исходного фильтра и пробразованного давало 1.
|
|
|
|
|
 |
Ответов
|
Jan 11 2014, 20:32
|
Частый гость
 
Группа: Свой
Сообщений: 151
Регистрация: 4-02-09
Из: Поволжье
Пользователь №: 44 403

|
Ну не знаю, читаю уже четвертую страницу, реальных советов ноль. ТС решил проблему еще на первой странице, а все остальное просто хлам на мой взгляд. А так хочется увидеть качественное решение, с исходными данными, конкретным результатом и списком налагаемых ограничений...
--------------------
Всеобщая дебилизация не повод наносить ущерб своему здоровью.
|
|
|
|
|
Jan 13 2014, 12:11
|
Профессионал
    
Группа: Свой
Сообщений: 1 351
Регистрация: 21-05-10
Пользователь №: 57 439

|
Цитата(bookd @ Jan 12 2014, 00:32)  Ну не знаю, читаю уже четвертую страницу, реальных советов ноль. ТС решил проблему еще на первой странице, а все остальное просто хлам на мой взгляд. А так хочется увидеть качественное решение, с исходными данными, конкретным результатом и списком налагаемых ограничений... Да ну. Я показал, что его решение можно использовать вопреки высказанному скептицизму и как произвести коррекцию ошибки измерения начального отсчета. Но это решение для скользящего среднего. Впоследствии выяснилось, что усреднение происходит физическими свойствами прибора и еще не доказано, что это скользящее среднее.
Сообщение отредактировал Tarbal - Jan 13 2014, 12:18
|
|
|
|
|
Jan 13 2014, 18:39
|

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

|
Понеслась опять нелёгкая, похоже... Цитата(Tarbal @ Jan 13 2014, 16:11)  Да ну. Я показал, что его решение можно использовать вопреки высказанному скептицизму и как произвести коррекцию ошибки измерения начального отсчета.
Но это решение для скользящего среднего. Впоследствии выяснилось, что усреднение происходит физическими свойствами прибора и еще не доказано, что это скользящее среднее. Вы полагаете, что совершили открытие, доказав мировую теорему о том, что:  , верно? Спешу Вас разочаровать - как раз для фильтра скользящего среднего обратный фильтр неустойчив. Т.е., прямая-обратная операция не является робастной. Цитата(Tanya @ Jan 13 2014, 19:50)  Удовлетворю Ваше любопытство - вот каждый день вижу хроматограф газовый. .......... Вот вспомнилось. Первое (почти) электронное устройство, вышедшее из моих рук, было для хроматографа. А вот и хымыки пожаловали. Жалобу на вас, что ли, накатать? Или предложить поискать площадку, где ваши таланты раскроются в полную силу? Пока попробую второе. Не соблаговолите ли сходить на ххымию, господа флудеры? http://www.ximuk.ru/
--------------------
Самонадеянность слепа. Сомнения - спутник разума. (с)
|
|
|
|
|
Jan 13 2014, 20:25
|
Местный
  
Группа: Участник
Сообщений: 336
Регистрация: 7-03-07
Из: Петербург
Пользователь №: 25 961

|
QUOTE (Stanislav @ Jan 13 2014, 22:39)  А вот и хымыки пожаловали Аппаратная функция (как пример интегрального преобразования) есть у всех измерительных приборов. Чем плох хроматограф в качестве примера?
|
|
|
|
|
Jan 13 2014, 20:43
|

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

|
Цитата(AndrewN @ Jan 14 2014, 00:25)  Аппаратная функция (как пример интегрального преобразования) есть у всех измерительных приборов. Чем плох хроматограф в качестве примера? Примера чего именно? И, самое интересное, какое этот пример имеет отношение к вопросу темы? ------------------------- Уважаемые модераторы. По-моему, тема не стоит и выеденного яйца, поскольку, по сути, посвящена "доказательству" (ну, или "опровержению") совершенно банального тождества (см. выше). Ответ получен топикстартером на первой же странице. Посему считаю, что её лучше закрыть, во избежание многостраничного флуда. Ну, или перенести в "курилку", на худой конец всё, что здесь написано после ответа на вопрос темы.
--------------------
Самонадеянность слепа. Сомнения - спутник разума. (с)
|
|
|
|
|
Jan 13 2014, 20:57
|
Местный
  
Группа: Участник
Сообщений: 336
Регистрация: 7-03-07
Из: Петербург
Пользователь №: 25 961

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

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

|
AndrewN, на мои вопросы ответа нет, я правильно понял? Хорошо, писать здесь больше не буду, если не попросят. И опять же, против поговорки: "мели, Емеля, - твоя неделя" ничего особо не имею. Только после этого называть Электроникс техническим форумом как-то язык не поворачивается... А о теории, уважаемый AndreyN, у нас настолько разные представления, что могу вам посоветовать держать свои советы при себе. Называть "интегрированием" фильтрацию скользящим средним можно только в понятийном пространстве, ортогональном к общепринятому. Желаю удачи на сей ниве, тем не менее.  Правила форума ещё бы не мешало почитать... У всех посетителей темы прошу прощения за оффтоп.
--------------------
Самонадеянность слепа. Сомнения - спутник разума. (с)
|
|
|
|
|
Jan 14 2014, 14:26
|
Знающий
   
Группа: Свой
Сообщений: 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). Если время его включения неизвестно, можно попытаться поиграться с корреляцией (по самим суммам) и максимально точно определив его время подогнать значения в буфере (уже по отсчетам). Можно сделать несколько итераций. Да, давайте считать все вычисления в целых, без деления и т.д. Понятно, что любая плавучка даст погрешности и они приведут к ошибкам, но это следующий уровень приближения. Если калибровочный сигнал велик, то отосительные погрешности восстановления буфера будут порядка погрешностей АЦП, а это уже неистребимые никакими цифровыми методами ошибки.
|
|
|
|
|
Jan 14 2014, 21:49
|
.
     
Группа: Участник
Сообщений: 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
--------------------
Заблуждаться - Ваше законное право :-)
|
|
|
|
Сообщений в этой теме
_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 страниц
1 2 >
2 чел. читают эту тему (гостей: 2, скрытых пользователей: 0)
Пользователей: 0
|
|
|