Синтез фильтра Калмана

Информация » Разработка системы автоматического управления углом тангажа легкого самолета » Синтез фильтра Калмана

Страница 2

(2.78)

где P(t) – ковариационная матрица ошибок оценивания.

Согласно (2.69) и (2.70) оказывается, что:

(2.79)

Начальное значение матрицы P(t):

(2.80)

Так как в соответствии с (2.76) справедливо , то согласно (2.67) получаем:

(2.81)

Ковариационная матрица ошибок является решением матричного дифференциального уравнения Риккати:

(2.82)

которое следует решать при начальных условиях (2.81).

Решим уравнение Риккати в численном виде. Для этого перейдем в дискретную область описания, используя следующие преобразования [13]:

(2.83)

(2.84)

где Т – период дискретизации,

или приближенно:

(2.85)

(2.86)

Принимая Т = 0.05 с и подставляя значения матриц А и В в формулы (2.85) и (2.86) соответственно получим:

Уравнение Риккати перепишем в дискретном виде [24] :

(2.87)

где к – номер шага.

Расчет проведем в пакете прикладных программ MatLab. Задаваясь начальным условием

R = 0.05, (2.88)

получим на первом шаге:

на восьмом шаге:

на пятнадцатом шаге:

на шестнадцатом шаге:

Таким образом, уравнение (2.87) сходится за 15 шагов.

Определим матрицу коэффициентов усиления фильтра K по формуле (2.78):

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

Пусть сигнал управления u(t) в системе (2.61) и (2.62) формируется по оценке состояния Z(t) в виде:

(2.89)

Подставляя это уравнение в уравнения (2.61) и (2.77), получаем:

(2.90)

В этом случае уравнение оптимального фильтра имеет форму (2.77) с учетом (2.89), т. е.

(2.91)

В пределах каждого из режимов работы будем считать систему стационарной. Тогда оценка Z(t) используется для замыкания стационарной системы

(2.92)

в виде:

(2.93)

где D – постоянная матрица, выбранная таким образом, чтобы система

(2.94)

была асимптотически устойчива.

Страницы: 1 2 

Другое по теме:

Горизонтальный анализ активов предприятия
Горизонтальный анализ внеоборотных активов (см. таблицу 14). Анализ внеоборотных активов предприятия показывает, что внеоборотные активы предприятия за период с 2000 года по отношению к 2001 году увеличились на 14,4 млн. руб. (на 22%), в ...

Определение коэффициента исспользования сцепной массы
Результирующая нагрузка на колёсную пару, кН , (5.1) где - изменение нагрузки от работы ТЭД, кН; - изменение нагрузки от действия силы тяги, кН, Изменение нагрузки от работы ТЭД, кН (5.2) где - расстояние между осью вращения кол ...

Предложения по совершенствованию практики ценообразования
Расчет цены на автомобиль ВАЗ может осуществляться смешанным методом, который учитывает с одной стороны конъюнктуру рынка, а с другой стороны необходимы уровень затрат и прибыли предприятия. Метод можно назвать «затраты-результаты». Рас ...

Навигация

Copyright © 2018 - All Rights Reserved - www.transportgood.ru