
Анализ временных характеристик дискретных и цифровых устройств
..pdf
|
|
|
|
|
|
|
|
|
|
|
|
|
|
61 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
c2,k |
1 |
|
|
c2,k |
|
k 1n b 1n 1 |
|
|
|
|
1 |
|
k 1n |
|
b 1n 1 |
|
|
|
b |
c2 , |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
n 1 d n (d 1) |
|
|
|
d 1 n 1 |
|
d n |
|
|
|
d (d 1) |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
где c1, c2 - новые постоянные суммирования. |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
Подставляя найденные значения c1,k |
и c2,k |
в предполагаемое общее |
||||||||||||||||||||||||||||
решение разностного уравнения, получаем его в виде |
|
|
|
|
|
|
|
|||||||||||||||||||||||
vk |
|
|
b 1k |
|
c1 |
1k |
b d k |
|
|
|
|
c2 |
d k |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
d |
1 |
|
d (d |
|
1) |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
b b d k 1 |
c |
c |
|
d k |
|
b (d k 1 |
1) |
|
c |
|
c |
|
d k . |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
d |
1 |
|
d |
1 |
|
1 |
2 |
|
|
|
|
|
d |
1 |
|
|
1 |
2 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Для определения постоянных суммирования c1 |
и |
c2 |
воспользуемся, |
|||||||||||||||||||||||||||
найденными ранее, |
начальными условиями v1 |
1 |
и |
v2 |
1 |
b d , так как |
решение разностного уравнения в виде импульсной характеристики определено при k 1. Так, приравнивая общее решение, при k 1 и k 2 , начальным условиям, находим
|
|
|
|
|
v1 |
1 |
0 |
c1 |
|
|
c2 |
d ; |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
v |
1 |
b |
d |
b |
|
|
c |
c |
d 2 |
|
|
|
|
|
|
|
||||
|
|
|
2 |
|
|
|
|
|
|
|
1 |
2 |
|
|
|
|
|
|
|
|
|
|
||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c1 |
d |
c2 |
1; |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
c |
|
|
d 2 |
c |
|
1 |
d . |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
1 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Определим значения постоянных суммирования c1 |
и |
c2 , |
используя |
|||||||||||||||||||||
правило Крамера для решения линейной системы |
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
1 |
d |
|
|
|
|
|
|
|
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
d |
|
d (d 1); |
c1 |
|
|
1 |
d |
d |
2 |
|
|
1 |
|
; c2 |
|
|
1 |
1 |
d |
|
1 |
. |
|
|
1 |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
d 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
1 |
|
|
|
|
|
|
|
|
|
d 1 |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d 1 |
Подставляя найденные значения постоянных суммирования в общее решение, получаем частное решение исходного неоднородного разностного уравнения
|
g |
|
v |
b (d k 1 |
1) |
|
1 |
|
d k |
1 d k |
|
|
b (1 d k 1) |
, |
|||||||
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
k |
d |
1 |
|
|
d 1 |
d 1 |
1 |
d |
1 |
d |
|
||||||||
|
|
|
|
|
|
|
|||||||||||||||
при k |
1 или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
gk |
|
vk |
|
|
b (1 |
d k ) |
1 |
d k |
1 |
, |
|
|
|||||
|
|
|
|
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
1 |
d |
1 |
d |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
при k |
0. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Полученное решение |
описывает |
импульсную |
характеристику, |
исследуемой дискретной системы, и, как видим, совпадает с выражениями, найденными операторным методом.

62
Решение в форме Коши (метод Коши). Рассматриваемый нами вариант метода Коши предполагает предварительное преобразование исходного неоднородного разностного уравнения
|
vk 2 |
(1 d ) vk 1 |
d vk |
11 |
b 10 |
fk |
||
в эквивалентную систему двух разностных уравнений первого порядка. |
||||||||
Так, |
вводя новые |
переменные |
x1,k |
vk ; |
x2,k x1,k 1 vk 1; |
|||
x3,k v2,k 1 |
vk 2 , получаем эквивалентную систему разностных уравнений |
|||||||
первого порядка |
|
|
|
|
|
|
|
|
|
x1,k |
1 |
0 |
1 |
x1,k |
0 |
|
|
|
x2,k 1 |
d 1 d |
x2,k |
11 b 10 |
||||
или |
|
|
|
|
|
|
|
|
|
|
|
X k 1 |
A X k |
Fk . |
|
|
Согласно методу Коши, частное решение неоднородной системы разностных уравнений первого порядка следует искать в виде
|
|
|
|
|
Ak X |
|
k |
|
|
|
|
|
|
X |
k |
0 |
Ak n F |
1 |
, |
||
|
|
|
|
|
|
n |
|
|||
|
|
|
|
|
|
|
|
n 1 |
|
|
где X0 |
x1,0 x2,0 |
t |
v0 |
v1 |
t - вектор начальных условий; Ak - степенная |
|||||
|
функция от матрицы коэффициентов системы.
Как известно, любая аналитическая функция от матрицы, имеющей различные и отличные от нуля собственные значения, на основании ее модального представления
A H H 1,
может быть определена в виде
F ( A) H F ( ) H 1,
где - диагональная матрица собственных значений; H - модальная матрица собственных векторов; F ( ) k диагональная матрица указанной функции
от каждого собственного значения.
Собственные значения матрицы коэффициентов системы определяются из характеристического уравнения
det([ A |
]) |
|
1 |
2 |
(1 d ) |
d 0 . |
|
|
|||||
d |
1 d |
|
||||
|
|
|
|
|
||
|
|
|
|
|
|
|
Как видим, данное уравнение полностью совпадает с характеристическими уравнениями, определяемыми либо знаменателем системной функции, либо левой (однородной) частью разностного уравнения. Таким образом, собственные значения матрицы коэффициентов системы или корни характеристического уравнения представляются в виде
1 |
0 |
1 |
0 |
|
|
d . |
|
0 |
2 |
0 |
|
|
|
|

|
63 |
|
|
Собственные вектора, как столбцы модальной |
матрицы H , по |
||
известным собственным значениям матрицы |
A, определяются из решения |
||
однородных систем уравнений |
|
|
|
A |
i hi 0 , |
|
|
где i - диагональная матрица, составленная из |
i . |
|
|
Доказывается, что модальная |
матрица |
H может |
быть определена |
алгебраическими дополнениями элементов одной из строк, например первой,
матрицы [ A |
i ] |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
H |
11( 1) |
11( 2 ) |
|
1 d 1 1 d d |
d 1 |
|||||||||||||||
12 ( 1) |
12 ( 2 ) |
|
d |
|
|
|
|
|
|
d |
d d . |
|||||||||
Определитель |
|
модальной |
матрицы равен |
|
|
|
H |
d (d |
1) . Используя |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
определитель, выразим матрицу обратную модальной |
|
|
||||||||||||||||||
|
H |
1 |
|
1 |
|
|
|
d |
1 |
|
|
1 |
|
|
|
1 |
|
1/ d |
. |
|
|
|
|
|
|
|
d |
d |
|
|
|
|
|
|
1 |
1 |
|||||
|
|
|
d (d |
1) |
|
|
|
d |
1 |
|||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||
В результате, получаем выражение степенной функции от матрицы |
||||||||||||||||||||
коэффициентов |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ak |
H |
|
k |
H 1 |
d |
|
1 |
1k |
0 |
|
|
|
1 |
|
1 |
1/ d |
||||
|
|
d |
|
d |
0 |
d |
k |
|
d |
1 |
1 |
1 |
||||||||
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
1 |
|
|
d 1k 1 d k |
1k |
|
d k |
|
|
|
. |
|
|
|
|||||
|
|
d |
1 |
k |
d |
d |
k |
k |
d d |
k |
|
|
|
|||||||
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
d 1 |
|
1 |
|
|
|
|
|
|
|
Заметим, что структура матрицы Ak n аналогична и отличается лишь показателем степени.
Теперь все подготовлено для представления решения в форме Коши. Предварительно отметим, что в нашем случае вектор начальных условий X 0
и вектор правой части эквивалентной системы разностных уравнений Fk имеют вид
X0 |
x1,0 |
v0 |
0 |
; Fk |
0 |
0 |
. |
|
x2,0 |
v1 |
1 |
fk |
11 b 10 |
||||
|
|
|
Кроме того, отметим, что нас интересует лишь первая компонента вектора решений x1,k vk .
Обратим внимание на тот факт, что, несмотря на определение импульсной характеристики при k 1, в отличие от метода Лагранжа, где в качестве начальных условий используются значения v1 и v2 , в методе Коши
вкачестве начальных условий используются значения v0 и v1.
Всвязи с отмеченными обстоятельствами, используя общее выражение для решения системы разностных уравнений первого порядка, и, учитывая структуры векторов X 0 и Fk , выразим выходное напряжение дискретной
системы

|
|
|
|
|
|
|
|
|
64 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
k |
|
|
|
|
||
vk |
|
|
( 1k |
d k ) 1 |
( 1k n d k n ) (1n |
b 1n 1) |
||||||||
|
|
|
|
|||||||||||
d |
1 |
|||||||||||||
|
|
|
|
n 1 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
1 |
|
|
(d k |
1 |
k |
(1 |
b 1 |
) d k k |
1n |
b 1n 1 |
. |
||
|
|
|
|
|
|
|
|
|||||||
|
|
d |
1 |
|
|
n |
n 1 |
|
|
d n |
||||
|
|
|
n 1 |
|
|
n 1 |
Используем для раскрытия сумм формулы арифметической и геометрической прогрессий. Учитывая тот факт, что входное воздействие в данном случае существует только при k 0, получаем значения сумм равные первым слагаемым, определяемым вторыми компонентами сумм.
В результате, получаем окончательное выражение для выходной реакции, исследуемой дискретной системы, как частное решение исходного разностного уравнения
|
g |
|
v |
|
1 |
|
|
(d k |
1) b b d k 1 |
1 d k |
|
b (1 d k 1) |
, |
|||||||||
|
k |
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
k |
d |
1 |
|
|
|
|
|
|
|
|
1 |
d |
1 |
d |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
при k |
1 или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
g |
|
|
v |
|
1 |
d k 1 |
b (1 |
d k ) |
, |
|
|
|
|||||
|
|
|
|
|
k |
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
k |
1 |
d |
1 |
d |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
при k |
0. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Полученные |
|
выражения |
|
совпадают |
с результатами, |
операторного |
метода и метода Лагранжа и описывают импульсную характеристику, исследуемой дискретной системы.
Пример F. Пусть задана системная функция дискретной системы второго порядка
S (z) Sz |
V (z) |
|
Vz |
|
1 |
, |
E(z) Ez |
|
(z d1) (z d2 ) |
||||
|
|
|
||||
где Ez - изображение входного |
воздействия; Vz - изображение выходной |
реакции и требуется определить частотную, переходную и импульсную характеристики системы.
Частотная характеристика дискретной системы определяется по
системной функции путем замены z e j |
|
T |
|
|
|
|
||||||||
S( ) |
V ( ) |
|
|
|
1 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
E( ) (e |
j T |
d1) (e j T |
d2 ) |
|||||||||||
|
||||||||||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
, |
||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
e 2 j T |
(d d |
2 |
) e j T |
d d |
2 |
|||||||
|
|
|
||||||||||||
|
|
|
|
|
|
1 |
|
1 |
|
|||||
где T - период дискретизации по времени. |
|
|
|
|
|
|||||||||
Амплитудно-частотная характеристика системы соответствует |
||||||||||||||
модулю комплексной частотной характеристики |
|
|
|
|
||||||||||
|
|
|
|
S( ) |
|
Abs(S( )) . |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|

65
Фазочастотная характеристика системы соответствует аргументу комплексной частотной характеристики
( |
) Arg(S( )) 180 / . |
|
|
|
|||||||
Изображение выходной реакции запишется |
|
|
|
||||||||
Vz |
|
|
Ez |
|
|
|
. |
|
|
|
|
|
|
d1) (z |
|
|
|
|
|
|
|||
|
|
|
(z |
|
d2 ) |
|
|
|
|||
Знаменатель системной (передаточной) функции, приравненный нулю, |
|||||||||||
определяет характеристическое уравнение |
|
|
|
|
|
|
|
||||
(z d ) (z d |
2 |
) z2 |
(d d |
2 |
) z d d |
2 |
0 , |
||||
1 |
|
|
1 |
|
|
|
1 |
|
|||
корни которого, соответственно, равны d1 |
d1; d2 |
d2 . |
|
|
|||||||
Переходная характеристика |
дискретной |
системы. Приступаем к |
определению переходной характеристики дискретной системы различными методами.
Как известно, переходная характеристика представляет собой реакцию дискретной системы, находящейся в исходном состоянии покоя, на входную последовательность 1k (единичных - импульсов при k 0 и периодом T ).
Под исходным состоянием покоя следует понимать полное установление реакции на предыдущие воздействия и отсутствие сторонних источников.
Операторный метод. Операторный метод определения выходной реакции дискретной системы основан на теории Z - преобразования дискретных функций, как оригиналов, в непрерывные функции комплексного аргумента z , называемых изображениями и наоборот.
Оригиналу входного воздействия 1k , согласно теории Z -
преобразования, соответствует изображение в плоскости комплексной переменной z вида
z ek 1k Ez z 1 .
Изображение выходной реакции дискретной системы будет иметь вид
z
Vz (z 1) (z d1) (z d2 ) .
Из таблиц обратного Z - преобразования находим оригинал выходной реакции, то есть переходную характеристику, исследуемой дискретной системы
Vz |
|
|
z |
vk |
A B1 e |
1 t |
B2 e |
2 t |
|
||||||
|
|
|
|
|
|
||||||||||
|
(z 1) (z d1) (z d2) |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
A |
B |
d k |
B |
d k |
|
|
||
|
|
|
|
|
|
|
|
1 |
1 |
|
2 |
2 |
|
|
|
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
d k |
|
|
|
|
d k |
|
|
|
hk |
vk |
|
|
|
|
1 |
|
|
|
|
2 |
|
, |
||
(1 d1) (1 d2 ) (1 d1) (d2 |
d1) (1 d2 ) (d1 d2) |
||||||||||||||
|
|
|
|

|
|
|
|
|
|
|
|
|
|
|
66 |
|
|
|
|
|
|
|
|
|
где |
|
A |
|
|
1 |
|
; |
|
B1 |
|
|
|
1 |
|
|
; B2 |
|
1 |
; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
(1 d1) (d2 |
d1) |
(1 d2 ) (d1 d2 ) |
|||||||||||||
|
|
|
(1 d1) (1 d2 ) |
|
|
|
|
|
||||||||||||
|
|
ln(d1) |
|
|
ln(d2 ) |
|
|
t |
|
|
1 t |
k |
|
|
2 t |
k |
|
|
||
1 |
|
T |
; |
2 |
T |
; |
k |
T |
; |
e |
|
d1 |
; |
e |
|
d2 ; |
T - период |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
входной последовательности. |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
Отметим, что при k |
0 имеем v0 |
|
h0 |
0 . |
|
|
|
|
|
|
|||||||||
|
Построение |
разностного |
уравнения |
|
дискретной |
системы. |
Построение разностного уравнения дискретной системы осуществляется по
системной функции путем |
замены |
|
|
переменной |
z |
на k , в левой части |
||||||||||||||
выражения и заменой оператора zn |
на оператор сдвига En , в правой части |
|||||||||||||||||||
выражения |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
Vz |
|
|
1 |
|
|
|
|
|
|
1 |
|
|
|
|
|
||||
|
Ez (z d1) (z d2) |
|
z2 |
(d d |
2 |
) z d d |
2 |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
1 |
|
||
|
|
|
|
vk |
|
|
|
|
|
1 |
|
|
|
. |
|
|
|
|
||
|
|
|
1k |
|
E2 |
(d d |
2 |
) E d d |
2 |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
1 |
|
|
1 |
|
|
|
|
|
|
|||
Преобразуя выражение, получаем неоднородное разностное уравнение |
||||||||||||||||||||
второго порядка |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
vk 2 |
(d1 |
d2 ) vk 1 |
d1 d2 vk |
1k |
fk |
|
|
|||||||||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
vk 2 |
|
(d1 d2 ) vk 1 |
d1 d2 vk |
|
1k . |
|
|
Отметим, что переход от системной функции к разностному уравнению осуществляется в предположении нулевых начальных значений, а истинные начальные значения учитываются позже при решении уравнения.
Определение начальных условий. Для однозначного определения решения разностного уравнения необходимы дополнительные независимые условия, в качестве которых удобно воспользоваться начальными условиями. Так как исходное разностное уравнение второго порядка, необходимо определить v0 и v1.
Начальные условия могут быть определены по изображению выходной переменной, в соответствии с теоремой теории Z - преобразования о начальном значении функции оригинала
|
v |
lim v |
lim V |
|
lim |
|
z |
|
|
|
|
|
1 |
|
|
|
0 . |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
0 |
k 0 k |
z |
|
z |
z |
|
z 1 (z d ) (z d |
2 |
) |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|||||
В соответствии с теоремой упреждения, |
значение |
функции vk 1 |
||||||||||||||||||||
определится выражением |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
vk |
|
z Vz |
|
z |
v0 |
|
|
|
|
z2 |
|
|
. |
|
|
|
|
|||
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
(z 1) (z d1) (z d2) |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Применяя повторно теорему о начальном значении функции, получаем |
||||||||||||||||||||||
v |
lim v |
|
lim (z V |
|
|
z v ) |
|
lim |
|
|
z |
z |
|
|
|
|
|
0 . |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
1 |
k 0 k 1 |
z |
|
z |
|
0 |
|
z |
|
|
z 1 (z d ) (z d |
2 |
) |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|

67
С другой стороны, для определения начальных условий можно воспользоваться исходным разностным уравнением, полагая соответствующим значение индекса k , и, учитывая, что входное воздействие и реакция системы в отрицательные моменты времени отсутствуют. Так, при
k |
2 и k |
1, последовательно получаем |
|
||||
|
v0 |
(d1 |
d2 ) v 1 |
d1 d2 v 2 |
1 2 |
(d1 |
d2 ) 0 d1 d2 0 0 0 ; |
|
v1 |
(d1 |
d2 ) v0 |
d1 d2 v 1 |
1 1 |
(d1 |
d2 ) 0 d1 d2 0 0 0 . |
|
Таким образом, получаем, что начальные значения нулевые, то есть |
||||||
v0 |
0 и v1 |
0 . |
|
|
|
|
|
Решение разностных уравнений. Приступаем к определению переходной характеристики дискретной системы путем решения разностного уравнения.
Метод вариации произвольных постоянных (метод Лагранжа).
Согласно методу Лагранжа, общее решение неоднородного разностного уравнения второго порядка
vk 2 (d1 d2 ) vk 1 d1 d2 vk 1k fk
следует искать в виде
|
|
|
|
v |
c |
v |
c |
v |
c |
d k |
c |
d k , |
|
|
|
|
|
|
|
k |
1,k |
1,k |
2,k |
2,k |
1,k |
1 |
2,k |
2 |
|
|
|
где |
d , |
d |
2 |
- корни |
характеристического |
уравнения; |
y |
d k , |
y |
d k - |
|||||
|
1 |
|
|
|
|
|
|
|
|
|
1,k |
1 |
2,k |
2 |
фундаментальная система решений соответствующего однородного уравнения; c1,k , c2,k - варьируемые постоянные – неизвестные пока функции.
Варьируемые постоянные находятся из определяющей системы
уравнений Лагранжа |
|
|
|
|
|
|
|
|
c |
d k 1 |
c |
|
d k |
1 |
0 ; |
|
|
1,k 1 |
2,k |
2 |
|
|
|
|
||
c |
d k 2 |
c |
d k 2 |
1 |
f |
k |
. |
|
1,k |
1 |
2,k |
|
2 |
k |
|
|
Напомним, что определяющая система уравнений Лагранжа образуется при подстановке предполагаемого общего решения в исходное разностное уравнение и наложении ограничения на сдвиг функций c1,k и c2,k . Первое
уравнение системы есть как раз данное ограничение, а второе уравнение есть результат подстановки предполагаемого решения в исходное разностное уравнение с учетом наложенного ограничения. Определитель системы уравнений есть определитель Касорати, построенный на основе фундаментальной системы решений и их сдвигов.
Выразим разности |
варьируемых постоянных c1,k и c2,k из |
|||||||||
определяющей системы уравнений, используя правило Крамера |
||||||||||
|
d k |
1 |
d k |
1 |
|
d k 1 |
d k 1 |
|
|
|
|
|
|
|
|
||||||
C |
1 |
|
2 |
|
|
(d |
2 |
d ) ; |
||
|
d k |
2 |
d k |
2 |
|
1 |
2 |
|
1 |
|
|
|
|
|
|
|
|
||||
|
1 |
|
2 |
|
|
|
|
|
|
|

|
|
|
|
|
|
|
|
|
68 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d2k |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
c1,k |
|
fk |
d2k 2 |
|
|
|
|
1k d2k 1 |
|
|
|
|
|
1k |
|
|
; |
||||||
|
|
|
|
|
|
d k 1 |
d k 1 |
(d |
|
|
d ) |
d k 1 |
(d |
|
|
|
|||||||
|
|
|
|
|
|
|
2 |
2 |
d ) |
||||||||||||||
|
|
|
|
1 |
2 |
|
1 |
|
|
1 |
|
|
1 |
|
|
||||||||
|
|
d k |
1 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c2,k |
d1k 2 |
fk |
|
|
1k d1k 1 |
|
|
|
|
|
1k |
|
|
|
|
. |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
d k 1 |
d k 1 |
(d |
d ) |
|
d k 1 |
(d |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
2 |
d ) |
|||||||||||||
|
|
|
|
1 |
2 |
|
|
2 |
1 |
|
|
2 |
|
|
1 |
|
|
||||||
Для определения |
варьируемых постоянных |
|
c1,k |
и c2,k |
применим |
обратный разностный оператор в виде суммы функциональной
последовательности, |
используя |
|
|
|
для |
|
раскрытия |
|
|
|
|
сумм |
формулы |
||||||||||||||||||||||||
арифметической либо геометрической прогрессий |
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
1n 1 |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
c1,k |
|
|
|
c1,k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
n 1 d1n (d2 |
d1) |
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
1 |
|
|
|
k 1 |
|
|
|
|
|
|
|
|
(d k |
1) |
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
n 1 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
c1 ; |
|
|||
|
|
|
d2 |
d1 n 1 d1n |
|
|
|
|
(d2 |
d1) d1k (d1 |
1) |
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
1n 1 |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
c2,k |
|
|
|
|
c2,k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
n 1 d2n (d2 |
d1) |
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
1 |
|
|
|
|
k |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
d k |
1 |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
n 1 |
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
c2 , |
|
||||
|
|
|
d 2 d1 n 1 d2n |
|
|
|
(d 2 d1) d2k (d2 |
1) |
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
где c1, c2 - новые постоянные суммирования. |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
Подставляя найденные значения c1,k |
и c2,k |
в предполагаемое общее |
|||||||||||||||||||||||||||||||||||
решение разностного уравнения, получаем его в виде |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
d1k |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
d2k |
1 |
|
|
|
|
k |
; |
|||||
vk |
|
|
|
|
|
|
|
|
|
|
|
|
c1 |
|
d1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
c2 d2 |
||||||||
|
(d 2 |
d1) (d1 |
1) |
|
|
|
(d 2 d1) (d2 |
1) |
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
d1k |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d2k |
1 |
|
|
|
|
|
k |
k |
|
|||||||
vk |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c1 |
d1 |
c2 d2 . |
||||
|
(d 2 |
d1) (d1 |
1) |
|
|
(d 2 d1) (d2 |
1) |
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
Для определения постоянных суммирования c1 |
и c2 воспользуемся, |
||||||||||||||||||||||||||||||||||||
найденными ранее, начальными условиями v0 |
0 и v1 |
0 . Так, приравнивая |
|||||||||||||||||||||||||||||||||||
общее решение, при k |
0 и k |
|
|
1, начальным условиям, находим |
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
v0 |
0 |
|
|
|
|
0 |
|
0 c1 |
c2 ; |
|
|
|
|
|
|
|
||||||||||||||
|
|
v1 |
0 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
c1 d1 |
c2 |
d2 |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
d2 |
|
|
|
|
d1 |
|
d2 |
d1 |
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c1 |
|
|
c2 |
0 ; |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
d1 |
c1 |
|
d2 |
c2 |
|
0 . |
|
|
|
|
|
|
|
|
|
|
|

69
Из полученной системы сразу следует, что c1 c2 0 .
Подставляя найденные значения постоянных суммирования в общее решение, получаем частное решение исходного неоднородного разностного уравнения
|
|
|
|
|
d k |
1 |
|
|
d k |
1 |
|
|
|||
|
|
vk |
|
|
1 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
(d 2 |
d1) (d2 |
1) (d 2 d1) (d2 1) |
|
|||||||||
|
|
|
|
|
|||||||||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
d k |
|
d k |
|
|||
hk |
vk |
|
|
|
|
|
|
|
1 |
|
2 |
. |
|||
(1 d1) (1 d2 ) |
|
(d 2 |
d1) (1 d1) |
(d 1 d2) (1 d2) |
|||||||||||
|
|
|
|
Полученное решение описывает переходную характеристику, исследуемой дискретной системы, и, как видим, совпадает с выражением, найденным операторным методом.
Решение в форме Коши (метод Коши). Рассматриваемый нами вариант метода Коши предполагает предварительное преобразование исходного неоднородного разностного уравнения
|
vk 2 |
(d1 |
d2 ) vk 1 |
d1 d2 vk |
1k |
fk |
|
в эквивалентную систему двух разностных уравнений первого порядка. |
|||||||
Так, |
вводя новые |
переменные |
x1,k |
vk ; |
x2,k x1,k 1 vk 1; |
||
x3,k v2,k 1 |
vk 2 , получаем эквивалентную систему разностных уравнений |
||||||
первого порядка |
|
|
|
|
|
|
|
|
x1,k |
1 |
0 |
1 |
x1,k |
|
0 |
|
x2,k 1 |
d1 d2 d1 |
d2 |
x2,k |
|
1k |
|
или |
|
|
|
|
|
|
|
X k 1 A X k Fk .
Согласно методу Коши, частное решение неоднородной системы разностных уравнений первого порядка следует искать в виде
|
|
|
|
|
Ak X |
|
k |
|
|
|
|
|
|
X |
k |
0 |
Ak n F |
1 |
, |
||
|
|
|
|
|
|
n |
|
|||
|
|
|
|
|
|
|
|
n 1 |
|
|
где X0 |
x1,0 x2,0 |
t |
v0 |
v1 |
t - вектор начальных условий; Ak - степенная |
|||||
|
функция от матрицы коэффициентов системы.
Как известно, любая аналитическая функция от матрицы, имеющей различные и отличные от нуля собственные значения, на основании ее модального представления
A H H 1,
может быть определена в виде
F ( A) H F ( ) H 1,
где - диагональная матрица собственных значений; H - модальная матрица собственных векторов; F ( ) k диагональная матрица указанной функции от каждого собственного значения.

70
Собственные значения матрицы коэффициентов системы определяются из характеристического уравнения
det([ A |
]) |
|
|
1 |
2 |
(d1 |
d2 ) |
d1 d2 0 . |
|
|
|
||||||
d1 |
d2 |
d1 d2 |
|
|||||
|
|
|
|
|
|
Как видим, данное уравнение полностью совпадает с характеристическими уравнениями, определяемыми либо знаменателем системной функции, либо левой (однородной) частью разностного уравнения. Таким образом, собственные значения матрицы коэффициентов системы или корни характеристического уравнения представляются в виде
1 |
0 |
d1 |
0 |
. |
0 |
|
0 |
d2 |
|
2 |
|
Собственные вектора, как столбцы модальной матрицы H , по известным собственным значениям матрицы A, определяются из решения однородных систем уравнений
A |
i hi 0 , |
|
где i - диагональная матрица, составленная из |
i . |
|
Доказывается, что модальная |
матрица |
H может быть определена |
алгебраическими дополнениями элементов одной из строк, например первой,
матрицы [ A |
i ] |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
H |
11( |
1) |
|
11( 2 ) |
|
|
|
d1 |
d2 |
|
d1 d1 |
d2 |
|
d2 |
|
|
d2 |
|
|
d1 |
. |
||||
12 ( 1) |
|
12 ( 2 ) |
|
|
|
d1 d2 |
|
d1 d2 |
|
|
|
d1 d2 |
d1 d2 |
||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
Определитель модальной матрицы равен |
H |
d1 |
d2 |
|
(d2 d1) . |
Используя |
|||||||||||||||||||
определитель, выразим матрицу обратную модальной |
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
1 |
|
|
|
d |
d |
2 |
d |
|
|
1 |
|
|
|
1 |
1/ d |
|
|
|
||
|
H 1 |
|
|
|
|
|
|
1 |
|
1 |
|
|
|
|
|
|
|
|
|
|
2 . |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
d1 d2 (d2 |
d1) |
d1 d2 |
d2 |
|
(d2 |
d1) |
1 |
1/ d1 |
|
|
||||||||||||
|
В результате, получаем выражение стеменной функции от матрицы |
||||||||||||||||||||||||
коэффициентов |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
d |
2 |
|
d |
|
d k |
0 |
|
|
1 |
|
|
1 |
|
1/ d |
2 |
|
|||
Ak H |
k H |
1 |
|
|
|
1 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
d1 d2 |
d1 d2 |
0 d2k |
|
d2 |
d1 |
1 1/ d1 |
|
||||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
1 |
|
|
d2 d1k |
|
d1 d2k |
|
|
d1k |
d2k |
|
|
. |
|
|
|
|
|
|
|
|
|||
|
|
d2 |
d1 d1 d2 (d1k |
d2k ) |
|
d1 d1k |
d2 d2k |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
Заметим, что структура матрицы Ak n аналогична и отличается лишь показателем степени.
Теперь все подготовлено для представления решения в форме Коши. Предварительно отметим, что в нашем случае вектор начальных условий X 0
и вектор правой части эквивалентной системы разностных уравнений Fk имеют вид