Стійкість системи лінійних алгебраїчних рівнянь

МІНІСТЕРСТВО ОСВІТИ І НАУКИ УКРАЇНИ

СУМСЬКИЙ ДЕРЖАВНИЙ УНІВЕРСИТЕТ

КАФЕДРА ІНФОРМАТИКИ

Курсова робота

по чисельним методам

на тему:”Стійкість СЛАР”

Суми 2005

Зміст

1. Постановка задачі.

2. Теоретична частина.

а) характеристичний многочлен

в) метод Левeр’є

б) критерій Калашнікова

Текст програми

Приклад

Список літератури

Постановка задачі

Дана система лінійних алгебраїчних рівнянь. Необхідно дослідити її на стійкість. Знайти характеристичний многочлен методом Левур’є. Зробити відповідні висновки щодо її стійкості.

Теоретична частина

Характеристичний многочлен

Нехай дана квадратна матриця А=[aij]. Розглянемо лінійне перетворення

у=Ах (1)

де х,у n-вимірні вектори (стовпові матриці) деякого, взагалі кажучи, комплексного n-вимірного простору.

Ненульовий вектор називається власним вектором даної матриці (або визначуваного нею лінійного перетворення), якщо в результаті відповідного лінійного перетворення цей вектор переходить в колінеарний йому, тобто якщо перетворений вектор відрізняється від початкового тільки скалярним множником.

Інакше кажучи, вектор х¹0 називається власним вектором матриці А, якщо ця матриця переводить вектор х у вектор

Ах=lx (2)

Число l, стоїть в рівності (2), називається власним значенням, або характеристичним числом, матриці А, відповідним даному власному вектору х.

Теорема 1. В комплексному векторному просторі кожне лінійне перетворення (матриця) має щонайменше один дійсний або комплексний власний вектор.

Доведення. нехай А ¾ матриця лінійного перетворення. Власні вектори матриці А є ненульовими розв'язками матричного рівняння

Ах=lх

або

(А- lЕ)х=0 (3)

де матриця (А- lЕ) називається характеристичною матрицею. Рівняння (3) представляє собою лінійну однорідну систему, яка має ненульові розв'язки тоді і лише тоді, коли визначник системи рівний нулю, тобто повинна виконуватися умова

det(А- lЕ)=0. (4)

Визначник (4) називається характеристичним (віковим) визначником матриці А, а рівняння (4) називається характеристичним (віковим) рівнянням матриці А. В розгорненому вигляді характеристичне рівняння (4) запишеться таким чином:

а11-l а12 ... а1n

а21 а22-l ... а2n =0

an1 an2 ann-l

або

ln- d₁ ln-1+ d₂ ln-2- ...+(-1)n-1 d_n-1 l+(-1)n d_n=0 . (5)

Поліном, що стоїть в лівій частині рівняння (5), називається характеристичним поліномом матриці А. Коефіцієнти його di(i=1,2,…,n) визначаються за наступними правилами. Коефіцієнт d₁₌ .

Це число називається услід матриці А і позначається так: d_1=Sp А. Коефіцієнт d₂ є сума всіх діагональних мінорів другого порядку матриці А. Взагалі, коефіцієнт d_k є сума всіх діагональних мінорів k-го порядку матриці А. Зрештою, вільний член d_n рівний визначнику матриці А:

d_n=det А.

Характеристичне рівняння (5) є алгебраїчне рівняння n-ої степені відносно l і, отже, як доводиться в алгебрі, має щонайменше один дійсний або комплексний корінь. Нехай l₁ l_2,… l_m(m £n) — різні корені рівняння (5). Ці корені називаються власними значеннями, або характеристичними числами, матриці А, а сукупність всіх власних значень називається спектром матриці А. Візьмемо який-небудь корінь l=l_j і підставимо його в рівняння (4). Тоді будемо мати (А-l_jЕ )х=0 або, в розгорненому вигляді

(а11-l_j )х₁ +а₁₂ х₂ +…+а_1n x_n =0

а21х1+(а22-l_j )х₂ + ...+а_2n x_n =0

an1х1+an2х2+ ...+(a_nn -l_j )x_n =0. . . . . . . . . . . . (6)

Оскільки визначник системи (6) det(А-l_jЕ )=0, то ця система явно має ненульові розв'язки, які і є власними векторами матриці А, відповідними власному значенню l_j . Якщо ранг матриці А-l_jЕ рівний r(r<n), то існує k=n r лінійно незалежних власних векторів

х(1j) , х(2j) ...,х(kj)

відповідаючих кореню l_j . Теорема доведена.

Метод Левер’є

Відомо багато інших способів одержання характеристичного многочлена.

Розглянемо метод Левер’є, що дозволяє вирішити проблему власних значень, в основу якого покладено обчислювання слідів степенів матриці А. Вказаний метод потребує більшої кількості операцій, ніж метод Данилевського, але зовсім не чутливий до частинних особливостей матриці, зокрема ”провалів” проміжних визначників.

Нехай характеристичний поліном матриці А записано у вигляді (5) де

l_1, l_2, l_{3, .........} l_n – його корені, серед яких деякі можуть бути рівні. Позначимо

(7)

Суми , k=1-n степенів коренів многочлена зв’язані з коефіцієнтами рівняння ( 5) формулами Ньютона

k= 1,…..,n (8)

Якщо обчислити сліди ,……., матриць , ….., ,то з (8) можна послідовно обчислити коефіцієнти

Внимание, отключите Adblock

Вы посетили наш сайт со включенным блокировщиком рекламы!
Ссылка для скачивания станет доступной сразу после отключения Adblock!

Скачать