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

книги из ГПНТБ / Видаль П. Нелинейные импульсные системы

.pdf
Скачиваний:
6
Добавлен:
24.10.2023
Размер:
10.96 Mб
Скачать

F4P)

(gp~ 1)

(ер — а)2

(gp— о (й?—я)3

(ер — О

(ер —а)п+1

(ер—1)

(еР— а) (ер Ь)

ер — 1

е2Р— 2cos (SeP-)-l

(ep — l)(ep— cosP) g2p — 2cos ре» -j- 1

1

(ep —- 1) (eP — a)

gap

(ep— 1)"

e ~ a p

(ep — 1)”

g - a p (gp — 1)

(ep — b)n+l

Продолжение табл. }-1

f(x) для х > О

Примечания

х ( х — \)ах ~г

2!

х (х — 1) ... (х — п + 1)ах~п

п\

ах Ь* а Ь

sin 0х

sin 8

а— действитель­ ная величина

То же

п — положитель­ ное целое число

а и Ь— дейст­ вительные ве­ личины

cos p

a* — 1

x

(a — l)2

a — 1

( х + а ) ( х + я —1)... (x-j-a—n + 1)

а и n— положи­

n\

тельные це­

 

лые , a ^.n

(x — а) (x a 1 )... X

 

nl

 

^ X - ( x - a - n + 1) ,;(x _д)

 

(x— a ) ( x — a— 1 )..-X n!

Х . . . ( х - а - л + 1)

^

X a ( x - i i ) fex-e-n

 

6— действитель­ ная величина

23

4. Изображение синусоидальной функции

Необходимо вначале найти преобразования двух функций sin рх и cos рх:

Д (sin рх) = sin р(х + 1) —sin fix=

='(cos р— 1) sinpx + sinp cos px;

A(cos px) = co s p(x + 1)—cos px =

= (cos 6— 1)cos px—sin p sin px,

Рис. 1-7. Ступенчатая функция f(x) = sin fk .

откуда

(еР— l)C(sin рх) =

=(cos р— 1) С (sin рх) +

+sin рС(cos рх),

т. е.

(еР—cos p) 0(sin рх) = = sin рС(cos px);

(ep—cos p) C (cos px) — = —sin pC (sin px) + e p—1.

Решая эти два уравнения, получаем:

C(sin px)

sin p (e? — 1)

e2p 2 cos p<?p + 1 *

 

с (и « и = - (^ _ с"с5 У .+ '| > •

Функция sin px изображена на рис. 1-7.

(1-43)

о - 44)

З а м е ч а н и е .

 

Если р=2л, т. е. период синусоиды равен перио­

ду дискретизации,

то

 

 

 

 

 

С (sin 2ях) = С (0) =0;

 

 

 

 

С (cos 2лх) = С (1) = 1.

 

 

 

 

5. Изображение

произведения

Найдем

преобразование

произведения /(х, п) —

=[х(х— 1)

... (х—п+ \)Уп\, когда п принимают разные

значения.

 

 

 

 

 

Для п= 2 находим:

 

 

. х ( х — 1) _ х(х + 1)

х ( х — 1) _

 

а

 

2

2

2

Л '

24

откуда

 

 

 

 

 

 

 

С(х) = (ер-

1) [ f * ( p ) - f ( 0) ]=

(ер _ l)F*(p),

т. е., учитывая уравнение (1-37),

 

 

 

 

 

С ' X ( х — I)

 

(ер — I)2

 

(1-45)

 

 

 

 

 

 

Для произвольного п

 

 

 

 

 

д х (х — 1 )... п 4 -1)

__ -f- 1) х ... (х — п +

2)

 

 

 

 

 

 

п\

 

у — 1) . .. (у - п +

1)

х ...(у — л + 2)

 

откуда

 

п\

 

 

( л - 1 ) !

 

 

 

 

 

 

 

 

С ]/ (*,

п)\ =

С If {х,

п -

1)].

(1-46)

Рекуррентное

соотношение

 

 

 

 

С

у (у — 1) ... (у — /г +

1) 1

__

1

(1-47)

 

п\

 

 

(еР— 1)п ‘

 

 

 

 

 

С некоторыми ограничениями на начальные условия доказано, что, зная дискретное преобразование Карсона некоторой функции, можно легко найти разностное уравнение, связывающее значения этой функции для разных моментов квантования.

Достаточно рассматривать ер как оператор сдвига. Преобразо­ вание Карсона дает тогда непоаредственио уравнение с конечными разностями.

Пример 1-2. Преобразование константы будет (1-36): С (а) —а.

Следовательно, если

у(п)=ах(п),

то при С(х) = ( е Р — I ) - 1 разностное уравнение будет:

(е?— \)у=х,

т. е.

у( п+1)—у(п)=х(п).

Справедливость выражения очевидна из рис. 1-5.

1-3. ОБРАТНОЕ ДИСКРЕТНОЕ ПРЕОБРАЗОВАНИЕ КАРСОНА

Изучение преобразований дискретных функций пока­ зывает, что все эти функции являются функциями ер, а не р. Следовательно, предпочтительнее во время вся­ ких алгебраических выкладок считать в качестве пере­ менной ер, а не р.

25

В пашем распоряжении два метода расчета обратно­ го преобразования Карсона применительно к дискретным функциям: разложение на элементарные дроби и разло­ жение в рядно степеням е*\

а) РАЗЛОЖЕНИЕ НА ЭЛЕМЕНТАРНЫЕ ДРОБИ

Метод состоит в разложении выражения на сумму элементарных дробей, в которых е? является перемен­ ной. Тогда с помощью таблицы дискретных 1преобразований Карсона некоторых функций можно найти ориги­ нал. Объясним этот метод на примере.

Пример 1-3. Нужно найти обратное преобразование функции

 

 

 

 

 

 

 

 

( е *

— 1)

 

С '48)

 

 

 

 

П Р ) =

{ePl a ) 4 J - b )

 

 

При а и 6 действительных числах это выражение можно напи­

сать следующим образом:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г

А

 

,

 

В

,

С

 

где

 

F{p) =

(a v - \ )y i^

zriа)2

еР а

е?

b

 

 

 

 

A

 

 

 

 

 

 

1

 

 

 

 

 

 

 

Г

_ep

1

1

 

 

 

 

 

 

 

 

 

 

 

a - b ’

 

 

 

 

 

 

i

 

— 6

eP-a

 

 

 

 

# =

d

 

 

 

 

1

 

 

— 1

 

 

dev

ep — b

 

eP=a

(ep b)2 JePa

 

 

 

 

 

 

 

 

f(a)

 

 

 

 

C =

Г

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[(ep — a)2 J eP-

(Ъ— а)2

 

 

 

 

 

откуда

( /7)] =

 

с -

 

 

 

 

 

 

 

 

 

 

 

 

(а — 6) (ер — а)2

 

 

 

 

 

 

 

 

 

 

 

 

 

ev

 

 

о 1

г з

ц

е

 

 

 

 

 

 

ер — 1

 

 

Рис. 1-8. Ступенча­

 

 

 

 

Ь)2 (ер а) ~

 

тая

функция

f

(х ) =

 

 

 

 

 

 

 

 

 

 

_ b x ах , хах ~1

 

 

 

+

_____ ^

_____ 1

 

(b — а)2 +

а — b

 

 

 

т

(6 — а)2 (ер — Ь) \

 

 

 

 

 

 

 

 

 

 

 

 

Согласно табл. 1-1 дискретных преобразований Карсона

 

С -, [^(р)] =

/(х) = ^ Г б хаХ~' ~

(a

- W

aX+ (Ь -а)2

Ь*'

 

 

С-'

[F* (/?)]

 

 

f a x

q x

 

XCLx ~1

для Х 3з 0.

( 1-49)

 

 

= (Й_

й)2

+

' а _ ь '

 

Эта

ступенчатая функция показана

на рис.

1-8

при а= 1,

6=0,5,

2(>

б) РАЗЛОЖЕНИЕ в РЯД ПО СТЕПЕНЯМ еР

Метод состоит в разложении в ряд по степеням ер. Наиболее трудная задача состоит в нахождении обрат­ ного преобразования бесконечного ряда в явном виде. Это суммирование составляет вообще трудную задачу, и, исключая отдельные частные случаи, применение ме­ тода затруднено.

Пример 1-4. Нужно найти обратное преобразование функции

 

 

 

? ( Р ) = (еР L а) (ер — ь) '

 

( ! ' 5 0 )

При а и b действительных величинах можно написать:

 

вР-- 1

 

 

1)1* '-*> + ае-*р +

...] \ер + Ье~*р + ...],

- ^ ДЯ=ау (еР — Ь ) ° ( * * -

откуда

F (р) = е ~ р (I е~р)[ 1 - f я<? ~ р -f- а2е ~2р + . •.] X

 

 

 

 

w

 

Ье-р +

Ь2е-*р +

...] =

e - v (1 — ер)

г

аг — Ь2

 

X [1 +

1 +

а'~1 Ге ' р +

 

 

 

Ь3

 

 

 

 

 

 

 

 

 

 

Т е

- 2 р

 

 

 

 

 

 

т. е.

 

 

 

 

 

 

 

 

 

 

1

(х - 1)

а 2 —

Ь2

1 (х ~

2) +

•••

 

 

 

f (х) =

+ ~—

ь

 

 

 

Суммируя ряд, получаем:

 

 

 

 

 

 

 

СО

ап Ьп

 

 

 

 

 

 

/ о о =

(х—п);

 

 

 

 

а Ь

Рис.

1-9.

Ступенча­

 

 

п = 1

 

 

 

 

 

 

 

 

 

 

тая

функция

f (х) =

f(x)

ах — Ьх

для 0 < х.

 

(1-51)

 

_ ах — Ьх

а — b

 

 

а Ь

 

 

 

 

 

 

 

Мы представили на рис. 1-9

f(x)

для частного случая,

когда а=

=0,5,

6=0,3.

 

 

 

 

 

 

 

 

1-4. РЕШЕНИЕ ЛИНЕЙНОГО РАЗНОСТНОГО УРАВНЕНИЯ

Для решения линейного разностного уравнения мож­ но использовать два метода: прямой метод и метод, ис­ пользующий преобразование Карсона.

27

а) ПРЯМОЙ МЕТОД

Доказано [Л. 1-1], что методы, позволяющие интегри­ ровать линейные уравнения, близки к методам решения линейных разностных уравнений. Например, общее ре­ шение линейного разностного уравнения третьего по­ рядка

f (х + 3) +/?/ (% + 2) + qf (х + 1) + г/ (х )— 0 (1-52)

имеет частные независимые решения: }(х)и /(хД, }(х)з, которые должны удовлетворять неравенству

/(0)!

/ (0)2

f(0),

(1-53)

/(1).

/ (1)2

f( Da Ф 0

f (2)i

f(2h

H2)3

 

для линейной комбинации с постоянными коэффициен­ тами.

Например, если р, q, г — постоянные и г ф 0, то корни характеристического уравнения а, р, у получаются из соотношений

/ -f- т) — хт;

(1-54)

х 3 рх2+ qx + г = 0.

Условие независимости решения ах, |3Х и ух запишет­ ся так:

1 1 1

а

Р

Y = (* — РХР — Т)(Т - а) 4= 0.

(1-55)

 

 

 

а2 pa

y 2

 

Величины а, у — различны, и общее решение сле­ дующее:

f (x )= C lax+ C2^ + C3yx,

(1-56)

причем предполагается, что величины f(0), /(1)

и /(2)

заданы как

 

f (0) — Ci + С2+ С3;

 

/ (1) = Сш + Cap + Cjy;

 

/(2 )= С ,а 2+ С 2р2+С3у2.

 

28

Исключая Ci, С2, C’3l получаем:

f U

) ах

fix

Ух

U 0 )

1

1

1

!

«р

У

f

2()

Р2

у2

Пример 1-5. Рассмотрим разностное уравнение

f ( x + l ) —af(x)=c,

где /( 0 )= р ; а, с, р — действительные числа. Общее решение имеет вид:

f(x) =А + Вах.

Определим А и В из двух условий:

с=А + Вах+1аАаВах;

Д0) = р=Л +В .

Получаем общее решение уравнения (1-58):

(1-57)

(1-58)

(1-59)

Схеме рис. 1-10 соответствует разностное уравнение первого по­ рядка. Импульсный элемент К управляется последовательностью импульсов с периодом Т. Напряжение на пластинах конденсатора в предположении а (0 )= 0 для л-го импульса

 

U (t) — \E U (л)][1 — e~th\+ V (л).

(1-60)

Постоянная

времени контура т = RC. Если t= T (n+ 1), то

можно

записать:

 

 

 

и

U (t)=U(n+l)

 

 

 

 

 

U ( n + l ) =

[ E - U (л)][1 -

+ U (л),

 

т. е.

 

 

 

 

U( л + 1)—ail (л) —с;

(1-61)

 

при

 

Рис. 1-10.

с — Е (1 — е ~ т^).

Выражение для напряжения на пластинах конденсатора на л-м периоде будет:

y ( « ) = T e r r (ап — ’ );

£ ( 1

— е ~ Т1' )

П (л) =

(1 — е~(Т>^п) = Е (1 — е ~ (Г1')п ). (1-62)

1

- е ~ т ! '

29

6) МЕТОД, ИСПОЛЬЗУЮЩИМ п р е о б р а з о в а н и е

КАРСОНА

Рассмотрим применение дискретного преобразования Карсона для решения линейного разностного уравнения

с

постоянными

коэффициентами. Решение получаем

в

три этапа:

разностное уравнение переписывается

в форме дискретного преобразования Карсона с учетом начальных условий; полученное уравнение решается алгебраически так, чтобы получить искомую функцию; с помощью обратного дискретного преобразования Кар­ сона находится решение разностного уравнения.

Заметим, что разностное уравнение может быть запи­ сано двумя различными способами: в обычном виде, если вводить величины f(x), f(x+l ) ■■., или с помощью разностей первого, второго и т. д. порядков f(x), Af(x),

Л В Д . . .

Как правило, эти два способа написания в действи­ тельности эквивалентны, но не одинаково удобны в раз­ ных случаях; первый способ решения быстрее. Важней­ шее преимущество по отношению к классическому мето­ ду решения конечно-разностных уравнений заключается в том, что нет необходимости вводить в рассмотрение произвольные постоянные. Для его использования необ­ ходимо иметь таблицу дискретных преобразований Кар­ сона (см. табл. 1-1).

Пример 1-6.

Рассмотрим еще раз

разностное уравнение

(1-58)

 

f { x + l ) —af(x)=c;

 

 

cf( 0)=Ц -

 

 

Предположим, что функция f(x)

имеет изображение

F*(p),

тогда вышеприведенное уравнение запишется:

 

 

Clf(x+l)-aO[f(x)] = C(c).

 

Уравнения (1-20) и (1-36), соответствующие преобразованию

Карсона f(x + l)

и константе, позволяют записать:

 

ePF*(p) — (eP1)/(0)—aF*(p)=c.

При начальных условиях /( 0 )= ц вышеприведенное уравнение запишется как

сер —-1

F*{P) ер — а ^ ер а

С * 1[F* (р)\ = С - '

с

 

ер~\

еР — а

ц .-------------- •

 

г

ер — а

30

Из табл. 1-1

а*

f (х) = с а

р.а*.

(1-63)

 

Это выражение было уже найдено прямым методом. Для этого примера оба метода тождественны по своей трудоемкости, но для разностных уравнений высшего порядка использование преобразова­ ний Карсона позволяет получить результат быстрее и с минимумом выкладок.

Пример 1-7. Решаем уравнение

при

f(x + 2)—5f(x+ 1) + 6f(x) = 4 Х

 

(1-64)

f (0) =

0;

1

 

 

 

 

 

 

 

Получаем:

HI) =

I-

I

 

 

 

e* —1

,

 

 

ер— 1

 

 

,

 

 

 

 

F

e2P5eP -f- 6

+

(c p 4)(e2? — 5e* +

6) ’

откуда

 

 

 

 

 

 

 

 

ev —- \ _ e p — 1

1

e p — l__ep— l

1

<?? — 1

t. e.

ep — 3 ep 2*~

2

ep — 2

<?? — 3~^" 2

e v — 4’

1

 

 

1

r 1

(

 

,

 

 

 

f (x) =

3* — 2* + -g - 2* — 3* +

-g - 4* = —

v4* — 2X). (1-65)

1-5. УСТОЙЧИВОСТЬ РЕШЕНИЯ РАЗНОСТНОГО УРАВНЕНИЯ

а) ВВОДНЫЕ ЗАМЕЧАНИЯ

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

lira /” (х) = const (или равен нулю),

(1-66)

*->00

 

 

где f(x) — ступенчатая функция.

 

 

Для разностного уравнения

 

 

a0f(x + n) +aif(x + п— 1) +

... an(x)=Lx

(1-67)

характеристическое уравнение запишется как

 

ctoXn+ йуХп~1+ ...

+ a n = 1.

(1-68)

Общее решение является линейной комбинацией ча­ стных решений для п корней г, уравнения (1-67):

f(x) = C0+ f i CirX.

(1-69)

М

Определим теперь условия устойчивости решения уравнения (1-67).

б) ОБОБЩЕНИЕ КРИТЕРИЯ РАУСА — ГУРВИЦА НА СЛУЧАИ РАЗНОСТНЫХ УРАВНЕНИИ

Решение уравнения (1-67) устойчиво, если корни г* характеристического уравнения (1-68) по модулю мень­ ше единицы. Это условие

М < 1

(1-70)

следует из общего выражения (1-69) решения уравнения

( 1- 68) .

Критерий устойчивости, изложенный А. Гурвицем и Е. Дж. Раусом, легко применим для рассматриваемого случая. Достаточно только произвести ^-преобразова­ ние, произведя подстановку

\ Wj +

1

(1-71)

П

1 '

'‘ Wi

 

тогда условие устойчивости (1-70) приобретает вид:

Re (гг;,-) <0. (1-72)

Проиллюстрируемвышеизложенное на нескольких примерах.

1. Уравнение первого порядка

Рассмотрим следующее разностное уравнение (нача­ ло координат соответствует точке равновесия):

х ( п + 1) +ах(п) = 0 .

(1-73)

Корень характеристического уравнения будет

г—а,

аусловие устойчивости запишется как

М < 1 ,

(1-74)

что ясно из рассмотрения уравнения (1-73).

 

2. Уравнение второго порядка

а0х(п + 2) +atx(n-f 1) + а 2х(п) =0,

(1-75)

т. е. характеристическое уравнение

(1-76)

aor2+air+az=Q,

22

Соседние файлы в папке книги из ГПНТБ