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

TAU_1 / Новая папка / Моделировние систем управления с помощью МАТЛАВ

.pdf
Скачиваний:
71
Добавлен:
10.02.2015
Размер:
1.22 Mб
Скачать

y – выход в момент t

t – моменты времени,

Начальные

T: – вектор состояния

в которые вычисляется

условия

x(t) – состояние в момент t

реакция

(по выбору)

 

u =вход

 

[y,T,x]=lsim(sys,u,t,x0)

 

Рис. 3.5 Функция lsim для вычисления состояния и выходной переменной

>>A=[0 -2;1 -3]; B=[2;0]; C=[1 0]; D=[0];

>>sys=ss(A,B,C,D);

>>x0=[1 1];

>>t=[0:0.01:1];

>>u=0*t;

>>[y,T,x]=lsim(sys,u,t,x0);

>>subplot(211), plot(T,x(:,1))

>>xlabel('Time (sec)'), ylabel('x_1')

>>subplot(212), plot(T,x(:,2))

>>xlabel('Time (sec)'), ylabel('x_2')

Модель в переменных состояния

Начальные условия

Нулевой входной сигнал

Рис. 3.6. Вычисление с помощью функции lsim временных характеристик при ненулевых начальных условиях и отсутствии входного сигнала

41

4. ЗАДАНИЯ ДЛЯ САМОСТОЯТЕЛЬНОЙ РАБОТЫ

1. Рассмотрите два полинома:

p(s) = s2 + 2s+1, q(s) = s+1.

Спомощью MATLAB вычислите следующие:

a)p(s)q(s)

б) полюсы и нули G(s)= qp((ss))

2. Рассмотрите систему с обратной связью, изображенную на рис. 3.1

 

 

Регулятор

Объект

Y (s)

+

1

X (s)

 

 

s+ 2

 

s+1

s+3

Рис .3.1 Система с отрицательной обратной связью

а) Вычислите передаточную функцию замкнутой системы с помощью функций

MATLAB series и feedback.

б) Определите реакцию системы на единичное ступенчатое воздействие с помощью функции step и убедитесь, что конечное значение выходной переменной равно 2/5.

3. Дано дифференциальное уравнение x+ 4x+ 4x = y, где x(0)= x(0)= 0 и y(t) есть единичная ступенчатая функция. Получите решение этого уравне-

ния аналитически и проверьте результат с помощью MATLAB, одновременно построив график x(t) с помощью функции step.

4. Система управления положением спутника по одной координате может быть представлена структурной схемой, изображенной на рис. 3.2. Переменные k, a иb являются параметрами регулятора, а J есть момент инерции спутника.

Примите следующие значения J =10,8E + 08, k =10,8E + 08, a =b =8.

42

Модель

вращения Регулятор спутника

θ

d

(t)

+

k(s+ a)

1

θ(t)

 

 

 

 

Заданное

s+b

Js2

Действительное

положение

 

 

 

положение

Рис. 3.2 Структурная схема системы управления положением спутника по одной координате

а) Напишите программу MATLAB, вычисляющую передаточную функцию замкнутой системы W (s)= θ(s)θd (s).

б) Вычислите и постройте график реакции системы на ступенчатое изменение

входного сигнала величиной 10o.

в) Точное значение момента инерции спутника в общем случае неизвестно и может медленно изменяться во времени. Сравните реакции системы в случаях, когда J уменьшается на 20% и на 50%. Параметры регулятора при этом остаются неизменными. Проанализируйте полученные результаты.

5. Рассмотрите структурную схему системы, представленную на рис. 3.3.

Y (s)

X (s)

Рис. 3.3 Структурная схема многоконтурной системы управления

а) С помощью MATLAB упростите структурную схему системы и определите передаточную функцию замкнутой системы.

б) С помощью функции pzmap определите положение полюсов и нулей передаточной функции замкнутой системы.

в) С помощью функций pole и zero вычислите точные значения полюсов и нулей передаточной функции замкнутой системы и сравните с результатами с данными, полученными в п.б).

6. Связь между выходом y(t) и входом x(t) нелинейной системы определяется уравнением

43

y(x) = x2 + xsinx.

Линейная аппроксимация этой зависимости имеет вид: y = ax,

где a – коэффициент подлежащий определению. Определите параметр a экспериментально, разработав программу MATLAB, которая вычисляет и строит график разности между y и y в зависимости от переменного параметра a. Па-

раметр a подбирайте та, чтобы максимальная разность между y и y не пре-

вышала значения 20 при 0x 10. Когда будет найдено соответствующие значение a, построите графики y и y в диапазоне 0x 10.

7. Система имеет передаточную функцию

X((s)) = (15z)(s + z)

Y s s2 +3s+15

Постройте график реакции системы на сигнал y(t) в виде единичной функции при значениях параметра z =3, 6 и 12.

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

а) W (s)=

1

, б) W (s)=

3s2 +10s+1

, в) W (s)=

s+14

.

 

 

 

 

 

s+10

s2 +8s+ 5

 

s3 +3s2 + 3s+1

 

 

 

 

9. С помощью функции tf определите передаточные функции для систем, модели которых в переменных состояния представлены следующими матрицами:

а) A =

0

1

 

 

1

 

 

 

,B =

,C=[1 0];

 

2

4

 

 

1

 

 

 

1

1

 

0

 

 

0

б) A = −2

0

 

4

,B = 0 ,C=[0 1 0];

 

 

6

2

 

 

 

 

 

 

 

10

 

1

в) A =

 

0

1

 

,B

0

,C=[1 1].

 

1

 

 

=

 

 

2

 

 

1

 

10. Рассмотрите систему, заданную уравнениями

44

 

0

1

0

0

X = 0

0

1 X+

0 u, y =[1 0 0]X.

 

3

2

 

 

 

5

1

а) С помощью функции tf определите передаточную функцию Y(s)U(s).

б) Постройте график реакции системы на начальные условия X(0)=[0 1 1]T

для 0t 10.

в) С помощью функции expm вычислите переходную матрицу состояния и определите X(t) при t =10 для начальных условий из п. (б). Сравните получен-

ные результат с данными их п. (б).

11. Рассмотрите две системы:

 

 

0

1

0

 

 

 

0

 

X

=

0

0

1

X

+ 0 u, y =[1 0 0]X

1

 

 

 

 

 

1

 

 

1

 

 

4

5

 

 

 

 

 

 

 

 

8

 

 

4

 

и

 

0,5

0,5

0,7071

 

0

X2 =

0,5

0,5

0,7071 X2

+ 0 u, y =[0,7071 0,7071 0]X2.

 

 

0,7071

8

 

 

 

6,354

 

 

4

а) С помощью функции tf определите передаточную функцию Y(s)U(s) для

системы (1).

б) Проделайте то же самое для системы (2)

в) Сравните результаты пп. (а) и (б) и прокомментируйте их.

12. Рассмотрите замкнутую систему, представленную на рис. 3.4.

 

 

Регулятор

Объект

Y (s)

+

1

X (s)

 

 

1

 

s+ 2

s2 + 2s+ 4

 

 

Рис. Замкнутая система управления

а) С помощью MATLAB получите модель в переменных состояния для регулятора.

б) Проделайте то же самое для объекта.

45

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

13. Рассмотрите систему

 

0

1

0

X =

2

X+ u, y =[1 0]X

 

3

1

при начальных условиях X(0)=[1 0]T.

С помощью функции isim вычислите и изобразите графически реакцию переменных x1(t) и x2 (t) при u(t) = 0.

ЛИТЕРАТУРА

1.Дьяков В.П., Абраменкова И.В. MATLAB 5.0/5.3. Система символьной математики. М.: Нолидж. – 1999 г. 640 с., ил.

2.Медведев В.С., Потемкин В.Г. Control System Toolbox. MATLAB для сту-

дентов / под общей ред. В.Г. Потемкина – М.: ДИАЛОГ-МИФИ. – 1999. – 287 с.

3.Методы классической и современной теории автоматического управления: Учебник в 3-х т. Т.1: Анализ и статистическая динамика систем автоматического управления / Под ред. Н.Д. Егупова. – М.: Изд-во МГТУ им. Н.Э. Баумана, 2000. – 748 с.

46