zoukankan      html  css  js  c++  java
  • 线性代数

    特征向量和特征值

    概念

    对于一个线性变换,若一个向量经过这个线性变换后,只是拉伸或者压缩,并没有离开该向量张成的空间,则称这个向量为该线性变换的特征向量,其中特征向量拉伸或者压缩的比例称为特征值。有些线性变换不存在特征向量,如旋转。

    考虑一个线性变换 (A),设 (v) 为它的某个特征向量,特征值为 (lambda),得 (Av=lambda v)。这里等式左边为矩阵乘法,右边为向量数乘,将右边也转化为矩阵乘法的形式,变形后得 ((A-lambda I)v=0),其中 (I) 为单位矩阵,也就是恒等变换。首先有 (v) 为零向量时,该等式成立。当 (v) 不是零向量时,因为有一个非零向量经过线性变换后变为了一个零向量,因此该线性变换一定降维,进一步得 ( ext{det}(A-Ilambda)=0),解该方程即可得特征值。

    应用

    [large A=egin{bmatrix} frac{1}{n} & 0 & 0 & 0 & cdots & 0 \ frac{1}{n} & frac{1}{n-1} & 0 & 0 & cdots & 0 \ frac{1}{n} & frac{1}{n-1} & frac{1}{n-2} & 0 & cdots & 0 \ frac{1}{n} & frac{1}{n-1} & frac{1}{n-2} & frac{1}{n-3} & cdots & 0 \ vdots & vdots & vdots & vdots & ddots & vdots \ frac{1}{n} & frac{1}{n-1} & frac{1}{n-2} & frac{1}{n-3} & cdots & 1 end{bmatrix} ]

    快速计算 (A^k)

    (E) 为该矩阵的特征向量组成的矩阵,即其第 (i) 列为第 (i) 个特征向量,(D) 为对角矩阵,对角线上的第 (i) 个元素为第 (i) 个特征值,得:

    [largeegin{aligned} AE&=ED\ A&=EDE^{-1}\ A^k&=left(EDE^{-1} ight)^k\ A^k&=ED^kE^{-1}\ end{aligned} ]

    这称为矩阵的对角化,当一个 (n) 阶矩阵有 (n) 个线性无关的特征向量时,可以使用对角化。那么得出 (A) 的特征向量和特征值后就能快速计算 (A^k) 了。对于 (A) 有:

    [large ext{det}(A-Ilambda)=(1-lambda)(frac{1}{2}-lambda)dots(frac{1}{n}-lambda)=0 ]

    (A)(n) 个特征向量,第 (i) 个特征向量的特征值为 (frac{1}{i})。设 (v_{i,j}) 表示该矩阵的第 (i) 个特征向量的倒数第 (j) 个元素,得:

    [large v_{i,j}=(-1)^{i+j}inom{i-1}{j-1} ]

    考虑证明,代入 (Av=lambda v) 得:

    [largeegin{aligned} frac{(-1)^{i+j}}{i}inom{i-1}{j-1}=&sum_{k=j}^ifrac{1}{k}(-1)^{i+k}inom{i-1}{k-1}\ =&frac{1}{i}sum_{k=j}^i(-1)^{i+k}inom{i}{k}\ =&frac{1}{i}sum_{k=0}^{i-j}(-1)^{k}inom{i}{k}\ =&frac{1}{i}sum_{k=0}^{i-j}inom{k-i-1}{k}\ =&frac{1}{i}inom{-j}{i-j}\ =&frac{(-1)^{i+j}}{i}inom{i-1}{j-1} end{aligned} ]

    通过卷积就能 (O(nlog n)) 计算 (A^k) 了。

    CF923E Perpetual Subtraction

    yhx 的题解

    CF 官方题解

  • 相关阅读:
    什么是JAVA?JAVA能做什么?
    vagrant 虚拟机搭建 以及virtualbox 安装 配置
    sourcetree跳过注册的方法
    layuiAdmin安装
    Navicat获取注册码
    laravel artisan命令大全
    vagrant命令
    安装python 3.7
    搭建golang webcron 定时任务管理平台
    mysql创建用户
  • 原文地址:https://www.cnblogs.com/lhm-/p/14347059.html
Copyright © 2011-2022 走看看