Применение многочленов от матриц для решения систем рекуррентных уравнений
Рассмотрим систему линейных рекуррентных уравнений
 (7.47)
где — коэффициенты системы, — заданные, а — неизвестные функции дискретного аргумента , . При описании дискретных динамических систем аргумент в (7.47) называют дискретным временем.
Систему (7.47) можно записать в матричном виде:
 (7.48)
где — квадратная матрица (n-го порядка) коэффициентов системы рекуррентных уравнений, — столбец заданных функций, a — столбец неизвестных.
Решением системы (7.48) называется последовательность столбцов , при подстановке которых в (7.48) получаются верные равенства для всех .
Поставим задачу нахождения решения системы (7.48), удовлетворяющего начальным условиям
 (7.49)
где — заданный столбец.
Рассмотрим сначала однородную систему с постоянными коэффициентами
 (7.50)
Записывая (7.50) для , последовательно получаем
 и т.д.
Следовательно, решение однородной системы (7.50), удовлетворяющее начальным условиям (7.49), имеет вид:
 (7.51)
Получим теперь решение системы (7.48). Учитывая (7.49), запишем (7.48) для 
и т.д. Следовательно, решение системы (7.48) имеет вид
 (7.52)
Первое слагаемое в (7.52) — решение однородной системы (7.50) с начальными условиями (7.49), второе слагаемое — решение системы (7.48) с нулевыми начальными условиями (т.е. при ).
Алгоритм решения системы рекуррентных уравнений
Для нахождения решения системы (7.48) с начальными условиями (7.49) требуется выполнить следующие действия.
1. Найти выражение для степени матрицы одним из способов, рассмотренных ранее.
2. Записать по формуле (7.52) искомое решение.
Замечания 7.9
1. Линейное рекуррентное уравнение с постоянными коэффициентами:
может быть сведено к эквивалентной системе рекуррентных уравнений вида (7.47). Действительно, используя обозначения получаем
или в матричной форме (7.48): с матрицами и 
Пример 7.19. Найти решение рекуррентного уравнения с начальными условиями 
Решение. В соответствии с пункту 1 замечаний 7.11 составим систему уравнений
где . Требуется найти решение этой системы, удовлетворяющее начальным условиям: .
1. Составим матрицу системы и найдем ее степень , используя второй способ. Характеристический многочлен имеет вторую степень и два простых корня: . Поэтому многочлен (7.44) — линейный: . Для его коэффициентов записываем систему (7.46):
Отсюда . Таким образом,
2. По формуле (7.52) записываем решение системы (учитывая, что ):
Нас интересует только первый элемент этого столбца:
Решение совпадает с найденным в примере 2.15.
Пример 7.20. Найти решение системы рекуррентных уравнений с начальными условиями 
![\begin{cases}x_1(k+1)=-2x_1(k)+5x_2(k)-3x_3(k)+k,\\[2pt] x_2(k+1)=-2x_1(k)+ 5x_2(k)-3x_3(k)+3^k,\\[2pt] x_3(k+1)=x_1(k)-x_2(k)-k. \end{cases}](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZ0AAABPBAMAAADCcAUBAAAAMFBMVEVHcEwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAlTPQ5AAAAD3RSTlMAQcFbqBCBMSGCcdDwkeCHEAbsAAAP6klEQVR42uxaf2wT1x3/nu/Ovtjxr0EpolUTOTBoCY1FItaNtD4gQEUAHxCmjTDiUejQGORgAbrCFmewqavaOlv2TzapeC20E0NLVgidBkEh/Fw7hsdAldahsLCObUxqnN9QEu/9ONv3nm3Mnwj1ifj87p7f+3y+P9/7HgAPZdvnf6joiB9q3J0v5Rg54b7nrNHMl9zNFiIXIZpvwo+Yyz1boIMnqIJ9fYbOVgNIORWpcMB/B/KWGWQmpsnkn6nt1mDSMnRdwE8osF17jEJaF+bksaE4Q0R3uTGwRwN7mYe7ueoYwnI+l4p7x5nxVhWUrV3oyzIO5ZktMzneMHEQXddxy630XWQk9IQBSWjnwNdn8HHc5u88hf6Cyc4KenH3DqPPn+Xg85hWGzf3X0MiKsQrRQxZiQafRGIhIwc0RhxAXwo4WMcgwKi2JAnphnHjJePaGMowt7O8qofQx61kb1LyLuZTp2fnU62L4+b+FbxSM/qQjOUsxiTzJjO/cyETco5iQ4yzdpnQI308pF787Tnjzk+NaxtvRxDkbBycSFbKWFY+jubsfMpCkDD3P0iuZL/O8uGkWYSMykGGdLIPqqCoy9SVMCRiRxGd5RPMwJKgsUWs2nGFjnWEsBHIb+mZfJzF2fns1GTE58DGCorXFicr7VsIwijH51H0+XqJvYT238Q6V0UUDT7mp6wPIUg/ukLt1YUhxeX1OhQakfCbxrBumMjJKEE1VvFqdwN1Njf6iTRUM1fN5GMlRuAljTU9cRAs/urn/0s72Dy6xQWnAAZYPtUb3gfb9MaSAMU5H9u7/289SJtcIJyEDOSdNWcaKAg3skupD0OyqAwf27D1F0OcqRI+SrFwp1ElcTeCbrjeX7i/meWD1S1gPkoraaxYXBegBXp+cgcgTGnbbj9rGQG4zfLpgOWaRXO3Xwor2JOmor+e118ap5Znjv9b/gA2VRgoUmEbetKIIU3HkCzFDB97fIG7i/V+ykfQrSOCZl2KOkXoL3JUw4Fthdd7wOt9NM0nzqUdoimMbKkOfjguexScEDAf+6Bf0Skfwet9DA2jiNuiqyDSvAvmbI0BPIlunFgMk+misI1MR8c1qkLYOipozvIbZkjE4pu83iV0WWnQY2PTjc2wN5CQXwp1Bp/ACRJXv9faOr+1tSUnH5FoCvmmiHOyPIYNBS1YgPhIn/0Lj8C/sre2vomGoXVlDZDEIYCI/gPPhXPaeDcGQPh8mUxH8TmQpVoQpIh/vhmSs52Om0qWhcI/H8mIBwYfN7ZMN+pgTsGv9fLxjfAZwbCp/zCBcq7mNDKJ1eBTqAaiwNvb2i4oQs+XkhA4TOxNHJBGk3zSQi4ByzCFZA1PpS4QVLrxRKy9NYZ69Kx8hBlFUTFK+GDn64VSMZTJBwMFmfqPecuF/KIQrqFM8jzhI/bhldyhxQB3GD5116EtVOPvhhpkq4i9D6cfaxw9DTCYpH6kH3FhwC/6wX7TDInj0+Yp8rcreiYfy0jQU6sRPijJoPTjk8KZ8a3gevZ4vcfne08eX9cnxAgfHK/bPHUb22nkTvNx6FCq1f9qEOYhzEj4m7FNWYfm0Mht2q4NQOQTy2jQsxLg2WkoFjTj9IMhuUIMnyAUeVXXQCYfa2fVsenU3pDLCQPw22uMvTnLEkiYFjVHPk0kBqGzpPOyRvjAaWxTlqkeoikTH1vngZvw+PpKvLWsoPnUrSoXp9FIZ2qXv31Sd3ZWfRFBwqkJcUfutvUadieGzz9h7fywYzCL/+yCJsN/hC6QPWDzsP5D7eVeJ6UVmqBT/0HpHV4BHLYkleEDq1o0UDzyI3jrjaJVoR9sGjRpqTyVcqB3PPj3TThV1yOxIEg6gVQJDJ8dgOeqYfiYMWI+qT0faXvNnefynz4In9eS2aTOkIo1Y1yD921kv0l9i9dzzVev12smSE9mO6J9JVt8I9P2LEKfs3OCvZGXjrL5nCe9DYWf5xx4KYFcS75g2pjmOBmWnzZBEs5mGzMtFx/DV3KdFgvU+zyeytMNszmbZ+D2vLxhAo5LkgGpNhtv0X9vPjj6ZG0T7/u8/UNqcLV6vrN+jIo9r6AW08tf8y+deLiqIUoW/XzO5wFqn+vngW7C53wetGZ75OHis2/o4eIjtufhsyvHD3fd9xJ2uj9Yk3djRJO5La9Mjf3nwWzPnCrDJ9Mei8E2KdsKK/X75bME7bs2oQ1fc55xaMamt9P7vdxebkDarWV56PBvMUlyPOP52jBIn/K7IHnJPbaXyvKjGleXtpX2GXVnM6wlE9bzvOfi09ieMHfqKGPfOLxhQLKr/MJnwlD3/aK0oKXhDHRVuCqWxKcYB7LSYeNJtvaGp/KCuY9fF7jwQehxjZlFSCRuMvxwTQmXGCQO5zJcHzW1RQYk+ZRxY3VyodrzUHl1lenIOZQhbFxvOcnVr4UmzCcSzs6nShP7zf3f4LExYgpsPXFWC6NHCzoGWEeSa5qM4W440sVDOo6/JdVr1K/bVOsYXDaLuSeUccpB4UIeBL5ehfm4cpyMykIy44b/IXVcYy6mfs3+rhHZiaUr+QtT+5NW1Gf2+HZa3UO/CDN8AmrBGBxxp18cCkngSrlWTmVniZFS2staJh8pxzllto7fL7zs2UujiELq8Ro6ZwlDWfiIv7Ztov3DpO6nvEsr2UwLqAgSlBsen4Lk9jN80MxxedS9L11pSp5Vt1ccPWAUv5txVcz5y1iWelVXzhzQD84FnSX0QE7eKtyyLTkFMlePv+K7CvJThzccpsLCLIqiFX8H2M8HtA81mFgxKwkJ2YVzCEPi6vHI4vxKTFycCs1/TFpshzAWaU/VMByffNSmZvIh5SfFRxqrKksc5kD9FBpccNlRGdteE8+ox5+H90DS3dc/pQHJhwse333mLmLFpYyZN/ErguFGWj7E9VFHx2IEiasnwuwzjEPaPkidwe3xr+pi9XRaeq3z6Tu1HPXerA15y1V42rZRnqJSPuJYTIzyfOQw7Pd/A+pih6DhaT+t1pS2wAt8vRdPGEWQRr6ui9cYSFIxZ28r2Xe6e1JilpAtNUC1UfzuJ3b6BZ/vpM83K82Hq8fbiaZQkrGROIp8MfV+QUqQKgTxOp+vFA2jKiH1eDT578VRUo+HxDksOMLnMpmOjnPFKaS1niozJCu2obd8vltkWawINn+m4gEU4oU8AZ0ot2fvRT6+ZeOTaj/QRcBV1XQ93tXRg9dn9YOsBdfjg4QoojoLp5/dZ4HXj/wtUo/HkFZ6ZtL3Pz2rxiCjHj9F4/drQRqClYONIcEvQCBMQshxKFP8mXzsOfxHvgkuKC8YgU2penyk2R16l+cT6YNAzK73ghN2yLQeXzBijTuRJBk3kPrBckeJNkYFPyinaZkTQbL5eT6XokKC3QJF6B5Giu/3WMJEYq5mnH58Tk8mH2uO+Fa76MWP5buOPkWl9fgRbFN1oXYaudN8CnUI6m2bBuEQmqsDGRhOPwWjEwG+w6bAMai7YBnp8UhhwO89C1OQHDGGz/884mf02wljyyNRCUpHZs5bQP4jBE4yyHKqKxh7EzePv5g7/wQTiX44t+H03HCqHh/UXIui9E1qmo945NUBqJy9fJJKt5YBbFPCwAw+/8jHJgf9UsnMvyygm6w0pGQF3eDzxDOV/zY0ZaRqO+UnH1o9B6lnjqLhWrKyEdZsZPjYvF4vqZ9nbdvw27+d/h8fNOrXywBlQuUF43Vamg80zNJB3CS2oAmnKSQ1/L996wuNowjj3+2f7Oau11xRRASTeqatGBqXWKKtqbcmFqGK2ZSrCTZwZwMHonKH1QqmmLXWJ0Hz6IPYK/hSLF6JkQpClbS2QkGOPuRJ7B/bgj5I2qRpfKhxvp3Zy/6b3YB7D2Lm6W53dnZ+MzvffN/v+02yAtsrNMPvKC/3fwry2N7jJgxBuYJpf7VqdcnOszI86sXd7MIHxYB4QdTS5LbLAh50/ukzIiMFC097w8iwAZD8vLTxALlcD+WlsZQzZdMhtLFdtWO+Db0aEM8dWl4OJZM3RzPX914h0yzY8/gUt+L3yMerNgyBGyiN7512dMnjtzo3dDMAz/tWLjuVWU1gy8GzYQNpWGGSLa6ACUUMyGLstHFzJ159Q3N0SSry6h23PQR/vC0/wnnmcPTnxgpj7oUobpwp0NRNUcPEKnzBrWGG8KMjnGc+XD2FRJEPR1akVlaNHKgRZ/X/E9+7hmcNzxqef4NHW5uf/xgeXn57WygrTh0oUQ/ZFbeEt98kPMSpkLqqATtauh7SDgG7H1U4A/wqxF2QujR/Gl/Mbmzm+iHecWLQg1K+7xyE6kPQtXwXqdtR/vfbgw2blOZ1PZvTY4MTkD/FXjd0wEy/I3QuQpgSxnL9k0j/8EUkFvncDQHjkovPJvnzPyIG1TOBfG+BP/YojEY+JsSlt8TYFv/ylefOJbOJ6wcpLuWvQDwCfwGdYBQRwB/c5aMzDqjBQtvldPz2oCX7xPP01/pJ1AGrz2p+PGlvANTzaJ6leJAAmKjmd/uSKunsO/2VBmvbuqSQf23u70tegP31mPEMiBf2UCtDxdtvtQfwvVQ/WrJK1TLOhVOX6H0U7NfefqYL4DN3+w+LZz+iDSOLJlze2T7pnWdxPv/gXLx4EnpifkKH3gp97bqH6rVqAJ+IdA3qYV8tjeFje6BQv0APeyDcF44ML3ipdbGYnO+oyztMyhm2DejlFRbNnsI7W9pOxotHNVqWFJDuv0FfOzENCpIhpVJ/qTTG5UdHiGUyYdtNRoT+/A0GXJZUt9eaQwMbluZUyN9zhjbccQNzdRaelxrTnLplqLHbAxRvF6q/0NfWXjyF16L08SwvyfAoC6dxXdQYg2VL3THx+KXBGj76K/oIFp6hhoy70HkufvtW0MFIGNOUXB2U/2QRuft7Q1us0vVDmYnvcLgRz++4/Ui3ZX8q55OibPQYZYqnW74Ghvd769BqFSMuPDSTLZ7PaUkNDm6iZmgGusWA9eM8X4LDcOjJu0h9IJ6tuP203j7io26v5zJIHffTLM4CXCWv8diDspnLnIw+E7i6Qs8zCT8Nam8CHCCfmVAEZdGnJxfvBtjro9sXkQxCPBfxhFnL3Cte6nb9j49lSMPKeesUD9l+plKGNWSO0g25A5Pr5uPBk/vB6u7041N43uF1A8+yJe7A01+79lMpmyWeSsrjy4xvfQ/TNYinj4x/Vb1KKrjXQjK749ouKsWQdEguQR9puMPtkJ6FfVOQ/jYePIJjXHrIxIPiyuC79eTB/g7iGbW7GOjvyLtEp8j68ybGC8rfK0vxhEnmxy3qGHL+4ZxvvOl0Q6UgfzQ1+5rTDf2tiXicJyPzs/imdiPMQfaX2U7dIQ0fDfIta8u3HO65uLGZeMSZlf4Ph4aYUpgJGrf9G36VVvb8x2Yz8cDhoufCMU7F3lXE2y1hoQwLtJ9rEpB/AGTZWy3ScxvFAAAAAElFTkSuQmCC) (7.53)
Решение. Запишем систему в матричной форме (7.48) и начальные условия 
 где  .
1. Выражение для степени матрицы найдено в примерах 7.17 и 7.18:
 при 
2. Выражение, полученное для , справедливо только при (см. пример 7.18). Поэтому для формулу (7.52) нельзя использовать. Найдем непосредственно из системы (7.53)
Далее записываем искомое решение для , выделяя в формуле (7.52) слагаемые с и 
В последней сумме слагаемые не зависят от индекса суммирования, т.е. это сумма одинаковых слагаемых. Поэтому, приводя подобные члены, получаем
Таким образом, решением заданной системы является последовательность столбцов:
 где  .
Если заметили ошибку, опечатку или есть предложения, напишите в комментариях.
|