51
комплексная переменная z , как оператор, заменяется оператором сдвига E , в предположении нулевых начальных значений
Vz |
|
1 |
|
vk |
|
1 |
. |
Ez |
|
z d |
10 |
|
E d |
||
Преобразуя полученное выражение, приходим к разностному уравнению вида
vk 1 d vk 10 fk
или
vk 1 d vk 10 d vk fk .
Определение начальных условий. Для однозначного определения решения разностного уравнения необходимы независимые условия, в качестве которых удобно воспользоваться начальными условиями.
Начальные условия могут быть определены по изображению выходной переменной, в соответствии с теоремой теории Z - преобразования о начальном значении функции
v |
lim v |
lim V |
|
lim |
1 |
0 . |
|
|
|||||
|
|
|||||
0 |
k 0 k |
z |
z |
z |
z d |
|
Так как импульсная характеристика для дискретных и цифровых систем определена при k 1, найдем значение v1, используя теорему Z - преобразования об упреждении на один такт
|
|
v1 |
z Vz |
z |
v0 . |
|
|
|
|
|
В соответствии |
с этим |
выражением и теоремой теории Z - |
||||||||
преобразования о начальном значении функции, получаем |
|
|||||||||
|
v1 |
lim z (Vz |
v0 ) |
lim |
|
z |
|
1 . |
|
|
|
|
|
|
|
||||||
|
z |
|
|
|
||||||
|
|
z |
|
z |
|
|
d |
|
||
С другой стороны для определения начальных условий можно |
||||||||||
воспользоваться |
исходным |
разностным |
|
|
уравнением, |
полагая |
||||
соответствующими значения индекса k . |
|
|
|
|
|
|
|
|||
Так при k |
1 получаем |
|
|
|
|
|
|
|
|
|
|
|
v0 |
d v 1 |
1 1. |
|
|
|
|
|
|
Учитывая, что входное воздействие и выходная реакция при отрицательных
значениях k |
равны нулю, окончательно получаем v0 0 . |
При k |
0 получаем |
|
v1 d v0 10 . |
Учитывая, что входное воздействие существует только в нулевой момент времени k 0, начальное значение v0 0 , окончательно получаем v1 1.
Метод вариации произвольных постоянных (метод Лагранжа). В
соответствии с методом Лагранжа, общее решение неоднородного разностного уравнения первого порядка записывается в виде
y |
c |
y |
c |
d k |
c |
d k , |
k |
1,k |
1,k |
1,k |
1 |
1,k |
|
52
где c1,k - варьируемая постоянная, то есть неизвестная пока дискретная функция; y1,k d1k - фундаментальное решение однородного уравнения; d1 d - корень характеристического уравнения.
Определяющее уравнение Лагранжа, позволяющее определить варьируемую постоянную, запишется в виде
c1,k |
fk |
|
1k |
. |
y1,k 1 |
|
d k 1 |
||
|
|
|
Применяя обратный разностный оператор, приходим к выражению
|
1 |
k |
fn 1 |
|
k |
1 |
fn |
|
k |
1n 1 |
|
|
c1,k |
c1,k |
c1 |
|
|
c1 |
|
c1 |
, |
||||
|
y1,n |
|
|
y1,n 1 |
n 1 d n |
|||||||
|
|
n 1 |
|
n |
0 |
|
|
|
||||
где c1- новая постоянная суммирования.
Учитывая то факт, что входное воздействие существует только в нулевой момент времени, получаем, что значение суммы соответствует первому слагаемому
c1,k |
1 |
c1. |
|
|
|
||
|
d |
||
|
|
|
|
Подставляя полученное выражение в общее решение, находим |
|||
v d k |
1 c d k . |
||
k |
|
|
1 |
Для определения частного |
решения неоднородного разностного |
||
уравнения необходимо из независимых условий, например, начальных условий, определить постоянную суммирования c1. Так как импульсная характеристика для дискретных и цифровых систем определена при k 1,
найдем постоянную суммирования c1 |
по значению v1 1 |
||||||||
|
v |
1 |
|
d 0 |
c |
d 1 |
c d , |
||
|
|
1 |
|
|
|
1 |
|
|
1 |
то есть c1 0 . |
|
|
|
|
|
|
|
|
|
|
Подставляя значение c1 |
в общее решение, находим частное решение |
|||||||
разностного уравнения |
|
|
|
|
|
|
|
|
|
|
g |
k |
v d k 1 |
0 d k |
d k 1, |
||||
|
|
k |
|
|
|
|
|
||
при k |
1 или |
|
|
|
|
|
|
|
|
|
|
|
|
g |
k 1 |
v |
1 |
d k , |
|
|
|
|
|
|
k |
|
|
||
при k |
0. |
|
|
|
|
|
|
|
|
Полученные решения разностного уравнения описывают импульсную характеристику дискретной системы и, как видим, совпадают с ранее полученными операторным методом выражениями.
Решение в форме Коши (метод Коши). Согласно методу Коши,
частное решение неоднородного разностного уравнения первого порядка vk 1 d vk 1k d vk fk
представляется в виде
|
|
|
53 |
|
|
|
|
|
|
|
d k |
|
k |
n f |
|
|
|
|
|
v |
v |
d k |
n |
1 |
, |
|
|
||
k |
|
0 |
|
|
|
|
|
||
|
|
|
n 1 |
|
|
|
|
|
|
где d - корень характеристического уравнения. |
|
|
|
|
|
||||
Обратим внимание, что решение в форме Коши определено через |
|||||||||
начальное значение функции при k |
0. |
|
|
|
|
|
|
||
Подставляя начальное значение и, раскрывая правую часть разностного |
|||||||||
уравнения, получаем решение в виде |
|
|
|
|
|
|
|
||
v d k 0 |
k |
d k n 1 |
d k |
k |
1n 1 |
. |
|||
|
|
||||||||
k |
|
|
n 1 |
|
|
n 1 d n |
|||
|
n 1 |
|
|
|
|||||
Учитывая то факт, что входное воздействие существует только в нулевой момент времени, получаем, что значение суммы соответствует первому слагаемому
|
g |
k |
|
v |
d |
k |
1 |
d |
k 1 |
, |
|
|
|
|
|
|
|
||||||
|
|
|
k |
|
|
|
d |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
при k |
1 или |
|
|
|
|
|
|
|
|
|
|
|
|
|
g |
k 1 |
v |
1 |
d k |
, |
|
||
|
|
|
|
|
k |
|
|
|
|||
при k |
0. |
|
|
|
|
|
|
|
|
|
|
Полученные решения разностного уравнения описывают импульсную характеристику дискретной системы и, как видим, совпадают с выражениями, ранее полученными операторным методом и методом Лагранжа.
5.7 Методы решения разностных уравнений высоких порядков
Общие предпосылки. Решить разностное уравнение, означает, из заданного уравнения связи неизвестной функции и ее сдвигов и дополнительных независимых условий найти функцию, обращающую уравнение в тождество.
В нашем случае разностные уравнения получаются либо по заданной функциональной схеме дискретного или цифрового устройства, либо путем перехода от заданной системной функции в предположении нулевых начальных условий.
Переход от системной функции к разностному уравнению осуществляется путем замены изображений входного воздействия и реакции системы оригиналами и формальной заменой в дробно-рациональном выражении комплексной переменной z оператором сдвига E .
Обратный переход от известного разностного уравнения к системной функции осуществляется путем замены оригиналов входного изображения и реакции системы на их изображения и заменой сдвигов функций оригиналов их изображениями с учетом начальных условий в соответствии теоремы о сдвиге функции оригинала на n - тактов. После приведения подобных, в числителе дробно-рационального представления системной функции
54
окажутся входное воздействие и начальные условия, а в знаменателе – характеристический полином, определяющий свободные колебания дискретной системы.
В качестве независимых дополнительных условий, необходимых для нахождения частного решения разностного уравнения, используются начальные значения функции и ее сдвигов. Начальные условия, в свою очередь определяются либо по системной функции на основании теоремы о начальном значении дискретной функции, либо непосредственно по разностному уравнению путем последовательной подстановки соответствующих значений индексов сдвига и учете того факта, что входное воздействие и реакция системы в отрицательные моменты времени равны нулю.
Системы разностных уравнений первого порядка. В связи с рассмотрением вопроса решения разностных уравнений n - го порядка и возможностью перехода к эквивалентной системе n разностных уравнений первого порядка уместно рассмотреть системы разностных уравнений первого порядка.
Нормированная система разностных уравнений первого порядка с постоянными коэффициентами может быть записана в виде
|
n |
|
y1,k 1 |
a1,i yi,k |
f1,k ; |
i |
1 |
|
|
n |
|
y2,k 1 |
a2,i yi,k |
f2,k ; |
i |
1 |
|
|
n |
|
yn,k 1 |
an,i yi,k |
fn,k . |
i |
1 |
|
Матричная форма записи системы разностных уравнений первого порядка имеет вид
X k 1 A X k Fk ,
где X k , X k 1, Fk - векторы неизвестных функций, сдвигов неизвестных
функций и правых частей системы.
Системы разностных уравнений, являются в общем случае математической моделью многоканальных дискретных и цифровых систем. Уравнения, входящие в систему могут быть любого порядка, однако системы уравнений первого порядка являются среди них наиболее простыми и проработанными в аналитическом плане. Зачастую системы n уравнений m - го порядка сводятся к системе (n m) уравнений первого порядка.
Как уже отмечалось, скалярное разностное уравнение n - го порядка yk n an 1 yk n 1 a1 yk 1 a0 yk fk
путем замены переменных может быть преобразовано в эквивалентную систему n уравнений первого порядка.
Так, вводя новые переменные
55 |
|
x1,k yk ; x2,k x1,k 1 yk 1; |
xn,k xn 1,k 1 yk n , |
получаем систему в виде
x1,k 1 x2,k ; x2,k 1 x3,k ;
xn 1,k 1 xn,k ;
xn,k 1 |
a0 x1,k a1 x2,k |
an xn,k fk . |
В матричной записи эквивалентная система имеет вид
x1,k 1 |
0 |
1 |
0 |
0 |
x1,k |
0 |
x2,k 1 |
0 |
0 |
1 |
0 |
x2,k |
0 |
xn 1,k 1 |
0 |
0 |
0 |
1 |
xn 1,k |
0 |
xn,k 1 |
a0 |
a1 |
a2 |
an |
xn,k |
fk |
или
X k 1 A X k Fk .
Таким образом, скалярное разностное уравнение n - го порядка эквивалентно системе n разностных уравнений первого порядка с матрицей коэффициентов определенного вида. Методы решения систем разностных уравнений первого порядка можно интерпретировать, как векторно-матричную интерпретацию методов решения скалярных разностных уравнений n - го порядка.
Забегая вперед, отметим, что решение систем разностных уравнений с использованием векторно-матричной символики, предполагает знакомство с функциями матричного аргумента, в частности с определением степенной матричной функцией. Понятие функции матричного аргумента тесно связано с проблемой собственных значений и векторов.
Проблема собственных значений и векторов. Как и в случае аналоговых систем, описываемых непрерывными функциями и обыкновенными дифференциальными уравнениями высокого порядка, аналитическое решение которых записывается через экспоненту матрицы коэффициентов эквивалентной системы дифференциальных уравнений первого порядка, решение разностных уравнений дискретных систем второго и более высокого порядка, описываемых дискретными (решетчатыми) функциями, записывается через степенную функцию матричного аргумента. Выражение функции матричного аргумента предполагает предварительное знакомство с проблемой собственных значений и векторов.
В разделе, посвященном аналитическим методам решения обыкновенных дифференциальных уравнений высокого порядка и систем дифференциальных уравнений первого порядка, проблема собственных значений и векторов уже была рассмотрена, поэтому здесь мы остановимся на ней кратко.
| 05_Холера |
| 1 |
| 10.4. Исследование регистров |
| 1112 |
| 12 |
| 1285 |
| 1560 |
| 1568 |
| 1604248853606027 |
| 1673 |