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

 
 
> Офлафн оценка скорости по дискретным отсчетам, Какие есть методы?
RHnd
сообщение Aug 15 2013, 20:43
Сообщение #1


Знающий
****

Группа: Свой
Сообщений: 518
Регистрация: 12-04-07
Из: Санкт-Петербург
Пользователь №: 26 997



Исходные данные - есть дискретные квантованные отсчеты некоторого сигнала. Мы знаем, что этот сигнал - выход некоторой механической системы, т.е. как минимум дважды дифференцируем. Более того, если надо, то мы можем априорно предположить некоторые ограничения на первую и вторую производные. Уровень квантования сигнала достаточно маленький по сравнению с диапазоном изменений сигнала. Частота дискретизации достаточно велика по сравнению со скоростью самого процесса. Скажем, дискретизация 1кГц при том, что в составе сигнала нет существенных компонент с частотами выше 50Гц.
Записанный процесс достаточно длинный - несколько сотен тысяч отсчетов, сигнал имеет совершенно произвольную форму (с учетом дифференцируемости), но преобладают низкие частоты.

Задача - как по имеющимся отсчетам оценить скорость? Какие есть оффлайновые методы? Если кто-то сможет дать ссылку на какой-то перечень или обзор, то было бы здорово. Рабочая среда - матлаб.

Я бы решал задачу так - прорядил бы отсчеты, взяв, например, каждый 20-й. Полученные отсчеты описал бы каким-то сплайном, а затем по найденному сплайну вычисли бы значения производной во всех исходных точках. Плюс - в матлабе это делается в две строчки, и если исходный сигнал достаточно гладкий, то получается приемлемая точность. Но при таком методе есть недостаток - выпадают промежуточные, прореженные точки, т.е. часть информации пренебрегается. Кто как решает, на какие методы стоит посмотреть?
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
Guest_TSerg_*
сообщение Aug 21 2013, 11:44
Сообщение #2





Guests






Если мы ведем речь о динамической, возможно упрощенно линейной системе, то методы идентификации ее параметров или уточнения их давным давно изучены.
Есс-но, нужна более правдоподобная модель, чем просто два интегратора, но даже и такой случай не критичен.
Что это дает?
Это дает правильное понимание методических и алгоритмических возможностей по извлечению полезных (внутренних) сигналов из выходных наблюдаемых координат.
Каким бы ни был вход, корреляционная функция выхода хорошо будет известна в рамках принятой модели, а ее уточнение - чисто техническая задача.
И опирается такой подход на теорию динамических систем, классический стат. аппарат и предельные теоремы.
Если говорить о теоретико-методическоих подходах, то это два основных: метод пространства состояний и метод ошибки предсказания.

P.S.
Заодно уж найдите и почитайте о "традиционном" подходе, предложенным Тихоновым Андреем Николаевичем в его трудах "Методы решения некорректных задач", Наука 1974-1979

P.P.S.
Задача ЦДС давным давно решена разными способами в рамках тех же ПВО и ПРО, иначе бы С-300 промахивалисьsm.gif
Go to the top of the page
 
+Quote Post
RHnd
сообщение Aug 21 2013, 12:56
Сообщение #3


Знающий
****

Группа: Свой
Сообщений: 518
Регистрация: 12-04-07
Из: Санкт-Петербург
Пользователь №: 26 997



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

Цитата(TSerg @ Aug 21 2013, 15:44) *
Задача ЦДС давным давно решена разными способами в рамках тех же ПВО и ПРО, иначе бы С-300 промахивалисьsm.gif

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

Цитата(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
 
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- 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


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

 


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


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