Loading [MathJax]/jax/output/HTML-CSS/jax.js
Skip to main content
LibreTexts - Ukrayinska

7.2: Діагоналізація

Результати
  1. Визначте, коли можна діагоналізувати матрицю.
  2. Коли це можливо, діагоналізуйте матрицю.

Подібність і діагоналізація

Починаємо цей розділ з згадки визначення подібних матриць. Нагадаємо, що якщоA,B двіn×n матриці, то вони схожі якщо і тільки в тому випадку, якщо існує оборотна матрицяP така, щоA=P1BP

В даному випадку пишемоAB. Поняття подібності є прикладом співвідношення еквівалентності.

Лемма7.2.1: Similarity is an Equivalence Relation

Подібність - відношення еквівалентності, тобто дляn×n матрицьA,B, іC,

  1. AA(рефлексивний)
  2. ЯкщоAB, тоBA (симетричний)
  3. ЯкщоAB іBC, тоAC (перехідний)
Доказ

ЗрозумілоAA, що, беручиP=I.

Тепер, якщоAB, тоді для деякихP обертається,A=P1BP і такPAP1=B Але тоді,(P1)1AP1=B який показує, щоBA.

Тепер припустимоAB іBC. Тоді існують оборотні матриціP,Q такі, щоA=P1BP, B=Q1CQ потім,A=P1(Q1CQ)P=(QP)1C(QP) показуючи,A що схоже наC.

Ще одним важливим поняттям, необхідним для цього розділу, є слід матриці. Розглянемо визначення.

Визначення7.2.1: Trace of a Matrix

ЯкщоA=[aij]n×n матриця, то слідA єtrace(A)=ni=1aii.

У словах слід матриці - це сума записів на головній діагоналі.

Лемма7.2.2: Properties of Trace

Дляn×n матрицьA іB, і будь-якихkR,

  1. trace(A+B)=trace(A)+trace(B)
  2. trace(kA)=ktrace(A)
  3. trace(AB)=trace(BA)

Наступна теорема включає посилання на характеристичний многочлен матриці. Нагадаємо, що для будь-якоїn×n матриціA характерний многочленA єcA(x)=det(xIA).

Теорема7.2.1: Properties of Similar Matrices

ЯкщоA іB єn×n матрицями іAB, то

  1. det(A)=det(B)
  2. rank(A)=rank(B)
  3. trace(A)=trace(B)
  4. cA(x)=cB(x)
  5. AіB мають однакові власні значення

Тепер перейдемо до основної концепції цього розділу. Коли матриця схожа на діагональну матрицю, матриця, як кажуть, є діагональною. Діагональну матрицю ми визначаємоD як матрицю, яка містить нуль у кожному записі, крім тих, що знаходяться на головній діагоналі. Точніше, якщоdij цеijth запис діагональної матриціD, тоdij=0 хібаi=j. Виглядають такі матриці наступним чином. D=[00]де число, яке не може бути нулем.

Далі наведено формальне визначення діагональної матриці.

Визначення7.2.2: Diagonalizable

AДозволяти бутиn×n матрицею. ТодіA, як кажуть, діагонально, якщо існує оборотна матрицяP така, щоP1AP=D деD діагональна матриця.

Зверніть увагу, що вищевказане рівняння можна переставити якA=PDP1. Припустимо, ми хотіли обчислитиA100. По діагоналіA спочатку досить потім обчислити(PDP1)100, що зводиться доPD100P1. Це останнє обчислення набагато простіше, ніжA100. Хоча цей процес докладно описаний далі, він забезпечує мотивацію для діагоналізації.

Діагоналізація матриці

Найважливішою теоремою про діагональності є наступний основний результат.

Теорема7.2.2: Eigenvectors and Diagonalizable Matrices

n×nAМатриця може бути діагональна тоді і тільки тоді, коли є оборотна матриця,P заданаP=[X1X2Xn] деXk є власні векториA.

Крім того,A якщо діагонально, відповідні власні значення єA діагональними записами діагональної матриціD.

Доказ

Припустимо,P наведено як вище, як обернена матриця, стовпці якої є власними векторамиA. ТодіP1 має форму,P1=[WT1WT2WTn] деWTkXj=δkj, є символом Кронекера, визначенимδij={1 if i=j0 if ij

ТодіP1AP=[WT1WT2WTn][AX1AX2AXn]=[WT1WT2WTn][λ1X1λ2X2λnXn]=[λ100λn]

І навпаки,A припустимо, діагонально так щоP1AP=D. НехайP=[X1X2Xn] де стовпці єXk іD=[λ100λn] ТодіAP=PD=[X1X2Xn][λ100λn] і так[AX1AX2AXn]=[λ1X1λ2X2λnXn] показуючиXk є власні векториA іλk є власними векторами.

Зверніть увагу, що оскільки матриця,P визначена вище, є оборотною, випливає, що набір власних векторівA{X1,X2,,Xn},, утворюють основуRn.

Ми демонструємо поняття, наведене у вищезгаданій теоремі в наступному прикладі. Зверніть увагу, що не тільки стовпці матриціP утворені власними векторами, алеP повинні бути оборотними, тому повинні складатися з найрізноманітніших власних векторів. Ми досягаємо цього за допомогою базових власних векторів для стовпцівP.

Приклад7.2.1: Diagonalize a Matrix

ДозвольтеA=[200141244] знайти оборотну матрицюP і діагональну матрицюD таку, щоP1AP=D.

Рішення

За теоремою7.2.2 ми використовуємо власні векториA як стовпчикиP, а відповідні власні значенняA як діагональні записиD.

Спочатку знайдемо власні значенняA. Для цього вирішуємоdet(λIA)=0 наступним чином. det(λ[100010001][200141244])=0

Це обчислення залишається як вправа, і ви повинні переконатися, що власні значенняλ1=2,λ2=2, іλ3=6.

Далі нам потрібно знайти власні вектори. Спочатку знайдемо власні вектори дляλ1,λ2=2. Розв'язуючи,(2IA)X=0 щоб знайти власні вектори, ми виявимо, що власніt,s векториt[210]+s[101] - це скаляри. Звідси є два основних власних вектора, які задаютьсяX1=[210],X2=[101]

Ви можете переконатися, що основним власним вектором дляλ3=6 єX3=[012]

Потім будуємо матрицюP наступним чином. P=[X1X2X3]=[210101012]Тобто стовпці єP основними власнимивекторамиA. Потім ви можете перевірити, щоP1=[14121412112141214] Таким чином,P1AP=[14121412112141214][200141244][210101012]=[200020006]

Ви можете бачити, що результатом тут є діагональна матриця, де записи на головній діагоналі є власними значеннямиA. Ми очікували цього на основі теореми7.2.2. Зверніть увагу, що власні значення на головній діагоналі повинні бути в тому ж порядку, що і відповідні власні вектори вP.

Розглянемо наступну важливу теорему.

Теорема7.2.3: Linearly Independent Eigenvectors

AДозволятиn×n матриця, і припустимо, щоA має різні власні значенняλ1,λ2,,λm. Для кожногоi, нехайXi будеλi -власним векторомA. Потім{X1,X2,,Xm} лінійно незалежний.

Наслідок, який випливає з цієї теореми, дає корисний інструмент дляA визначення діагоналізації.

Слідство7.2.1: Distinct Eigenvalues

AДозволятиn×n матриця і припустимо, що вона маєn різні власні значення. Потім випливає,A що діагонально.

Можливо, що матрицяA не може бути діагоналізована. Іншими словами, ми не можемо знайти оборотну матрицюP такP1AP=D.

Розглянемо наступний приклад.

Приклад7.2.2: A Matrix which cannot be Diagonalized

НехайA=[1101] Якщо можливо, знайдіть оборотну матрицюP і діагональну матрицюD так, щобP1AP=D.

Рішення

За допомогою звичайної процедури ми знаходимо, що власні значенняA єλ1=1,λ2=1. Щоб знайти власні вектори, вирішуємо рівняння(λIA)X=0. (λIA)Матриця задається[λ110λ1]

Підставляючи вλ=1, ми маємо матрицю[111011]=[0100]

Потім рішення рівняння(λIA)X=0 передбачає перенесення наступної доповненої матриці до її зменшеної рядково-ешелонної форми. [010000][010000]

Тоді власні вектори мають вигляд,t[10] а основний власний вектор -X1=[10]

У цьому випадку матрицяA має одне власне значення кратності два, але тільки один основний власний вектор. Для того щоб діагоналізуватиA, нам потрібно побудувати оборотну2×2 матрицюP. Однак, оскільки маєA лише один базовий власні вектори, ми не можемо побудувати цеP. Зверніть увагу, що якби ми повинні були використовуватиX1 як обидва стовпціP, неP буде обертається. З цієї причини ми не можемо повторити власні вектори вP.

Отже, ця матриця не може бути діагоналізована.

Ідея про те, що матриця не може бути діагональною, свідчить про те, що існують умови, щоб визначити, коли можна діагоналізувати матрицю. Ми бачили раніше в Слідство,7.2.1 щоn×n матриця зn різними власними значеннями діагонально. Виявляється, є й інші корисні тести на діагональність.

Для початку нам знадобиться наступне визначення.

Визначення7.2.3: Eigenspace

AДозволятиn×n матриці іλR. У власному просторіA відповідногоλ, записаногоEλ(A) є множиною всіх власних векторів, що відповідаютьλ.

Іншими словами, власне простірEλ(A) - це всеX таке, щоAX=λX. Зверніть увагу, що цей набір може бути записанийEλ(A)=null(λIA), показуючи, щоEλ(A) є підпросторомRn.

Нагадаємо, що кратність власноївеличиниλ - це кількість разів, яке воно зустрічається як корінь характеристичного многочлена.

Розглянемо тепер наступну лему.

Лемма7.2.3: Dimension of the Eigenspace

ЯкщоAn×n матриця,λ тоdim(Eλ(A))m де - власне значенняA кратностіm.

Цей результат говорить нам, що якщоλ є власним значеннямA, то кількість лінійно незалежнихλ -власних векторів ніколи не перевищує кратністьλ. Тепер ми використовуємо цей факт, щоб забезпечити корисну умову діагональності.

Теорема7.2.4: Diagonalizability Condition

AДозволяти бутиn×n матрицеюA. ПотімA діагонально, якщо і тільки якщо для кожного власне значенняλA,dim(Eλ(A)) дорівнює кратностіλ.

Складні власні значення

У деяких додатках матриця може мати власні значення, які є комплексними числами. Наприклад, це часто відбувається в диференціальних рівняннях. До цих питань підходять так само, як і вище.

Розглянемо наступний приклад.

Приклад7.2.3: A Real Matrix with Complex Eigenvalues

ДозвольтеA=[100021012] знайти власні значення та власні векториA.

Рішення

Спочатку ми знайдемо власні значення, як зазвичай, вирішивши наступне рівняння.

det(λ[100010001][100021012])=0Це зводиться(λ1)(λ24λ+5)=0. до Рішення єλ1=1,λ2=2+i іλ3=2i.

Немає нічого нового в пошуку власних векторів,λ1=1 тому це залишається як вправа.

Розглянемо тепер власне значення.λ2=2+i. Як завжди, ми вирішуємо рівняння,(λIA)X=0 як задано.((2+i)[100010001][100021012])X=[000] Іншими словами, нам потрібно вирішити систему, представлену доповненою матрицею.[1+i0000i1001i0]

Тепер ми використовуємо наші операції рядків для вирішення системи. Розділіть перший ряд на,(1+i) а потім візьмітьi раз другий ряд і додайте до третього ряду. Це дає[10000i100000] Тепер множимо другий рядок на,i щоб отримати зменшену форму рядка-ешелону, задану[100001i00000] Отже, власні вектори мають вигляд,t[0i1] а основний власний вектор задаєтьсяX2=[0i1]

В якості вправи переконайтеся, що власніλ3=2i вектори для мають виглядt[0i1] Отже, основний власний вектор задаєтьсяX3=[0i1]

Як завжди, обов'язково перевірте свої відповіді! Щоб перевірити, ми перевіряємо цеAX3=(2i)X3 наступним чином. [100021012][0i1]=[012i2i]=(2i)[0i1]

Тому ми знаємо, що це власне вектор і власне значення є правильними.

Зверніть увагу7.2.3, що в прикладі два власних значення були заданіλ2=2+i іλ3=2i. Ви можете згадати, що ці два комплексних числа є сполученими. Виявляється, що всякий раз, коли матриця, що містить реальні записиλ, має складне власне значення, вона також має власне значення, яке дорівнює¯λ, сполучений зλ.