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

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

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

ляемую натуральными числами, можно считать f(n) не­ прерывной функцией п. В таком случае имеем уравнение первой степени с конечными разностями. Для соотноше­ ния вида

f(n + p)=F^f(n), f ( n + 1),

.... f(n + p— 1),

я],

(1-1)

или

 

. . .+ fp f(n )= g

(1 - 2 )

f(n + p)+hf(n + p—l ) +

знание р значений последовательности f(ri)

позволяет

определить следующее значение.

 

 

 

 

Эти соотношения являются рекуррентными последо­

вательностями р-го порядка.

Это

же f(n)

можно

рас­

сматривать как непрерывную

функцию от

п.

Это

соот­

ношение является разностным уравнением р-го порядка. Решить рекуррентное соотношение — это значит опре­

делить функцию |f(n)

по р

начальным

значениям /о,

fi> •••>fp-i-

 

 

 

 

 

 

б)

 

 

РАЗНОСТИ В

По определению первая разность, обозначаемая Д,

определяется так:

 

 

 

 

 

A f(n )= /( n + l) - f ( n ) .

 

(1-3)

Последующее применение этого определения приво­

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

 

 

A2 {f (л) ] = А [А/ (п) ] = A {f (п+ 1) —/ (п) ]=

 

 

= f (н + 2 )—2f(n+1) 4-f(n);

 

(1-4)

Д’ [f (п)\=

Д [Д2f (п)} =

Д [f(n +

2) - 2/ (я +

1)+ "I

(1-5)

+ f(n)] =

f(n+3)- Zf(n4 -2) +

Sf(n+ \) - f(nj l

 

Amlf (п)} = Д”» '1 / (п + 1) - Д”»-1/ (я).

(1-6)

Примем обозначения

 

 

 

 

/( л + 1 )= :(1 + Д )/( я );

 

(1-7)

f ^ + 2)= ( l + A ) f ( / i + l ) = (l+ A )*f(*);

(1-8)

/ ( « + я » ) ; = ( 1 - к д ) да/ ( « ) -

(1-9)

Пример 1-1. Для линейного разностного уравнения второго по­ рядка с постоянными коэффициентами

aof(n) + a if(n + l)+ a 2f(n+2) = g (« ),

(ЫО)

где f(n) — неизвестная функция, п — независимая переменная,

коэф­

фициенты do, fli, а2 и функция g заданы.

 

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

Уравнение (1-10), полностью соответствующее вышеназванным определениям, может быть представлено в другой форме, более ясно иллюстрирующей аналогию с дифференциальным уравнением.

Подстановка уравнений (1-7) и (1-8) в (1-10)

позволяет на­

писать:

 

W V (n ) +biAf(n) +bof(n) =g(n).

(1-11)

Коэффициенты b ввязаны с коэффициентами а следующими со­

отношениями:

 

Ьо й0-f- я, -f- Яг!

 

6, = я, + 2аг\

( 1- 12)

Ьг= яг.

Аналогия между уравнением (1-11) и линейным дифференциаль­ ным уравнением второго порядка с постоянными коэффициентами ясна. Др соответствует dP/dnP. Различие заключается в конечной величине независимой переменной п в разностных уравнениях, стре­ мящейся к нулю в дифференциальном уравнении.

в) ВЕКТОРНО-МАТРИЧНЫЕ РАЗНОСТНЫЕ УРАВНЕНИЯ

Рассмотрим вектор F(n), определенный через

f{n)

f ( n + p -

1)

и обозначим через А матрицу ( п Х р )

 

0

0

01

- и

14 «

й через В(п) вектор-столбец

О

В (л )= ;

о

е

Если для данного п функция fp удовлетворяет соот­ ношению (1-2), то F (п) удовлетворяет рекуррентному соотношению

F (/i+ l)= A (n )F (n ) + B(n)

(1-13)

и соответственно уравнение (1-13) называется неодно­ родным уравнением,а уравнение

F(n+l) =A{n) F( n)

(1-14)

называется однородным уравнением, соответствующим уравнению (1-13).

1-2. ПРЕОБРАЗОВАНИЕ КАРСОНА

Дискретное преобразование Карсона F*(p)

функции

J f(n) определяется как

 

СО

 

с [f(x)] = F*(p ) = p f lf{x)\e-v*dx,

(1-15)

9

 

где p= a+ jсо — комплексная величина (параметр преоб­ разования).

Выражение (1-15) можно сравнить с классическим преобразованием Карсона

ПР) = Р \ e-*4{t)dt.

(1-16)

о

 

Сходство между этими двумя выражениями очевид­ но [непрерывная переменная t, соответствует дискретной переменной х, и непрерывная функция f (t) — ступенча­ той функции f(x)].

Интеграл, определяемый выражением (1-15), сходит­ ся к конечному пределу, когда р>оо, где сг0 — радиус сходимости.

15

Приведем

примеры

использования

преобразования Карсона.

1) fW = 1(*),

 

 

 

 

 

 

 

 

где

 

 

 

(

1 при х

0;

 

 

 

1

(х) =

 

 

 

<

0 при

х <

(рис. 1-3)

 

 

 

 

(

0.

 

 

Используя (1-16),

получаем:

 

 

 

 

 

 

 

 

00

 

 

 

 

 

 

F* (р) =

р J 1 e~v*dx — — [<?“ рх]§° = 1•

 

е**,

 

 

о

 

 

 

 

 

2) } (х) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1f ( X

х '

 

 

 

 

 

 

 

 

 

7,1

 

 

 

 

 

 

ч Г Г "

 

 

 

 

 

 

 

J

b

<х<0

 

 

 

 

 

X

 

1Д С

 

 

 

 

 

 

 

0

1 г з ь s

 

Рис. 1-3.

 

 

 

 

 

Рис. 1-4.

где а — действительное

число (рис.

1-4).

 

 

 

+00

 

 

 

 

00

п+\

F * ( p ) = p |

[ J

вах \ e - v * d x = ^ p

j

еапe~Pxd x =

оо

О

я+ 1

 

-|

оо

/1=0

u

 

р

 

 

 

 

= 2 i

р

j ' e~p*dx |=

 

[ c - nP — <?-("+»)?] ==

 

я

 

 

-I

n=.0

 

 

 

ti—О

 

 

 

 

 

 

 

 

00

 

 

 

1—<?-p

 

e? —■1

(1 — <?-J>)^' a — p ) n ==

 

 

(P

^

ev ea

 

 

n=0

 

1 — e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а)

ОСНОВНЫЕ ТЕОРЕМЫ

1. Теорема суперпозиции

J Дискретное преобразование Карсона обладает свой­ ством суперпозиции, т. е. преобразование произведения функции на константу можно записать:

C{af(x)] = aC[f(x)]=aF*(p),

(1-17)

а если

(1-18)

f(x)— ag(x)+bh(x),

то получим:

C{f(x)]=F*(p)=aG*(p)+bH*(p),

16

где

б* (Р) = С [£(.V)]; И* (Р) =С{1г (х)].

2. Теорема смещения

Дискретное преобразование Карсона функции J f(x) обозначается так:

С [(/ (х)] = F* (р) = р J J / (х) e~Pxdx,

О

откуда следует, что

Clf(x + k)]= p Jj(f(x + k)e-p*dx =

О

00

= р ( J f (г)е~р{г~к) dr.

г—к

Считая r = x-\-k, получаем:

+ 00

n—k—1

C[f{x + k)) = p0* j

J /О') e~Prdr—pevh J J f(r)e~Prdr.

Интегрирование от 0 до бесконечности не зависит от обозначения, принятого для переменной. Учитывая это, получаем:

С ( /( * + k)] = er>* F*(p) — 2

J/OO^P j

e-p’-dr'j

 

k-\

 

 

 

■.ePhF *(p )-eP h'£l (l ~e~p)e-Pr l

f (r);

 

C [/ {x + k)\= eP*F* (p) -

1))£ e "

J / (r). (1-19)

 

 

о

 

 

Исходя из этого соотношения следует, что

 

C[f(x+l )]=ePF* (р) - ( е р - 1)/ (0);

(1-20)

C{f(x+2)]=e*PF*(p) — (e2P—eP)f(Q) — (еР— 1)/(1);

 

 

 

 

( 1-21)

С (/(х + 3)]

еЗРр* (р) — (е3р _е2р)f (0) _

 

— (е2Р—

е Р )/(1 > —

( е Р 1)/(

И&учно -

 

2—352

библиотек* СС СР

 

сКссГУоП/иц*3

3.Изображение приращений (разност

C W ( x ) ] = C [ f ( x + l ) - f ( x ) ] = C[f(x+ \)] -

—С [/(х )]= (е?— 1) {F*(p)—/(0)];

 

СШ( х)]= ( е р - 1) {F*(p)- f( 0)]-

(1-23)

В частности,

 

F * (p )= t^ r= T -+ f(°)-

О '24)

Аналогичным образом докажем, что изображение второй разности будет:

С[Д2Д*)]=(еР— 1 )2F*{p)-f(0) (е*— 1]*—Д/(0) (еР— 1) (1-25)

A /( 0 )= f( l) - f (0 ) .

Таким же образом

С,[Д2/(х) ]= ( е Р - 1)3F* (р) —f (0) [еР— 1Р -

—А/(0) (еР— I)2—А2/(0) (еР— 1),

(1-26)

А2/(0) = /(2 ) —2f (1) + /(0 ).

 

 

 

4. Теорема свертки

Можно показать, что если F*i(p)

и F*z(p)

дискрет­

ные преобразования Карсона от fi(x)

и /г(я)> т0

 

F*(p) = F\ (p)F\(p) =

С S fi (m)f2 (x — m)\,

 

m=.0

 

 

ИЛИ

 

 

 

F*(p) = C S

f, (■« — да) f. (да)

(1-27)

m=0

 

 

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

5.Теорема о конечных и начальных значениях

Применение выражения (1-24) позволяет показать,

)

Пт / (*) =

lim e?F* (р)

(1-28)

х~>оо

Р+о

 

18

и что

 

lim f(x) = f (0) =

lim

CV

 

(1-29)

 

ер — 1 F(p)-

 

 

 

р-*00

 

' Г

 

 

 

 

 

 

 

 

 

б. Дискретное преобразование Лапласа

Можно определить дискретное преобразование Лап­

ласа ![Л. 1-6] по формуле

 

 

 

 

 

 

00

• •

 

 

 

 

 

F\ (Р) = V e - p * f (х) = L [f (х)}.

(1-30)

 

 

jc=0

 

 

 

 

Дискретное преобразование

Карсона

имеет вид:

 

 

 

 

+ 00

 

 

 

 

F*(p) = F\(p) =

p f

Jf(x)e~r>*dx =

 

 

оо

лс+1

о

оо

х + \

 

 

 

 

=

S p

[ I f {x)e-v*dx=Yi If(x)p

^e--Pxdx =

 

x=zO

x

*= 0

x

 

00

 

■X

 

00

 

 

 

4

 

 

 

= — S

I /Чх)[г~р(*+1) — e-T>x\— — S i

f (x)e~r>x(er>~ 1).

*= 0

 

 

 

x = 0

 

 

Отсюда соотношение, связывающее эти два преоб­

разования, будет:

 

 

 

 

 

F*(p) = - ( e P - l ) f L(p) = - ? £ ± F * L(p).

(1-31)

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

е р(е Р— 1 ) - \

Дискретное преобразование Карсона имеет то пре­ имущество, что изображение константы есть константа. Это соотношение важно при изучении разностных урав­ нений системы (переходная функция системы есть: пре­ образование Карсона реакции на единичный входной импульс).

Доказательства, касающиеся свойств этого преобра­ зования, очевидно, одинаковы. Отметим, однако, что тесь рема смещения записывается как

L {f (•* + 1)] = ev [Fl(р) — 1(0)},

(1-32)

2*

и

fl(р) =

L- ^ r + f (0)

(1-33)

L If (x +

2)] =

g2p [F*L (p ) -

f (0) -

f (1) e-p]: (1-34)

L [f (x +

3)] =

e*p [P\ (p) -

f (0) -

f (0) -

 

- f ( l ) e - p '- f ( 2 ) e ~ 4

(1-35)

 

б) ИЗОБРАЖЕНИЯ НЕКОТОРЫХ ФУНКЦИИ

Преобразования

выполняются с

использованием

уравнений (1-3) и (1-23):

Af ( * ) = / ( * + ! ) —f(x);

C [ A f ( x ) l = ( e P - l ) [ F * ( p ) - f m

Рассмотрим несколько примеров. Используем при этом табл. 1-1, в которой приведены преобразования

обычных функций.

 

1.

И

Найдем С (а), где а — действительная величина,

f(x)=a и A f{x )= a —а = 0.

Применив уравнение (1-23), получим:

(еР— 1)С (а)—а(еР— 1) = 0 ,

С (а)=а. (1-36)

2.Изображение равномерно нарастающей функции

Пусть f(x) — х (рис. 1-5),

А(х) = х + 1—х = 1.

Согласно уравнению (1-36) (п = 1 )

(ер— l)C (x) = 1,

откуда

0-37)

Таким же образом можно показать, что если f(x) — =х?, то

Д (*2) = + 1) 2 _ Х2= + 1;

С(х2)= 2 С ( х )+ С (1 ),

20

и с помощью вышенайденных результатов

О Д

ер + 1

(1-38)

(ер — I)2 ‘

 

 

Рис. 1-5. Ступенчатая

Рис. 1-6. Ступенчатая

функция f ( x ) —x.

функция f(x) = С Х.

Аналогичным способом показывается, что

О Д

е*р + 4ер + 1

 

(1-39)

(ev— I)3

 

 

еьР +

\\е*р+ \\ер + 1

(1-40)

О Д =

(ер — 1)*

 

3. Изображение экспоненциальной нарастающей функции

Найдем преобразование \{x)— dx, где d — действи­ тельная величина. Если d= 0, то мы получаем преоб­ разование Лапласа (рис. 1-6) от единичной ступенчатой функции

A (dx) = dx+l—dx(d— 1) dx.

Отметим, что разность Ad* имеет тот же характер, что и функция d* (так же, как и для непрерывной экс­ поненциально-показательной функции):

С [Adx]= (d— l)C{d*];

(er>— 1) С [dx]—eP+1 = (d— 1) С (dx) ,

откуда

С (d*) =

ер — 1

(1-41)

еРd

Если d— ea, то

ер — 1

 

С (еах) =

(1-42)

ер е{

 

 

21

Таблица 1-1

Дискретные преобразования Карсона

F*(P)

- у - и ( Р )

ер _ 1 e - av----—

ер

1

1

еР — 1

(^р + 1)

(ер — I)2

(62р + 4й>р + 1)

(е е — I)3

е+ 11е2Р+

(ер — I)4

+11ер+ 1

---------------

1 (ер — I)2

1

( е р - 1 ) 3

1 (еР — 1)п

ер — 1 ер — а

f(х) ДЛЯ X5: 0

и ( х ) — и(х — 1),

Л

или и (х) и (1 — х)

)1 а (х)

1 (х — а)

 

1

 

X

 

X2

 

Примечания

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

X3 —

X4 —

X 1)

2!

 

х (х 1) — 2)

 

3!

 

х — 1) ...'(х — п + 1)

п —положитель-

п\

ное целое число

ах

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

 

ная величина

22

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