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

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

|
>Да, если бы решалась онлайн.
Не имеет значения - офф или он. Если параметры модели известны - на основе обратной модели вычисляете вход, затем на основе входа и прямой модели - ненаблюдаемые переменные состояния.
|
|
|
|
|
Aug 21 2013, 08:13
|
Знающий
   
Группа: Свой
Сообщений: 518
Регистрация: 12-04-07
Из: Санкт-Петербург
Пользователь №: 26 997

|
Цитата(TSerg @ Aug 21 2013, 11:38)  Если параметры модели известны - на основе обратной модели вычисляете вход, затем на основе входа и прямой модели - ненаблюдаемые переменные состояния. Модель - два интегратора с неизвестным сигналом на входе. Какая такая обратная модель? Дифференциатор второго порядка? Так задача и сформулирована - как по набору измерений провести дифференцирование хотя бы первого порядка. А вы предлагаете сразу восстанавливать оценку входа (ускорения). Как? Цитата(Tanya @ Aug 21 2013, 11:13)  Вот если бы у Вас была модель, то Вы бы могли обучить нейронную сеть ... Если бы у меня была какая-то более информативная модель, чем два интегратора, а, еще лучше, какие-то входные сигналы этой какой-то модели, то я бы, как правильно пишет TSerg, рассматривал бы задачу в контексте идентификации и построения некоторого наблюдателя. Но у меня нет ни того, ни другого. Давайте оставаться в рамках постановки задачи. Цитата(Tanya @ Aug 21 2013, 11:13)  Это если бы Вы знали какие-то зависимости и ограничения на функцию и ее производные. А если не знаем, то может быть любое, построение которого я и предложила. Если результат Вы отвергаете, значит знаете нечто априори, что отвергает такое решение. Пожалуйста, прочитайте тему с начала. Там есть про ограничения на значения скорости и усорения, про ограничения на частотные компоненты (что будет порушено при взятии первой разности). И даже немного про выбор кривой из семейства за счет штрафа на ускорение. Более того, уже было предложение использовать первую разность. И был мой ответ. Цитата(Tanya @ Aug 21 2013, 11:13)  какую-нибудь искаженную (женскую) логику. Давайте воздерживаться от таких высказываний.
|
|
|
|
|
Aug 21 2013, 08:50
|
Гуру
     
Группа: Модераторы
Сообщений: 8 752
Регистрация: 6-01-06
Пользователь №: 12 883

|
Цитата(RHnd @ Aug 21 2013, 12:13)  Пожалуйста, прочитайте тему с начала. Там есть про ограничения на значения скорости и усорения, про ограничения на частотные компоненты (что будет порушено при взятии первой разности). И даже немного про выбор кривой из семейства за счет штрафа на ускорение. Более того, уже было предложение использовать первую разность. И был мой ответ.
Давайте воздерживаться от таких высказываний. Если бы Вы тоже внимательно читали, то воздержались бы от подобных высказываний. Я предлагала сглаживать ступеньки с учетом ограничений. Вот тогда и симулируйте - случайное значение ускорения (или еще какие условия от положения или скорости) интегрируем (знаем идеальную скорость), интегрируем второй раз - идеальное положение, квантуем и получаем наблюдаемый сигнал... И тренируем женскую логику или имитатор кусочка мозга. Вот что я и предлагала. А чем Вам не нравятся мои высказывания?
|
|
|
|
|
Aug 21 2013, 09:03
|
Знающий
   
Группа: Свой
Сообщений: 518
Регистрация: 12-04-07
Из: Санкт-Петербург
Пользователь №: 26 997

|
Цитата(Tanya @ Aug 21 2013, 12:50)  Вот тогда и симулируйте - случайное значение ускорения (или еще какие условия от положения или скорости) интегрируем (знаем идеальную скорость), интегрируем второй раз - идеальное положение, квантуем и получаем наблюдаемый сигнал... И тренируем женскую логику или имитатор кусочка мозга. Вот что я и предлагала. Т.е. вы прелагаете построить и обучить некоторую нейронную сеть, которая по набору квантованных измерений положения будет выдавать хорошую и не сдвинутую оценку скорости, так? И будет хорошо работать во всем диапазоне входных сигналов? Я с нейронными сетями знаком очень поверхностно, но в области управления не припомню никаких работ, где бы использовались дифференциаторы, построенные на нейронных сетях. Можете указать такие статьи/книги? Или сами что-то подобное делали? Пожалуйста, не надо писать про женскую логику. Мне неприятно.
|
|
|
|
|
Aug 21 2013, 09:21
|
Гуру
     
Группа: Модераторы
Сообщений: 8 752
Регистрация: 6-01-06
Пользователь №: 12 883

|
Цитата(RHnd @ Aug 21 2013, 13:03)  будет выдавать хорошую и не сдвинутую оценку скорости, так? И будет хорошо работать во всем диапазоне входных сигналов? Пожалуйста, не надо писать про женскую логику. Мне неприятно. Вот у меня нет критерия "хорошести". Ведь квантование искажает картинку. Даже если бы квантование по ординате было более мелкое, квантование по времени внесет всегда ошибку в вычисление скорости и ускорения. Если нет ограничений, то невозможно получить ничего хорошего. Чем теснее ограничения, тем лучше будет получаться. Но ведь нет другого пути? Хорошо, буду использовать alias (или эвфемизм) - FUZZY. Сети очень давно не трогала... Когда-то были free пакеты.
|
|
|
|
|
Aug 21 2013, 09:35
|
Знающий
   
Группа: Свой
Сообщений: 518
Регистрация: 12-04-07
Из: Санкт-Петербург
Пользователь №: 26 997

|
Цитата(Tanya @ Aug 21 2013, 13:21)  Вот у меня нет критерия "хорошести". Ведь квантование искажает картинку. Если нет ограничений, то невозможно получить ничего хорошего. Чем теснее ограничения. тем лучше будет получаться. Но ведь нет другого пути? Видимо, я раньше не достаточно отчетливо высказал мысль, что из множества кривых, описывающий квантованные данные, я хочу выбрать ту, где наименьший квадрат ускорения. Или где наименьший вклад компонентов с б ольшими частотами. Т.е. если две кривые описывают квантованные данные, но у одной существенно преобладают низкие частоты, а у другой заметно присутствуют среднечастотные, то первая кривая лучше. Вот такой у меня критерий "хорошести", который проистекает из физическиой природы сигнала. Цитата(Tanya @ Aug 21 2013, 13:21)  Сети очень давно не трогала... Когда-то были free пакеты. Я же не про пакеты спрашивал, их бесплатность меня на данном этапе вообще не волнует. Я спрашивал, известны ли вам успешные примеры использования нейронных сетей для решения поставленной задачи. Из литературы или из личного опыта. Мне - не известны.
|
|
|
|
Сообщений в этой теме
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      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 TSerg Если мы ведем речь о динамической, возможно упроще... Aug 21 2013, 11:44 RHnd В целом, все ваши рассуждения об известной динамич... Aug 21 2013, 12:56 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
|
|
|