|
Офлафн оценка скорости по дискретным отсчетам, Какие есть методы? |
|
|
|
Aug 15 2013, 20:43
|
Знающий
   
Группа: Свой
Сообщений: 518
Регистрация: 12-04-07
Из: Санкт-Петербург
Пользователь №: 26 997

|
Исходные данные - есть дискретные квантованные отсчеты некоторого сигнала. Мы знаем, что этот сигнал - выход некоторой механической системы, т.е. как минимум дважды дифференцируем. Более того, если надо, то мы можем априорно предположить некоторые ограничения на первую и вторую производные. Уровень квантования сигнала достаточно маленький по сравнению с диапазоном изменений сигнала. Частота дискретизации достаточно велика по сравнению со скоростью самого процесса. Скажем, дискретизация 1кГц при том, что в составе сигнала нет существенных компонент с частотами выше 50Гц. Записанный процесс достаточно длинный - несколько сотен тысяч отсчетов, сигнал имеет совершенно произвольную форму (с учетом дифференцируемости), но преобладают низкие частоты.
Задача - как по имеющимся отсчетам оценить скорость? Какие есть оффлайновые методы? Если кто-то сможет дать ссылку на какой-то перечень или обзор, то было бы здорово. Рабочая среда - матлаб.
Я бы решал задачу так - прорядил бы отсчеты, взяв, например, каждый 20-й. Полученные отсчеты описал бы каким-то сплайном, а затем по найденному сплайну вычисли бы значения производной во всех исходных точках. Плюс - в матлабе это делается в две строчки, и если исходный сигнал достаточно гладкий, то получается приемлемая точность. Но при таком методе есть недостаток - выпадают промежуточные, прореженные точки, т.е. часть информации пренебрегается. Кто как решает, на какие методы стоит посмотреть?
|
|
|
|
|
 |
Ответов
Guest_TSerg_*
|
Aug 21 2013, 11:44
|
Guests

|
Если мы ведем речь о динамической, возможно упрощенно линейной системе, то методы идентификации ее параметров или уточнения их давным давно изучены. Есс-но, нужна более правдоподобная модель, чем просто два интегратора, но даже и такой случай не критичен. Что это дает? Это дает правильное понимание методических и алгоритмических возможностей по извлечению полезных (внутренних) сигналов из выходных наблюдаемых координат. Каким бы ни был вход, корреляционная функция выхода хорошо будет известна в рамках принятой модели, а ее уточнение - чисто техническая задача. И опирается такой подход на теорию динамических систем, классический стат. аппарат и предельные теоремы. Если говорить о теоретико-методическоих подходах, то это два основных: метод пространства состояний и метод ошибки предсказания. P.S. Заодно уж найдите и почитайте о "традиционном" подходе, предложенным Тихоновым Андреем Николаевичем в его трудах "Методы решения некорректных задач", Наука 1974-1979 P.P.S. Задача ЦДС давным давно решена разными способами в рамках тех же ПВО и ПРО, иначе бы С-300 промахивались
|
|
|
|
|
Aug 21 2013, 12:56
|
Знающий
   
Группа: Свой
Сообщений: 518
Регистрация: 12-04-07
Из: Санкт-Петербург
Пользователь №: 26 997

|
В целом, все ваши рассуждения об известной динамической моедли и тому подобном - справедливы. Но это все - общие слова. Я испытываю некоторые сомнения, что они существенно пригодятся для решения данной задачи. Какие методические и алгоритмические возможности можно вытащить из этой информации, которых нет в изначальной постановке задачи? Цитата(TSerg @ Aug 21 2013, 15:44)  Задача ЦДС давным давно решена разными способами в рамках тех же ПВО и ПРО, иначе бы С-300 промахивались  Я интересуюсь оффлайновыми методами, которые, по моему убеждению, могут дать более точную оценку скорости, так как используют больший объем информации, в том числе информацию из будущего, которая онлайновым методам недоступна. Цитата(TSerg @ Aug 21 2013, 15:44)  Заодно уж найдите и почитайте о "традиционном" подходе, предложенным Тихоновым Андреем Николаевичем в его трудах "Методы решения некорректных задач", Наука 1974-1979 Постараюсь посмотреть. Что-то конкретно, или всю книгу вообще? Выкладываю данные. Это модельные данные, так что исходная скорость известна. Тут t - время, x1, w1 - положение и скорость в первом эксперименте, x2, w2 - положение и скорость во втором эксперименте, quant_size - размер кванта. Отличаются два эксперимента наличием среднечастотных составляющих, до 50 Гц. Я ищу метод, который в обоих случаях дает оценку, близкую к модельной скорости. Пока лучше всего себя показывает подход, который я описал в сообщении #27. Полученные этим методом скорости в переменных est_w1, est_w2.
Прикрепленные файлы
data.zip ( 371.82 килобайт )
Кол-во скачиваний: 29
|
|
|
|
Сообщений в этой теме
RHnd Офлафн оценка скорости по дискретным отсчетам Aug 15 2013, 20:43 _pv тогда сразу натягивать на все точки кубический спл... Aug 15 2013, 20:56 RHnd Спасибо за наводку на сглаживающий сплайн. Не очен... Aug 15 2013, 22:09 _pv Цитата(RHnd @ Aug 16 2013, 05:09) А что в... Aug 16 2013, 06:04 Tanya Цитата(RHnd @ Aug 16 2013, 02:09) Дополню... Aug 19 2013, 13:14 RHnd Это работало бы, если бы пренебрежимо малое кванто... Aug 16 2013, 06:25 thermit Сначала ФНЧ, затем дифференцировать. ФНЧ ограничит... Aug 16 2013, 16:11 RHnd Использование ФНЧ приведет к появлению задержки, ч... Aug 16 2013, 16:33 _pv ЦитатаСначала ФНЧ, затем дифференцировать. ФНЧ огр... Aug 16 2013, 17:05 Alex11 Если Вы работаете offline, то что Вам до задержки.... Aug 16 2013, 17:03 RHnd Цитата(Alex11 @ Aug 16 2013, 21:03) Если ... Aug 16 2013, 17:10 thermit Цитата_pv:
операции линейные, от смены очерёдности... Aug 16 2013, 20:16 RHnd Цитата(thermit @ Aug 17 2013, 00:16) Но т... Aug 16 2013, 21:57  _pv Цитата(RHnd @ Aug 17 2013, 03:57) Вспомни... Aug 17 2013, 07:37   RHnd Цитата(_pv @ Aug 17 2013, 11:37) у БИХ фа... Aug 17 2013, 08:29 thermit ЦитатаRHnd:
Почему?
Да зачем прореживать, чтобы п... Aug 17 2013, 09:38 RHnd 1) Если в самом сигнале шумов нет и все проблемы с... Aug 17 2013, 10:43 thermit Вот скрипт иллюстрирующий оценку производной при п... Aug 18 2013, 19:52 RHnd Спасибо за код.
Смотрите, убираем вашу строчку
К... Aug 19 2013, 06:17 thermit Это была иллюстрация. Ваши поправки требуют измене... Aug 19 2013, 12:09 RHnd Цитата(thermit @ Aug 19 2013, 16:09) Это ... Aug 19 2013, 14:34  Tanya Цитата(RHnd @ Aug 19 2013, 18:34) В общем... Aug 19 2013, 15:59 RHnd Tanya, простите, а в чем смысл вашего коментария? ... Aug 19 2013, 17:17 Tanya Цитата(RHnd @ Aug 19 2013, 21:17) Tanya, ... Aug 20 2013, 14:59 thermit ЦитатаRHnd:
1) Из каких соображений L=300 заменили... Aug 20 2013, 12:55 RHnd Но если взять такой фнч, то при наличии синусоидал... Aug 20 2013, 14:06 RHnd Я не пойму, к чему вы ведете. Вы считаете, что рас... Aug 20 2013, 16:20 Tanya Цитата(RHnd @ Aug 20 2013, 19:59) Я не по... Aug 20 2013, 16:27  RHnd Это вы мне первую разность предложили использовать... Aug 20 2013, 16:37   Tanya Цитата(RHnd @ Aug 20 2013, 20:37) Это вы ... Aug 21 2013, 07:13 TSerg Это не задача ЦОС, это задача идентификации динами... Aug 21 2013, 06:57 RHnd Цитата(TSerg @ Aug 21 2013, 10:57) Это не... Aug 21 2013, 07:10 TSerg >Да, если бы решалась онлайн.
Не имеет значени... Aug 21 2013, 07:38 RHnd Цитата(TSerg @ Aug 21 2013, 11:38) Если п... Aug 21 2013, 08:13  Tanya Цитата(RHnd @ Aug 21 2013, 12:13) Пожалуй... Aug 21 2013, 08:50   RHnd Цитата(Tanya @ Aug 21 2013, 12:50) Вот то... Aug 21 2013, 09:03    Tanya Цитата(RHnd @ Aug 21 2013, 13:03) будет в... Aug 21 2013, 09:21     RHnd Цитата(Tanya @ Aug 21 2013, 13:21) Вот у ... Aug 21 2013, 09:35      Tanya Цитата(RHnd @ Aug 21 2013, 13:35) Видимо,... Aug 21 2013, 09:46       RHnd Цитата(Tanya @ Aug 21 2013, 13:46) Что зн... Aug 21 2013, 10:27        Tanya Цитата(RHnd @ Aug 21 2013, 14:27) Минимиз... Aug 21 2013, 10:41 TSerg Рискну предположить, что речь идет, к примеру, об ... Aug 21 2013, 11:20 RHnd Цитата(TSerg @ Aug 21 2013, 15:20) Рискну... Aug 21 2013, 11:30 thermit ЦитатаПостараюсь посмотреть. Что-то конкретно, или... Aug 21 2013, 13:13 RHnd Цитата(thermit @ Aug 21 2013, 17:13) 4-я ... Aug 21 2013, 13:40 TSerg Из новомодных - Hilbert–Huang transform (HH... Aug 21 2013, 13:55 RHnd Цитата(TSerg @ Aug 21 2013, 17:55) Из нов... Aug 21 2013, 15:43 thermit ЦитатаRHnd:
Спасибо, посмотрел. На сколько я с нас... Aug 21 2013, 14:21 TSerg HHT
http://www.mathworks.com/matlabcentral/fil...m... Aug 21 2013, 16:13 RHnd Цитата(TSerg @ Aug 21 2013, 20:13) HHT
Сп... Aug 21 2013, 17:21 _Ivana Цитата(RHnd @ Aug 21 2013, 21:21) Какое о... Aug 21 2013, 22:50 RHnd Если вы для разных наборов использовали разные пар... Aug 21 2013, 23:08 iDiode Пользовался вот этим:
http://en.wikipedia.org/wiki... Oct 28 2013, 05:19
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0
|
|
|