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

 
 
> Фильтр второго порядка
Yevdokimenko
сообщение Jan 8 2013, 19:50
Сообщение #1


Участник
*

Группа: Участник
Сообщений: 44
Регистрация: 23-05-12
Пользователь №: 71 972



Добрый день.
Я не знаю в этом ли разделе должна быть эта тема, но всё же...
В ходе разработки устройства, одной из функций которого является определение ускорения, столкнулся с проблемой вычисления оного.
Итак, на автомобиле (на ступице колеса) установлен зубчатый венец, который вращается "вместе" с колесом.
Этот венец своими зубьям создает "наводки" на датчик, который "генерирует" сигнал с частотой, пропорциональной частоте вращения колеса (скорости движения).
Если венец и возможно изготовить в промышленных масштабах с хорошей точностью (но все же неидеально), то колесо как таковое не является абсолютно упругим телом, и имеет переменный радиус обката, что приводит к переменной частоте сигнала даже при условии равномерного движения (с постоянной скоростью).
Судя по проведенным тестам разброс частоты относительно некой средней подчиняется нормальному закону распределения случайной величины.
В случае же неравномерного (и даже не равноускоренного) движения разброс частоты становится еще большим.
Если для вычисления скорости (первой производной) такая ситуация вполне допустима, то при вычислении ускорения (второй производной) расчетный параметр ОЧЕНЬ зашумлен.
Посоветуйте, пожалуйста, с какой стороны подойти к решению данной задачи?
Пока что в голову пришла лишь мысль об аппроксимации МНК (методом наименьших квадратов) для нахождения линейной функции по последним N наблюдениям... Наклон данной функции и будет искомым ускорением.
Достаточно ли функции линейного вида, либо нужны функции больших порядков (квадратичные и т.п.)?
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
Yevdokimenko
сообщение Jan 12 2013, 21:28
Сообщение #2


Участник
*

Группа: Участник
Сообщений: 44
Регистрация: 23-05-12
Пользователь №: 71 972



Вы правы. Тогда что за фильтр, который я назвал скользящим средним?
Go to the top of the page
 
+Quote Post
Xenia
сообщение Jan 12 2013, 22:18
Сообщение #3


Гуру
******

Группа: Модератор FTP
Сообщений: 4 479
Регистрация: 20-02-08
Из: Москва
Пользователь №: 35 237



Цитата(Yevdokimenko @ Jan 13 2013, 01:28) *
Тогда что за фильтр, который я назвал скользящим средним?


Один из вариантов рекурсивного. Это потому что он использует уже отфильтрованные (исправленные на прошлых шагах) значения.

Тест здесь довольно простой. Представьте себе, что ваше колесо остановилось sm.gif, и с некоторого момента вы начинаете получать исключительно одни нулевые значения.
В этой ситуации нерекурсивный фильтр после n шагов (n - это предел суммирования тех самых сумм) полностью забудет предысторию. А на протяжении этих n шагов плавно опустит среднее до нуля. А ваш вариант рекурсивен, т.к. он не забудет "былое величие" никогда. Точнее говоря, оно будет рассасываться теоретически бесконечно долго (подобно тому как сходятся к нулю обратные величины натурального ряда).

Второй тест на ту же тему. Нерекурсивный фильтр дает одно и тоже значение среднего, независимо от того, двигается ли окно слева направо или справа налево. А для рекурсивных это не одно и тоже, т.е. слева и справа разная предыстория.
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- Yevdokimenko   Фильтр второго порядка   Jan 8 2013, 19:50
- - diwil   Цитата(Yevdokimenko @ Jan 8 2013, 23:50) ...   Jan 9 2013, 07:01
|- - Yevdokimenko   Цитата(diwil @ Jan 9 2013, 10:01) И то и ...   Jan 9 2013, 07:37
- - Xenia   Цитата(Yevdokimenko @ Jan 8 2013, 23:50) ...   Jan 9 2013, 13:04
- - Xenia   Предположим, что мы пока остановились на N=7 (± 3 ...   Jan 9 2013, 15:37
- - Yevdokimenko   Спасибо, суть расчета я если честно не понял (как ...   Jan 11 2013, 17:55
|- - Xenia   Цитата(Yevdokimenko @ Jan 11 2013, 21:55)...   Jan 11 2013, 18:15
- - Yevdokimenko   Равномерный шаг нагадит в точность, а это то, о че...   Jan 11 2013, 19:21
|- - Xenia   Цитата(Yevdokimenko @ Jan 11 2013, 23:21)...   Jan 11 2013, 20:20
|- - Yevdokimenko   Цитата(Xenia @ Jan 12 2013, 00:20) Ваш оп...   Jan 11 2013, 21:43
- - TSerg   Откуда Вы взяли пары x,y? У Вас измеряется лишь де...   Jan 12 2013, 09:48
|- - _Anatoliy   Цитата(TSerg @ Jan 12 2013, 11:48) Длител...   Jan 12 2013, 10:53
|- - Yevdokimenko   Цитата(_Anatoliy @ Jan 12 2013, 14:53) На...   Jan 12 2013, 14:08
|- - TSerg   Цитата(_Anatoliy @ Jan 12 2013, 14:53) На...   Jan 12 2013, 14:16
|- - Yevdokimenko   Цитата(TSerg @ Jan 12 2013, 18:16) То, чт...   Jan 12 2013, 18:14
|- - TSerg   Цитата(Yevdokimenko @ Jan 12 2013, 22:14)...   Jan 13 2013, 08:24
- - TSerg   http://ru.wikipedia.org/wiki/%D0%A1%D0%BA%...%BD%D...   Jan 12 2013, 20:23
|- - Yevdokimenko   Цитата(Xenia @ Jan 13 2013, 02:18) Один и...   Jan 13 2013, 09:23
- - TSerg   Вам надо решение найти или поспорить тут? Если спо...   Jan 13 2013, 10:33


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

 


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


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