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

 
 
 
Reply to this topicStart new topic
> Ковариационные матрицы. Как вычислить.
Ugedey
сообщение Aug 20 2008, 07:58
Сообщение #1





Группа: Validating
Сообщений: 3
Регистрация: 20-08-08
Пользователь №: 39 705



Столкнулся с необходимостью разработки фильтра Калмана для INS. Задача для меня новая, посему плаваю во многих аспектах.
С алгоритмом работы разобрался. С матрицами перехода и измерений разобрался. Проблема возникла с составлением ковариационных матриц шума процесса и шума наблюдений. Не могу понять, откуда брать компоненты.
Можете помочь чем-то? Объяснить, что и откуда берется.
Go to the top of the page
 
+Quote Post
Builder
сообщение Aug 21 2008, 16:24
Сообщение #2


iBuilder©
****

Группа: Свой
Сообщений: 519
Регистрация: 14-07-04
Из: Минск
Пользователь №: 322



Цитата(Ugedey @ Aug 20 2008, 10:58) *
Столкнулся с необходимостью разработки фильтра Калмана для INS. Задача для меня новая, посему плаваю во многих аспектах.
С алгоритмом работы разобрался. С матрицами перехода и измерений разобрался. Проблема возникла с составлением ковариационных матриц шума процесса и шума наблюдений. Не могу понять, откуда брать компоненты.
Можете помочь чем-то? Объяснить, что и откуда берется.


Как вариант - методом тыка подобрать, я как-то использовол этот метод - инфы по шумам для моей задачи небыло в принципе, так подбирал эксперементально, по наилучшему результату. Хотя в какой-то книге что-то по этому поводу было, но там всё тоже не так просто получалось, мне оказалось проще побобрать тыком (экспериментом). Книгу с ходу не помню, нет с собой архива - если не подскажут - напишите в личку, посмотрю как смогу.
Go to the top of the page
 
+Quote Post
Ugedey
сообщение Aug 26 2008, 11:56
Сообщение #3





Группа: Validating
Сообщений: 3
Регистрация: 20-08-08
Пользователь №: 39 705



Цитата(Builder @ Aug 21 2008, 19:24) *
Как вариант - методом тыка подобрать, я как-то использовол этот метод - инфы по шумам для моей задачи небыло в принципе, так подбирал эксперементально, по наилучшему результату. Хотя в какой-то книге что-то по этому поводу было, но там всё тоже не так просто получалось, мне оказалось проще побобрать тыком (экспериментом). Книгу с ходу не помню, нет с собой архива - если не подскажут - напишите в личку, посмотрю как смогу.


Ну шум будет стандартный - гироскоп+одометр, в принципе, взаимонезависимые и некореллируемые, поэтому, насколько я понял, первая ковариационная матрица процесса будет диагональной с главными элементами в виде квадрата дисперсий.
Методом ты - как именно?
Go to the top of the page
 
+Quote Post
Builder
сообщение Aug 27 2008, 12:03
Сообщение #4


iBuilder©
****

Группа: Свой
Сообщений: 519
Регистрация: 14-07-04
Из: Минск
Пользователь №: 322



Цитата(Ugedey @ Aug 26 2008, 14:56) *
Ну шум будет стандартный - гироскоп+одометр, в принципе, взаимонезависимые и некореллируемые, поэтому, насколько я понял, первая ковариационная матрица процесса будет диагональной с главными элементами в виде квадрата дисперсий.
Методом ты - как именно?

у меня было 2 сигнала - сигнал с датчика и сигнал управления системой.
Соответственно 2 шума - шум измерения и шум системы (управления), эти 2 значения я и подбирал.
Если речь идёт о этих данных.
Всё остальное впринципе просто, если есть модель системы и данные по шумам.
Go to the top of the page
 
+Quote Post
diwil
сообщение Sep 4 2008, 08:15
Сообщение #5


Местный
***

Группа: Свой
Сообщений: 366
Регистрация: 5-09-06
Из: Санкт-Петербург
Пользователь №: 20 107



я не большой специалист в кальмановской фильтрации, но сталкиватся с этим приходилось часто.
В моем случае измеряемая величина была одна - растояние. Тогда:

a = [1 dt dt^2/2; 0 1 dt; 0 0 1]; % transition matrix
c = [1 0 0]; % measurement matrix

xhat = [yym(1); 0; 0]; % initial state vector
E = ([dt^3/6 dt^2/2 dt]);
Q = accelnoise^2 * E' * E;

где Q - искомая ковариационная матрица
Go to the top of the page
 
+Quote Post
Ugedey
сообщение Sep 22 2008, 14:23
Сообщение #6





Группа: Validating
Сообщений: 3
Регистрация: 20-08-08
Пользователь №: 39 705



diwil, можно описать, что за формулы Вы привели? И что за операторы применяли? Я не совсем понял суть.

Builder, можете описать алгоритм своего примера? Данные по шумам гироскопа и одометра берутся из характеристик изделия, так что можно полагать, что шумовые характеристики известны.
Go to the top of the page
 
+Quote Post

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

 


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


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