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

7 сем / Vse_lektsii_TAU

.pdf
Скачиваний:
0
Добавлен:
26.01.2024
Размер:
13.15 Mб
Скачать

Способы определения значений временных сигналов

1) По описанию ИСАУ во временной области (т.е. по (***)).

Например, для определения выходного сигнала и сигнала ошибки используются рекуррентные соотношения2:

 

l

l

 

y[lT ] wп (l s)T g[sT ] wп (l s)T y[sT ]

 

 

s 0

s 0

(***)

 

 

l

 

 

x[lT ] g[lT ] wп (l s)T x[sT ]

 

 

 

s 0

 

2) На основании разностных уравнений, которые составляются относительно решетчатой функции.

Опр. Разностное уравнение – это соотношение, связывающее решетчатую функцию и ее разности до некоторого порядка n.

y lT

y (l 1)T y lT прямая разность 1-го порядка (прямая первая разность)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 y lT

 

y (l 1)T

y lT

y (l 2)T 2 y

(l 1)T

y lT

прямая

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

разность 2-го порядка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y lT

y lT

y (l 1)T обратная разность 1-го порядка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общий вид разностного уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

y lT , y lT

, 2 y lT ,

, n y lT

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Линейное разностное уравнение:

q n y lT

q

n 1 y lT

q y lT

x lT

n

 

 

n 1

 

 

0

 

 

 

Оно может быть также представлено в виде:

a y (l n)T a

n 1

y (l n 1)T

a y lT

x lT

n

 

 

 

0

 

 

 

Порядок разностного уравнения равен n, если оно содержит решетчатые функции

ylT

иy (l n)T (порядок РУ может отличаться от порядка старшей разности).

Способы получения разностных уравнений для ИСАУ

I способ: описание реальной ИСАУ разностными уравнениями (когда физические процессы в системе так описываются).

II способ: по известной дискретной передаточной функции ИСАУ

Например, по

 

 

 

 

 

 

 

W

*

(p)=

Y* (p)

 

b e pmT b

e p(m 1)T ... b e pT b

 

 

m

m 1

1

0

– дробно-рациональной.

 

 

 

 

 

 

 

р

 

X * (p)

 

ane pnT an 1e p(n 1)T ... a1e pT a0

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

6

Делим в выражении для Wр* (p) числитель и знаменатель на e pnT :

W* (p)

b e p(m n)T ... b e p(1 n)T b e pnT

 

Y* (p)

m

 

 

1

 

0

 

an an 1e pT an 2e 2 pT ... a0e pnT

 

р

 

X * (p)

Составим уравнение относительно Y* (p) и

X * (p) :

 

 

 

 

 

 

A* (p)Y* (p) B* (p)X * (p)

 

 

Применим обратное D -преобразование:

 

 

 

 

 

 

D 1

 

 

D 1

 

 

 

 

A* (p)Y* (p)

 

B* (p)X * (p)

 

 

(согласно свойству 20 (сдвиг аргумента РФ во временной области) D -преобразования:

D{f l m T } e pmT F* (p) )

a y lT

a

y (l 1)T

n

 

n 1

 

 

a0 y (l

n)T

b

x (l n m)T

... b x (l n 1)T

 

m

 

 

1

 

 

 

 

 

 

 

b x (l n)T

 

 

 

 

 

0

 

(учтено, что e p(m n)T e p(n m)T )

 

далее находят

y lT

по

x lT

и значениям y

в предыдущие моменты времени.

 

 

 

 

 

 

 

 

Пример.

Дискретная передаточная функция замкнутой ИСАУ имеет вид:

W

*

(p)=

Y* (p)

 

2e pT

3

 

 

e2 pT 5e pT 1

з

 

G* (p)

 

Определить значения выходного сигнала для l 0,2 , если на входе действует сигнал g(t) g010 (t) .

Решение:

Делим в выражении для W * (p) числитель и знаменатель на e2 pT :

з

 

 

 

 

2e pT 3e 2 pT

 

Y* (p)

 

1 5e pT e 2 pT

G* (p)

 

Составим уравнение относительно Y* (p) и G* (p) :

Y* (p) 1 5e pT e 2 pT G* (p) 2e pT 3e 2 pT

Переходим к оригиналам и составляем разностное уравнение:

y[lT ] 5 y[(l 1)T ] y[(l 2)T ] 2g[(l 1)T ] 3g[(l 2)T ] .

Отсюда значение выходной переменной определяется по рекуррентному соотношению

7

y[lT ] 2g[(l 1)T ] 3g[(l 2)T ] 5 y[(l 1)T ] y[(l 2)T ]

Вычисляем y[lT ] по тактам:

при l 0 : y[0] 0

при l 1: y[T ] 2g[0] 5 y[0] 2g0

при l 2 : y[2T ] 2g[T ]+3g[0] 5 y[T ] y[0] 2g0 3g0 5 2g0 5g0

3) (3-й способ определения значений временных сигналов) По формуле обратного дискретного

преобразования Лапласа:

а) в общем виде:

 

 

 

 

 

 

c j 0

 

 

 

 

 

f l T D 1{F* (p)}

 

1

 

2 F* (p)ep l T dp

 

 

 

 

 

2 j

 

 

 

 

 

 

 

 

c j 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

б) для случая, когда

 

дискретное изображение F* (p)

представляет собой дробно-

рациональную функцию

 

( F* (p)

B* (p)

), знаменатель

которой

имеет простые и

 

*

 

 

 

 

 

 

 

A (p)

 

 

 

 

 

действительные корни p :

 

 

 

 

 

 

 

 

 

 

 

 

f l T

 

 

 

 

 

 

*

(p ) ep (l 1)T

,

(5)

 

 

D 1{F* (p)} B*

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

1 A (p )

 

 

где A* (p)

d

A* (epT ) , n – степень полинома A* (p) .

 

 

depT

 

 

(Для кратных корней формула имеет более сложный вид.)

 

 

в) также можно сначала разложить F* (p)

на простые дроби, а затем по таблицам

искать оригиналы этих дробей.

 

 

 

 

 

8

Лекция № 5 (4 октября 2021)

2.3. Описание ИСАУ в области изображений. Дискретная передаточная функция.

Рассмотрим описание импульсной системы в области изображений на примере типовой структуры:

g(t)

x(t)

x* (t)

u(t)

y(t)

 

 

 

WФ(p)

WН(p)

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

Описание замкнутой ИСАУ в области изображений:

Y (p) W (p) X * (p)

п (*)

X (p) G(p) Y (p)

Данная система уравнений содержит одновременно изображения непрерывных и дискретных сигналов (т.е. она неоднородна) → нельзя разрешить ее относительно одного из изображений. Для того чтобы избавиться от этой трудности, подвергнем эту систему D преобразованию:

По свойству 2) D -преобразования:

преобразование произведения изображений непрерывной и решетчатой функций равно произведению изображений решетчатых функций: D{F1 (p)F2* (p)} F1* (p)F2* (p)

Y* (p) W* (p) X * (p)

 

 

п

(**)

X * (p) G* (p) Y* (p)

 

 

 

 

 

 

 

 

Формальная математическая

операция

D

-преобразования означает, что при описании

системы используются только дискретные значения непрерывных сигналов g(t) и y(t) .

Теоретически это соответствует наличию ИИЭ на входе и выходе системы (показаны на рисунке пунктиром):

1

По системе (**) можно определить дискретные1 передаточные функции:

Опр. Дискретная передаточная функция замкнутой ИСАУ – это отношение изображения дискретного выходного сигнала Y* (p) к изображению дискретного входного сигнала G* (p) при нулевых начальных условиях:

W* (p)=

Y* (p)

 

 

 

W* (p)

 

 

 

 

 

 

п

 

дискретная передаточная функция замкнутой ИСАУ

G* (p)

 

 

W* (p)

з

ННУ

1

 

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

(часто обозначается как K* (p) )

Из (**) Y * (p) Wп* (p) (G* (p) Y * (p))

Y* (p) 1 Wп* (p) Wп* (p) G* (p)

Дискретная передаточная функция разомкнутой ИСАУ (которая состоит только из ИИЭ и ПНЧ):

 

 

Y* (p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W * (p)=

 

W * (p)=D

{W (p)}=D{W (p) W (p)}

 

 

 

*

 

 

 

р

 

п

 

 

 

п

ф

н

 

 

 

 

 

G (p)

 

ННУ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

нет ОС

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P* (p)

*

 

*

 

 

pT

 

– в общем случае это дробь

 

, где P (p)

и Q (p) – многочлены e

 

.

Q* (p)

 

Дискретная передаточная ошибки по задающему воздействию:

 

 

 

 

 

 

 

 

W * (p)=

X * (p)

 

 

 

 

1

 

G* (p)

 

 

W* (p)

 

x

ННУ

1

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

Если ИСАУ имеет более сложную структуру, то для получения передаточных функций W * (p) используются правила структурных преобразований (те же, что и для непрерывных систем).

Т.о., в отличие от непрерывных систем, для которых передаточные функции составляются непосредственно по уравнениям элементов, в ИСАУ дискретные передаточные функции определяются по передаточным функциям непрерывной части или по временным характеристикам.

Способы нахождения W * (p)

Для определенности рассмотрим нахождение дискретной передаточной функции

разомкнутой ИСАУ.

 

 

 

 

 

 

I способ:

 

 

 

 

 

 

 

 

 

 

 

 

, где T 2

 

Wр* (p)=D{Wп (p)} 1

Wп (p j 0l)

 

wп (0)

 

 

 

 

 

 

 

 

 

0

 

 

 

 

Т l

2

 

 

 

 

 

 

 

 

если wп (0) 0

 

 

 

 

1 Слово «дискретный» далее иногда будем опускать

 

 

 

 

 

 

 

 

 

 

 

2

wп (0)=0 , если степень числителя Wн (p) меньше степени знаменателя и импульс s(t)

на входе непрерывной части ИСАУ конечной высоты (оба условия практически всегда имеют место).

II способ:

L 1

 

t lT

D

 

 

 

 

 

W (p) w (t) w [ lT ] W * (p)

 

 

 

 

 

п

п

п

р

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W* (p)= w [ lT ]e plT

 

 

 

 

 

 

 

р

0

п

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

III способ:

 

 

 

 

 

 

 

 

Если W (p) представляет собой дробно-рациональную функцию: W (p)

B(p)

 

, то

 

н

 

 

 

 

н

A(p)

 

 

 

 

 

 

 

 

 

 

 

 

 

A(p) 0

 

 

 

1) в случае, когда ХУ

имеет простые действительные корни pi , i 1,n

, среди

которых нет нулевого, и

а) на выходе импульсного элемента имеет место сигнал типа δ-функции ( Wф (p)=1), можно воспользоваться формулой разложения:

 

n

B(p )

 

epT

 

 

 

W * (p)

i

 

 

 

,

(4)

 

epT epiT

р

i 1

A(p )

 

 

 

 

 

 

i

 

 

 

 

 

где A(p ) dA

p pi

i

dp

 

б) на выходе импульсного элемента – прямоугольные импульсы ( W (p)=1 e pT

):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ф

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

B(pi )

 

e

piT

1

 

 

 

 

 

 

 

 

 

 

W* (p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

epT epiT

 

 

 

 

 

 

 

 

 

р

 

i 1

p A(p )

 

 

 

 

 

 

 

 

 

 

 

 

 

i

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Легко запомнить, т.к.

1 e pT

 

epT

1

 

 

 

 

 

 

p

pepT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) в случае, когда один из полюсов Wн (p) равен нулю ( A(p) pA1 (p) ) и

 

 

а) Wф (p)=1 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W * (p) B(0)

e

pT

B(pi )

 

 

e

pT

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

р

 

A (0)

epT 1

i 2

p A (p )

 

epT epiT

 

 

 

1

 

 

 

 

 

 

 

 

i 1 i

 

 

 

 

 

 

 

 

 

см. (4): A(p) pA1 (p)+pA1 (p)

 

 

 

 

 

 

 

 

 

 

 

 

 

=A1 (0) , A(pi 0) pi A1 (pi )

 

 

 

 

 

A(p1

0) pA1 (p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

Wп (p)

б) W (p)=1 e pT

,

 

 

 

 

 

 

ф

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W * (p) T

B(0)

1

 

 

B(pi )

e

piT

1

 

 

 

 

 

 

 

n

 

 

 

 

 

р

A(0)

epT 1

 

 

i 2

p A(p )

epT epiT

 

 

 

 

 

 

 

 

i i

 

 

 

 

IV способ:

представляют в виде суммы простых дробей (если она дробно-рациональная)

находят Wр* (p) по таблицам, где представлены соотношения между W * (p) и W (p) (см. пример на предыдущей лекции).

Дискретная передаточная функция является периодической функцией вдоль мнимой оси

комплексной плоскости Р с периодом

j 0 :

 

 

 

 

 

W* (p jk )=

 

w lT e p jk 0 lT

 

w lT e plT

 

e j 2

 

kl W* (p)

0

 

 

 

 

 

 

 

 

 

l 0

 

 

l 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

Полоса

0

Im p

0

называется основной.

 

2

 

2

 

Задание W * (p) в основной полосе в силу ее периодичности полностью определяет W * (p).

Полюсы и нули дискретной передаточной функции, лежащие в основной полосе, называют основными.

Если полюсы p1 ,

..., pn

передаточной функции непрерывной части системы таковы, что

0

Im p

0

(т.е.

лежат в основной полосе), то они совпадают с основными

2

 

2

 

 

 

 

 

полюсами Wр* (p) .

2.4. Описание ИСАУ во временной области. Определение значений временных сигналов.

Пусть имеется импульсная система, заданная структурной схемой:

4

Т.к. сигнал на выходе линейного звена с известной весовой функцией определяется на основании интеграла свертки

 

 

 

 

 

 

y( t ) w( )x( t )d ,

 

 

0

 

 

 

можно получить описание замкнутой ИСАУ:

 

 

 

 

 

 

 

 

 

 

 

 

y( t ) wп ( t )x*

( )d

(*)

 

0

 

 

 

 

 

 

 

 

 

 

x( t ) g( t ) y( t )

 

 

 

Подставим в первое уравнение системы (*)

выражение для сигнала на выходе ИИЭ

( x* (t)=x(t) (t sT ) ), поменяв порядок операций интегрирования и суммирования: s 0

y( t ) wп ( t )x( ) ( sT )d

s 0 0

 

 

теорема :

 

=

 

 

 

 

(t) (t )dt ( )

wп ( t sT )x[sT ]

 

 

s 0

0непрер .функ .

формула свертки для дискретных сигналов.

По (**) находим y( t ) в дискретные моменты времени t lT , l 0,1,2... :

 

l

 

y[lT ] w (l s)T x[sT ]

 

 

п

 

 

s 0

 

x[lT ] g[lT ] y[lT ]

 

– описание замкнутой ИСАУ во временной области.

 

(Пишем, что верхний предел сумы равен l , а не , т.к. при s l

wп (l s)T 0 )

(**)

(***)

Применяя к уравнениям (***) D -преобразование, получим описание ИСАУ в области изображений:

Y

 

(p) Wп

(p) X

 

(p)

 

Y

(p)

 

 

W

(p)

 

*

 

 

 

*

 

 

 

*

 

 

 

*

 

 

 

 

 

*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W* (p)=

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

 

 

G* (p)

 

 

W* (p)

X

*

(p)

*

 

 

*

(p)

з

ННУ

1

 

 

 

G (p)

 

Y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Использование дискретных значений непрерывных сигналов g(t) и y(t) эквивалентно тому, что мы ставим ИИЭ на входе и выходе системы (показаны на рисунке пунктиром):

5

Способы определения значений временных сигналов

1) По описанию ИСАУ во временной области (т.е. по (***)).

Например, для определения выходного сигнала и сигнала ошибки используются рекуррентные соотношения2:

 

l

l

 

y[lT ] wп (l s)T g[sT ] wп (l s)T y[sT ]

 

 

s 0

s 0

(***)

 

 

l

 

 

x[lT ] g[lT ] wп (l s)T x[sT ]

 

 

 

s 0

 

2) На основании разностных уравнений, которые составляются относительно решетчатой функции.

Опр. Разностное уравнение – это соотношение, связывающее решетчатую функцию и ее разности до некоторого порядка n.

y lT

y (l 1)T y lT прямая разность 1-го порядка (прямая первая разность)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 y lT

 

y (l 1)T

y lT

y (l 2)T 2 y

(l 1)T

y lT

прямая

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

разность 2-го порядка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y lT

y lT

y (l 1)T обратная разность 1-го порядка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общий вид разностного уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

y lT , y lT

, 2 y lT ,

, n y lT

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Линейное разностное уравнение:

q n y lT

q

n 1 y lT

q y lT

x lT

n

 

 

n 1

 

 

0

 

 

 

Оно может быть также представлено в виде:

a y (l n)T a

n 1

y (l n 1)T

a y lT

x lT

n

 

 

 

0

 

 

 

Порядок разностного уравнения равен n, если оно содержит решетчатые функции

ylT

иy (l n)T (порядок РУ может отличаться от порядка старшей разности).

Способы получения разностных уравнений для ИСАУ

I способ: описание реальной ИСАУ разностными уравнениями (когда физические процессы в системе так описываются).

II способ: по известной дискретной передаточной функции ИСАУ

Например, по

 

 

 

 

 

 

 

W

*

(p)=

Y* (p)

 

b e pmT b

e p(m 1)T ... b e pT b

 

 

m

m 1

1

0

– дробно-рациональной.

 

 

 

 

 

 

 

р

 

X * (p)

 

ane pnT an 1e p(n 1)T ... a1e pT a0

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

6

Делим в выражении для Wр* (p) числитель и знаменатель на e pnT :

W* (p)

b e p(m n)T ... b e p(1 n)T b e pnT

 

Y* (p)

m

 

 

1

 

0

 

an an 1e pT an 2e 2 pT ... a0e pnT

 

р

 

X * (p)

Составим уравнение относительно Y* (p) и

X * (p) :

 

 

 

 

 

 

A* (p)Y* (p) B* (p)X * (p)

 

 

Применим обратное D -преобразование:

 

 

 

 

 

 

D 1

 

 

D 1

 

 

 

 

A* (p)Y* (p)

 

B* (p)X * (p)

 

 

(согласно свойству 20 (сдвиг аргумента РФ во временной области) D -преобразования:

D{f l m T } e pmT F* (p) )

a y lT

a

y (l 1)T

n

 

n 1

 

 

a0 y (l

n)T

b

x (l n m)T

... b x (l n 1)T

 

m

 

 

1

 

 

 

 

 

 

 

b x (l n)T

 

 

 

 

 

0

 

(учтено, что e p(m n)T e p(n m)T )

 

далее находят

y lT

по

x lT

и значениям y

в предыдущие моменты времени.

 

 

 

 

 

 

 

 

Пример.

Дискретная передаточная функция замкнутой ИСАУ имеет вид:

W

*

(p)=

Y* (p)

 

2e pT

3

 

 

e2 pT 5e pT 1

з

 

G* (p)

 

Определить значения выходного сигнала для l 0,2 , если на входе действует сигнал g(t) g010 (t) .

Решение:

Делим в выражении для W * (p) числитель и знаменатель на e2 pT :

з

 

 

 

 

2e pT 3e 2 pT

 

Y* (p)

 

1 5e pT e 2 pT

G* (p)

 

Составим уравнение относительно Y* (p) и G* (p) :

Y* (p) 1 5e pT e 2 pT G* (p) 2e pT 3e 2 pT

Переходим к оригиналам и составляем разностное уравнение:

y[lT ] 5 y[(l 1)T ] y[(l 2)T ] 2g[(l 1)T ] 3g[(l 2)T ] .

Отсюда значение выходной переменной определяется по рекуррентному соотношению

7

Соседние файлы в папке 7 сем