- •Министерство общего и профессионального образования российской федерации
- •Передаточная функция звена.
- •Весовая функция звена.
- •Переходная функция звена.
- •1. Типы звеньев систем автоматического управления и их характеристики
- •1.1 Позиционные звенья и их характеристики
- •Апериодическое звено второго порядка.
- •Колебательное звено.
- •1.2 Интегрирующие и дифференцирующие звенья и их характеристики
- •2. Описание пакета matlab V. 3.5g
- •3. Содержание лабораторной работы
- •4. Пример выполнения лабораторной работы
- •Приложение Варианты заданий к лабораторной работе.
- •Литература
3. Содержание лабораторной работы
1. Решить дифференциальное уравнение, соответствующее данному звену.
2. Построить переходную функцию звена.
3. Построить частотную характеристику звена.
4. Пример выполнения лабораторной работы
Пусть дано колебательное звено с параметрами к=2,T=1.5 и =0.7.
1. Решение дифференциального уравнения, соответствующего данному звену.
Дифференциальное уравнение колебательного звена имеет вид:
(Т2р2+2Тр+l)X2(t)=kXl(t) .
Принимая x1(t)= [1 + sign(t)] /2 решаем уравнение относительноx”(t)
и сделав подстановку y2 =x2’ , y1=x2 ,получим систему обыкновенных дифференциальных уравнений:
Решение у2 и есть искомая переходная функция звена.
Получим выражения для построения частотной характеристики:
2. Построение графика зависимости выходной функции звена от входной
Создаем файл описания полученной системы уравнений MYFUN.M:
function yprime = myfun(t,y);
k=-2;
Т=1.5;
z=0.7;
yprime = [y(2); (k-2 . * z. * Т. * у (2) - y(1))./T./T];
Во второй строке видно задание правых частей системы дифференциальных уравнений. Создадим файл LAB1 .М:
t0-= 0;
tfinl = 15;
y0=[0 0] '
tol= 1.e-3;
trace = 1;
k=2;
T=1,5;
z=0.7;
[t,v] =ode23('mvfun'.t0.tfinal,y0,tol,trace);
plot(t,y(.;1));
title('Time history');
xlabel(‘time');
pause
Теперь, запустив MATLAB, вводим lab1 и получаем график выходной
функции звена (рис. 4.1)
3. Построение частотной характеристики звена
Дописав в конец файла LAB 1 .М следующий текст MATLAB построит два заданных графика (рис. 4.1 и рис. 4 2)
for 1=1:70,
w=(i-l)./20;
znam=(1-T^2.*w.^ 2).^ 2+(2.*z *T *w)^2;
u(i)=(k-k.*T^2.*w.^2)./znam;
v(i)=((-k).*z.*T.*w)./znam;
end;
plot(u,v);
xlabel('U(w)');
ylabel('V(w)');
На данном годографе передаточной функции U(0)=2 V(0)=0
U(+)=0, V(+)=0 '~
5. Рекомендации по выполнению лабораторной работы
Некоторые типы звеньев нельзя описать в MATLAB, т. к. их решением является 8-функция (рис. 5.2). Идеальное дифференцирующее звено - пример такого звена. Его решение 5-функция такая, что
Другое подобное звено - идеальное звено с введением производной. Его
pешение (рис. 5.3) есть функция Хевисайда (рис. 5.1) плюс-функция.
Для описания инерциального дифференцирующего звена можно вместо - функции (эквивалентной в данном случае X1 ‘) подставить функцию
Достаточная степень точности
обеспечится импульсом длительностью t=0.2 единицы времени и амплитудой А =5 единиц.
Приложение Варианты заданий к лабораторной работе.
|
ВАРИАНТЫ
| ||||||||||||
I
|
2
|
3
|
4
|
5
|
6
|
7
|
8
|
9
|
10
|
11
|
12
|
13
| |
k
|
2
|
0.5
|
1
|
3
|
2
|
3
|
1
|
5
|
4
|
3
|
2
|
4
|
1
|
Т
|
0.5
|
2
|
4
|
1
|
1
|
0.8
|
1.2
|
1.6
|
0.6
|
1
|
0.7
|
1.2
|
0.9
|
|
|
|
|
|
10
|
5
|
3
|
6
|
6
|
0.1
|
0.5
|
0.3
|
0.2
|
|
Апериодическое
|
Апериодическое 2 порядка |
Колебательное
|
|
В А РИА Н Т Ы
| ||||||||||||
|
14 -15
|
16
|
17
|
18
|
19
|
20
|
21
|
22
|
23
|
24
|
25
|
26
| |
k
|
1
|
0.5
|
10
|
3
|
2
|
3
|
3
|
1
|
4
|
2
|
5
|
11
|
6.5
|
Т
|
|
|
|
|
1
|
3
|
0.5
|
5
|
3.5
|
5
|
10
|
10
|
8
|
|
Идеальное интегрирующее
|
Инерциальное интегрирующее
|
Инерциальное дифференцирующ
|