Материал: Фарфоровская Ю. Б. Разностные уравнения и Z-преобразование

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

nqn

 

 

z / q

 

 

 

 

 

 

zq

 

 

 

,

 

 

 

 

 

 

 

 

 

z / q 1

2

 

z

q

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zq2

 

 

 

 

 

 

 

n(n 1)

q

n

 

 

 

 

z / q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

z / q

1 3

z

q 3 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n(n 1)...(n m 1)

q

n

 

 

 

 

 

 

z / q

 

 

zqm

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m!

 

 

 

 

 

 

 

 

 

 

 

 

m 1

z q

m 1

 

 

 

 

 

 

 

 

 

 

z

/ q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

Заметим, что с помощью линейных комбинаций можно из этих фор-

мул получать z-преобразования для последовательностей вида (пm) и (птqn). Напримeр,

 

 

 

n(n 1)

 

 

 

 

 

 

 

 

 

2z

 

 

 

 

 

 

 

 

z

 

 

 

 

 

z(z 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n2 2

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

z

1

3

 

z 1

2

 

z

1

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n2qn

z / q z / q 1

 

 

z z q q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z / q

 

1

 

 

 

 

z q 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим

теперь

 

z-преобразования

последовательностей

 

(qncos(nβ)

 

(qnsin(nβ)). Для этого достаточно сначала найти z-преобразования последо-

 

вательностей (cos(nβ)) и (sin(nβ)).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим

 

 

сначала последовательность {cos(nβ)}. Применим

 

к членам этой

 

 

 

 

последовательности

 

 

 

формулы

 

 

Эйлера,

 

а

 

именно

 

 

ein

e in

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iпβ

 

 

 

-iпβ

 

 

 

 

 

 

 

 

cos(n )

 

 

 

 

 

 

 

 

. Так как последовательности e

 

 

 

и e

являются гео-

 

 

2

 

 

 

 

 

 

 

 

 

метрическими прогрессиями со знаменателями еiβ и e-iβ соответственно, то

 

члены нашей последовательности равны сумме двух геометрических про-

 

грессий. Тогда мы получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

in

e

n

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

z

 

 

 

 

z

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos(n )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

2

 

 

 

 

2

1

 

(ei / z)

1 (e i / z)

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z ei

 

 

z e i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее, приводя в последней формуле слагаемые к общему знаменате-

 

лю и снова применяя формулы Эйлера, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos(n

1

z2

ze i z2 zei

 

1

(

 

 

2z z ei e i / 2

 

 

 

z z cos

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

i

 

 

i

 

 

 

 

 

 

 

 

 

2

 

 

2z e

i

 

i

/ 2 1

 

 

2

 

 

 

 

 

 

2

 

 

z

ze

ze

1

 

 

 

 

2

 

z

 

e

 

z

2z cos 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее, применяя свойство 5, после несложных преобразований, полу-

 

чим

 

 

 

 

 

 

 

 

 

qn cos(n )

 

 

 

z(z q cos )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z2

2qz cos q2

 

 

 

 

 

 

 

 

 

 

 

 

 

16

Аналогичным образом выводятся z-преобразования для последова-

тельностей {sin(nβ)} и {qnsin(nβ)}. Эти z-преобразования приведены в таблице, содержащей все основные применяющиеся на практике z- преобразования.

Напомним, что в математическом анализе рациональной функцией (или рациональной дробью) называется функция, равная отношению двух полиномов. Рациональная дробь называется правильной, если степень ее числителя строго меньше степени ее знаменателя.

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

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

Если правильная рациональная дробь F(z) не равнa нулю или бесконечности при z = 0, то ее можно умножить и разделить на z и искать ориги-

нал от F1(z) = F(z) · z. Из свойств преобразований следует, что нужный нам

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

корень z = 0 с кратностью k то домножить и разделить ее надо на zk + 1 и ис-

кать оригинал для F1(z) = F(z)/zk + 1. Оригинал же данной функции будет таков: первые k членов будут равны нулю, а далее будет располагаться

оригинал для F1(z) с заменой п на п + k. Подробнее все это мы рассмотрим на конкретных примерах в следующем разделе.

Таблица соответствий последовательностей и их z-преобразований

 

{an}

 

 

 

 

 

a

f(z) = nn

 

 

 

 

 

 

k 0 z

1.

{1}

 

 

 

 

z

 

 

 

 

 

 

 

z 1

 

 

 

 

 

 

 

 

 

2.

{qn}

 

 

 

 

z

 

 

 

 

 

z

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.

{n}

 

 

 

 

z

 

 

 

 

(z

1)2

 

 

 

17

4.

 

 

{nqn}

 

 

 

 

 

 

 

 

 

zq

 

 

 

 

 

 

 

 

 

 

 

(z q)2

 

 

 

 

 

5.

 

n(n 1)

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

(z 1)3

 

 

 

 

 

 

 

 

 

 

 

 

 

6.

n(n 1)

q

n

 

 

 

 

 

 

 

 

zq2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

(z q)3

 

 

 

 

 

 

 

 

 

 

 

 

7.

n(n 1)(n 2)...(n m 1)

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m!

 

 

 

 

 

 

 

(z 1)m 1

 

 

 

 

 

 

 

 

 

 

 

8.

n(n 1)(n 2)...(n m 1)

q

n

 

 

 

 

 

 

zqm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m!

 

 

 

 

 

(z q)m 1

 

 

 

 

 

 

 

 

 

 

9.

 

 

{n2}

 

 

 

 

 

 

 

 

 

 

z(z 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

(z 1)3

 

 

 

10.

 

 

{n2qn}

 

 

 

 

 

 

 

qz(z q)

 

 

 

 

 

 

 

 

 

 

 

(z q)3

11.

 

{cos(nβ)}

 

 

 

 

z(z cos )

 

 

 

 

 

 

 

z2 2z cos 1

 

 

 

 

 

 

 

 

 

 

12.

 

{sin(nβ)}

 

 

 

 

 

 

 

 

z sin

 

 

 

 

 

 

 

z2 2z cos 1

 

 

 

 

 

 

 

 

 

 

13.

{qncos(nβ)}

 

 

 

 

z(z q cos )

 

 

 

 

z2 2qz cos q2

 

14.

{qnsin(nβ)}

 

 

 

 

 

 

 

 

zq sin

 

 

 

 

z2 2qz cos q2

 

6. РЕШЕНИЕ ТИПОВЫХ ПРИМЕРОВ

Авторы считают, что усвоение некоторой математической теории непосредственно связано с умением решать (по крайней мере, в процессе этого усвоения) типовые задачи, характерные для этой теории. Мы предлагаем (в случае разностных уравнений) 6 задач трех типов. Первые 2 задачи – это решение однородных разностных уравнений методом Эйлера (задачи 1–20 а) и б)). Задачи 21–40 а) и б) – это задачи, связанные с нахождением оригинала по данному z-преобразованию. В задачах 41–60 а), б) требуется решить неоднородные разностные уравнения с помощью z-преобразования. Во всех заданиях возможные последовательности составлены из вещественных чисел.

6.1. Решение задач 120 методом Эйлера. В этих задачах могут уча-

ствовать разностные уравнения 1, 2 и 3-го порядков. Сначала мы составля-

18

ем характеристическое уравнение, (которое является алгебраическим) 1, 2 или 3-го порядка соответственно. Во всех случаях решение характеристического уравнения не должно вызывать затруднений. Среди корней этого уравнения могут быть кратные или комплексные корни. Рассмотрим несколько вариантов возможных заданий.

Решение задач 1–20, т. е. решение однородных линейных разностных уравнений методом Эйлера. Здесь могут встретиться разностные уравнения 1, 2 и 3-го порядков. Во всех случаях решение характеристического уравнения не должно вызывать трудностей. Корни характеристического уравнения могут быть простыми или кратными, вещественными или комплексными.

Пример 1. Корни характеристического уравнения вещественны и различны (уравнение 2-го порядка).

ап + 2 + 2ап + 1 – 8ап = 0, а0 = 1, а1 = 4.

Запишем характеристическое уравнение

q2 + 2q – 8 = 0.

Решаем его

q1,2 2 24 32 22 6 1 3 , т. е. q1 = 2, q2 = –4.

Значит, общее решение данного разностного уравнения имеет вид: an

= C1q1n + C2q2n, т. е. an = C12n + C2(–4)n = C12n + (–1)nC24n. Используем начальные условия для нахождения постоянных:

п = 0: а0 = С1 + С2 = 1;

п = 1: а1 = 2С1 – 4С2 = 4.

Решаем систему

 

 

 

 

C C

 

1

 

2

 

 

1

2

 

 

 

2C1 4C2 4

 

 

Получаем C1 43 , C2 13 . Итак, an 43 2n 13 ( 1)n 4n .

Рассмотрим теперь самый простой случай. Пример 2. Уравнение 1-й степени.

ап + 1 + 2ап = 0, а0 = 1.

Характеристическое уравнение имеет вид:

q + 2 = 0; q = –2, an = C(–2)n, n = 0: a0 = C = 1, C = 1, an = (–2)n = (–1)n2n.

Пример 3. Характеристическое уравнение (2-го порядка) имеет кратные, т. е. равные, вещественные корни.

19

 

 

ап + 2

– 6ап + 1

+ 9ап = 0, а0 = 1, а1 = 2.

 

 

 

Характеристическое уравнение:

 

 

 

 

 

 

 

 

 

 

 

 

q2 – 6q + 9 = 0, q

1,2

= 3. a

n

= C 3n

+ C n3n.

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

Подставляем начальные условия.

 

 

 

 

 

 

 

 

 

п = 0: а

0

= С = 1, С = 1; п = 1: а

1

= С 31 + С 31 = (С

1

+ С )3 = 2,

 

1

1

 

 

 

 

 

 

 

1

2

 

2

 

 

 

С + С = 2 ,

 

 

C

2

1 .

 

 

 

 

 

 

 

1

 

2

 

3

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общее решение имеет вид:

a

3n

1 n 3n 3n n 3n 1

3n 1(3 n).

 

 

 

 

 

n

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Наконец, самый «неприятный» случай уравнения 2-го порядка, когда характеристические уравнение имеет комплексно сопряженные корни. Заметим, что теоретически этот случай отличается от предыдущих только тем, что ответ получится в комплексной форме (несмотря на то что сама последовательность будет вещественной при всех п. Но мы покажем как получается (через тригонометрическую форму комплексных корней) вещественный ответ.

Пример 4. Характеристическое уравнение имеет комплексно сопряженные корни.

an + 2 + 2an + 1 + 10an = 0, a0 = 0, a1 = 1. Тогда q1,2 = –1± 3i, q1 = –1 + 3i, q2 = –1 – 3i.

В соответствии с теорией найдем показательную форму числа q1 (для этого надо найти модуль и аргумент этого числа): q1 q2 10 ;

q1 = –1 + 3i лежит во второй четверти; tg φ = –3, φ = π – arctg3. Тогда

q1 10 ei ,

q2 10 e i ; an = C1q1n + C2q2n = 10 n C1ein C2e in .

Найдем С1 и С2. Подставим п = 0. Тогда а0 = 0 = С1 + С2, то есть С2 = С1.

inφ

-inφ

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

Отсюда ап = С1(е

-e

 

n

 

ein e in

 

n

2i, получим a 2i 102 C

, то есть a 2 102 i C sin n .

 

n

1

2i

n

1

 

 

 

 

 

Подставим теперь п = 1. Получим уравнение 2

10 iC1 sin 1. Так как

sin φ положителен (угол φ лежит во второй четверти ), то

sin

 

 

tg 3

 

 

3

 

3

.

 

 

 

 

1 tg

2

1 9

10

 

 

 

 

Итак 2 10

iC

3

1;

iC

 

1

, и получаем окончательный ответ:

10

6

 

1

 

1

 

 

20

Источник: https://studfile.net/preview/15049180/