Пример. Разложение на множители Холецкого действительных матриц

Функции > Векторы и матрицы > Разложение матриц на множители > Пример. Разложение на множители Холецкого действительных матриц

Используйте функцию Cholesky, чтобы выполнить разложение Cholesky для вещественной матрицы.

Во избежание логических несоответствий при выполнении логических сравнений включите Приблизительное равенство (Approximate Equality) в выпадающем списке Параметры расчета (Calculation Options).

1. Определите вещественную положительно определенную квадратную матрицу M.

2. Примените функцию eigenvals, чтобы обеспечить положительную определенность матрицы.

3. Примените функцию rank, чтобы обеспечить, что M является матрицей полного ранга.

4. Задайте аргументы p и u, чтобы управлять включением и отключением выбора главного элемента и нижним или верхним разложением на множители.

5. Используйте функцию Cholesky, чтобы выполнить разложение на множители по умолчанию матрицыM - с выбором главного элемента и нижним разложением на множители.

Щелкните для копирования этого выражения

Функция по умолчанию Cholesky(M) эквивалентна Cholesky(M,1,0).

<region id="ID0EOZIFB" actualWidth="284.50666666666672" actualHeight="27.791200000000003" top="844.80000000000018" left="76.800000000000011" xmlns="http://schemas.mathsoft.com/worksheet50">
 <math resultRef="80">
 <eval xmlns="http://schemas.mathsoft.com/math50">
 <apply>
 <equal />
 <apply>
 <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">Cholesky</id>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">M</id>
 </apply>
 <apply>
 <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">Cholesky</id>
 <sequence>
 <id xml:space="preserve">M</id>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">
 p<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript>
 </id>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">
 u<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript>
 </id>
 </sequence>
 </apply>
 </apply>
 <unitOverride>
 <placeholder />
 </unitOverride>
 </eval>
 </math>
</region>

6. Покажите, что P10T x M x P10 = L10 x L10T.

<region id="ID0EN3IFB" actualWidth="165.17000000000002" actualHeight="62.800000000000004" top="940.80000000000018" left="38.400000000000006" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="82"> <define xmlns="http://schemas.mathsoft.com/math50"> <id labels="VARIABLE" xml:space="preserve">P10</id> <eval> <apply> <indexer /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">m10</id> <real>0</real> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> </define> <resultFormat> <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" /> </resultFormat> </math> </region>	<region id="ID0E23IFB" actualWidth="260.32666666666665" actualHeight="62.800000000000004" top="1017.6000000000001" left="38.400000000000006" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="84"> <define xmlns="http://schemas.mathsoft.com/math50"> <id labels="VARIABLE" xml:space="preserve">L10</id> <eval> <apply> <indexer /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">m10</id> <real>1</real> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> </define> <resultFormat> <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="true"> <nestedMatrixFormats> <nestedMatrixFormat size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="true" position="0,0" /> </nestedMatrixFormats> </matrix> </resultFormat> </math> </region>
<region id="ID0EO4IFB" actualWidth="202.72666666666666" actualHeight="62.800000000000004" top="1094.4" left="38.400000000000006" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="86"> <eval xmlns="http://schemas.mathsoft.com/math50"> <apply> <mult /> <apply> <mult /> <apply> <transpose /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">P10</id> </apply> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">M</id> </apply> <id xml:space="preserve">P10</id> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> </math> </region>	<region id="ID0E34IFB" actualWidth="177.81666666666666" actualHeight="62.800000000000004" top="1094.4" left="307.20000000000005" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="88"> <eval xmlns="http://schemas.mathsoft.com/math50"> <apply> <mult /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">L10</id> <apply> <transpose /> <id xml:space="preserve">L10</id> </apply> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> </math> </region>
<region id="ID0EP5IFB" actualWidth="218.39000000000002" actualHeight="27.752000000000002" top="1152.0000000000002" left="38.400000000000006" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="90"> <eval xmlns="http://schemas.mathsoft.com/math50"> <apply> <equal /> <apply> <mult /> <apply> <mult /> <apply> <transpose /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">P10</id> </apply> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">M</id> </apply> <id xml:space="preserve">P10</id> </apply> <apply> <mult /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">L10</id> <apply> <transpose /> <id xml:space="preserve">L10</id> </apply> </apply> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> </math> </region>

Отношение логически истинно.

7. Используйте функцию Cholesky, чтобы выполнить разложение на множители матрицы M - без выбора главного элемента и нижним разложением на множители (по умолчанию).

<region id="ID0EUAJFB" actualWidth="344.06666666666672" actualHeight="62.800000000000004" top="1305.6000000000001" left="38.400000000000006" xmlns="http://schemas.mathsoft.com/worksheet50">
 <math resultRef="92">
 <define xmlns="http://schemas.mathsoft.com/math50">
 <id labels="VARIABLE" xml:space="preserve">m00</id>
 <eval>
 <apply>
 <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">Cholesky</id>
 <sequence>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">M</id>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">
 p<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript>
 </id>
 </sequence>
 </apply>
 <unitOverride>
 <placeholder />
 </unitOverride>
 </eval>
 </define>
 <resultFormat>
 <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="true" />
 </resultFormat>
 </math>
</region>

Если не указывать аргумент u, как в Cholesky(M, 0), это эквивалентно заданию для него значения 0, как в Cholesky(M, 0, 0).

<region id="ID0ENCJFB" actualWidth="307.08666666666676" actualHeight="27.791200000000003" top="1401.6000000000001" left="76.800000000000011" xmlns="http://schemas.mathsoft.com/worksheet50">
 <math resultRef="94">
 <eval xmlns="http://schemas.mathsoft.com/math50">
 <apply>
 <equal />
 <apply>
 <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">Cholesky</id>
 <sequence>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">M</id>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">
 p<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript>
 </id>
 </sequence>
 </apply>
 <apply>
 <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">Cholesky</id>
 <sequence>
 <id xml:space="preserve">M</id>
 <id xml:space="preserve">
 p<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript>
 </id>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">
 u<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript>
 </id>
 </sequence>
 </apply>
 </apply>
 <unitOverride>
 <placeholder />
 </unitOverride>
 </eval>
 <resultFormat>
 <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" />
 </resultFormat>
 </math>
</region>

Возвращенная нижняя матрица L10, когда выбор главного элемента включен, не равна возвращенной нижней матрице L00, когда выбор главного элемента отключен.

Отношение логически ложно.

8. Покажите, что M = L00 x L00T.

Отношение логически истинно.

9. Используйте функцию Cholesky, чтобы выполнить разложение на множители матрицы M - с выбором главного элемента и верхним разложением на множители.

<region id="ID0EQLJFB" actualWidth="353.51333333333338" actualHeight="118.60000000000001" top="1881.6000000000004" left="38.400000000000006" xmlns="http://schemas.mathsoft.com/worksheet50">
 <math resultRef="106">
 <define xmlns="http://schemas.mathsoft.com/math50">
 <id labels="VARIABLE" xml:space="preserve">m11</id>
 <eval>
 <apply>
 <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">Cholesky</id>
 <sequence>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">M</id>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">
 p<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript>
 </id>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">
 u<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript>
 </id>
 </sequence>
 </apply>
 <unitOverride>
 <placeholder />
 </unitOverride>
 </eval>
 </define>
 <resultFormat>
 <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="true" />
 </resultFormat>
 </math>
</region>

10. Покажите, что P11T x M x P11 = U11T x U11.

Отношение логически истинно.

11. Используйте функцию Cholesky, чтобы выполнить разложение на множители матрицы M - без выбора главного элемента и с верхним разложением на множители.

<region id="ID0EZSJFB" actualWidth="366.78000000000003" actualHeight="62.800000000000004" top="2304.0000000000005" left="38.400000000000006" xmlns="http://schemas.mathsoft.com/worksheet50">
 <math resultRef="118">
 <define xmlns="http://schemas.mathsoft.com/math50">
 <id labels="VARIABLE" xml:space="preserve">m01</id>
 <eval>
 <apply>
 <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">Cholesky</id>
 <sequence>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">M</id>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">
 p<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript>
 </id>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">
 u<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript>
 </id>
 </sequence>
 </apply>
 <unitOverride>
 <placeholder />
 </unitOverride>
 </eval>
 </define>
 <resultFormat>
 <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="true" />
 </resultFormat>
 </math>
</region>

12. Покажите, что M = U01T x U01.

Отношение логически истинно.

Было ли это полезно?