Example: Finding the Left and Right Eigenvectors

Functions > Vector and Matrix > Array Characteristics > Example: Finding the Left and Right Eigenvectors

To avoid logical mismatches when performing boolean comparisons, enable Approximate Equality in the Calculation Options drop-down list.

1. Define a square real matrix.

2. Use the basic definition of eigenval to find the lambdas of matrix A.

<region id="ID0E3XSFB" actualWidth="220.30666666666667" actualHeight="39.2" top="211.20000000000002" left="0" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="26">
    <symEval xmlns="http://schemas.mathsoft.com/math50">
      <apply>
        <absval />
        <apply>
          <minus />
          <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">A</id>
          <apply>
            <mult />
            <id xml:space="preserve">λ</id>
            <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">I</id>
          </apply>
        </apply>
      </apply>
      <command>
        <id xml:space="preserve">factor</id>
      </command>
      <symResult>
        <apply>
          <mult />
          <parens>
            <apply>
              <minus />
              <id labels="VARIABLE" xml:space="preserve">λ</id>
              <real>4</real>
            </apply>
          </parens>
          <parens>
            <apply>
              <plus />
              <id labels="VARIABLE" xml:space="preserve">λ</id>
              <real>3</real>
            </apply>
          </parens>
        </apply>
      </symResult>
    </symEval>
    <formattingOverride BackgroundColor="#00FFFFFF" />
  </math>
</region>

3. Alternatively, use the function eigenvals to find the lambdas of matrix A.

Finding the Right eigenvectors

1. Use the function eigenvecs to find the right (default) eigenvectors of matrix A.

2. Check the deviation between A x Rn and λn x Rn

<region id="ID0EWBTFB" actualWidth="120.21666666666667" actualHeight="44.2" top="921.60000000000014" left="0" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="33"> <eval xmlns="http://schemas.mathsoft.com/math50"> <apply> <mult /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">A</id> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">R<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript></Span> </id> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> </math> </region>	<region id="ID0EGCTFB" actualWidth="123.25" actualHeight="44.2" top="921.60000000000014" left="172.8" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="34"> <eval xmlns="http://schemas.mathsoft.com/math50"> <apply> <mult /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">λ<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript></Span> </id> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">R<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript></Span> </id> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> </math> </region>	<region id="ID0EXCTFB" actualWidth="146.17666666666668" actualHeight="44.2" top="921.60000000000014" left="345.6" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="35"> <eval xmlns="http://schemas.mathsoft.com/math50"> <apply> <minus /> <apply> <mult /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">A</id> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">R<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript></Span> </id> </apply> <apply> <mult /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">λ<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript></Span> </id> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">R<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript></Span> </id> </apply> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> <resultFormat> <decimal precision="3" show-trailing-zeros="false" radix="dec" imaginary-value="i" /> </resultFormat> </math> </region>
<region id="ID0EKDTFB" actualWidth="130.96666666666664" actualHeight="44.2" top="979.20000000000016" left="0" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="36"> <eval xmlns="http://schemas.mathsoft.com/math50"> <apply> <mult /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">A</id> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">R<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript></Span> </id> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> <resultFormat> <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" /> </resultFormat> </math> </region>	<region id="ID0EZDTFB" actualWidth="134" actualHeight="44.2" top="979.20000000000016" left="172.8" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="37"> <eval xmlns="http://schemas.mathsoft.com/math50"> <apply> <mult /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">λ<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript></Span> </id> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">R<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript></Span> </id> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> <resultFormat> <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" /> </resultFormat> </math> </region>	<region id="ID0EJETFB" actualWidth="146.17666666666668" actualHeight="44.2" top="979.20000000000016" left="345.6" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="38"> <eval xmlns="http://schemas.mathsoft.com/math50"> <apply> <minus /> <apply> <mult /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">A</id> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">R<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript></Span> </id> </apply> <apply> <mult /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">λ<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript></Span> </id> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">R<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript></Span> </id> </apply> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> </math> </region>

The deviation is zero.

Finding the Left eigenvectors

1. Use function eigenvecs to find the left eigenvectors of matrix A.

2. Check the deviation between LnT x A and λn x LnT.

<region id="ID0EJLTFB" actualWidth="169.52333333333331" actualHeight="29.943200000000004" top="1382.4" left="0" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="42"> <eval xmlns="http://schemas.mathsoft.com/math50"> <apply> <mult /> <apply> <transpose /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">L<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript></Span> </id> </apply> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">A</id> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> </math> </region>	<region id="ID0EZLTFB" actualWidth="184.89" actualHeight="29.3712" top="1382.4" left="172.8" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="43"> <eval xmlns="http://schemas.mathsoft.com/math50"> <apply> <mult /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">λ<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript></Span> </id> <apply> <transpose /> <parens> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">L<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript></Span> </id> </parens> </apply> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> </math> </region>	<region id="ID0EKMTFB" actualWidth="174.46333333333337" actualHeight="29.943200000000004" top="1382.4" left="403.20000000000005" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="44"> <eval xmlns="http://schemas.mathsoft.com/math50"> <apply> <minus /> <apply> <mult /> <apply> <transpose /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">L<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript></Span> </id> </apply> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">A</id> </apply> <apply> <mult /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">λ<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript></Span> </id> <apply> <transpose /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">L<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript></Span> </id> </apply> </apply> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> </math> </region>
<region id="ID0E4MTFB" actualWidth="180.27333333333331" actualHeight="29.943200000000004" top="1440.0000000000002" left="0" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="45"> <eval xmlns="http://schemas.mathsoft.com/math50"> <apply> <mult /> <apply> <transpose /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">L<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript></Span> </id> </apply> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">A</id> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> <resultFormat> <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" /> </resultFormat> </math> </region>	<region id="ID0EMNTFB" actualWidth="195.64" actualHeight="29.3712" top="1440.0000000000002" left="172.8" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="46"> <eval xmlns="http://schemas.mathsoft.com/math50"> <apply> <mult /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">λ<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript></Span> </id> <apply> <transpose /> <parens> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">L<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript></Span> </id> </parens> </apply> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> <resultFormat> <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" /> </resultFormat> </math> </region>	<region id="ID0E3NTFB" actualWidth="174.46333333333337" actualHeight="29.943200000000004" top="1440.0000000000002" left="403.20000000000005" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="47"> <eval xmlns="http://schemas.mathsoft.com/math50"> <apply> <minus /> <apply> <mult /> <apply> <transpose /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">L<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript></Span> </id> </apply> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">A</id> </apply> <apply> <mult /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">λ<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript></Span> </id> <apply> <transpose /> <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">L<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript></Span> </id> </apply> </apply> </apply> <unitOverride> <placeholder /> </unitOverride> </eval> <resultFormat> <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" /> </resultFormat> </math> </region>

The deviation is zero.