
Образец контрольной работы 2 по ТАУ
.pdf
3.7 |
Ф(j w) |
амплитудно-частотной характеристики (АЧХ) замкнутой системы: |
|
|||||||||||||||
Ф |
|
(w) |
= |
Ф |
вещ |
(w)2 |
+ Ф |
мним |
(w)2 |
|
|
|
|
|
|
|
|
|
мод |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
K K K |
éK K K + 1 − |
(T T + T T + T T ) |
w2ù2 |
+ é(T + T + T ) |
w − T T T w3ù2 |
|||||||||
Фмод(w) |
= |
1 2 3 |
ë 1 2 3 |
1 2 |
1 3 |
2 3 |
û |
ë 1 |
2 |
3 |
1 2 3 |
û |
||||||
|
|
|
|
|
|
A6 w6 + A4 w4 + A2 w2 + A0 |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
упростим выражение для числителя: |
|
|
|
|
|
|
|
|
||||||||||
B |
= T 2 |
T |
2 T |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
|
1 |
2 |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
B4 = (T1 T2 + T1 T3 + T2 T3)2 − 2 T1 T2 T3 (T1 + T2 + T3) |
|
|
|
|
|
|
|
|||||||||||
B2 = (T1 + T2 + T3)2 − 2 (K1 K2 K3 + 1) (T1 T2 + T1 T3 + T2 T3) |
|
|
|
|
|
|
||||||||||||
B0 = (K1 K2 K3 + 1)2 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Фмод(w) = |
K1 K2 K3 |
B6 w6 + B4 w4 + B2 w2 + B0 |
|
|
|
|
|
|
|
|||||||||
|
A6 w6 + A4 w4 + A2 w2 + A0 |
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Построим график амплитудно-частотной характеристики (АЧХ) замкнутой системы на |
||||||||||||||||||
интервале 0-20 рад/с (из соображений наглядности графика), по вертикальной оси |
|
|||||||||||||||||
будем откладывать амплитуду, а по горизонтальной оси - частоту: |
|
|
|
|
||||||||||||||
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3.5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2.5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Фмод(w) |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
1.5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
4 |
|
|
8 |
|
12 |
|
|
16 |
|
20 |
|
|
|
|
|
|
|
|
|
|
|
|
w рад/с |
|
|
|
|
|
|

4. Используя полученные характеристики и построенные графики, найдем и оценим |
|
||||||
следующие показатели качества системы: |
|
|
|||||
4.1 eст - статическую ошибку при подаче на вход единичного ступеньчатого |
|
||||||
воздействия: |
|
|
|
|
|
||
eст = |
|
1 |
= |
1 |
= 0.055 |
|
|
1 |
+ K1 K2 K3 |
1 + 12 0.749 1.9 |
|
|
|||
|
|
|
|
|
|||
4.2 wc - частоту среза системы , запасы устойчивости системы по амплитуде Lзап и фазе |
|||||||
фзап: |
|
|
|
|
|
|
|
частота среза системы определяется по графику ЛАЧХ как частота, на которой |
|
||||||
коэффициент усиления равен 0 дБ (корень уравнения ЛАЧХ(w) = 0): |
|
||||||
|
|
|
|
|
30 |
|
|
|
|
|
|
|
24 |
|
|
|
|
|
|
|
18 |
|
|
|
|
|
|
|
12 |
|
|
ЛАЧХ(w) |
|
|
6 |
|
|
||
|
|
|
|
|
|||
дБ |
|
0.1 |
|
|
1 |
10 |
100 |
|
|
|
|
|
− 6 |
|
|
|
|
|
|
|
− 12 |
|
|
|
|
|
|
|
− 18 |
|
|
|
|
|
|
|
− 24 |
|
|
|
|
|
|
|
− 30 |
|
|
wc = 6.591 рад/с |
|
|
|
w,рад/с |
|
||
|
|
|
|
|

запас устойчивости системы по фазе фзап показывает, на какое значение ФЧХ разомкнутой системы на частоте среза отличается от п (180 град):
180
108
36
ф(w)
град |
0.1 |
1 |
10 |
100 |
− 36
− 108
− 180
w,рад/с
фзап = 180 − ф(wc) = 180 − −163.643 = 16.357 град
критическая частота определяется из графика ФЧХ как частота, на которой фаза разомкнутой системы равна п (180 град):
wкр = 9.393 рад/с
запас устойчивости системы по амплитуде (усилению) Lзап показывает, во сколько раз
нужно увеличить коэффициент усиления, чтобы система оказалась на границе
устойчивости
Lзап |
= |
1 |
= |
1 |
= 2.008 |
|
Wмод(wкр) |
0.498 |
|||||
|
|
|
|
20 log(Lзап) = 6.056 дБ
4.3 M - показатель колебательности системы:
показатель колебательности системы определяется как максимальное значение АЧХ замкнутой системы (обычно стремятся, стобы показатель колебательности не превышал двух)
M = max(Фмод(w)) = 3.648

4.4 tp - время регулирования и переругулирование G:
определим переходную характеристику системы как обратное преобразование Лапласа от:
|
|
|
æ W |
|
(p) ö |
|
|
|
|
|
|
|
|
||||
h(t) = |
− 1ç |
|
замк |
÷ |
|
|
|
|
|
|
|
|
|||||
L |
è |
|
|
|
p |
|
ø |
|
|
|
|
|
|
|
|
||
|
|
− 1é1 |
|
|
|
|
|
|
|
K1 K2 K3 |
|
ù |
|||||
h(t) |
= |
L |
ê |
|
|
|
|
|
|
|
|
|
|
|
|
|
ú |
|
|
|
|
|
|
3 |
+ (T1 |
|
T3 + T2 T3) (p) |
2 |
+ (T1 + T2 + T3) p + K1 K2 K3 |
|
|||||
|
|
|
êp |
|
T1 T2 T3 (p) |
T2 + T1 |
ú |
||||||||||
|
|
|
ë |
|
|
|
|
|
+ 1û |
||||||||
h(t) |
= −0.089 e− 22.5 t + −0.856 cos(6.92 t) |
e− 0.871 t + −0.397 sin(6.92 t) e− 0.871 t + 0.945 |
|
|
определим установившееся значение на выходе системы yуст при воздействии на нее ступеньчатой функции:
yуст = lim h(t) = 0.945 t → oo
построим график переходного процесса:
1.6
1.2
h(t)
l(t) 0.8
0.4
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
|
|
|
|
t |
|
|
|
|
|
длительность переходного процесса определим как время, с момента подачи сигнала до момента времени, когда выходной сигнал не будет отличаться от его установившегося значения не более чем на 5%:
tp = 3.505 с
Определим по графику переходного процесса максимальное значение выходного сигнала:
hmax = 1.55
Определим перерегулирование:
G = |
|
hmax − yуст |
|
= |
|
1.55 − 0.945 |
|
= 64.077 % |
|
|
|
|
|||||
|
yуст |
|
|
0.945 |
|
|||
|
|
|
|
|
|
|
5. Найдем дифференциальное уравнение замкнутой системы, связывающее координаты y и v (полагая f=0):
Воспользуемся передаточной функцией замкнутой системы, для этого представим ее в виде полинома, где символ s эквивалентен операции дифференцирования:
Ф(s) = |
|
|
|
|
|
K1 K2 K3 |
|||
|
T1 T2 T3 s3 + (T1 T2 + T1 T3 + T2 T3) s2 + (T1 + T2 + T3) s + K1 K2 K3 + 1 |
||||||||
|
|
|
|
|
|||||
a3 |
= T1 T2 T3 = 0.05 1.1 0.3 = 0.017 |
||||||||
a2 |
= T1 T2 + T1 T3 + T2 T3 = 0.05 1.1 + 0.05 0.3 + 1.1 0.3 = 0.4 |
||||||||
a1 |
= T1 + T2 + T3 = 0.05 + 1.1 + 0.3 = 1.45 |
||||||||
a0 |
= K1 K2 K3 + 1 = 12 0.749 1.9 + 1 = 18.076 |
||||||||
K = K1 K2 K3 = 12 0.749 1.9 = 17.076 |
|||||||||
Ф(s) = |
|
|
|
|
K |
|
|||
|
a3 s3 + a2 s2 + a1 s + a0 |
||||||||
|
|
|
|
|
|||||
æa s3 + a s2 |
+ a s + a ö y = K v |
||||||||
è 3 |
|
|
2 |
1 |
0ø |
|
|||
a |
|
d3 |
|
y + a |
d2 |
y + a d y + y = K v |
|||
|
dt3 |
|
|||||||
3 |
|
2 |
dt2 |
1 dt |

6. Найдем уравнение состояния замкнутой системы в векторно-матричном виде, в нормальной форме, связывающее координаты y и v (полагая f=0):
6.1 Найдем уравнение состояния замкнутой системы в векторно-матричном виде, связывающее координаты y и v (полагая f=0):
v+ |
|
y(3) |
1 |
y(2) |
2 |
y(1) |
3 |
y |
|
|
. |
K |
. |
K |
. |
K |
|
|
|
T1s+1 |
T2s+1 |
T3s+1 |
x1 |
|||
|
- |
x3 |
|
x3=x2 |
|
x2=x1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
æ |
|
|
|
1 |
|
|
|
K3 |
|
|
|
|
|
|
ö |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
ç |
− |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
÷ |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
T3 |
|
|
|
T3 |
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
ç |
|
|
|
|
|
|
|
|
|
|
|
|
÷ |
|
æ 0 |
ö |
|
|
|
|
|
|||||||
|
|
|
|
ç |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
K |
÷ |
|
|
æ |
1 |
ö |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
ç |
÷ |
|
|
|
||||||||
d X = |
ç |
|
|
|
0 − |
|
|
|
|
2 |
|
÷ |
X + ç 0 |
÷ V |
|
Y = ç |
0 |
÷ |
X |
|||||||||||||
|
|
|
T2 |
|
|
T2 |
|
|
||||||||||||||||||||||||
dt |
|
ç |
|
|
|
|
|
|
|
|
|
|
|
|
÷ |
|
ç |
÷ |
|
ç ÷ |
|
|||||||||||
|
|
|
|
ç |
|
|
|
K1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
÷ |
|
è K1 |
ø |
|
è |
0 |
ø |
|
||
|
|
|
|
ç |
− |
0 |
|
|
|
− |
|
1 |
|
÷ |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
ç |
|
|
|
|
|
|
|
÷ |
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
T |
|
T |
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
è |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
1 |
ø |
|
|
|
|
|
|
|
|
|||
æ |
|
|
1 |
|
|
K3 |
|
|
|
|
|
|
|
ö |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
ç |
− |
|
|
|
|
|
|
|
|
0 |
|
|
÷ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
T3 |
|
|
T3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
ç |
|
|
|
|
|
|
|
|
|
|
|
÷ |
|
|
|
æ |
−3.33 |
|
|
0 ö |
|
|
|
|||||||||
ç |
|
|
|
|
|
|
|
1 |
|
|
|
K2 |
|
÷ |
|
|
|
6.33 |
|
|
|
|||||||||||
ç |
|
|
0 |
|
− |
|
|
|
|
÷ |
= |
ç |
|
0 |
−0.909 |
0.681 |
÷ |
|
|
|
||||||||||||
ç |
|
|
|
|
T |
|
T |
÷ |
ç |
|
÷ |
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
−240 |
|
|
|
|
|
|
|||||||||||||||
ç |
|
|
K1 |
|
|
|
|
2 |
|
2 |
|
|
÷ |
|
|
|
è |
0 |
|
−20 ø |
|
|
|
|||||||||
ç |
|
|
|
|
|
|
|
|
|
1 |
|
|
÷ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
ç |
− |
|
|
|
|
|
0 |
|
|
− |
|
|
÷ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
T |
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
è |
|
|
1 |
|
|
|
|
|
|
|
1 |
|
ø |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6.2 Найдем уравнение состояния замкнутой системы в нормальной форме, связывающее координаты y и v (полагая f=0):
vK1K2K3 T1T2T3
+ |
|
y (3) |
y (2) |
y (1) |
y |
|
. |
. |
. |
|
|
|
|
x1 |
|||
|
- |
x3 |
x3=x2 |
x2=x1 |
|
|
|
T1T2+T1T3+T2T3 |
|
|
|
|
|
|
|
|
|
|
|
|
T1T2T3 |
T1+T2+T3 |
|
|
|
|
|
|
|
|
|
|
|
T1T2T3 |
|
|
|
|
|
|
K1K2K3+1 |
|
|
|
|
|
T1T2T3 |
|
|
æ |
|
|
|
0 |
|
|
|
ç |
|
|
|
0 |
|
d |
X = |
ç |
|
|
|
|
|
|
K |
K |
|
K + 1 |
|||
dt |
ç |
|
2 |
||||
|
ç |
− |
1 |
|
3 |
||
|
|
T1 T2 T3 |
|||||
|
|
è |
|
1
0
− T1 + T2 + T3
T1 T2 T3
|
|
|
0 |
|
|
|
ö |
|
|
|
|
1 |
|
|
|
÷ |
|
|
|
|
|
|
|
÷ |
X + |
|
|
T T + T |
|
T + T T |
|||||
|
1 |
÷ |
||||||
− |
1 |
2 |
3 |
2 |
3 |
÷ |
|
|
|
T1 T2 T3 |
|
|
|
||||
|
|
|
|
ø |
|
æ |
|
0 |
|
ö |
|
æ 1 |
ö |
|
ç |
|
0 |
|
÷ |
|
|
||
ç |
|
|
÷ |
V |
Y = ç 0 |
÷ |
X |
|
|
|
|
||||||
ç K1 K2 K3 |
÷ |
|
ç |
÷ |
|
|||
ç |
|
|
|
÷ |
|
è 0 |
ø |
|
T |
T |
T |
|
|
||||
è |
1 |
2 |
3 |
ø |
|
|
|
|
æ |
|
|
|
|
0 |
|
|
|
ç |
|
|
|
|
0 |
|
|
|
ç |
|
|
|
|
|
|
|
|
|
|
K |
K |
K |
|
+ 1 |
||
ç |
|
|
|
|||||
ç |
− |
|
1 2 3 |
|
||||
|
T1 T2 T3 |
|||||||
è |
|
|
|
|||||
æ |
|
|
|
0 |
|
|
ö |
|
ç |
|
|
|
0 |
|
|
÷ |
|
ç |
|
|
|
|
|
÷ |
= |
|
|
K1 K2 K3 |
|||||||
ç |
|
÷ |
|
|||||
ç |
|
|
|
|
|
÷ |
|
|
|
T |
T |
T |
|
||||
è |
1 |
2 |
3 |
|
ø |
|
|
|
|
1 |
|
|
|
0 |
|
|
ö |
|
æ |
|
|
|
ö |
|
|
|
0 |
|
|
|
1 |
|
|
÷ |
|
0 |
1 |
0 |
||
|
|
|
|
|
|
|
|
÷ |
= |
ç |
0 |
0 |
1 |
÷ |
||
|
|
T1 + T2 + T3 |
|
T1 T2 + T1 T3 + T2 |
T3 |
|||||||||||
|
|
|
÷ |
|
ç |
|
|
|
÷ |
|||||||
|
− |
|
|
|
− |
|
|
|
|
÷ |
|
è −1095 |
−87.9 |
−24.2 ø |
||
|
|
T1 T2 T3 |
T1 |
T2 T3 |
|
|
|
|||||||||
|
|
|
|
|
|
ø |
|
|
|
|
|
|
||||
æ |
0 |
ö |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ç |
0 |
÷ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ç |
÷ |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
è 1034 ø |
|
|
|
|
|
|
|
|
|
|
|
|
|
7. Обоснуем, при каком значение периода дискретизации рассматриваемая система эвивалентна импульсной системе:
Рассмотрим разомкнутую импульсную систему, приведенную на передаточная функция ЭЛНЧ которой имеет вид W (s) = WФ (s)W0 (s).
Передаточная функция импульсной системы будет W (z), а частотная характеристика W ( jw). Существует связь между частотными характеристиками ЭЛНЧ W ( jw) и частотной характеристикой импульсной системы W ( jw), которая имеет вид:
W ( jw) = |
1 ∞ |
W é j(w+ rw )ù , |
(1.68) |
|||
|
|
|||||
T rå=−∞ |
||||||
|
ë |
0 û |
|
где, как и ранее, w0 = 2pT - частота дискретизации.
Итак, характеристика W ( jw) получается суммированием смещенных относительно друг друга вдоль оси w на частоту повторения w0 характеристик W ( jw) ЭЛНЧ, умноженных на 1T . Из (1.68) вещественные и
мнимые части частотных характеристик
Re W ( jw) = |
1 ∞ |
Re W é j (w+ rw |
)ù , |
|||
|
|
|||||
T rå=−∞ |
||||||
|
ë |
0 |
û |
|||
|
1 ∞ |
|
|
(1.69) |
||
Im W ( jw) = |
Im W é j (w+ rw )ù . |
|||||
|
|
|||||
T rå=−∞ |
||||||
|
ë |
0 |
û |

Предположим, что вещественная и мнимая частотные характеристики |
|||||||||||
заданы на интервале частоты ω от -ωn до ωn |
( ωn − полоса пропускания) и |
||||||||||
вне этого интервала равны нулю, а спектральные характеристики входного |
|||||||||||
сигнала V ( jω) |
определены на некотором интервале |
−ωc < ω< ωc |
и равны |
||||||||
нулю вне этого интервала. В этом случае, |
если ω |
c |
< ω < ω0 , справедливо |
||||||||
|
|
|
|
|
|
|
|
n |
2 |
|
|
следующее соотношение |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||
|
Y ( jω) = W ( |
jω)V ( jω) = |
1 W ( jω)V |
( jω). |
|
(1.70) |
|||||
|
|
|
T |
|
|
|
|
|
|
|
|
Итак, при |
выполнении |
условия ω |
c |
< ω < ω0 |
|
импульсная система с |
|||||
|
|
|
|
n |
2 |
|
|
|
|
|
|
передаточной функцией W ( |
|
|
|
|
|
|
|
|
|
||
jω) преобразует входной сигнал точно так же, |
|||||||||||
как некоторая непрерывная система с передаточной функцией 1 W ( jω). |
|||||||||||
|
|
|
|
|
|
|
|
|
|
T |
|
Фактически сформулирован аналог известной теоремы Котельникова: |
|||||||||||
если спектр частот входного воздействия ограничен и лежит в диапазоне |
|||||||||||
частот −ωn < ω< ωn , то свойство системы |
с |
АИМ, |
у |
которой |
ω0 > 2ωn |
||||||
тождественны свойствам эквивалентной непрерывной системы |
с АФЧХ |
||||||||||
1 W ( jω). |
|
|
|
|
|
|
|
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
3.5 |
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
2.5 |
|
|
|
|
|
|
|
|
|
|
|
Фмод(w) 2 |
|
|
|
|
|
|
|
|
|
|
|
1.5 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
0.5 |
|
|
|
|
|
|
|
|
|
|
|
0 |
5 |
10 |
|
|
|
|
|
|
15 |
20 |
|
|
|
w, рад/с |
|
|
|
|
|
|
|

Пос графику амплитудно-частотной характеристики (АЧХ) замкнутой системы
определим полосу пропускания как диапазон частот в котором Фмод(w) ³ Фмод(0) :
2
wп = 10.415 рад/с
тогда минимальное значение w0 равно:
w0 = 2 wп = 2 10.415 = 20.831 рад/с
Определим максимальное значение периода дискретизации как:
tд = 2 п = 2 3.14 = 0.301 с w0 20.831
8. Оценим влияние насыщения первого из звеньев системы на ее процессы, для этого воспользуемся специальным программным продуктом - Matlab:
где НЭ - нелинейность первого из элементов системы с характеристикой типа "насыщение":
u(e) = (-c) if e £ -a
ca e if -a < e < a c if e ³ a

|
|
|
20 |
|
|
|
|
|
10 |
|
|
u(e) |
− 20 |
− 10 |
0 |
10 |
20 |
|
|||||
|
|
|
− 10 |
|
|
|
|
|
− 20 |
|
|
|
|
|
e |
|
|
При подаче на вход гармонического сигнала A sin (0.25 t) с амплитудой А=(10, 200, |
|
||||
500) выходной сигнал тем сильнее будет отличаться от входного, чем больше будет |
|
||||
амплитуда: |
|
|
|
|
|
|
|
|
- линейная область. |
|
|