|
|
|
56 |
|
|
|
|
|
|
|
|
Известно, что система линейных алгебраических уравнений n - го |
|||||||||
порядка в векторно-матричной форме записывается в виде |
|
|
|
|
||||||
|
|
|
A X |
Y , |
|
|
|
|
|
|
где |
A- матрица |
коэффициентов системы; X , Y - |
входной |
и выходной |
||||||
векторы. Естественно, что результирующий вектор Y |
зависит, |
как |
от |
|||||||
матрицы коэффициентов A, так и от входного вектора X . |
|
|
|
|
||||||
|
Возникает вопрос, существуют ли при данной матрице коэффициентов |
|||||||||
A |
такие векторы xi , преобразование |
которых |
сводится |
лишь |
к |
их |
||||
умножению на |
масштабный |
коэффициент i , |
то |
есть |
выполняется |
|||||
соотношение |
|
|
|
|
|
|
|
|
|
|
|
|
A |
xi yi |
i |
xi . |
|
|
|
|
|
Оказывается, что в общем случае невырожденной матрицы имеется ровно n таких векторов, где n - размерность системы векторов. При этом, векторы, удовлетворяющие указанному условию, называются собственными векторами матрицы, а масштабные коэффициенты при векторах называются собственными значениями матрицы. Собственные векторы xi принято
обозначать в литературе как hi .
Записывая данные соотношения, определяющие собственные вектора и значения в векторно-матричной форме, получим
A [hi ] A H H ,
где H - модальная матрица, составленная из собственных векторов; - в случае невырожденной матрицы коэффициентов это диагональная матрица собственных значений.
В общем случае, вырожденной матрицы, вместо диагональной матрицы имеем Жорданову матрицу J , на главной диагонали которой располагаются так называемые “ящики или клетки Жордана”, соответствующие кратным собственным значениям. В данной дисциплине мы не будем подробно останавливаться на вырожденных случаях соответствующих нулевым либо кратным собственным значениям, так как в аналитическом случае всегда можно предположить, что все собственные значения различны и отличны от нуля, довести аналитическое решение до конца, а затем использовать предельный переход к истинным значениям.
Собственные вектора образуют собственный канонический базис линейно независимых векторов. Каноническое представление матрицы в собственном базисе запишется на основании предыдущего соотношения в виде
H |
1 A H . |
Наоборот любая невырожденная матрица на основании этого |
|
соотношения может быть представлена в виде |
|
A H |
H 1. |
На основании канонического разложения невырожденной матрицы вводится понятие аналитической функции от матричного аргумента
57
F ( A) H F ( ) H 1,
где F ( ) - понимается как диагональная матрица с элементами f ( i ) .
На основании векторно-матричного соотношения, определяющего собственные вектора и значения, учитывая перестановочность диагональной матрицы [ ], можно записать основное определяющее однородное уравнение
|
A |
H 0 , |
где A |
- характеристическая матрица. |
|
Далее, учитывая, что в общем |
случае модальная матрица H не |
|
вырождена, получаем, что характеристическая матрица равна нулю, то есть вырождена
A 0 .
Приравнивая нулю, определитель характеристической матрицы, получаем характеристическое уравнение
det( A |
) |
A |
0 . |
|
|
|
|
Раскрывая определитель характеристической матрицы, и, приводя подобные, получаем характеристический степенной полином относительно переменной
n |
an 1 |
n 1 |
a1 |
a0 ( |
1) ( |
2 ) |
( |
n ) 0 . |
|
|
|||||||
Корни |
характеристического |
полинома |
i |
соответствуют |
собственным |
|||
|
|
|
|
|
|
|
|
|
значениям матрицы коэффициентов системы. |
|
|
|
|||||
Собственные вектора hi , соответствующие собственным значениям i |
||||||||
находятся из систем однородных уравнений вида |
|
|
||||||
|
|
|
|
A |
i hi B( i ) hi 0 , |
где |
i |
[ i ]- диагональная |
матрица, содержащая одинаковое i - тое |
||
собственное значение по диагонали. |
|||||
|
Так как матрица |
A |
i |
вырождена, то есть ее ранг равен, по крайней |
|
|
|
|
|
|
|
мере, |
n |
1, получаем, что однородная система, определяющая собственный |
|||
вектор hi , должна иметь нетривиальное решение. Аналитические
соотношения, определяющие собственные вектора из решения однородных систем, можно получить, записывая определяющую систему в виде
n |
|
|
|
|
bij ( |
i ) hij |
0 . |
|
|
i 1 |
|
|
|
|
Далее, раскрывая определитель матрицы |
A |
i |
B( i ) , по какой либо |
|
строке i , и, учитывая, что он равен нулю, получаем |
|
|||
n |
|
|
|
|
det(B( i )) |
bij ( i ) |
ij ( |
i ) |
0 . |
i 1 |
|
|
|
|
Теперь, сравнивая два последних выражения, получаем
|
58 |
|
|
hij k ij ( |
i ) , |
где j 1, 2, |
, n ; k - некоторая константа; |
ij ( i ) - алгебраическое дополнение |
матрицы A |
i . |
|
Напомним, что матрица, полученная из транспонированной исходной матрицы, заменой элементов на их алгебраические дополнения, называется присоединенной матрицей
C [cij ] Adj(B( i )) [ ij ] .
Из полученного соотношения следует, что собственный вектор hi , с точностью до константы, определяется алгебраическими дополнениями элементов i - той строки матрицы A i , то есть соответствующим
столбцом присоединенной матрицы C .
Разностные уравнения второго порядка. Прежде чем перейти к разностным уравнениям высокого порядка остановимся подробнее на уравнениях второго порядка.
Как уже отмечалось, в изучаемой дисциплине предпочтение отдано универсальным аналитическим методам – операторному, Лагранжа и Коши.
В общем случае нормированное неоднородное разностное уравнение
второго порядка имеет вид |
|
|
|
|
|
|
|
|
|
yk 2 |
a1 yk 1 |
a0 yk |
|
|
fk |
||
или |
|
|
|
|
|
|
|
|
|
(E2 |
a E |
a ) |
y |
f |
k |
, |
|
|
|
1 |
|
0 |
k |
|
|
|
где yk - неизвестная |
дискретная |
функция; |
|
yk n - сдвиг неизвестной |
||||
дискретной функции |
на n - |
тактов; |
ai - |
коэффициенты уравнения; fk - |
||||
произвольная известная дискретная функция, определяемая внешним воздействием.
Операторный метод. Операторный метод решения разностных уравнений основан на переходе от оригиналов правой и левой частей уравнений к изображениям с учетом начальных условий, используя теорему о смещении дискретной функции на n - тактов
|
zn |
|
n 1 |
m |
|
y |
Y |
y z |
, |
||
k n |
|
z |
m |
|
|
|
|
|
m 0 |
|
|
или частный случай при n |
1 |
|
|
|
|
|
yk 1 |
z Yz |
z y0 . |
|
|
В результате преобразований приходим дробно-рациональному выражению изображения выходной реакции системы
Yz |
Fz |
Nz |
, |
|
z2 a |
z a |
|||
|
|
|||
|
1 |
0 |
|
|
где Fz - изображение входного |
воздействия; N z - изображение, |
|||
обусловленное начальными условиями. |
|
|
||
59
Далее, применяя обратное Z - преобразование к дробно-рациональному представлению выходной переменной, находим оригинал выходной реакции системы, то есть решение разностного уравнения.
В том случае, если задана системная функция дискретной системы дробно-рациональным отношением общего вида
|
Y |
b z2 |
b |
z |
b |
||
Sz |
z |
|
2 |
1 |
|
0 |
, |
X z |
z2 |
|
|
|
|||
|
a z |
a |
|||||
|
|
|
|
1 |
|
0 |
|
где X z - изображение входного |
воздействия, |
то, выражая изображение |
|||||
выходной реакции |
|
|
|
|
|
|
|
|
|
|
|
|
(b z2 |
b |
z |
b ) |
X |
z |
|
Y |
S |
z |
X |
z |
2 |
1 |
|
0 |
|
, |
|
|
|
|
|
|
|
||||||
z |
|
|
z2 |
a |
z |
a |
|
|
|
||
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
1 |
|
0 |
|
|
|
и, применяя обратное Z - преобразование, находим оригинал выходной реакции системы, то есть решение разностного уравнения.
Операторный метод легко обобщается на решение разностных уравнений произвольного порядка и основная трудность заключается в выполнении обратного Z - преобразования для сложных дробнорациональных функций. Для выполнения обратного Z - преобразования используются, либо готовые таблицы, либо теория вычетов функций комплексного переменного, позволяющая заменить интеграл по замкнутому контуру на комплексной плоскости суммой предельных значений функции в особых точках, называемых вычетами.
Метод Лагранжа для разностных уравнений второго порядка.
Рассмотрим аналитическое решение нормированного неоднородного разностного уравнения с постоянными коэффициентами второго порядка
E2 a E a |
y |
k |
y |
a y |
k 1 |
a y |
k |
f |
k |
|
1 |
0 |
|
|
k 2 1 |
0 |
|
||||
методом вариации произвольных постоянных (методом Лагранжа).
Согласно методу Лагранжа, общее решение неоднородного разностного уравнения записывается в том же виде, что и однородного уравнения, только вместо постоянных используются варьируемые постоянные c1,k , c2,k , то есть неизвестные функции
|
|
y |
c |
|
y |
c |
y |
c |
d k |
c |
d k , |
|
|
k |
1,k |
1,k |
2,k |
2,k |
1,k |
1 |
2,k |
2 |
|
где y |
d k , |
y |
d k |
- фундаментальные решения однородного уравнения; |
|||||||
1,k |
1 |
2,k |
2 |
|
|
|
|
|
|
|
|
d1, d2 - корни характеристического уравнения.
С целью определения варьируемых постоянных подставим общее решение в исходное разностное уравнение. Для этого последовательно подействуем на предполагаемое решение оператором сдвига. После первого действия получим
yk 1 c1,k 1 y1,k 1 c2,k 1 y2,k 1.
Добавляя и вычитая дополнительные слагаемые, преобразуем выражение к виду
|
60 |
|
yk 1 c1,k 1 y1,k 1 |
c1,k y1,k 1 c2,k 1 y2,k 1 |
c2,k y2,k 1 |
c1,k y1,k 1 c2,k y2,k 1 |
|
|
c1,k y1,k 1 |
c2,k y2,k 1 c1,k y1,k 1 |
c2,k y2,k 1 . |
Повторное действие оператора сдвига, может привести к появлению разности второго порядка от варьируемых постоянных, что исключает возможность их определения. Для ограничения роста порядка разности варьируемых постоянных наложим ограничение, полагая
|
c1,k y1,k 1 |
c2,k y2,k 1 0 . |
|
||||
В результате наложения ограничения имеем |
|
||||||
|
yk 1 |
c1,k y1,k 1 |
c2,k y2,k 1. |
|
|||
Теперь повторное действие оператора сдвига дает выражение |
|||||||
yk 2 |
c1,k 1 y1,k 2 |
c2,k 1 y2,k 2 |
|
|
|
||
|
c1,k y1,k 2 |
|
c2,k y2,k 2 |
c1,k y1,k 2 |
c2,k y2,k 2 . |
||
Подстановка полученных выражений в исходное разностное уравнение |
|||||||
и перегруппировка дают |
|
|
|
|
|
|
|
|
c1,k y1,k 2 |
|
c2,k y2,k 2 |
|
|||
|
c1,k [ y1,k 2 |
a1 y1,k 1 |
a0 y1,k ] |
|
|||
|
c2,k [ y2,k 2 |
a1 y2,k 1 |
a0 y2,k ] |
fk . |
|||
Заметим, что выражения в скобках равны нулю, так как удовлетворяют |
|||||||
однородному |
разностному |
уравнению |
в |
силу того, |
что y1,k , y2,k - его |
||
фундаментальные решения.
В результате подстановка общего предполагаемого решения в исходное разностное уравнения приводит к уравнению
c1,k y1,k 2 c2,k y2,k 2 fk .
Таким образом, соотношение ограничивающее рост порядка разности варьируемых постоянных и результат подстановки предполагаемого общего решения в исходное неоднородное разностное уравнение дают, так называемую определяющую систему уравнений Лагранжа, позволяющую определить разности варьируемых постоянных
c1,k y1,k 1 |
c2,k y2,k 1 |
0 ; |
c1,k y1,k 2 |
c2,k y2,k 2 |
fk . |
Определяющая система уравнений Лагранжа линейна относительно разностей варьируемых постоянных
c1,k |
y1,k 1 |
y2,k 1 |
0 |
. |
|
c2,k |
y1,k 2 |
y2,k 2 |
fk |
||
|
Определитель системы, построенный из набора фундаментальных решений однородного уравнения и их сдвигов, называется определителем Касорати. Независимость фундаментальной системы решений, соответствует не равному нулю определителю Касорати
| 05_Холера |
| 1 |
| 10.4. Исследование регистров |
| 1112 |
| 12 |
| 1285 |
| 1560 |
| 1568 |
| 1604248853606027 |
| 1673 |