Logo GenDocs.ru

Поиск по сайту:  


Загрузка...

Контрольная работа №2. Диагностика и идентификация систем. Вариант 14 - файл 1.doc


Контрольная работа №2. Диагностика и идентификация систем. Вариант 14
скачать (3199 kb.)

Доступные файлы (1):

1.doc3199kb.02.12.2011 08:21скачать

содержание
Загрузка...

1.doc

Реклама MarketGid:
Загрузка...
Федеральное агентство по образованию

Московский государственный открытый университет

Чебоксарский политехнический институт

Кафедра

Управления и информатики в технических системах

Специальность 220201

Контрольная работа №2
по курсу «Идентификация и диагностика систем»
Вариант № 14


Дата проверки: Выполнила студентка:

Цветкова Н.В.


Результат проверки: Учебный шифр: 607081
Курс: 3 (сокращ.)


Замечания: Проверила: Изосимова Т.А.

2009 год


Содержание.


1. Анализ методов аппроксимации кривых разгона. 2

1.1. Аппроксимация кривых разгона передаточными функциями более высокого порядка (метод Шварца). 5

1.2. Определение передаточных функций объектов регулирования по кривым разгона методом площадей. 7

1.3. Определение передаточных функций объектов регулирования по кривым разгона методом Ротача. 10

2. Практическая часть. Определение по кривой разгона передаточной функции объекта методом М.П. Симою (метод площадей). 11

^

1. Анализ методов аппроксимации кривых разгона.



Аппроксимации графиков переходных характеристик объектов регулирования обычно включает в себя два этапа:

  1. Выбор общей аналитической формулы;

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

В практике чаще всего встречаются следующие объекты:

1) объекты с самовыравниванием без транспортного запаздывания;



Рис. 1.1.
2) объекты без самовыравнивания и без транспортного запаздывания;


Рис. 1.2.

3) объекты обоих видов, но с транспортным запаздыванием.


Рис. 1.3.
Имея экспериментально снятую кривую разгона, можно вычислить передаточную функций объекта регулирования. Известны несколько способов нахождения передаточных функций. Простейший из них основан на аппроксимации переходной функций объекта некоторой кривой, вид передаточной функции которой известен.







Рис. 1.4 Рис. 1.5. Рис. 1.6.
Рассмотрим типовую кривую разгона объекта с самовыравниванием (рис. 1.5.). Проведем к ней через точку перегиба касательную и обозначим отрезок отсекаемый касательной на оси абсцисс, буквой , а отрезок от точки пересечения касательной с линией нового установившегося состояния до буквой , как это показано на рис. 1.5. Тогда кривую разгона можно приближенно заменить экспонентой (рис. 1.7.) с постоянной времени , чистым запаздыванием и установившимся значением равным установившемуся значению кривой разгона. Передаточную функцию такой экспоненты:

, (1.1)

где - коэффициент усиления объекта, приближенно можно считать передаточной фикцией объекта регулирования, имеющего кривую разгона (рис. 1.5).




Рис. 1.7. Рис. 1.8.


Для приближенной аппроксимации кривой разгона объекта без самовыравнивания (рис. 1.6) к кривой разгона проводится касательная. Отрезок, отсекаемый касательной на оси абсцисс, обозначается (рис.1.8). Угол наклона касательной определится из формулы:

, (1.2)
Тогда кривую разгона объекта без самовыравнивания можно приближенно заменить прямой (рис. 1.8), имеющей с осью абсцисс угол наклона и отсекающей на оси абсцисс отрезок . Этой прямой соответствует передаточная функция:

, (1.3)

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


^

1.1. Аппроксимация кривых разгона передаточными функциями более высокого порядка (метод Шварца).



Более точное совпадение кривой разгона и аппроксимирующей кривой дает способ, предложенный Г. Шварцем. Рассматриваются несколько вариантов аппроксимации. Первый случай - объект с самовыравниванием. Передаточная функция представляется в виде:

, (1.4)

то есть модель объекта составляется из одинаковых апериодических звеньев, соединенных последовательно. Значения коэффициента усиления объекта, постоянной времени и показатели степени определяются с помощью графиков (рис. 1.9).


Рис. 1.9.
Для этого берется кривая разгона, приведенная к единичному возмущению. Проводится касательная к ней в точке перегиба и отмечаются отрезки , , (рис. 1.9, а), а также определяется коэффициент усиления . Далее на графике (рис. 1.9, 6) по известному отношению кривой разгона находится показатель степени . Затем по графику (рис. 1.9, в) определяется отношение и вычисляется постоянная времени .

В другом варианте передаточная функция объекта с самовыравниванием представляется в виде:

, (1.5)

где - постоянный коэффициент.

По кривой разгона определяются величины , , в соответствии с рис. 1.10, а. Далее, воспользовавшись графиками (рис. 1.10, б и 1.10, в), находят значение и, а затем постоянную времени . Коэффициент усиления , в предположении, что кривая разгона снята при единичном возмущении, определяется непосредственно по кривой разгона.


Рис. 1.10.
Передаточную функции объекта без самовыравнивания можно представить в виде:

, (1.6)

Значения величин , находятся по графикам (рис. 1.11,б и 1.11, в). Для этого следует по кривой разгона предварительно найти вспомогательные величины , . Коэффициент усиления объекта при единичном скачкообразном возмущении определится угла (рис. 1.11, а) наклона касательной к оси абсцисс.

Другой разновидностью приближенной передаточной функции объекта без самовыравнивания является:

, (1.7)

Постоянная времени и коэффициент в этом случае определяются по графикам (рис. 1.12, а, б).

Рис. 1.11. Рис. 1.12.
Если исследуемый объект регулирования обладает дифференцирующими свойствами и кривая разгона имеет вид (рис. 1.13, а), его передаточная функция может быть приближенно представлена выражением:

, (1.8)

Значения вспомогательных величин , , , , а также , , находятся с помощью графиков (рис. 1.13, б, в, г).


Рис. 1.13.

^

1.2. Определение передаточных функций объектов регулирования по кривым разгона методом площадей.



В ряде случаев точность представления передаточной функции, определяемой приближенными методами, сказывается недостаточной. Более точный способ вычисления передаточных функции по экспериментально снятым кривым разгона был предложен М.Д. Симою и получил название метода площадей. Теоретически этот метод может дать любую точность. Но реально эта точность не может быть выше точности исходной информации, т.е. точности экспериментального определения кривой разгона.

Рассмотрим кривую разгона изучаемого объекта (рис. 1.14).


Рис. 1.14.

Обозначим звездочкой входные и выходные величины, записанные в размерном виде, и представим кривую разгона в безразмерной форме, приняв обозначения:

; , (1.9)

где - выходная величина в безразмерной форме,

- входная безразмерная величина.

Прочерк в квадратных скобках - символ размерности - означает, что данная величина представлена в безразмерной форме. Перестроенная кривая разгона приведена на рис. 1.15.



Рис. 1.15.

В основе метода лежит предположение, что исследуемый объект регулирования может быть описан линейным дифференциальным уравнением с постоянными коэффициентами:

, (1.10)

где , - постоянные коэффициенты.
Передаточная функция объекта, описываемого уравнением (1.10) может быть представлена как:

, (1.11)

или в размерной форме:

, (1.12)
Задача состоит в том, чтобы определить неизвестные коэффициенты , используя для этого систему уравнений (1.13). В этой системе уравнений для всех значений , , а для , .

, (1.13)
Входящие в данную систему уравнений коэффициенты вычисляются по следующим формулам:

. (1.14)
Эти коэффициенты получили наименование «площадей». Для - это действительно геометрическая площадь (рис. 1.15), а для остальных коэффициентов это название условно. В формулах (1.14) введена новая переменная:

, (1.15)

В практике чаще всего встречаются следующие объекты:

1) объекты с самовыравниванием без транспортного запаздывания;

2) объекты без самовыравнивания и без транспортного запаздывания;

3) объекты обоих видов, но с транспортным запаздыванием.


^

1.3. Определение передаточных функций объектов регулирования по кривым разгона методом Ротача.



Надежный и простой метод аппроксимации переходных характеристик объекта (отсутствуют колебательные составляющие).

, (1.16)



Рис. 1.16.
Цель аппроксимации в выборе оптимальных коэффициентов , , , , .

Условимся считать, что значения этих коэффициентов оптимальны, если они обеспечивают совпадение исходные и аппроксимирующие переходные характеристики.

Критерий может быть записан:

,





Простейшее выражение запаздывающее звено:

, (1.17)

ее переходная характеристика:









Геометрический смысл производной











При второе выражение можно заменить:



Получившийся график показан на рисунке пунктиром (рис. 1.16).


^

2. Практическая часть. Определение по кривой разгона передаточной функции объекта методом М.П. Симою (метод площадей).




Задание.
Определить по кривой разгона (рис. 2.1) передаточную функцию объекта.

Рис. 2.1.


Решение.

Для расчета и построения системы автоматического управления технологическим объектом необходимо знать его математическую модель. В инженерной практике в качестве моделей технологических объектов чаще используются передаточные функции. Передаточная функция конкретного объекта управления находится, как правило, из кривой разгона объекта. Кривая разгона, представляющая собой график изменения выходной (регулируемой) величины во времени при подаче на вход объекта ступенчатого воздействия, может быть легко получена опытным путем.

Определение характеристик объектов регулирования по данным экспериментальных исследований называется идентификацией. В данной работе рассмотрен метод М.П. Симою, предложенный в 1956 году, который позволяет определить параметры передаточной функции по кривой разгона объекта.

График (рис. 2.1) представлен в виде переходной характеристики , где . Кривая разгона объекта с самовыравниванием без транспортного запаздывания. Передаточная функция объекта определяется как произведение двух передаточных функций и .

1. На кривой разгона ось абсцисс разбиваем на равные отрезки с интервалом времени . При выборе величины отрезка следует учесть:

а) кривая разгона на участке должна мало отличаться от прямой;

б) чем меньше , тем более точным будет конечный результат. Однако с уменьшением интервала разбиения резко возрастает объем вычислений, а излишне высокая точность вычисления может не дать желаемого результата.




Рис. 2.2.

Примем , тогда









  1. Значения в конце каждого интервала делим на и получившиеся значения заносим в таблицу 2.1.


Таблица 2.1.









1

2

3

4

0

0,000

1,000

0,000

7,5

0,17

0,83

0,32

15

0,43

0,57

0,64

22,5

0,60

0,40

0,96

30

0,73

0,27

1,28

37,5

0,80

0,20

1,60

45

0,87

0,13

1,92

52,5

0,90

0,10

2,24

60

0,94

0,06

2,56

67,5

0,96

0,04

2,87

75

0,98

0,03

3,19

82,5

0,99

0,01

3,51

90

1,00

0

3,83

97,5

1,00

0

4,15


3. Вычисляем значение 3-го столбца и записываем в таблицу 2.1.

  1. Определяем площадь по формуле:

, (2.1)

где - сумма 3-го столбца (1). Сумма 3-го столбца равна 3,63.


5. Функцию перестраиваем в другом масштабе времени, за независимую переменную берем время , значения которой также записываем в таблицу 2.1.
6. Заполняем таблицу 2.2 (столбцы 1 и 2 переписываем из таблицы 2.1).
Таблица 2.2













1

2

3

4

5

6

0,000

1,000

1,000

1,000

1,00

1,00

0,32

0,83

0,680

0,563

0,41

0,34

0,64

0,57

0,360

0,205

-0,08

-0,04

0,96

0,40

0,040

0,016

-0,46

-0,18

1,28

0,27

-0,280

-0,075

-0,74

-0,20

1,60

0,20

-0,600

-0,120

-0,92

-0,18

1,92

0,13

-0,920

-0,121

-1,00

-0,13

2,24

0,10

-1,240

-0,124

-0,97

-0,10

2,56

0,06

-1,560

-0,098

-0,84

-0,05

2,87

0,04

-1,870

-0,082

-0,62

-0,03

3,19

0,03

-2,190

-0,055

-0,29

-0,01

3,51

0,01

-2,510

-0,016

0,14

0,00

3,83

0

-2,830

0,000

0,67

0,00

4,15

0

-3,150

0,000

1,31

0,00

7. Определяем площади и по формулам:

, (2.2)

где - сумма 4-го столбца.

, (2.3)

где - сумма 6-го столбца.
Сумма 4-го столбца равна 1,09, а сумма 6-го столбца равна 0,42.


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

, (2.4)

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

, (2.5)

Практически в этом случае можно выбирать передаточную функцию в более простом виде:

(2.6)

При расчете ограничимся всего тремя площадями: , , . В этом случае имеем:

.
Если в этом случае передаточная функция принимает вид (2.6), то получим

. (2.7)
9. Отрицательные коэффициенты не учитываются, поэтому получаем передаточную функцию вида
k = 8 (для нашего примера)
Для проверки результата снимем для полученной передаточной функции при помощи математического пакета Maple 7 (рис. 2.3).

Рис. 2.3.

Полученная временная характеристика полностью соответствует заданной.
Использованная литература.
Яковлев Ю.С. Локальные системы автоматики: Текст лекций. - Чебоксары: Изд-во Чуваш. ун-та, 1993
Зиатдинова А.А. Идентификация объекта управления с самовыравниванием, без запаздывания методом площадей Симою М.П. в пакете Maple. - Нижнекамск: Нижнекамский хим.-технол. институт.
Комиссарчик В.Ф. Автоматическое регулирование технологических процессов. Учебное пособие. – Тверь.: Тверской госуд-й технич-й университет, 2001


Скачать файл (3199 kb.)

Поиск по сайту:  

© gendocs.ru
При копировании укажите ссылку.
обратиться к администрации
Рейтинг@Mail.ru