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

СММиФ заоч

.pdf
Скачиваний:
44
Добавлен:
21.03.2016
Размер:
746.52 Кб
Скачать

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

z +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z 1)2

 

 

 

 

 

 

 

 

 

 

4

 

 

n2

 

 

 

 

 

 

 

 

 

 

 

 

z(z + 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z 1)3

 

 

 

 

 

 

 

 

 

5

 

 

an

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

nan1

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

(z a)2

 

 

 

 

 

 

 

 

7

 

Cnk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

(z 1)k +1

 

 

 

 

 

 

 

8

 

an sin nτ

 

 

 

 

 

 

 

 

 

 

 

azsin τ

 

 

 

 

 

 

 

 

 

z2 2az cosτ+a2

 

 

 

9

 

an cos nτ

 

 

 

 

 

 

 

z(z a cos τ)

 

 

 

 

 

 

 

z2 2az cos τ+ a2

 

 

10

 

an sh nτ

 

 

 

 

 

 

 

 

 

 

 

az sh τ

 

 

 

 

 

 

 

 

 

z2 2az ch τ+ a2

 

 

11

 

an ch nτ

 

 

 

 

 

 

 

 

z(z a ch τ)

 

 

 

 

 

 

 

 

 

z2 2az ch τ+ a2

 

 

12

T (x)

= cos(n arccos x)

 

 

 

 

 

 

 

 

z(z x)

 

 

 

 

 

 

 

 

z2 2xz +1

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

13

 

 

an

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ea / z

 

 

 

n!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Некоторые свойства z-преобразования

I. Линейность.

Если

fn

 

F* (z) и gn

 

G* (z), то для любых ком-

 

 

плексных постоянных α и β

αf (n) g(n) αF* (n)G* (n).

II. Теорема запаздывания (первая теорема смещения).

Если fn

 

F* (z), то для любого целого k > 0 справедливо пре-

 

образование

 

 

 

31

fnk zk F* (z), k = 1,2,...

III. Теорема опережения (вторая теорема смещения).

Если fn F* (z), то для любого целого k > 0 справедливо пре-

образование

 

 

 

k 1

 

 

 

 

fn+k

 

zk

F* (z)fm zm

.

 

 

 

 

 

m=0

 

IV. Изображение суммы.

 

 

n1

1

 

 

Если fn

 

F* (z), то fm

 

 

 

F * (z).

 

 

 

 

 

z 1

 

 

m=0

 

Для того, чтобы по известному изображению F* (z) найти оригинал fn можно:

1)воспользоваться таблицей, после разбиения дроби F* (z) на сумму простейших дробей;

2)в случае, когда F* (z) есть правильная рациональная дробь отно-

сительно z ,

функцию fn можно найти по формуле

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

fn = res(F * (z)zn1 ),

 

 

 

 

 

 

i=1 zi

 

 

 

 

 

 

 

где сумма берется по всем полюсам функции F* (z).

 

 

 

 

 

Решение разностных уравнений

 

 

 

 

Разностью первого порядка решетчатой функции fn

называется

величина, обозначаемая как

fn , и равная

fn = fn+1 fn .

 

 

 

Разностью второго порядка

2 fn называется величина, опреде-

ляемая как

2 fn = fn+1 fn .

k fn называется величина

 

 

 

Разностью k -го порядка

 

 

 

 

k fn = k 1 fn+1 k 1 fn .

 

 

 

 

 

 

k

 

 

 

 

k!

 

 

 

k fn = (1)m Ckm fn+k m ,

где

Ckm =

 

 

 

.

m!(k m)!

 

m=0

 

 

 

 

32

ТЕОРЕМА z -преобразовании разности). Пусть fn F* (z). Тогда z -преобразование разностей равны

fn = fn+1 (z 1)F* (z)zf0 ,

( ) ( ) k1 ( ) − − ( )

k fn = fn+k z 1 k F* z z z 1 k m 1 m f0 .

m=0

Уравнение вида

F(n, fn , fn+1,..., fn+ k ) = 0,

где f (n) fn – решетчатая функция, называется разностным уравнением k -го порядка.

Замечание. Далее в задачах для решетчатой функции будет использоваться обозначение xn ( fn xn ).

Рассмотрим процедуру решения линейного неоднородного раз-

ностного уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

a0 xn+k

+ a1xn+k1 +... + ak xn =ϕ(n) .

1) Применяем к обеим частям уравнения дискретное преобразо-

вание Лапласа

 

 

 

 

 

 

 

 

 

 

 

 

xn

 

 

 

X * (z),

ϕ(n)

 

Φ* (z)

 

 

 

 

и, учитывая, что

 

xn+1

 

 

z(X * (z)x0 )z ,

 

 

 

 

 

x

n+ 2

 

 

 

 

z2 (X * (z) x

0

x z1 ),

 

 

 

 

 

 

 

 

 

 

 

 

1

. . .

получим линейное алгебраическое уравнение относительно изображения X * (z).

2) Разрешив полученное уравнение относительно X * (z), возвращаемся назад к оригиналу – последовательности. Общее решение будет содержать неопределенные константы x0 , x1,xk1 , которые фиксируются, исходя из начальных условий.

Решение типовых задач

33

Пример 1. Пользуясь определением, найти изображение F* (z) для функции fn = eαn .

Решение.

F * (z) = eαn zn = (eα z1 )n =

1 = z

при z > eReα .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=0

n=0

1

eα

 

z eα

 

 

 

 

 

 

 

 

 

 

z

найти изображение F* (z)

 

 

 

 

 

Пример 2. Пользуясь определением,

для функции fn

= sin n.

Решение.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin n = ein ein

zi

Тогда

sin n

 

 

1

 

 

z

 

 

 

 

 

z

 

 

; ein

 

(z1ei )n =

 

=

 

;

ein

 

 

 

.

 

 

 

 

 

 

1 z

1

i

z e

i

 

z e

i

 

 

n=0

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

1

z

 

 

 

 

 

 

 

i

2i z e

 

z

 

 

=

 

 

 

z e

i

 

 

 

 

z

ei ei

 

z sin1

 

 

 

(z ei )(z ei ) =

 

 

.

2i

z2 + 2z cos1+1

Ответ: sin n

 

 

 

z sin1

. ▲

 

 

z2 + 2z cos1+1

 

 

Пример 3. Найти оригинал для изображения F* (z) =

3z2 z

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение.

 

 

 

 

 

 

 

 

 

 

 

 

(z +1)4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Точка z0 = −1 является полюсом четвертого порядка для функ-

ции F* (z). Тогда оригинал

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fn

= res

F (z)zn1 = res

 

(3z2 z)zn1

= res

 

 

3zn+1 zn

=

 

 

 

(z +1)4

 

 

 

 

 

 

 

 

 

 

 

z0 =−1

 

 

 

z0 =−1

 

 

 

 

z0 =−1

 

 

(z +1)4

 

 

 

 

 

 

 

1

 

d 3

3zn+1 zn

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

n+1

 

n

 

′′′

 

 

 

 

=

 

 

lim

 

3

 

4 (z +1)

 

lim (3z

z

)

=

 

 

 

 

 

 

(z +1)

 

=

 

 

 

 

 

 

 

 

 

1

 

 

3! z →−1 dz

 

 

 

 

1

 

 

 

6 z →−1

 

 

 

 

 

 

 

 

 

 

 

=

lim (3(n +1)zn nzn1)′′ =

lim (3(n +1)nzn1 n(n 1)zn2 )′ =

 

 

 

 

 

6 z →−1

 

 

 

 

 

 

 

6 z →−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

1

lim n(n 1) (3(n +1)zn2 (n 2)zn3 )=

(1)n

 

n(n 1) (4n +1).

 

 

 

 

 

6 z →−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ:

fn

=

(1)n

n(n 1)(4n +1).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

34

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 4. Найти оригинал для изображения F * (z) =

z

.

(z 1)(z 2)2

 

 

Решение.

 

 

 

 

 

 

 

 

 

Вынесем z в числителе изображения F* (z) за скобки:

 

 

 

 

1

 

 

 

 

F

*

 

 

 

 

 

 

 

 

 

 

 

 

(z) = z

(z 1)(z 2)

2

 

 

 

 

 

 

 

 

 

и разобьем дробь в скобках на сумму простейших дробей:

 

1

 

=

A

 

+

B

 

+

C

 

=

(z 1)(z

2)2

z

1

z

2

(z 2)2

 

 

 

 

 

 

=

A(z 2)

2 + B(z 1)(z 2) + C(z 1)

.

 

(z 1)(z 2)2

 

 

 

 

 

 

 

 

 

 

 

 

После чего приравняем числители начальной и конечной дроби: 1 = A(z 2)2 + B(z 1)(z 2) +C(z 1).

Подставляя в полученное равенство различные значения переменной

z , найдем неизвестные коэффициенты A, B, C :

 

 

 

 

 

 

 

 

 

 

 

z = 2 :

1 = C;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z =1:

 

1 = A;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z = 0 : 1 = 4A + 2B C, B = −1.

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*

1

 

 

 

1

 

1

 

 

z

 

 

z

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

(z) = z

 

 

 

 

+

 

 

 

=

 

 

 

 

+

 

 

 

.

 

 

 

z 2

(z 2)

2

z 1

z 2

(z

2)

2

 

 

z 1

 

 

 

 

 

 

 

 

 

Оригинал находим из таблицы:

fn =1 2n + n 2n1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ:

fn =12n +n 2n1.

Пример 5. С помощью дискретного преобразования Лапласа решить линейное разностное уравнение

xn+2 6xn+1 +9xn = 0; x0 =1, x1 = −1.

Решение.

Способ 1. Пусть xn X * (z), тогда

xn+1 z(X * (z)x0 )= z(X * (z) 1)= zX * (z)z ,

35

xn+2 z2 (X * (z) x0 x1z1 )= z2 X * (z)z2 + z .

Применяя к обеим частям уравнения дискретное преобразование Лапласа, получим операторное уравнение

z2 X * (z)z2 + z 6zX * (z) + 6z +9X * (z) = 0 ,

откуда

 

 

 

 

 

 

 

(z 3)2 X * (z) z2 + 7z = 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выразим функцию X * (z):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X * (z) =

z(z 3) 4z

 

=

z

 

4

 

 

 

z

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z 3

(z

3)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z 3)2

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как

 

z

 

 

 

3n ,

 

 

z

 

 

 

 

n3n ;

то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z 3

 

 

(z 3)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn = 3n 4n3n1 = 3n1 (3 4n) .

 

 

 

 

 

 

 

 

Способ 2. Используем для определения функции xn

формулу:

 

k

 

 

*

 

n1

 

 

 

z

2

 

7z n1

 

d

z

n+1

7z

n

2

 

 

xn = res(X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z)z

 

 

 

 

 

 

 

2 z

 

 

 

 

 

 

 

2

 

(z 3)

 

=

 

 

)= res

(z

3)

 

= lim

dz

 

(z 3)

 

 

 

i=1

zi

 

 

 

 

 

 

z=3

 

 

 

 

 

z3

 

 

 

 

 

 

 

lim ((n +1)zn 7nzn1 )= n3n + 3n 7n3n1 = 3n1 (3 4n) .

z3

Ответ: xn = 3n1 (3 4n) .

Пример 6. С помощью дискретного преобразования Лапласа решить линейное разностное уравнение

x

n+2

4x

n+1

+ 4x

n

= 5n ,

x = 0, x = 1.

 

 

 

 

0

1

Решение. Способ 1. Пусть xn X * (z), тогда

x

n+1

 

z(X * (z)x

),

x

n+2

 

z2 (X * (z)x x z1 ), 5n

 

 

z

.

 

 

 

 

 

 

 

 

 

 

0

 

0

1

 

 

 

 

z 5

Имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

z2 X * (z)z 4zX * (z)+ 4X * (z) =

 

 

 

 

 

 

z

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

36

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

*

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

*

(z) =

 

 

 

 

z2 4z

 

 

 

 

 

 

 

 

 

 

 

 

 

(z 2) X

 

(z)z =

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

z 5

 

 

 

 

 

 

(z 5)(z 2)2

 

 

 

Разложим функцию X * (z) на слагаемые следующим образом

 

 

 

 

 

 

z2 4z

 

 

= z

 

 

 

z 4

 

 

 

 

=

 

1 z

 

 

 

 

 

1 z

+

2

 

 

 

 

 

 

z

 

 

.

 

 

(z 5)(z 2)2

(z 5)(z 2)2

 

9

 

z 5

 

9

 

 

z 2

 

 

3

 

(z 2)2

Так как

 

 

z

 

 

 

 

 

2n ,

 

 

 

z

 

 

 

 

 

 

 

 

 

n 2n1

,

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

5n , то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

2

 

 

(z

− 2)2

 

 

 

z 5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn =

 

 

1

(5n 2n )+

2

n 2n1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Способ 2. Используем для определения функции xn

формулу

 

 

 

k

 

 

 

 

 

*

 

 

n1

 

 

 

 

 

 

z

n+1

− 4z

n

 

 

 

 

 

 

 

 

 

 

 

z

n+1

− 4z

n

 

 

xn = res(X

 

(z)z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

)= res

(z − 5)(z

2)

 

 

+ res

(z − 5)(z

2)

=

 

 

 

i=1

zi

 

 

 

 

 

 

 

 

 

 

 

z=2

 

 

 

 

 

 

z=5

 

 

 

 

 

 

 

d

 

 

zn+1 − 4zn

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zn+1 − 4zn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=lim

 

 

 

 

 

 

 

 

 

 

 

 

 

(z − 2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 (z

5)

dz

 

(z

− 5)(z − 2)

 

 

 

 

+ lim

(z

5)(z − 2)

 

 

 

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=lim

((n +1)zn 4nzn1 )(z 5)

(zn+1 4zn )

 

 

5n+1 4

5n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z 5)

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

 

 

 

 

 

 

 

 

 

 

 

 

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

((n + 1)2n − 4n2n1 )(−3) − (2n+1 − 4 2n )+

5n

 

=

 

1

 

(5n − 2n )

+

 

2

 

n 2n1.

 

 

 

9

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ: xn =

 

1

(5n 2n )

+

 

 

2

n 2n1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

Пример 7. С помощью дискретного преобразования Лапласа решить линейное разностное уравнение

x

n+2

5x

n+1

+ 6x

n

= 3n ,

x = 1, x = 0.

 

 

 

 

0

1

Решение.

Способ 1. Пусть xn X * (z), тогда

37

 

 

 

 

 

 

 

 

 

x

n+1

 

 

 

z(X *

(z)x

0

)= zX * (z)z ,

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

z2 (X *

(z)x x z1 )= z2 X * (z)

z2 ,

3n

 

 

 

 

 

.

 

 

 

 

n+2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z 3

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z2 X * (z)z2 5zX

* (z)+ 5z + 6 X * (z) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z2 5z + 6)X * (z) =

 

 

 

z

+ z2 5z

 

 

 

X * (z) =

z3 8z2 +16z

.

 

z

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z 2)(z 3)2

Разложим функцию X * (z) на слагаемые следующим образом

 

 

 

 

 

 

 

z3 8z2 +16z

= z

 

 

 

z2 8z +16

 

 

 

=

3z

+

 

 

 

z

+

 

4z

 

 

.

 

 

 

 

(z 2)(z 3)2

 

(z 2)(z 3)2

 

 

z 3

(z 3)2

 

z 2

 

 

 

Из таблиц находим оригинал:

xn

= −3n +1 + n 3n 1 + 4 2n.

 

 

 

 

 

 

 

 

 

 

Способ 2. Используем для определения функции xn

 

формулу

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z 4)

2

z

n

 

 

 

 

 

 

 

(z

4)

2

z

n

 

 

xn = res(X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*

(z)z

 

n1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

)= res

 

 

 

 

 

 

 

3)

+ res

(z

2)(z

3)

=

 

i =1 zi

 

 

 

 

 

 

 

 

 

 

 

z =3 (z 2)(z

 

 

 

z =2

 

 

 

 

 

d

 

 

 

(z 4)2 zn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z 4)2 zn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= lim

 

(z 2)(z 3)

2 (z 3)

+ lim

 

 

2)(z 3)

2 (z 2)

=

 

z3 dz

 

 

 

 

 

 

 

 

z2

(z

 

 

 

 

 

 

 

 

 

 

 

 

=lim (2(z 4)zn + (z 4)2 nzn1)(z 2) (z 4)2 zn + 4 2n =

3 2)2(z 1z

=

2(1)3n + n3n1

3n

+

4

2n

= −3n+1

+ n 3n1 + 4 2n.

1

 

 

1

 

 

 

 

 

 

 

 

 

 

Ответ: xn

= −3n +1 + n 3n 1 + 4 2n.

ЗАДАНИЯ.

38

1. Пользуясь определением, найти изображение F* (z) для сле-

дующих функций:

a) fn = en . б) fn = eαn .

в) fn = n2 .

2. Найти оригиналы для следующих изображений:

а) F (z) =

z

.

б) F (z) =

z

.

(z 3)2

 

 

 

z2 +1

3. С помощью дискретного преобразования Лапласа решить линейное разностное уравнение

1.

x

n+2

8x

n+1

+16x

n

= 2n , x

0

=1, x =1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2.

x

n+2

2x

n+1

+ x

n

= 4n , x

0

=1, x = 0.

3.

 

 

 

 

 

= 3n

 

 

 

 

 

1

x

n+2

6x

n+1

+ 5x

n

, x

0

=1, x =1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

4.

xn+2 + 4xn+1 + 3xn = n, x0 = 1, x1 = 1.

5.

x

n+2

5x

n+1

+ 6x

n

= 8n

, x

0

= 0, x =1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

6.

x

n+2

2x

n+1

= 5n , x

0

=1, x = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

7.

x

n+2

6x

n+1

+ 5x

n

= (1)n , x

0

= 0, x = −1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

8.

x

n+2

4x

n

= 5n , x

0

= −1, x = 2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

9.

xn+2 +8xn+1 +16xn = 3, x0 = 1, x1 = 0.

10.

x

n+2

2x

n+1

+ x

n

= 6n ,

x

0

= 0,

x = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

11.

x

n+2

+ x

n

= (1)n , x

0

=1,

 

 

x

=1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

12.

xn+2 4xn+1 + 3xn

= n,

x0

= 0,

x1 = 0.

13.

xn+2 6xn+1 + 8xn

= 2,

x0

= 0,

x1 = 1.

8. Вариационное исчисление

Определение. Если каждой функции y(x) из некоторого множе-

ства поставлено в соответствие некоторое число J, то говорят, что на этом множестве задан функционал J ( y) = J[ y] .

Приведем примеры функционалов:

39

1. Длина плоской кривой, заданной уравнением y = y(x) , a x b :

b

J[ y] = 1 + y2 (x)dx .

a

2.Стоимость проезда по дорогам, имеющим вид кривых, соединяющих пункты А и В.

Простейшая задача вариационного исчисления ставится так:

среди всех функций

y(x) , найти ту,

которая доставляет экстремум

функционалу J[ y] . Эта кривая удовлетворяет уравнению Эйлера

 

 

(a, A)

(b, B)

Fy

d

 

Fy= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

y(a) = A,

y(b) = B.

 

 

 

 

 

F

 

 

F

 

 

a

b

(Здесь Fy =

;

Fy=

).

 

 

 

 

 

 

 

 

y

 

y

Данная краевая задача может иметь единственное решение, может иметь множество решений или не иметь ни одного.

Решение типовых задач

Пример 1. Найти экстремали функционала

0

 

 

 

 

 

 

 

 

 

J ( y) = ( y2 + 2 yy′+ y2 )dx ,

y(1) =1, y(2) = 0 .

1

 

 

 

 

 

 

 

 

 

Решение.

 

 

 

 

 

2

+

2 yy

+ y

2

.

F (x, y, y ) = y

 

 

 

Найдем частные производные функции F (x, y, y ) :

 

 

 

 

 

 

 

 

 

Fy = 2 y+ 2 y , Fy= 2 y+ 2 y .

Уравнение Эйлера примет вид

 

 

 

 

 

 

 

 

2 y′+ 2 y

d

(2 y′+ 2 y) = 0 ,

dx

2 y′+ 2 y 2 y′′− 2 y′ = 0

или

y′′− y = 0 .

40