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

Образец контрольной работы 2 по ТАУ

.pdf
Скачиваний:
71
Добавлен:
15.06.2014
Размер:
411.73 Кб
Скачать

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) выходной сигнал тем сильнее будет отличаться от входного, чем больше будет

 

амплитуда:

 

 

 

 

 

 

 

- линейная область.