【连续介质力学】特征值和特征向量问题

news2024/11/24 12:39:32

特征值和特征向量问题

二阶张量和一个向量(单位向量 n ^ ′ \hat n' n^)的点积会得到一个向量,也就是说,将一个二阶张量投影到某个方向所得到的向量的方向实际上与 n ^ ′ \hat n' n^ 的方向不一样

特征值和特征向量问题的目标是找到一个方向 n ^ \hat n n^,与投影向量 t ⃗ n ^ = T ⋅ n ^ \vec t^{\hat n} = T \cdot \hat n t n^=Tn^的方向一致

在这里插入图片描述
n ^ \hat n n^是张量 T T T特征向量,如果存在一个标量 λ \lambda λ,即特征值,使得:
T ⋅ n ^ = λ n ^ T \cdot \hat n=\lambda \hat n Tn^=λn^

指标形式:
在这里插入图片描述

以上齐次方程只有非平凡解,例如: n ^ ≠ 0 \hat n \neq 0 n^=0 , 当且仅当:
det ⁡ ( T − λ 1 ) = 0 ; ∣ T i j − λ δ i j ∣ = 0 \det (T - \lambda 1)=0; \quad |T_{ij} - \lambda \delta_{ij}| = 0 det(Tλ1)=0;Tijλδij=0
以上称为张量的特征行列式,显式表示为:
在这里插入图片描述
展开行列式,得到特征多项式:
λ 3 − λ 2 I T + λ I I T − I I I T = 0 \lambda^3 - \lambda^2I_T +\lambda II_T - III_T=0 λ3λ2IT+λIITIIIT=0
其中, I t , I I T , I I I T I_t, II_T, III_T It,IIT,IIIT是张量T的不变量:
在这里插入图片描述
其中 M i i M_{ii} Mii是余子式矩阵的迹: M i i = M 11 + M 22 + M 33 M_{ii} = M_{11} + M_{22} + M_{33} Mii=M11+M22+M33
显式表示:

在这里插入图片描述

如果 T T T 是对称张量,主不变量如下:
在这里插入图片描述
通过求解特征多项式可以得到特征值,一旦得到特征值,可以应用以下等式求出特征向量:
( T i j − λ 1 δ i j ) n ^ j ( 1 ) = 0 i (T_{ij}-\lambda_1 \delta_{ij})\hat n_j^{(1)}=0_i (Tijλ1δij)n^j(1)=0i
( T i j − λ 2 δ i j ) n ^ j ( 2 ) = 0 i (T_{ij}-\lambda_2 \delta_{ij})\hat n_j^{(2)}=0_i (Tijλ2δij)n^j(2)=0i
( T i j − λ 3 δ i j ) n ^ j ( 3 ) = 0 i (T_{ij}-\lambda_3 \delta_{ij})\hat n_j^{(3)}=0_i (Tijλ3δij)n^j(3)=0i
这些特征向量构成了一个新空间,称为主空间

如果 T T T是一个对称张量,那么主空间由正交基定义且所有的特征值都是实数

如果三个特征值都互不相同 λ 1 ≠ λ 2 ≠ λ 3 \lambda_1 \neq \lambda_2 \neq \lambda_3 λ1=λ2=λ3,那么主方向是唯一的

如果其中有两个相等, λ 1 = λ 2 ≠ λ 3 \lambda_1 = \lambda_2 \neq \lambda_3 λ1=λ2=λ3,可以说明特征值 λ 3 \lambda_3 λ3的主方向 n ^ ( 3 ) \hat n^{(3)} n^(3)是唯一的,并且任何与 n ^ ( 3 ) \hat n^{(3)} n^(3)垂直的平面是主平面,正交性是确定 n ^ ( 1 ) \hat n^{(1)} n^(1) n ^ ( 2 ) \hat n^{(2)} n^(2)的唯一约束

如果 λ 1 = λ 2 = λ 3 \lambda_1 = \lambda_2 =\lambda_3 λ1=λ2=λ3,那么任何方向都是主方向

如果一个张量有3个相等特征值,那么这个张量称为球形张量

在主空间的张量分量只由主分量构成
在这里插入图片描述
所以,主不变量可以表示为:
I T = T 1 + T 2 + T 3 I_T = T_1 + T_2 + T_3 IT=T1+T2+T3
I I T = T 1 T 2 + T 2 T 3 + T 1 T 3 II_T = T_1T_2 + T_2T_3 + T_1T_3 IIT=T1T2+T2T3+T1T3
I I I T = T 1 T 2 T 3 III_T = T_1 T_2 T_3 IIIT=T1T2T3

如果 T T T是球形张量,有 T 1 = T 2 = T 3 = T , T_1=T_2 =T_3=T, T1=T2=T3=T那么主不变量为
I T 2 = 3 I I T I_T^2 = 3II_T IT2=3IIT
I I I T = T 3 III_T=T^3 IIIT=T3

如果 W W W 是一个反对称张量, W W W 的主不变量为

在这里插入图片描述
其中, w 2 = ∣ ∣ w ⃗ ∣ ∣ 2 = w ⃗ ⋅ w ⃗ = W 23 2 + W 13 2 + W 12 2 w^2 = ||\vec w||^2=\vec w \cdot \vec w = W_{23}^2+W_{13}^2+W_{12}^2 w2=∣∣w 2=w w =W232+W132+W122,然后,反对称张量的特征方程变为:
λ 3 − λ 2 I W + λ I I W − I I I W = 0    ⟹    λ 3 + w 2 λ = 0    ⟹    λ ( λ 2 + w 2 ) = 0 \lambda^3 - \lambda^2I_W+\lambda II_W-III_W=0 \implies \lambda^3 + w^2\lambda = 0 \implies \lambda(\lambda^2+w^2)=0 λ3λ2IW+λIIWIIIW=0λ3+w2λ=0λ(λ2+w2)=0

在这种情况,有一个特征值是实数0,而其他的特征值是虚根:
λ 2 + w 2 = 0    ⟹    λ 2 = − w 2 = 0    ⟹    λ ( 1 , 2 ) = ± w − 1 = ± w i \lambda^2 + w^2=0 \implies \lambda^2=-w^2 = 0 \implies \lambda_{(1, 2)}= \pm w\sqrt{-1}=\pm wi λ2+w2=0λ2=w2=0λ(1,2)=±w1 =±wi

特征向量的正交性

λ 1 , λ 2 , λ 3 \lambda_1, \lambda_2, \lambda_3 λ1,λ2,λ3是T的特征值,那么:
T ⋅ n ^ ( 1 ) = λ 1 n ^ ( 1 ) T\cdot \hat n^{(1)} = \lambda_1\hat n^{(1)} Tn^(1)=λ1n^(1)
T ⋅ n ^ ( 2 ) = λ 2 n ^ ( 2 ) T\cdot \hat n^{(2)} = \lambda_2\hat n^{(2)} Tn^(2)=λ2n^(2)
T ⋅ n ^ ( 3 ) = λ 3 n ^ ( 3 ) T\cdot \hat n^{(3)} = \lambda_3\hat n^{(3)} Tn^(3)=λ3n^(3)

n ^ ( 2 ) \hat n^{(2)} n^(2) T ⋅ n ^ ( 1 ) = λ 1 n ^ ( 1 ) T\cdot \hat n^{(1)} = \lambda_1\hat n^{(1)} Tn^(1)=λ1n^(1) 进行点积, 并且将 n ^ ( 1 ) \hat n^{(1)} n^(1) T ⋅ n ^ ( 2 ) = λ 2 n ^ ( 2 ) T\cdot \hat n^{(2)} = \lambda_2\hat n^{(2)} Tn^(2)=λ2n^(2) 进行点积:

n ^ 2 ⋅ T ⋅ n ^ ( 1 ) = λ 1 n ^ ( 2 ) ⋅ n ^ ( 1 ) n ^ 1 ⋅ T ⋅ n ^ ( 2 ) = λ 2 n ^ ( 1 ) ⋅ n ^ ( 2 ) \hat n^{2} \cdot T \cdot \hat n^{(1)} = \lambda_1 \hat n^{(2)} \cdot \hat n^{(1)} \\ \hat n^{1} \cdot T \cdot \hat n^{(2)} = \lambda_2 \hat n^{(1)} \cdot \hat n^{(2)} n^2Tn^(1)=λ1n^(2)n^(1)n^1Tn^(2)=λ2n^(1)n^(2)

由于 T T T 是对称的,所以 n ^ 2 ⋅ T ⋅ n ^ ( 1 ) = n ^ 1 ⋅ T ⋅ n ^ ( 2 ) \hat n^{2} \cdot T \cdot \hat n^{(1)} = \hat n^{1} \cdot T \cdot \hat n^{(2)} n^2Tn^(1)=n^1Tn^(2), 所以:
λ 1 n ^ ( 2 ) ⋅ n ^ ( 1 ) = λ 2 n ^ ( 1 ) ⋅ n ^ ( 2 ) = λ 2 n ^ ( 2 ) ⋅ n ^ ( 1 )    ⟹    ( λ 1 − λ 2 ) n ^ ( 1 ) ⋅ n ^ ( 2 ) = 0 \lambda_1 \hat n^{(2)} \cdot \hat n^{(1)} = \lambda_2 \hat n^{(1)} \cdot \hat n^{(2)} = \lambda_2 \hat n^{(2)} \cdot \hat n^{(1)} \\ \implies (\lambda_1 - \lambda_2)\hat n^{(1)} \cdot \hat n^{(2)} = 0 λ1n^(2)n^(1)=λ2n^(1)n^(2)=λ2n^(2)n^(1)(λ1λ2)n^(1)n^(2)=0
为了在 λ 1 ≠ λ 2 ≠ 0 \lambda_1 \neq \lambda_2 \neq 0 λ1=λ2=0情况下使上面等式成立,需要:
n ^ ( 1 ) ⋅ n ^ ( 2 ) = 0 \hat n^{(1)} \cdot \hat n^{(2)} =0 n^(1)n^(2)=0

同理,可得: n ^ ( 2 ) ⋅ n ^ ( 3 ) = 0 n ^ ( 1 ) ⋅ n ^ ( 3 ) = 0 \hat n^{(2)} \cdot \hat n^{(3)} =0\quad \hat n^{(1)} \cdot \hat n^{(3)} =0 n^(2)n^(3)=0n^(1)n^(3)=0,则可以下结论特征向量是正交的,且构成正交基,其中坐标系之间的变换矩阵是:
在这里插入图片描述

对角化
在这里插入图片描述

问题1.31 证明以下为不变量: C 1 2 + C 2 2 + C 3 2 C 1 3 + C 2 3 + C 3 3 C 1 4 + C 2 4 + C 3 4 C_1^2+C_2^2+C_3^2 \quad C_1^3+C_2^3+C_3^3 \quad C_1^4+C_2^4+C_3^4 C12+C22+C32C13+C23+C33C14+C24+C34,其中 C 1 , C 2 , C 3 C_1, C_2, C_3 C1,C2,C3是二阶张量 C C C 的特征值

在这里插入图片描述

问题1.32 Q是正交张量,E是任意二阶张量,证明:E的特征值不会影响以下正交变换: E ∗ = Q ⋅ E ⋅ Q T E^*=Q\cdot E \cdot Q^T E=QEQT

在这里插入图片描述
结论:任意张量经过正交变换不改变特征值

三次方程的解

T T T是对称二阶张量,其特征方程 λ 3 − λ 2 I T + λ I I T − I I I T = 0 \lambda^3 - \lambda^2I_T +\lambda II_T - III_T = 0 λ3λ2IT+λIITIIIT=0的根是都是实数,表示成:
在这里插入图片描述
其中:
在这里插入图片描述
矩阵形式:
在这里插入图片描述
其中,我们清楚地将张量在主空间分解为球部分和偏部分
注意到,当 T T T 是一个球形张量,满足: I T 2 = 3 I I T I_T^2 =3II_T IT2=3IIT, 因此 S = 0 S = 0 S=0

问题1.33 求T的主值和主方向

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
总结:张量的特征向量构成原坐标系 ( x 1 , x 2 , x 3 ) (x_1, x_2, x_3) (x1,x2,x3)到主空间 ( x 1 ′ , x 2 , ′ x 3 ′ ) (x_1', x_2,' x_3') (x1,x2,x3)的变换矩阵 A A A

参考教材:
Eduardo W.V. Chaves, Notes On Continuum Mechanics

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/526574.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

IDEA添加.gitignore忽略不需要提交的文件

问题 git上传的时候,我们已经将 xxx 文件添加到了.gitignore 中,但是在push 后,远程仓库还是会显示此文件,比如我们在.gitignore文件当中添加了不需要提交的target目录,但是提交的时候,还是会被提交。 原因…

2023.5.14总结

这周平时在刷蓝桥杯的题目,周天打了一场2021年陕西省的省赛的重现赛。 重现赛我们没打满,打了三个小时,A了四个,不过应该也差不多了。 登录—专业IT笔试面试备考平台_牛客网 以前没学过数论分块,今天学了学。对于一些向…

Hadoop上传及下载数据流程

网络拓扑及机架感知 网络拓扑 节点距离:两个节点到达共同父节点的距离和 机架感知 ( 副本节点的选择 ) 例如:500个节点,上传数据my.tar.gz,副本数为3, 根据机架感知,副本数据…

防止攻击者对您使用合法工具

恶意行为者越来越多地利用合法工具来实现其目标,其中包括禁用安全措施、横向移动和传输文件。使用常用工具可以让攻击者逃避检测。 虽然端点产品可以将定制工具或恶意软件标记为恶意软件,但商业上可用的工具通常被组织标记为干净或列入允许列表。 这让…

MacBook Pro合上盖子不休眠的问题简单分析

15年款的MacBook Pro每次不用的时候都是直接合上盖子(开着一堆程序)系统会自动休眠,但是升级了新系统Sierra之后就发现合上盖子竟然没有休眠(第二次打开盖子后发现掉了50%多的电,而且温度比较高)&#xff0…

软考A计划-真题-分类精讲汇总-第十一章(多媒体基础)

点击跳转专栏>Unity3D特效百例点击跳转专栏>案例项目实战源码点击跳转专栏>游戏脚本-辅助自动化点击跳转专栏>Android控件全解手册点击跳转专栏>Scratch编程案例 👉关于作者 专注于Android/Unity和各种游戏开发技巧,以及各种资源分享&am…

MySQL只有information_schema一个数据库

背景:centos新安装的mysql数据库。使用DbEaver连接mysql库时,发现左边不显示表。使用命令框mysql -uroot回车登录时,发现只能看到information_schema一个数据库了。 原因:   因为mysql数据库的user表里,存在用户名为…

Julia入门-1、使用C++调用Julia脚本语言

文章目录 0、开发环境1、测试Julia环境2、调用Julia脚本语言准备3、使用C++调用Julia脚本语言(1)使用C++调用简单的Julia脚本语言(2)使用C++调用复杂的Julia脚本语言0、开发环境 操作系统: ①Windows 10 开发编译器: ①VS 2015 Professional ②VS Code + julia-vscode插件(…

《计算机网络—自顶向下方法》 第六章Wireshark实验:IP 协议分析

IP 协议(Internet Protocol),又译为网际协议或互联网协议,是用在 TCP/IP 协议簇中的网络层协议。主要功能是无连接数据报传送、数据报路由选择和差错控制。IP 协议是 TCP/IP 协议族的核心协议,其主要包含两个方面&…

Spring Boot 配置文件总结

前言 Spring Boot 中提供一个全局的配置文件:application.properties,这个配置文件的作用就是,允许我们通过这个配置文件去修改 Spring Boot 自动配置的默认值。 Spring Boot 支持两种格式的配置文件:application.properties 和…

C/C++每日一练(20230515) 区间和的个数、BST最近公共祖先、最接近元素

目录 1. 区间和的个数 🌟🌟🌟 2. 二叉搜索树的最近公共祖先 🌟 3. 找最接近元素 🌟🌟 🌟 每日一练刷题专栏 🌟 Golang每日一练 专栏 Python每日一练 专栏 C/C每日一练 专栏…

前端路由、vue-router常见用法、路由重定向、动态路由匹配、声明式导航 编程式导航 、导航守卫

前端路由、vue-router常见用法、路由重定向、动态路由匹配、声明式导航 & 编程式导航 、导航守卫 前端路由的概念与原理前端路由 vue-router 的基本使用vue-router 的常见用法路由重定向动态路由匹配声明式导航 & 编程式导航导航守卫 后台管理案例 前端路由的概念与原理…

mysqld之mha高可用

1.MHA的相关知识 1.1 什么是 MHA MHA(MasterHigh Availability)是一套优秀的MySQL高可用环境下故障切换和主从复制的软件。 MHA 的出现就是解决MySQL 单点故障的问题。 MySQL故障切换过程中,MHA能做到0-30秒内自动完成故障切换操作。 MHA能在…

C语言设计三子棋

引入 谈到三子棋,大家应该都不陌生,学生时代我们大多人都爱拿作文本有事没事就跟同桌下两把,只要任意一方三点连成一线,就可以胜利。今天我作为一个计算机方面的博主,将会用C语言实现这个简单的小游戏(人机…

sort命令 uniq命令 tr命令 cut命令

sort命令 ——以行为单位对文件内容进行排序,也可以根据不同的数据类型来排序 比较原则是从首字符向后,依次按ASCII码值进行比较,最后将他们按升序输出 语法格式: sort [选项] 参数 cat file | sort 选项 -n按照数字进行排序…

MySQL高可用之MHA集群

一、MHA概述 1.1 什么是 MHA MHA(MasterHigh Availability)是一套优秀的MySQL高可用环境下故障切换和主从复制的软件。 MHA 的出现就是解决MySQL 单点故障的问题。 MySQL故障切换过程中,MHA能做到0-30秒内自动完成故障切换操作。 MHA能在…

【模式识别9】python计算目标检测IoU、TP、FP、FN、Precision、Recall指标

python计算目标检测IoU、TP、FP、FN、Precision、Recall指标 1. 基础概念1.1 TP、TN、FP、FN1.2 IoU1.3 Precision(P)、Recall(R)、F1-score 2. python代码3. 总结 代码资源:IoU_P_R.py 1. 基础概念 1.1 TP、TN、FP、…

2023/5/14总结

哈夫曼树 哈夫曼树:给定n个权值作为n个叶子结点,构造一棵二叉树,若该树的带权路径长度(WPL)达到最小,则称该二叉树为哈夫曼树,也被称为最优二叉树。 怎样才能使带权路径长度最短:根据…

CTF-PHP反序列化漏洞5-反序列化字符逃逸

作者:Eason_LYC 悲观者预言失败,十言九中。 乐观者创造奇迹,一次即可。 一个人的价值,在于他所拥有的。可以不学无术,但不能一无所有! 技术领域:WEB安全、网络攻防 关注WEB安全、网络攻防。我的…

简单聊聊微前端

简单聊聊微前端 介绍微前端的优点应用间相互独立,互不依赖可以同时使用不同的技术栈可拓展性高可维护性更强,减少代码量提高开发和部署的效率团队的高度自主权错误隔离 微前端的缺点依赖项冗余CSS样式冲突和重叠性能比较差应用间的通信不够便捷 实现微前…