【考研数学】线性代数第五章 —— 特征值和特征向量(1,理论背景与基本概念)


引言

特征值与特征向量是线性代数的重要内容,它的主要应用有:求矩阵的幂、矩阵的对角化、二次型的标准型。


一、特征值和特征向量的基本概念

1.1 理论背景

在一个多项式中,未知数的个数为任意多个,且每一项的次数都是 2 的多项式称为二次型,其有两种类型:非标准二次型和标准二次型。

什么是标准和非标准呢?我们看两个例子。

f ( x 1 , x 2 , x 3 ) = 2 x 1 2 − x 2 2 + 3 x 3 2 + 2 x 1 x 2 − 4 x 2 x 3 f(x_1,x_2,x_3)=2x_1^2-x_2^2+3x_3^2+2x_1x_2-4x_2x_3 f(x1,x2,x3)=2x12x22+3x32+2x1x24x2x3 为非标准二次型,而 g ( x 1 , x 2 , x 3 ) = 2 x 1 2 − x 2 2 + 3 x 3 2 g(x_1,x_2,x_3)=2x_1^2-x_2^2+3x_3^2 g(x1,x2,x3)=2x12x22+3x32 为标准二次型。多项式 f f f 尽管每项也都是 2 次,但含有杂项,所以不是标准型。

f f f 进行矩阵化,即将其写成几个矩阵相乘的形式,有 f = X T A X f=X^TAX f=XTAX ,其中 A = [ 2 1 0 1 − 1 − 2 0 − 2 3 ] , X = [ x 1 x 2 x 3 ] . A=\begin{bmatrix} 2 & 1 & 0\\ 1 & -1 & -2 \\0 & -2 & 3 \end{bmatrix},X=\begin{bmatrix} x_1 \\ x_2 \\ x_3 \end{bmatrix}. A= 210112023 ,X= x1x2x3 . 同理,可以将 g g g 进行矩阵化,有 g = X T A X g=X^TAX g=XTAX ,其中 A = [ 2 0 0 0 1 0 0 0 − 3 ] , X = [ x 1 x 2 x 3 ] . A=\begin{bmatrix} 2 & 0 & 0\\ 0 & 1 & 0 \\0 & 0 & -3 \end{bmatrix},X=\begin{bmatrix} x_1 \\ x_2 \\ x_3 \end{bmatrix}. A= 200010003 ,X= x1x2x3 . 从中,我们有以下结论:

  1. 二次型 X T A X X^TAX XTAX 为非标准型的充要条件为 A A A 为非对角矩阵但 A T = A A^T=A AT=A ;二次型 X T A X X^TAX XTAX 为标准型的充要条件为 A A A对角矩阵
  2. 将非标准二次型化为标准型,等价于将非对角矩阵化为对角矩阵,称这一过程为对角化。特征值与特征向量的理论即矩阵对角化理论。

1.2 基本概念

特征值与特征向量 —— 设 A A A n n n 阶矩阵,若存在常数 λ \lambda λ n n n 维非零列向量 α \pmb{\alpha} α ,使得 A α = λ α A\pmb{\alpha}=\lambda\pmb{\alpha} Aα=λα λ \lambda λ 为矩阵 A A A 的特征值, α \pmb{\alpha} α 为矩阵 A A A 的属于特征值 λ \lambda λ 的特征向量。

A α = λ α A\pmb{\alpha}=\lambda\pmb{\alpha} Aα=λα 等价与 ( λ E − A ) α = 0 (\lambda E-A)\pmb{\alpha}=0 (λEA)α=0 ,由 α \pmb{\alpha} α 为非零向量,得到方程组 ( λ E − A ) X = 0 (\lambda E-A)X=0 (λEA)X=0 有非零解。根据前面方程组的知识,齐次方程组有非零解,等价于系数矩阵非满秩,即 r ( λ E − A ) < n r(\lambda E-A)<n r(λEA)<n ,或 ∣ λ E − A ∣ = 0 |\lambda E-A|=0 λEA=0

λ \lambda λ 为矩阵 A A A 的特征值,那么一定有 ∣ λ E − A ∣ = 0 |\lambda E-A|=0 λEA=0 ;反之,若 ∣ λ E − A ∣ = 0 |\lambda E-A|=0 λEA=0 ,则齐次线性方程组 ( λ E − A ) X = 0 (\lambda E-A)X=0 (λEA)X=0 有非零解,从而存在非零向量 α \pmb{\alpha} α ,使得 ( λ E − A ) α = 0 (\lambda E-A)\pmb{\alpha}=0 (λEA)α=0 ,即 A α = λ α A\pmb{\alpha}=\lambda\pmb{\alpha} Aα=λα ,故 λ \lambda λ 为矩阵 A A A 的特征值。

特征方程 —— 设 A = [ a 11 a 12 ⋯ a 1 n a 21 a 22 ⋯ a 2 n ⋮ ⋮ ⋮ a n 1 a n 2 ⋯ a n n ] A=\begin{bmatrix} a_{11} & a_{12} & \cdots & a_{1n}\\ a_{21} & a_{22} & \cdots & a_{2n} \\ \vdots & \vdots & & \vdots \\ a_{n1} & a_{n2} & \cdots & a_{nn} \end{bmatrix} A= a11a21an1a12a22an2a1na2nann ,称 ∣ λ E − A ∣ = 0 |\lambda E-A|=0 λEA=0 为矩阵 A A A 的特征方程。

对角线元素之和,即 a 11 + a 22 + ⋯ + a n n a_{11}+a_{22}+\cdots+a_{nn} a11+a22++ann 称为矩阵 A A A 的迹,记作 t r ( A ) tr(A) tr(A)

矩阵相似 —— 设 A , B A,B A,B n n n 阶矩阵,若存在可逆矩阵 P P P ,使得 P − 1 A P = B P^{-1}AP=B P1AP=B ,称矩阵 A A A 和矩阵 B B B 相似,记为 A ∼ B A\sim B AB 。若存在可逆矩阵 P P P ,使得 P − 1 A P = Λ P^{-1}AP=\Lambda P1AP=Λ ,其中 Λ \Lambda Λ 为对角矩阵,则称 A A A 可以相似对角化。

一个矩阵和本身是相似的;若 A ∼ B A \sim B AB ,则有 B ∼ A B\sim A BA ,且相似具有传递性。

f ( x ) = a n x n + a n − 1 x n − 1 + ⋯ + a 1 x + a 0 f(x)=a_nx^n+a_{n-1}x^{n-1}+\cdots+a_1x+a_0 f(x)=anxn+an1xn1++a1x+a0 ,若 A ∼ B A \sim B AB ,有 A T ∼ B T , f ( A ) ∼ f ( B ) . A^T \sim B^T,f(A) \sim f(B). ATBT,f(A)f(B). 特别地,当 A , B A,B A,B 可逆时,有 A − 1 ∼ B − 1 , A ∗ ∼ B ∗ . A^{-1} \sim B^{-1},A^* \sim B^*. A1B1,AB.

A ∼ B A \sim B AB ,则有 r ( A ) = r ( B ) r(A)=r(B) r(A)=r(B) ,但反过来不一定成立。即秩相等是相似的必要条件,若两个矩阵秩不相等,肯定不相似。

A ∼ B A \sim B AB ,则有 ∣ λ E − A ∣ = ∣ λ E − B ∣ |\lambda E-A|=|\lambda E-B| λEA=λEB ,从而 A , B A,B A,B 的特征值相同,但反过来不一定成立。

A ∼ B A \sim B AB ,有 ∣ A ∣ = ∣ B ∣ , t r ( A ) = t r ( B ) . |A|=|B|,tr(A)=tr(B). A=B,tr(A)=tr(B).

施密特正交化 —— 把一组线性无关的向量组转化为一组两两正交且规范的向量组的过程称为施密特正交化。

α 1 , α 2 , ⋯   , α n \pmb{\alpha_1,\alpha_2,\cdots,\alpha_n} α1,α2,,αn 线性无关,其正交化过程如下:

(1)正交化:

β 1 = α 1 , β 2 = α 2 − ( α 2 , β 1 ) ( β 1 , β 1 ) β 1 , ⋯   , \pmb{\beta_1=\alpha_1,\beta_2=\alpha_2-\frac{(\alpha_2,\beta_1)}{(\beta_1,\beta_1)}\beta_1,\cdots,} β1=α1,β2=α2(β1,β1)(α2,β1)β1,, β n = α n − ( α n , β 1 ) ( β 1 , β 1 ) β 1 − ( α n , β 2 ) ( β 2 , β 2 ) β 2 − ⋯ − ( α n , β n − 1 ) ( β n − 1 , β n − 1 ) β n − 1 . \pmb{\beta_n=\alpha_n-\frac{(\alpha_n,\beta_1)}{(\beta_1,\beta_1)}\beta_1-\frac{(\alpha_n,\beta_2)}{(\beta_2,\beta_2)}\beta_2-\cdots-\frac{(\alpha_n,\beta_{n-1})}{(\beta_{n-1},\beta_{n-1})}\beta_{n-1}}. βn=αn(β1,β1)(αn,β1)β1(β2,β2)(αn,β2)β2(βn1,βn1)(αn,βn1)βn1. 则向量组 β 1 , β 2 , ⋯   , β n \pmb{\beta_1,\beta_2,\cdots,\beta_n} β1,β2,,βn 两两正交。

(2)规范化:

γ 1 = β 1 / ∣ β 1 ∣ , γ 1 = β 2 / ∣ β 2 ∣ , ⋯   , γ n = β n / ∣ β n ∣ . \pmb{\gamma_1=\beta_1/|\beta_1|,\gamma_1=\beta_2/|\beta_2|,\cdots,\gamma_n=\beta_n/|\beta_n|.} γ1=β1/∣β1,γ1=β2/∣β2,,γn=βn/∣βn∣. ,则向量组 γ 1 , γ 2 , ⋯   , γ n \pmb{\gamma_1,\gamma_2,\cdots,\gamma_n} γ1,γ2,,γn 为两两正交且规范的向量组。

正交矩阵 —— 设 Q Q Q n n n 阶矩阵,若 Q T Q = E Q^TQ=E QTQ=E ,称 Q Q Q 为正交矩阵。

显然, Q Q Q 可逆,且 Q − 1 = Q T Q^{-1}=Q^T Q1=QT

Q = ( γ 1 , γ 2 , ⋯   , γ n ) Q=(\pmb{\gamma_1,\gamma_2,\cdots,\gamma_n}) Q=(γ1,γ2,,γn) ,矩阵 Q Q Q 为正交矩阵的充要条件是 γ 1 , γ 2 , ⋯   , γ n \pmb{\gamma_1,\gamma_2,\cdots,\gamma_n} γ1,γ2,,γn 两两正交且规范。

Q Q Q 为正交矩阵,有 ∣ Q ∣ = ± 1 |Q|=\pm1 Q=±1 ,根据定义式两边取行列式即可得到。

Q Q Q 为正交矩阵,则 Q Q Q 的特征值为 ± 1 \pm1 ±1

证明: Q Q Q 为正交矩阵,根据定义,有 Q T Q = E Q^TQ=E QTQ=E 。设矩阵 Q Q Q 的属于特征值 λ \lambda λ 的特征向量为 α \alpha α ,则有 Q α = λ α Q\alpha=\lambda\alpha Qα=λα 。等式两边同时进行转置运算,有 α T Q T = λ α T \alpha^TQ^T=\lambda\alpha^T αTQT=λαT ,等式两边同时乘以 Q α Q\alpha Qα ,有 α T Q T Q ^ α = λ α T Q α ^ \alpha^T\widehat{Q^TQ}\alpha=\lambda\alpha^T\widehat{Q\alpha} αTQTQ α=λαTQα ,即 α T E α = λ α T λ α \alpha^TE\alpha=\lambda\alpha^T\lambda\alpha αTEα=λαTλα ,可得到 α T α = λ 2 α T α \alpha^T\alpha=\lambda^2\alpha^T\alpha αTα=λ2αTα 。由 α \alpha α 为特征向量,有 α ≠ 0 \alpha \ne 0 α=0 ,所以 λ 2 = 1 \lambda^2=1 λ2=1 ,原命题得证。


写在最后

就先到这里吧,下一篇文章我们来学习特征值与特征向量的性质。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
LAPACK是一种线性代数库,用于解决各种数值线性代数问题,包括矩阵特征值特征向量的求解。下面我们来介绍一下LAPACK库中求解矩阵特征值特征向量的函数。 1. DSYEV函数 DSYEV函数用于求解实对称矩阵的特征值特征向量。该函数的原型如下: ``` void dsyev(char jobz, char uplo, int n, double* a, int lda, double* w, double* work, int lwork, int* info); ``` 参数说明: - jobz:指定计算特征值还是特征向量,取值为'N'(只计算特征值)或'V'(同时计算特征值特征向量)。 - uplo:指定矩阵的上三角部分还是下三角部分存储在数组a中,取值为'U'(上三角部分)或'L'(下三角部分)。 - n:矩阵的维数。 - a:存储矩阵的一维数组。 - lda:指定a数组中每个列向量的存储长度(通常为n)。 - w:存储特征值的一维数组。 - work:工作空间数组。 - lwork:指定work数组的长度(通常为3n)。 - info:返回求解结果,取值为0表示成功,其他值表示出错。 2. ZGEEV函数 ZGEEV函数用于求解复矩阵的特征值特征向量。该函数的原型如下: ``` void zgeev(char jobvl, char jobvr, int n, std::complex<double>* a, int lda, std::complex<double>* w, std::complex<double>* vl, int ldvl, std::complex<double>* vr, int ldvr, std::complex<double>* work, int lwork, double* rwork, int* info); ``` 参数说明: - jobvl:指定是否计算左特征向量,取值为'N'(不计算)或'V'(计算)。 - jobvr:指定是否计算右特征向量,取值为'N'或'V'。 - n:矩阵的维数。 - a:存储矩阵的一维数组。 - lda:指定a数组中每个列向量的存储长度(通常为n)。 - w:存储特征值的一维数组。 - vl:存储左特征向量的一维数组。 - ldvl:指定vl数组中每个列向量的存储长度(通常为n)。 - vr:存储右特征向量的一维数组。 - ldvr:指定vr数组中每个列向量的存储长度(通常为n)。 - work:工作空间数组。 - lwork:指定work数组的长度(通常为2n)。 - rwork:实数数组,长度为2n(用于存储中间计算结果)。 - info:返回求解结果,取值为0表示成功,其他值表示出错。 以上就是LAPACK库中求解矩阵特征值特征向量的函数介绍。需要注意的是,在调用这些函数之前,需要先将矩阵按列存储方式存储在一维数组中,并传入一些参数,如矩阵的维数、存储方式等。具体的参数可以参考LAPACK库的文档。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Douglassssssss

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值