示例：协方差和主分量分析

函数 > 向量和矩阵 > 数组特征 > 示例：协方差和主分量分析

示例：协方差和主分量分析

计算数据集的协方差，并使用 eigenvals 和 eigenvec 函数来执行主分量分析。

1. 定义以下矩阵 D。

D 的每一行代表一个观察值，每一列代表一个测量的特征。

2. 计算样本的协方差矩阵。

3. 计算和排序特征值。

<region id="ID0EIBUFB" top="921.6" left="38.4" width="233.816666666667" height="25.6" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="5">
    <define xmlns="http://schemas.mathsoft.com/math50">
      <id label="None" xml:space="preserve">V</id>
      <apply>
        <id label="Function" xml:space="preserve" label-is-contextual="true">reverse</id>
        <apply>
          <id label="Function" xml:space="preserve" label-is-contextual="true">sort</id>
          <apply>
            <id label="Function" xml:space="preserve" label-is-contextual="true">eigenvals</id>
            <id label="Variable" xml:space="preserve" label-is-contextual="true">S1</id>
          </apply>
        </apply>
      </apply>
    </define>
    <resultFormat>
      <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
      <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" />
    </resultFormat>
  </math>
</region>

4. 计算变换矩阵。

<region id="ID0EGDUFB" top="1075.2" left="38.4" width="150.513333333333" height="36.8" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="7">
    <define xmlns="http://schemas.mathsoft.com/math50">
      <apply>
        <indexer />
        <id label="Variable" xml:space="preserve">x</id>
        <id label="None" xml:space="preserve">j</id>
      </apply>
      <apply>
        <id label="Function" xml:space="preserve" label-is-contextual="true">eigenvec</id>
        <sequence>
          <id label="Variable" xml:space="preserve" label-is-contextual="true">S1</id>
          <apply>
            <indexer />
            <id label="Variable" xml:space="preserve" label-is-contextual="true">V</id>
            <id xml:space="preserve">j</id>
          </apply>
        </sequence>
      </apply>
    </define>
    <resultFormat>
      <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
      <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" />
    </resultFormat>
  </math>
</region>

<region id="ID0EUDUFB" top="1132.8" left="38.4" width="166.673333333333" height="36.8" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="8">
    <define xmlns="http://schemas.mathsoft.com/math50">
      <id label="None" xml:space="preserve">T</id>
      <apply>
        <id label="Function" xml:space="preserve" label-is-contextual="true">augment</id>
        <sequence>
          <apply>
            <indexer />
            <id label="Variable" xml:space="preserve" label-is-contextual="true">x</id>
            <real>0</real>
          </apply>
          <apply>
            <indexer />
            <id label="Variable" xml:space="preserve" label-is-contextual="true">x</id>
            <real>1</real>
          </apply>
          <apply>
            <indexer />
            <id label="Variable" xml:space="preserve" label-is-contextual="true">x</id>
            <real>2</real>
          </apply>
        </sequence>
      </apply>
    </define>
    <resultFormat>
      <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
      <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" />
    </resultFormat>
  </math>
</region>

5. 变换原始数据。

6. 计算变换数据的协方差矩阵。

数据的主分量是矩阵 S2 的对角元素。