Example: Special Characteristics of Matrices

Vectors, Matrices, and Tables > Working With Arrays > Example: Special Characteristics of Matrices

Example: Special Characteristics of Matrices

Find the trace, rank, generalized inverse, norms, and condition numbers of a square matrix.

Click to copy this expression

The Trace, Rank and Generalized Inverse of a Matrix

1. Use the tr function to find the trace, or the sum of the diagonal elements, of M.

Click to copy this expression

2. Use the rank function to find the rank of the real-valued matrix M.

Click to copy this expression

3. Use the geninv function to find the generalized inverse of matrix M.

Click to copy this expression

The different Norms of a Matrix

1. Find the L1 norm of M, and compare the result with the output of function norm1

The L1 norm is the maximum of the absolute column sums (max taken over j= 0, 1, 2).

Click to copy this expression

Click to copy this expression

Click to copy this expression

Click to copy this expression

<region id="ID0EWBNFB" top="1056.0000000000002" left="19.200000000000003" height="136.00000000000003" width="191.66000000000003" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="53">
    <eval xmlns="http://schemas.mathsoft.com/math50">
      <apply>
        <id xml:space="preserve">max</id>
        <matrix rows="3" cols="1">
          <apply>
            <summation />
            <lambda>
              <boundVars>
                <placeholder />
              </boundVars>
              <apply>
                <matcol />
                <parens>
                  <apply>
                    <vectorize />
                    <apply>
                      <absval />
                      <id xml:space="preserve">M</id>
                    </apply>
                  </apply>
                </parens>
                <real>0</real>
              </apply>
            </lambda>
            <upperBound>
              <placeholder />
            </upperBound>
          </apply>
          <apply>
            <summation />
            <lambda>
              <boundVars>
                <placeholder />
              </boundVars>
              <apply>
                <matcol />
                <parens>
                  <apply>
                    <vectorize />
                    <apply>
                      <absval />
                      <id xml:space="preserve">M</id>
                    </apply>
                  </apply>
                </parens>
                <real>1</real>
              </apply>
            </lambda>
            <upperBound>
              <placeholder />
            </upperBound>
          </apply>
          <apply>
            <summation />
            <lambda>
              <boundVars>
                <placeholder />
              </boundVars>
              <apply>
                <matcol />
                <parens>
                  <apply>
                    <vectorize />
                    <apply>
                      <absval />
                      <id xml:space="preserve">M</id>
                    </apply>
                  </apply>
                </parens>
                <real>2</real>
              </apply>
            </lambda>
            <upperBound>
              <placeholder />
            </upperBound>
          </apply>
        </matrix>
      </apply>
      <unitOverride>
        <placeholder />
      </unitOverride>
    </eval>
  </math>
</region>

Click to copy this expression

2. Use the norm2 function to find the L2 norm of M.

Click to copy this expression

3. Use the norme function to find the Euclidean norm of M:

The Euclidean norm for a matrix is analogous to that for a vector:

Click to copy this expression

<region id="ID0EAGNFB" top="1440.0000000000002" left="38.400000000000006" height="60.441333557128914" width="176.66000000000003" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="59">
    <eval xmlns="http://schemas.mathsoft.com/math50">
      <apply>
        <pow />
        <parens>
          <apply>
            <summation />
            <lambda>
              <boundVars>
                <id xml:space="preserve">i</id>
              </boundVars>
              <apply>
                <summation />
                <lambda>
                  <boundVars>
                    <id xml:space="preserve">j</id>
                  </boundVars>
                  <apply>
                    <pow />
                    <parens>
                      <apply>
                        <indexer />
                        <id xml:space="preserve">M</id>
                        <sequence>
                          <id xml:space="preserve">i</id>
                          <id xml:space="preserve">j</id>
                        </sequence>
                      </apply>
                    </parens>
                    <real>2</real>
                  </apply>
                </lambda>
                <lowerBound>
                  <real>0</real>
                </lowerBound>
                <upperBound>
                  <real>2</real>
                </upperBound>
              </apply>
            </lambda>
            <lowerBound>
              <real>0</real>
            </lowerBound>
            <upperBound>
              <real>2</real>
            </upperBound>
          </apply>
        </parens>
        <apply>
          <div />
          <real>1</real>
          <real>2</real>
        </apply>
      </apply>
      <unitOverride>
        <placeholder />
      </unitOverride>
    </eval>
  </math>
</region>

Click to copy this expression

4. Find the Infinity norm of M, and compare the result with the output of function normi.

The Infinity norm is the maximum of the absolute row sums (max taken over i=0, 1, 2)

Click to copy this expression

Click to copy this expression

<region id="ID0E4INFB" top="1766.4000000000003" left="19.200000000000003" height="143.056" width="200.52000000000004" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="65">
    <eval xmlns="http://schemas.mathsoft.com/math50">
      <apply>
        <id xml:space="preserve">max</id>
        <matrix rows="3" cols="1">
          <apply>
            <summation />
            <lambda>
              <boundVars>
                <placeholder />
              </boundVars>
              <apply>
                <matcol />
                <parens>
                  <apply>
                    <vectorize />
                    <apply>
                      <absval />
                      <apply>
                        <transpose />
                        <id xml:space="preserve">M</id>
                      </apply>
                    </apply>
                  </apply>
                </parens>
                <real>0</real>
              </apply>
            </lambda>
            <upperBound>
              <placeholder />
            </upperBound>
          </apply>
          <apply>
            <summation />
            <lambda>
              <boundVars>
                <placeholder />
              </boundVars>
              <apply>
                <matcol />
                <parens>
                  <apply>
                    <vectorize />
                    <apply>
                      <absval />
                      <apply>
                        <transpose />
                        <id xml:space="preserve">M</id>
                      </apply>
                    </apply>
                  </apply>
                </parens>
                <real>1</real>
              </apply>
            </lambda>
            <upperBound>
              <placeholder />
            </upperBound>
          </apply>
          <apply>
            <summation />
            <lambda>
              <boundVars>
                <placeholder />
              </boundVars>
              <apply>
                <matcol />
                <parens>
                  <apply>
                    <vectorize />
                    <apply>
                      <absval />
                      <apply>
                        <transpose />
                        <id xml:space="preserve">M</id>
                      </apply>
                    </apply>
                  </apply>
                </parens>
                <real>2</real>
              </apply>
            </lambda>
            <upperBound>
              <placeholder />
            </upperBound>
          </apply>
        </matrix>
      </apply>
      <unitOverride>
        <placeholder />
      </unitOverride>
    </eval>
  </math>
</region>

Click to copy this expression

The Different Condition Numbers of a Matrix

The Condition number of a matrix is the product of two matrix norms. It measures the sensitivity of a linear system solution to errors in the input vector:

Click to copy this expression

1. Use the cond1 function to find the L1 condition number of M.

Click to copy this expression

2. Use the cond2 function to find the L2 condition number of M.

Click to copy this expression

3. Use the conde function to find the Euclidean condition number of M.

Click to copy this expression

4. Use the condi function to find the Infinity condition number of M.

Click to copy this expression