Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Численные расчеты в среде MatLab. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
15.08.2026
Размер:
1 Мб
Скачать

 

 

Окончание таблицы

 

 

 

Кнопка

 

Назначение

 

 

 

 

 

Режим создания флажка

 

 

 

 

 

Режим создания поля для ввода текстовой инфор-

 

 

мации

 

 

 

 

 

Режим создания поля для ввода текстовой инфор-

 

 

мации без возможности редактирования

 

 

 

 

 

Режим создания списка

 

 

 

23. Установите соответствие между названием команды и ее результатом (с помощью стрелки):

Команда

 

Действие

 

 

 

GUIDE

 

Конструктор графического интерфейса

 

 

 

Property Inspec-

 

Управление сеткой и линейками

tor

 

 

Object Browser

 

Создание и редактирование меню

 

 

 

M-file Editor

 

Инспектор свойств

 

 

 

Component

 

Просмотр объектов

Callbacks

 

 

Figure Callbacks

 

Редактор М-файлов

 

 

 

Align Objects

 

Создание функций обработки событий для

 

 

компонентов

Grid and Rulers

 

Создание функций обработки событий для

 

 

окон

Menu Editor

 

Выравнивание объектов

 

 

 

51

Контрольные вопросы и задания

1.Перечислите компоненты, используемые при построении графического интерфейса.

2.С какой целью используется функция uicontrol при построении графического интерфейса?

3.Каким образом можно изменить свойства у существующего интерфейсного объекта?

4.Каким образом можно создать командную кнопку типа Push But-

ton?

5.Что позволяют задать свойства Backgroundcolor и ForegroundColor?

6.Чем отличается кнопка типа Toggle Button от командной кнопки?

7.В каком качестве используются компоненты типа panel и Button Group?

8.Каким образом используются кнопки типа checkbox?

9.Какой объект используется для вывода символьной строки (или нескольких строк) в выделенном прямоугольнике?

10.В каких случаях может использоваться интерфейсный элемент типа Edit Text?

11.Что собой представляет интерфейсный компонент типа Listbox?

12.Каким образом можно использовать интерфейсный компонент типа Slider?

13.Как установить полосу прокрутки в горизонтальном или вертикальном положении?

14.В каких случаях происходит вызов события callback?

15.Перечислите сервисные средства среды MatLab, упрощающие проектирование графического интерфейса приложения.

16.Какие варианты компоновки, которые можно положить в основу проектируемого приложения, предлагает список шаблонов GUIDE templace?

17.Каким образом осуществляется управление конструктором графического интерфейса?

18.В чем заключается основное назначение Property Inspector?

19.Каким образом можно осуществить выравнивание объектов на

форме?

20.Какое преобразование выполняет функция str2num ()?

21.Какая функция выполняет преобразование строкового значения

вчисловое значение?

52

3. Установившийся режим в задаче оценивания параметров динамических объектов

3.1.Формирование матрицы Фишера

вустановившемся режиме

Рассматривается динамический объект в следующей форме: xk 1 xk uk wk ,

yk 1 H xk 1 vk 1 ,

где xk 1 вектор состояния; uk вектор управления; wk вектор возмущения; yk 1 вектор наблюдения, vk 1 вектор ошибки измерения,

– матрица состояния; – матрица возмущения; – матрица управления; H – матрица наблюдения. При этом обозначим через Q неотри-

цательно определенную ковариационную матрицу вектора возмущения; через R – положительно определенную ковариационную матрицу вектора ошибки измерения. Предполагается, что все переходные процессы закончились, то естьрассматриваетсяустановившийсярежим.

В задаче активной идентификации при оценивании параметров линейных динамических дискретных объектов во многих работах используется фильтр Калмана с обновленной последовательностью и информационная матрица Фишера [25–28]. В данном разделе рассматривается установившийся режим для вычисления фильтра Калмана и матрицы Фишера в среде Simulink.

Приведем процедуру вычислений оценок параметров динамических дискретных объектов для случая, когда неизвестные параметры находятся в матрице состояния. Вывод формул, используемых в данном разделе, подробно изложен в работе [29]. В среде Simulink построена модель в виде блоков для вычисления оценок параметров динами-

53

ческих объектов на основе метода наименьших квадратов. Каждый блок соответствует определенной группе формул. Ниже дан набор формул, необходимый для проведения расчетов.

Первая группа:

 

 

 

 

 

 

 

,

 

(1)

Pk 1,k Pk,k

 

Q

 

 

 

k 1

HP

 

H R 1/2

,

 

(2)

 

k 1,k

 

 

 

 

 

 

 

 

K1k 1 Pk 1,k H k11 ,

 

 

 

(3)

P

 

I K1

1

1 H P

 

.

(4)

k 1,k 1

 

k

 

k 1

k 1,k

 

 

Первая группа формул соответствует вычислению фильтра Калмана с обновленной последовательностью в установившемся режиме [31]. На рис. 15 приведены результаты построения уравнения (1). Матрицы

, , Q являются входными данными для уравнения (1).

Рис. 15. Вычисление по уравнению (1)

На рис. 16 приведены результаты формирования уравнения (2). Матрицы H и R являются входными данными для уравнения (2).

54

Рис. 16. Вычисление по уравнению (2)

На рис. 17 даны результаты моделирования уравнения (3).

Рис. 17. Вычисление по уравнению (3)

Схема вычисления уравнения (4) показана на рис. 18. Уравнения второй группы формул:

P

 

 

 

P

T

 

1

 

 

P

T

0 T P

 

 

,

 

 

 

i

i

0

 

0

i

 

 

 

 

i

 

 

1

 

 

1

 

 

 

P

 

 

 

 

 

2

 

3H

1

HT ,

 

 

 

i

 

 

 

 

 

i

 

 

K1

P

HT 1

P HT

1

 

 

 

1

 

,

 

 

i

i

 

 

 

 

 

1

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

55

 

 

 

 

 

(5)

(6)

(7)

P0

i

 

 

K1

 

1

 

 

 

 

 

H

 

 

 

i

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

P1

 

 

 

 

 

 

 

 

 

 

K1

 

H P

 

I K1

 

 

H

 

, (8)

 

 

 

i

 

1

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 18. Вычисление по уравнению (4)

На рис. 19 приведена схема моделирования уравнения (5).

Рис. 19. Вычисление по уравнению (5)

56

На рис. 20 показан порядок формирования уравнения (6).

Рис. 20. Вычисление по уравнению (6)

На рис. 21 приводится схема вычислений уравнения (7).

Рис. 21. Вычисление по уравнению (7)

57

Третья группа формул:

 

 

K 2

 

K1

,

 

K 2

 

K1

K1 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

i

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A,1

 

 

K 2 1H

 

 

 

 

0

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

K 2 1H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K 2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

K 2

1

 

 

,

A,1

 

0

 

 

,

 

 

 

K A,1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K 2

 

 

 

1

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

A,1

 

A,1

x

A,1

 

A

u ,

 

A,1

 

 

 

 

 

T

K

A,1

KT

.

 

 

 

 

 

 

 

 

 

 

A,1 A,1 A,1

 

 

A,1

 

Exˆk 1,ki

 

xˆ

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

k 1,k

 

 

 

 

 

ci A,k 1,k

xA,k 1,k xAT,k 1,k

cTj ,

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ci

 

0 I 0 ,

 

 

 

 

 

 

 

 

N

1

 

 

 

 

 

 

 

 

 

 

 

M

 

 

 

 

 

 

 

k 1 2

k 1

 

 

 

 

Sp

 

 

 

 

 

ij

 

k 0

 

 

 

 

i

k 1

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

58

 

 

 

 

 

xˆk 1,k

 

xˆk 1,k

 

T

 

Sp HT 2 HE

 

 

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

 

 

 

 

 

 

 

 

i

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

k 1

1

 

k 1

 

 

 

 

 

 

 

 

 

 

 

 

Sp k 1

 

i

k 1

j

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

3.2. Моделирование итерационной процедуры оценивания динамических параметров в задаче активной идентификации

Для оценивания параметров динамических объектов используется рекуррентная схема метода наименьших квадратов [33]. Как известно, метод наименьших квадратов заключается в минимизации следующей функции квадратичного отклонения [34]:

s( )

N

y xT 2 .

 

 

i i

 

i 1

 

Тогда оценку неизвестных параметров на основе метода наименьших квадратов можно получить исходя из формулы

 

ˆ

 

 

T

 

 

1

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N X N X N

 

 

X NYN ,

 

 

 

где приняты обозначения:

 

 

 

 

 

 

 

 

 

 

 

 

 

xT

 

x1

x2

 

 

 

y

 

 

 

 

1

 

 

1

1

 

 

 

 

1

 

X N

xT

 

x1

x2

 

, YN

 

y

2

 

 

2

 

 

2

2

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

1

2

 

 

 

 

 

 

 

 

xN

xN xN

 

yN

59

Для получения рекуррентной оценки неизвестных параметров ˆ

на основе N 1 измерений необходимо воспользоваться нижеприведенной группой формул [34]:

ˆ

 

 

ˆ

 

 

 

 

 

 

 

T

 

 

ˆ

 

,

N 1

N KN

1 yN 1 xN 1 N

K

N 1

P x

N

1

1 xT

P x

N 1

,

 

 

N

 

N 1 N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

P

 

I P

 

 

xN 1xN 1

 

 

P

 

,

 

1

 

 

 

 

N 1

 

 

N

xT

 

P x

N

 

 

 

N

 

 

 

 

 

 

 

 

N 11

 

N

1

 

 

 

где KN 1 , PN 1 вычисляются по результатам N 1 измерений.

Результаты вычислений в значительной мере зависят как от задания начального значения 0 , так и от начальных значений K0 и P0 .

Заметим, что при нулевом значении P0 0 рекуррентный процесс не

сходится и получить удовлетворительную оценку ˆ N не удается.

Моделирование итерационной процедуры вычисления оценки параметров объекта выполнено в среде Simulink. Схема вычисления параметров объекта представлена на рис. 22, где введены следующие блоки: Input_Blok – генератор входного сигнала; Simulation – блок моделирования системы; Estimation – блок вычисления неизвестных параметров.

Рис. 22. Итерационная процедура вычисления оцениваемых параметров

60

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]