线代:1.6矩阵的特征值和特征向量
文章目录
- 任务详解:
- 1.向量的内积和范数
- 向量的内积以及正交性
- 定义1:
- 定义2
- 定义3
- 定义4
- 判定矩阵A可逆的小结
- 2.特征值特征向量以及矩阵的相似
- 方阵的特征值与特征向量
- 定义6
本课程来自 深度之眼,部分截图来自课程视频。
【第一章 线性代数】1.6矩阵的特征值和特征向量
在线LaTeX公式编辑器
任务详解:
1、掌握内积,正交,线性相关,线性无关的概念
2、掌握规范正交基,正交矩阵
3、掌握特征值特征向量的几何意义与算法
1.向量的内积和范数
向量的内积以及正交性
定义1:
设有n维向量(如果不做特殊说明,n维向量都是指列向量)
[x,y]称为向量x与y的内积(或者叫点积,elementwise).
内积是两个向量之间的一种运算,其结果是一个实数,用矩阵记号表示,当x与y都是列向量时,有
[x,y]=xTy=yTx[x,y]=x^Ty=y^Tx[x,y]=xTy=yTx
有时候也记做:<x,y>
还有一个重要性质:柯西不等式
[x,y]2≤[x,x][y,y][x,y]^2≤[x,x][y,y][x,y]2≤[x,x][y,y]
由以上性质加上我们中学在二维空间里面向量夹角的概念,我们可以推广到高维空间,也可以用来衡量高维空间中两个样本的相似度的一种度量(不同于欧式距离)。
定义2
令
∣∣x∣∣=[x,x]=x12+x22+...+xn2||x||=\sqrt{[x,x]}=\sqrt{x_1^2+x_2^2+...+x_n^2}∣∣x∣∣=[x,x]=x12+x22+...+xn2
∣∣x∣∣||x||∣∣x∣∣称为n维向量x的长度或者范数或者模长
当∣∣x∣∣=1||x||=1∣∣x∣∣=1时,称x为单位向量。
向量的长度具有下述性质:
(i)非负性:当x≠0时,∣∣x∣∣>0||x||>0∣∣x∣∣>0;当x=0时,∣∣x∣∣=0||x||=0∣∣x∣∣=0;
(i)齐次性:∣λx∣=∣λ∣∣∣x∣∣|\lambda x|=|\lambda|||x||∣λx∣=∣λ∣∣∣x∣∣;右边的实数外面是绝对值
(ii)三角不等式:∣∣x+y∣∣≤∣∣x∣∣+∣∣y∣∣||x+y||≤||x||+||y||∣∣x+y∣∣≤∣∣x∣∣+∣∣y∣∣。
当[x,y]=0时,称向量x与y正交(二维上看就在垂直关系).显然,若x=0,则x与任何向量都正交。
定理1:若n维向量a1,a2,…,ana_1,a_2,…,a_na1,a2,…,an是一组两两正交的非零向量([ai,aj]=0,i≠j[a_i,a_j]=0,i\neq j[ai,aj]=0,i=j),则a1,a2,…,ana_1,a_2,…,a_na1,a2,…,an线性无关.
以下是百度百科中的线性无关定义:
在向量空间V的一组向量A: a1,a2,⋅⋅⋅,ama_1, a_2, ···,a_ma1,a2,⋅⋅⋅,am如果存在不全为零的数 k1,k2,⋅⋅⋅,kmk_1, k_2, ···,k_mk1,k2,⋅⋅⋅,km , 使
k1a1+k2a2+...+kmam=0k_1a_1+k_2a_2+...+k_ma_m=0k1a1+k2a2+...+kmam=0
则称向量组A是线性相关的 ,否则数 k1,k2,⋅⋅⋅,kmk_1, k_2, ···,k_mk1,k2,⋅⋅⋅,km全为0时,称它是线性无关。
由此定义看出 是否线性相关,就看是否存在一组不全为零的数 k1,k2,⋅⋅⋅,kmk_1, k_2, ···,k_mk1,k2,⋅⋅⋅,km使得上式成立。
定理1证明:
在式子k1a1+k2a2+...+kmam=0k_1a_1+k_2a_2+...+k_ma_m=0k1a1+k2a2+...+kmam=0的左右两边同时点乘a1a_1a1得
k1[a1,a1]+k2[a2,a1]+...+km[am,a1]=0k_1[a_1,a_1]+k_2[a_2,a_1]+...+k_m[a_m,a_1]=0k1[a1,a1]+k2[a2,a1]+...+km[am,a1]=0
由于a1,a2,…,ama_1,a_2,…,a_ma1,a2,…,am两两正交,因此:[a2,a1]=0,...[am,a1]=0[a_2,a_1]=0,...[a_m,a_1]=0[a2,a1]=0,...[am,a1]=0
k1[a1,a1]=0k_1[a_1,a_1]=0k1[a1,a1]=0,由条件可知a1a_1a1是非零向量,[a_1,a_1]≠0,
因此k1=0k_1=0k1=0,同理k2=0,⋅⋅⋅,km=0k_2=0, ···,k_m=0k2=0,⋅⋅⋅,km=0
a1,a2,…,ana_1,a_2,…,a_na1,a2,…,an线性无关.得证。
定义3
设n维向量e1,e2,…,ere_1,e_2,…,e_re1,e2,…,er,是向量空间V(V⊂Rn)V(V\subset R^n)V(V⊂Rn)的一个基,如果e1,e2,…,ere_1,e_2,…,e_re1,e2,…,er两两正交,且都是单位向量,则称e1,e2,…,ere_1,e_2,…,e_re1,e2,…,er是V的一个规范正交基。例如:
就是R4R^4R4的一个规范正交基.
若e1,e2,…,ere_1,e_2,…,e_re1,e2,…,er是V的一个规范正交基,那么V中任一向量a应能由e1,e2,…,ere_1,e_2,…,e_re1,e2,…,er线性表示,设表示式为
a=λ1e1+λ2e2+,…,+λrera=\lambda_1 e_1+\lambda_2e_2+,…,+\lambda_re_ra=λ1e1+λ2e2+,…,+λrer
λr=[a,er]\lambda_r=[a,e_r]λr=[a,er]
定义4
如果n阶矩阵A满足
ATA=E,即A−1=ATA^TA=E,即A^{-1}=A^TATA=E,即A−1=AT
那么称A为正交矩阵,简称正交阵。
上式用列向量表示,即是
[a1Ta2T⋮anT](a1,a2,⋯,an)=E\begin{bmatrix} a_1^T\\a_2^T \\ \vdots \\a_n^T \end{bmatrix}(a_1,a_2,\cdots,a_n)=E⎣⎢⎢⎢⎡a1Ta2T⋮anT⎦⎥⎥⎥⎤(a1,a2,⋯,an)=E
因为ATA=EA^TA=EATA=E与AAT=EAA^T=EAAT=E等价,所以上述结论对A的行向量也成立。
由此可见,你、阶正交阵A的n个列(或者行)向量构成的向量空间ℜn\real^nℜn的一个规范正交基。
判定矩阵A可逆的小结
1、A的行列式不等于0
2、A的秩等于A的维度n
3、a1,a2,…,ana_1,a_2,…,a_na1,a2,…,an线性无关
2.特征值特征向量以及矩阵的相似
方阵的特征值与特征向量
定义6
设A是n阶矩阵,如果数λ和μ维非零列向量x使下面关系式成立,
Ax=λxAx=\lambda xAx=λx
那么,这样的数λ称为矩阵A的特征值,非零向量x称为A的对应于特征值λ的特征向量。
人话版本(物理意义):刚开始讲矩阵的时候,讲过矩阵的本质是对应线性变换,如果从线性变换的角度看待这个问题,那么就是:现在我们有一个可以做线性变换的矩阵A,如果有一个向量x(注意不是变量),通过这个矩阵进行线性变换(就是乘上A)后的到x~\tilde xx~相对于原来的x方向不变,仅仅是大小变化而已(变大了λ倍),(说明这个x还蛮特殊的,一般的向量经过线性变换后大小方向都会变化)那么就把这个特殊的x叫做A的特征向量,变大的倍数λ称为特征值。
如果给我们一个A,如何来求特征值λ和特征向量x呢?就是把上面的公式Ax=λxAx=\lambda xAx=λx解方程,把x提取出来,x向量提取出来后,还剩下单位向量E,变成下面的公式:
(A−λE)x=0(A-\lambda E)x=0(A−λE)x=0
根据之前学过的克莱姆法则(如果Ax=0Ax=0Ax=0有非零解,则|A|=0,如果是|A|≠0则方程只有唯一解,那么x只能=0),则要使得上面的式子要有非零解的充分必要条件是∣A−λE∣=0|A-\lambda E|=0∣A−λE∣=0:
∣a11−λa12⋯a1na21a22−λ⋯a2n⋮⋮⋮an1an2⋯ann−λ∣=0\begin{vmatrix} a_{11}-\lambda& a_{12}&\cdots & a_{1n}\\ a_{21}& a_{22}-\lambda&\cdots & a_{2n}\\ \vdots& \vdots&& \vdots\\ a_{n1}& a_{n2}&\cdots & a_{nn}-\lambda \end{vmatrix}=0∣∣∣∣∣∣∣∣∣a11−λa21⋮an1a12a22−λ⋮an2⋯⋯⋯a1na2n⋮ann−λ∣∣∣∣∣∣∣∣∣=0
把上面的式子看做是关于λ的方程f(λ)=0f(\lambda)=0f(λ)=0
(i)λ1+λ2+⋯+λn=a11+a22+⋯+ann\lambda_1+\lambda_2+\cdots+\lambda_n=a_{11}+a_{22}+\cdots+a_{nn}λ1+λ2+⋯+λn=a11+a22+⋯+ann
(II)λ1λ2⋯λn=∣A∣\lambda_1\lambda_2\cdots\lambda_n=|A|λ1λ2⋯λn=∣A∣
设λ=λi\lambda=\lambda_iλ=λi为矩阵A的一个特征值,则由方程
(A−λiE)x=0(A-\lambda_iE)x=0(A−λiE)x=0
可求得非零解x=pix=p_ix=pi,那么pip_ipi便是A的对应于特征值λi\lambda_iλi的特征向量。
例子:求矩阵A=[3−1−13]A=\begin{bmatrix} 3 & -1\\ -1& 3 \end{bmatrix}A=[3−1−13]的特征值和特征向量。
解:先求∣A−λE∣=∣3−λ−1−13−λ∣=(3−λ)2−1=0|A-\lambda E|=\begin{vmatrix} 3-\lambda & -1\\ -1& 3-\lambda \end{vmatrix}=(3-\lambda)^2-1=0∣A−λE∣=∣∣∣∣3−λ−1−13−λ∣∣∣∣=(3−λ)2−1=0
3−λ=±13-\lambda=\pm 13−λ=±1求得两个特征值:λ1=2,λ2=4\lambda_1=2,\lambda_2=4λ1=2,λ2=4
分两步
第一步求λ1=2\lambda_1=2λ1=2对应的特征向量,解下面方程
(A−λ1E)x1=0(A-\lambda_1 E)x_1=0(A−λ1E)x1=0
[3−λ1−1−13−λ1]x1=0\begin{bmatrix} 3-\lambda_1 & -1\\ -1& 3-\lambda_1 \end{bmatrix}x_1=0[3−λ1−1−13−λ1]x1=0
[1−1−11][x11x12]=0\begin{bmatrix} 1 & -1\\ -1& 1 \end{bmatrix}\begin{bmatrix}x_{11}\\x_{12}\end{bmatrix}=0[1−1−11][x11x12]=0
解得:x1=[11]x_1=\begin{bmatrix}1\\1\end{bmatrix}x1=[11]归一化后得:x1=[2222]x_1=\begin{bmatrix}\cfrac{\sqrt{2}}{2}\\\cfrac{\sqrt{2}}{2}\end{bmatrix}x1=⎣⎢⎢⎡2222⎦⎥⎥⎤
第一步求λ1=4\lambda_1=4λ1=4对应的特征向量,解下面方程
(A−λ1E)x2=0(A-\lambda_1 E)x_2=0(A−λ1E)x2=0
[3−λ2−1−13−λ2]x2=0\begin{bmatrix} 3-\lambda_2 & -1\\ -1& 3-\lambda_2 \end{bmatrix}x_2=0[3−λ2−1−13−λ2]x2=0
[−1−1−1−1][x21x22]=0\begin{bmatrix} -1 & -1\\ -1& -1 \end{bmatrix}\begin{bmatrix}x_{21}\\x_{22}\end{bmatrix}=0[−1−1−1−1][x21x22]=0
解得:x2=[1−1]x_2=\begin{bmatrix}1\\-1\end{bmatrix}x2=[1−1]归一化后得:x2=[22−22]x_2=\begin{bmatrix}\cfrac{\sqrt{2}}{2}\\-\cfrac{\sqrt{2}}{2}\end{bmatrix}x2=⎣⎢⎢⎡22−22⎦⎥⎥⎤
再看一例:
求矩阵A=[−110−430102]A=\begin{bmatrix} -1& 1&0\\ -4& 3&0\\ 1 &0 &2 \end{bmatrix}A=⎣⎡−1−41130002⎦⎤的特征值和特征向量。
解:A的特征多项式为
∣A−λE∣=∣−1−λ10−43−λ0102−λ∣=(3−λ)2−1=(2−λ)(1−λ)2|A-\lambda E|=\begin{vmatrix} -1-\lambda & 1&0\\ -4& 3-\lambda&0\\ 1 &0&2-\lambda \end{vmatrix}=(3-\lambda)^2-1=(2-\lambda)(1-\lambda)^2∣A−λE∣=∣∣∣∣∣∣−1−λ−4113−λ0002−λ∣∣∣∣∣∣=(3−λ)2−1=(2−λ)(1−λ)2
所以A的特征值为λ1=2,λ2=λ3=1\lambda_1=2,\lambda_2=\lambda_3=1λ1=2,λ2=λ3=1
当λ1=2\lambda_1=2λ1=2时,解方程(A−2E)x=0(A-2E)x=0(A−2E)x=0.由
得基础解系:p1=[001]p_1=\begin{bmatrix}0\\0\\1\end{bmatrix}p1=⎣⎡001⎦⎤
所以kp1(k≠0)kp_1(k\neq0)kp1(k=0)是对应于λ1=2\lambda_1=2λ1=2的全部特征向量。
另外一组解:
当λ2=λ3=1\lambda_2=\lambda_3=1λ2=λ3=1时,解方程(A−2E)x=0(A-2E)x=0(A−2E)x=0.由
得基础解系:p2=[−1−21]p_2=\begin{bmatrix}-1\\-2\\1\end{bmatrix}p2=⎣⎡−1−21⎦⎤
所以kp2(k≠0)kp_2(k\neq0)kp2(k=0)是对应于λ2=λ3=1\lambda_2=\lambda_3=1λ2=λ3=1的全部特征向量。
由于有重根,所以只要两个特征向量
例8设λ\lambdaλ是方阵A的特征值,证明
(1)λ2\lambda^2λ2是A2A^2A2的特征值;
(2)当A可逆时,1λ\frac{1}{\lambda}λ1是A−1A^{-1}A−1的特征值.
证明(1):由λ\lambdaλ是方阵A的特征值可知:Ax=λxAx=\lambda xAx=λx
A2x=λAx=λ2xA^2x=\lambda Ax=\lambda^2xA2x=λAx=λ2x
以此类推:AnA^nAn的特征值为λn\lambda^nλn,特征向量为x
n可以为负数,例如A−2A^{-2}A−2的特征值为λ−2\lambda^{-2}λ−2
证明(2):由λ\lambdaλ是方阵A的特征值可知:Ax=λxAx=\lambda xAx=λx,两边同时乘以A的逆矩阵得:
x=λA−1xx=\lambda A^{-1}xx=λA−1x,两边同时除以λ\lambdaλ得
1λx=A−1x\frac{1}{\lambda}x=A^{-1}xλ1x=A−1x,即A−1x=1λxA^{-1}x=\frac{1}{\lambda}xA−1x=λ1x,根据特征值的定义可知:
1λ\frac{1}{\lambda}λ1是A−1A^{-1}A−1的特征值,特征向量为x
再推广:如果λ\lambdaλ是方阵A的特征值,那么f(λ)f(\lambda)f(λ)是方阵f(A)f(A)f(A)的特征值。
例子:设3阶矩阵A的特征值为1,-1,2,求A2+3A−2EA^2+3A-2EA2+3A−2E的特征值。
解:把A的特征值1,-1,2分别代入上式
12+3∗1−2=21^2+3*1-2=212+3∗1−2=2
(−1)2+3(−1)−2=−4(-1)^2+3(-1)-2=-4(−1)2+3(−1)−2=−4
22+3∗2−2=82^2+3*2-2=822+3∗2−2=8
A2+3A−2EA^2+3A-2EA2+3A−2E的特征值为2,-4,8
线代:1.6矩阵的特征值和特征向量相关推荐
- 【线代NumPy】第八章 - 特征值和特征向量 | Eigenvalue and Eigenvector | 简述并提供代码
- 如何用计算机求特征值特征向量,利用QR算法求解矩阵的特征值和特征向量
利用QR算法求解矩阵的特征值和特征向量 为了求解一般矩阵(不是那种幼稚到shi的2 x 2矩阵)的特征值. 根据定义的话,很可能需要求解高阶方程... 这明显是个坑...高阶方程你肿么破... 折腾了 ...
- 利用QR算法求解矩阵的特征值和特征向量
利用QR算法求解矩阵的特征值和特征向量 为了求解一般矩阵(不是那种幼稚到shi的2 x 2矩阵)的特征值. 根据定义的话,很可能需要求解高阶方程... 这明显是个坑...高阶方程你肿么破... 折腾了 ...
- 特征值与特征向量_矩阵的特征值和特征向量
不少学习过线性代数的同学可能都有这样的疑惑,就是线性代数到底是什么?我们算的这些东西究竟有什么用?回忆起这门课来可能仅有的印象也就是矩阵.向量.还有一个特征什么来着? 线性代数是一门相对较为年轻的学科 ...
- 矩阵的特征值和特征向量的雅克比算法C/C++实现
矩阵的特征值和特征向量是线性代数以及矩阵论中非常重要的一个概念.在遥感领域也是经常用到,比如多光谱以及高光谱图像的主成分分析要求解波段间协方差矩阵或者相关系数矩阵的特征值和特征向量. 根据普通线性代数 ...
- 矩阵的特征值、特征向量
这部分包括:正交矩阵.矩阵的特征值.特征向量.相似矩阵.实对称矩阵对角化.
- numpy求解矩阵的特征值和特征向量
python2.7代码如下: #-*- encoding:utf-8 -*- import sys reload(sys) sys.setdefaultencoding('utf-8') import ...
- 雅可比旋转求解对称二维矩阵的特征值和特征向量
问题描述: 给定一个矩阵,如下: A=[a11a21a12a22] A=\begin{bmatrix} a_{11}&a_{12}\\ a_{21}& a_{22} \end{bmat ...
- python中向量长度_线性代数精华——矩阵的特征值与特征向量
点击上方蓝字,和我一起学技术. 今天和大家聊一个非常重要,在机器学习领域也广泛使用的一个概念--矩阵的特征值与特征向量. 我们先来看它的定义,定义本身很简单,假设我们有一个n阶的矩阵A以及一个实数λ, ...
最新文章
- 柳叶刀发布陈薇团队新冠疫苗试验结果:安全,能诱导免疫反应
- 在组策略中用户策略仅对特定计算机生效,如何对本地组策略设置使之不对特定用户生效?...
- Oracle ebs(E-Business Suite) 电子商务套件 简介
- powerdesigner连接db2生成模型步骤
- MySQL索引背后的数据结构及算法原理-转
- 设计模式(五)行为型模式
- 解决Eclipse 项目报错:Unbound classpath container: ‘JRE System Library [JavaSE-1.7]
- 实训作业 4(界面2)
- 仿英雄联盟网页HTML代码 学生网页设计与制作期末作业下载 大学生网页设计与制作成品下载 DW游戏介绍网页作业代码下载
- 泛微协同办公系统移动服务器,泛微协同办公平台Ecology系统重装迁移指导手册.pdf...
- 计算机无法删除tf卡的内容,SD卡无法格式化数据可以读取,但是删除不掉
- 软件开发人员是一种很棒的职业选择的五大理由
- Android5.0,IDA远程调试 The debugger could not attach to the selected process. irs_recv 等待的操作过时
- 最简明扼要的 Systemd 教程,只需十分钟
- 基于WEB多媒体电子贺卡平台
- 2023年湖北武汉建设厅特种工架子工证怎么报名?报考条件?启程别
- 不同大小硬盘对拷oracle,大小不一样的硬盘怎么实现对拷?
- 以太坊解析之二——POA共识过程与一些可能的修改方案
- 2013年中国android智能手机用户调查研究报告,ZDC:2013年7月中国智能手机市场分析报告...
- 正点原子STM32F429核心板的插座型号