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

лр4

.pdf
Скачиваний:
6
Добавлен:
22.05.2015
Размер:
642.89 Кб
Скачать

Решение будем искать в матричном виде:

B ( X

T

X )

1

( X

T

Y )

 

 

 

Введите матрицу входных сигналов. Для этого растяните окно Mathcad за правый нижний угол. В меню Вид в строке Панели инструментов установите флажок в строке Математика (рис. 10).

Рис. 10. Выбор палитры Математика.

Появтся сводная палитра математичеких знаков Математика. С ее помощью выведите палитры Оценка и Матрица, щелкнув соответственно

на символы и . Появятся указанные палитры рис. 11.

Рис. 11. Палитры математических знаков.

Установите в удобном длля вас месте окна Mathcad курсор в виде красного крестика и нажмите на клавиатуре символ Х. введите оператор присваивания :=, взяв его в палитре Оценка (щелкнув по нему).

11

Введите шаблон матрицы, щелкнув в палитре Матрица по символу . Появится окно установки размеров (параметров) матрицы рис. 12.

Рис. 12. Окно установки размеров (параметров) матрицы.

Установите необходимые параметры, в данном случае Ряды – 10, Столбцы – 3, и щелкните ОК. Появится шаблон матрицы (рис. 13).

Рис. 13. Шаблон матрицы.

Введите из таблицы 1 входные сигналы. Ниже запишите выражение Х Т . Для ввода символа Т в степени Х щелкните в панели Матрица по

символу . Далее нажмите на кнопку клавиатуры равно . Появится транспонированная матрица рис. 14.

12

Рис.14. Транспонированная матрица.

Ниже транспонированной матрицы запишите возьмите из палитры Оценка. Появится результат :

А : Х

Т

Х

 

. Стрелку

Стрелку в этом выражении возьмите из палитры Оценка, а после

символа Х

Т

не забудьте ввести знак умножения, нажав на клавиатуре

звездочку * над цифрой 8. После символа Х

Т

появится точка. В противном

 

 

 

 

 

 

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

Далее введите матрицу-столбец У и ниже выражение D : Х

Т

Y .

 

 

 

 

 

 

Появится результат:

Потом ниже введите выражение окончательный ответ: оценка коэффициентов:

B :

A

1

D

 

.

Появится

Для введения символа A 1 введите A и в палитре Матрица щелкните

по знаку .

Следовательно, коэффициенты равны b1 2,b2 3,b3 5.

Эксперимент 2. Оценка параметров модели методом наименьших квадратов при наличии ошибок измерения входных сигналов.

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

13

индикаторов dispay позволяют измерять соответствующие входные сигналы в момент окончания моделирования. Блок gaussian находится в библиотеке Random Generator (рис. 16), генерирует псевдослучайный шумовой сигнал с нормальным распределением и позволяет сформировать входные сигналы с нормальным распределением.

Рис. 15. Модель объекта управления.

Рис. 16. Ввод блоки задания Гауссова шума gaussian.

14

Для получения входных сигналов, содержащих случайные ошибки измерения, установите дисперсию этих сигналов Dх = 0,01 и математические ожидания, равные значениям входных сигналов, указанных в таблице 1. Для этого щелчком правой кнопки по блоку gaussian, вызовите его диалоговое окно (рис. 17) и установите математическое ожидание Mean и дисперсию Standart Deviation.

Рис. 17. Диалоговое окно блока gaussian.

Однако блок dispay фиксирует выходной сигнал в момент окончания моделирования. По умолчанию время моделирования равно 20 секунд. Если не менять время моделирования, то при неизменном математическом ожидании входной сигнал не будет меняться случайным образом. следовательно необходимо изменять время окончания моделирования.

Соберите модель рис. 18.

Рис. 18. Модель для оценки случайного сигнала.

15

Включите процесс моделирования рис. 19 (установите на осциллографе масштабную сетку).

Рис. 19. Выходной сигнал блока gaussian.

Из осциллограммы следует, что процесс моделирования необходимо завершать в точках излома осциллограммы (в вершинах). Следовательно, эти точки надо найти (в данном случае их достаточно получить 10). Это можно сделать, растянув (увеличив) осциллограф и установив время окончания процесса моделирования 5 с (рис. 20).

.

Рис. 20. Определение времени окончания моделирования.

Возьмите эти значения времени из этой осциллограммы и поместите в таблицу 2.

16

Таблица 2. Входные и выходные сигналы модели объекта.

1

2

3

4

5

6

7

8

9

10

t, с

0,2

0,4

0,6

1

1,2

1,4

1,6

1,8

2

2,6

х1

1

0

0

1

1

0

1

2

1

1

х2

0

1

0

1

0

1

1

1

2

1

х3

0

0

1

0

1

1

1

1

1

2

 

 

 

 

 

 

 

 

 

 

 

х11

0,992

-0,004

0,016

0,997

1,019

-0,020

1,002

2,017

0,983

0,996

х21

0,006

1,002

-0,003

1,011

0,007

1,005

0,999

1,021

1,986

1,008

х31

0,002

-0,011

0,994

0,014

1,002

1,005

0,999

1.021

1,006

1.996

у

2

-3

5

-1

7

2

4

6

1

9

 

 

 

 

 

 

 

 

 

 

 

у

1,973

-3,072

5,009

-

7,027

1,969

3,998

6,075

1,040

8,952

1

 

 

 

0,964

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Запустите процесс моделирования. Устанавливайте поочередно соответствующее время окончания (End) моделирования и входные сигналы в соответствии с таблицей 2 и в нее же занесите полученные

входные сигналы

х

, х

21

, х

31

11

 

 

(это измеренные с ошибкой входные

сигналы), а также измеренный с ошибкой выходной сигнал Как и ранее, решение найдите в матричном виде:

у1

.

B ( X

T

X )

1

( X

T

Y )

 

 

 

Находя для этого соответствующие матрицы:

Здесь У1 - матрица выходных сигналов при их измерении с ошибкой. Далее:

17

Следовательно, при наличии ошибок измерения входных сигналов коэффициенты имеют вид:

b 1,96, b

2,989, b

5,015.

1

2

3

 

При наличии ошибок измерения входных и выходного сигналов коэффициенты имеют вид:

b 1,99987, b

2,9994, b

4,9999.

1

2

3

 

5. Сохранение работы.

На диске D создайте папку «Идентификация… Лабораторная работа № 4». В меню «File» выберите пункт «Save as…». В открывшемся окне наберите имя диаграммы «Лабораторная работа № 4. Группа № .Студент Иванов И.И.». Указать место сохранения диаграммы на диске D созданная вами папка. После того, как все сделано, нажать кнопку «ОК». Диаграмма сохранится на диске в виде файла с расширением .vsm.

6. Отчет и защита работы.

1.Отчет должен содержать:

титульный лист;

цель и задачи работы;

полученные данные;

выводы.

2.Защита работы включает доклад студента и его ответы на вопросы по теме лабораторной работы.

18

Контрольные вопросы

 

1. В чем заключается сущность

метода активного

эксперимента ?

 

2.Что понимается под оцениванием параметров модели исследуемого объекта ?

3.Какие основные требования предъявляются к процедуре оценивания?

4.Какие свойства оценок параметров, получаемых на основе статистических данных, принято выделять?

5.Какие условия необходимо соблюсти, чтобы МНК-оценки были несмещенными, эффективными и состоятельными?

6.В чем математическая сущность метода наименьших квадратов?

7.Запишите формулу определения коэффициентов при оценки модели методом наименьших квадратов в матричной форме.

19

ОТЧЕТ

Лабораторная работа № 4

Идентификация статических объектов методами активного эксперимента. Оценка параметров модели

методом наименьших квадратов

студента гр.__________________

__________________________________

__________________________________

__________________________________

«___»_____________20___г.

Результаты экспериментов

Эксперимент 1. Оценка параметров модели методом наименьших квадратов при отсутствии ошибок измерения

1.Приведите модель исследования.

2.Приведите таблицу с входными и выходными сигналами.

3.Запишите все расчетные матрицы.

4.Сделайте вывод об оценках полученных параметров.

Эксперимент 2. Оценка параметров модели методом наименьших квадратов при наличии ошибок измерения входных и выходного сигналов.

1.Приведите модель исследования.

2.Приведите таблицу с входными и выходными сигналами.

3.Запишите все расчетные матрицы.

4.Сделайте вывод об оценках полученных параметров.

Преподаватель:

______________________________

«_____»___________________20____г.

20