核—幂零分解
若向量空间 V \mathcal V V存在子空间 X \mathcal X X与 Y \mathcal Y Y,当
X + Y = V X ∩ Y = 0 \mathcal {X\text{+}Y\text{=}V}\\ \mathcal {X}\cap \mathcal {Y}=0 X+Y=VX∩Y=0
时称子空间 X \mathcal X X与 Y \mathcal Y Y是完备的,其中记为 X ⊕ Y = V \mathcal X \oplus \mathcal Y = \mathcal V X⊕Y=V
若存在 X ⊕ Y = V \mathcal X \oplus \mathcal Y = \mathcal V X⊕Y=V, 对于 x ∈ X , y ∈ Y , v ∈ V x\in \mathcal X,y \in \mathcal Y,v \in \mathcal V x∈X,y∈Y,v∈V,满足 v = x + y v=x+y v=x+y,则向量 x x x被称为向量 v v v沿着 Y \mathcal Y Y到 X \mathcal X X 空间的投影,向量 y y y被称为向量 v v v沿着 X \mathcal X X到 Y \mathcal Y Y 空间的投影,若存在 P v = x Pv=x Pv=x, P P P被称为沿着 Y \mathcal Y Y到 X \mathcal X X 空间的投影算子,其中
- P 2 = P P^2=P P2=P
- 1 − P 1-P 1−P沿着 X \mathcal X X到 Y \mathcal Y Y 空间的投影算子
- R ( P ) = N ( 1 − P ) = X R(P)=N(1-P)=\mathcal X R(P)=N(1−P)=X
- N ( P ) = R ( 1 − P ) = Y N(P)=R(1-P)=\mathcal Y N(P)=R(1−P)=Y
若 V = R n V=\mathfrak R^n V=Rn,则 P [ X ∣ Y ] = [ X ∣ 0 ] P[\mathbf X|\mathbf Y]=[\mathbf X|\mathbf 0] P[X∣Y]=[X∣0],即 P = [ X ∣ 0 ] [ X ∣ Y ] − 1 = [ X ∣ 0 ] ( I 0 0 0 ) [ X ∣ Y ] − 1 P=[\mathbf X|\mathbf 0][\mathbf X|\mathbf Y]^{-1}=[\mathbf X|\mathbf 0]\begin{pmatrix}\mathbf I&\mathbf0\\\mathbf 0&\mathbf 0\end{pmatrix}[\mathbf X|\mathbf Y]^{-1} P=[X∣0][X∣Y]−1=[X∣0](I000)[X∣Y]−1,其中 X , Y \mathbf X,\mathbf Y X,Y分别表示 X , Y \mathcal X ,\mathcal Y X,Y的一组基
值域零空间分解
若存在一个k,满足 rank ( A k ) = rank ( A k + 1 ) \text{rank}(A^k)=\text{rank}(A^{k+1}) rank(Ak)=rank(Ak+1),则将最小的那个k值称为index,其中非奇异矩阵的index为0
对于奇异矩阵 A n × n A_{n\times n} An×n,存在一个index k,使得$R(A^k)\oplus N(A^k)=\mathfrak R^n $
若存在一个矩阵 A k = 0 A^k=0 Ak=0,其中index(A)=0,则矩阵A被称为幂零矩阵
核—幂零分解
如果A是一个 n × n n\times n n×n 的index为k的奇异矩阵,其中 rank ( A k ) = r \text{rank}(A^k)=r rank(Ak)=r,则存在一个非奇异矩阵 Q Q Q, 满足
Q − 1 A Q = ( C r × r 0 0 N ) \left.\mathbf{Q}^{-1}\mathbf{A}\mathbf{Q}=\left(\begin{array}{cc}\mathbf{C}_{r\times r}&\mathbf{0}\\\mathbf{0}&\mathbf{N}\end{array}\right.\right) Q−1AQ=(Cr×r00N)
其中 C C C是非奇异矩阵, N N N是index为k的幂零矩阵,其中 Q Q Q为矩阵 A k A^k Ak的值域空间和零空间的基的组合
若存在 A = Q ( C 0 0 N ) Q − 1 \left.\mathbf{A}=\mathbf{Q}\left(\begin{array}{ll}\mathbf{C}&0\\0&\mathbf{N}\end{array}\right.\right)\mathbf{Q}^{-1} A=Q(C00N)Q−1,则 A D = Q ( C − 1 0 0 0 ) Q − 1 \left.\mathbf{A}^D=\mathbf{Q}\left(\begin{array}{ll}\mathbf{C}^{-1}&0\\0&0\end{array}\right.\right)\mathbf{Q}^{-1} AD=Q(C−1000)Q−1,其中 A D A^D AD被称为A的广义逆
对于矩阵 A = ( − 2 0 − 4 4 2 4 3 2 2 ) \left.\textbf{A}=\left(\begin{array}{rrr}-2&0&-4\\4&2&4\\3&2&2\end{array}\right.\right) A= −243022−442 ,计算出 core-nilpoten 的分解形式,并给出对应的 Drazin 逆的形式。
直接计算可得 : r a n k ( A ) = 2 , r a n k ( A 2 ) = 1 , r a n k ( A 3 ) = 1 :\:rank(\mathbf{A})=2,\:rank(\mathbf{A^2})=1,\:rank(\mathbf{A^3})=1 :rank(A)=2,rank(A2)=1,rank(A3)=1, 由此可知 : i n d e x ( A ) = 2. :index(\mathbf{A})=2. :index(A)=2. 由 core-nilpotent 分解可知,矩阵 Q = [ X ∣ Y ] \mathbf{Q}=[\mathbf{X}|\mathbf{Y}] Q=[X∣Y], 这里 X \mathbf{X} X 和 Y 分别为 R ( A 2 ) R(\mathbf{A}^2) R(A2) 和 N ( A 2 ) N(\mathbf{A}^2) N(A2) 的一组基。从而直接计算可得,
X = ( − 8 12 8 ) , Y = ( − 1 0 1 0 0 1 ) , \left.\mathbf{X}=\left(\begin{array}{rr}-8\\12\\8\end{array}\right.\right),\quad\mathbf{Y}=\left(\begin{array}{rr}-1&0\\1&0\\0&1\end{array}\right), X= −8128 ,Y= −110001 ,
可得
Q = ( − 8 − 1 0 12 1 0 8 0 1 ) \left.\mathbf{Q}=\left(\begin{array}{rrr}-8&-1&0\\12&1&0\\8&0&1\end{array}\right.\right) Q= −8128−110001
所以
Q − 1 A Q = ( 1 4 1 4 0 − 3 − 2 0 − 2 − 2 1 ) ( − 2 0 − 4 4 2 4 3 2 2 ) ( − 8 − 1 0 12 1 0 8 0 1 ) = ( 2 0 0 0 − 2 4 0 − 1 2 ) \left.\mathbf{Q}^{-1}\mathbf{A}\mathbf{Q}=\left(\begin{array}{rrr}\frac{1}{4}&\frac{1}{4}&0\\-3&-2&0\\-2&-2&1\end{array}\right.\right)\left(\begin{array}{rrr}-2&0&-4\\4&2&4\\3&2&2\end{array}\right)\left(\begin{array}{rrrr}-8&-1&0\\12&1&0\\8&0&1\end{array}\right)=\left(\begin{array}{rrr}2&0&0\\0&-2&4\\0&-1&2\end{array}\right) Q−1AQ= 41−3−241−2−2001 −243022−442 −8128−110001 = 2000−2−1042
因为 Q − 1 A Q = ( C 0 0 N ) , C = ( 2 ) , N = ( − 2 4 − 1 2 ) \left.\mathbf{Q}^{-1}\mathbf{AQ}=\left(\begin{array}{ll}\mathbf{C}&\mathbf{0}\\\mathbf{0}&\mathbf{N}\end{array}\right.\right),\left.\mathbf{C}=(2),\mathbf{N}=\left(\begin{array}{cc}-2&4\\-1&2\end{array}\right.\right) Q−1AQ=(C00N),C=(2),N=(−2−142)
所以
相关文章:
核—幂零分解
若向量空间 V \mathcal V V存在子空间 X \mathcal X X与 Y \mathcal Y Y,当 X Y V X ∩ Y 0 \mathcal {X\text{}Y\text{}V}\\ \mathcal {X}\cap \mathcal {Y}0 XYVX∩Y0 时称子空间 X \mathcal X X与 Y \mathcal Y Y是完备的,其中记为 X ⊕ Y V \ma…...
轻松掌控财务,分析账户花销,明细记录支出情况
随着科技的发展,我们的生活变得越来越智能化。然而,对于许多忙碌的现代人来说,管理财务可能是一件令人头疼的事情。复杂的账单、花销、收入,这些可能会让你感到无从下手。但现在,我们有一个全新的解决方案——一款全新…...
竞赛 题目:基于机器视觉opencv的手势检测 手势识别 算法 - 深度学习 卷积神经网络 opencv python
文章目录 1 简介2 传统机器视觉的手势检测2.1 轮廓检测法2.2 算法结果2.3 整体代码实现2.3.1 算法流程 3 深度学习方法做手势识别3.1 经典的卷积神经网络3.2 YOLO系列3.3 SSD3.4 实现步骤3.4.1 数据集3.4.2 图像预处理3.4.3 构建卷积神经网络结构3.4.4 实验训练过程及结果 3.5 …...
11. Spring源码篇之实例化前的后置处理器
简介 spring在创建Bean的过程中,提供了很多个生命周期,实例化前就是比较早的一个生命周期,顾名思义就是在Bean被实例化之前的处理,这个时候还没实例化,只能拿到该Bean的Class对象,如果在这个时候直接返回一…...
Python-Python高阶技巧:HTTP协议、静态Web服务器程序开发、循环接收客户端的连接请求
版本说明 当前版本号[20231114]。 版本修改说明20231114初版 目录 文章目录 版本说明目录HTTP协议1、网址1.1 网址的概念1.2 URL的组成1.3 知识要点 2、HTTP协议的介绍2.1 HTTP协议的概念及作用2.2 HTTP协议的概念及作用2.3 浏览器访问Web服务器的过程 3、HTTP请求报文3.1 H…...
P1304 哥德巴赫猜想
题目描述 输入一个偶数 N,验证 4∼N 所有偶数是否符合哥德巴赫猜想:任一大于 22 的偶数都可写成两个质数之和。如果一个数不止一种分法,则输出第一个加数相比其他分法最小的方案。例如 1010,10=3+7=5+510=3+7=5+5,则 10=5+510=5+5 是错误答案。 输入格式 第一行输入一个…...
CSDN每日一题学习训练——Python版(搜索插入位置、最大子序和)
版本说明 当前版本号[20231118]。 版本修改说明20231118初版 目录 文章目录 版本说明目录搜索插入位置题目解题思路代码思路参考代码 最大子序和题目解题思路代码思路参考代码 搜索插入位置 题目 给定一个排序数组和一个目标值,在数组中找到目标值,…...
Java在物联网中的重要性
【点我-这里送书】 本人详解 作者:王文峰,参加过 CSDN 2020年度博客之星,《Java王大师王天师》 公众号:JAVA开发王大师,专注于天道酬勤的 Java 开发问题中国国学、传统文化和代码爱好者的程序人生,期待你的关注和支持!本人外号:神秘小峯 山峯 转载说明:务必注明来源(…...
动态规划解背包问题
题目 题解 def knapsac(W: int, N: int, wt: List[int], val: List[int]) -> int:# 定义状态动作价值函数: dp[i][j],对于前i个物品,当前背包容量为j,最大的可装载价值dp [[0 for j in range(W1)] for i in range(N1)]# 状态动作转移for…...
PCL内置点云类型
PCL内置了许多点云类型供我们使用,下面先介绍PLC内置的点云数据类型 PCL中的点云类型为PointT;至于为什么是PointT类型需要追随到原来的ros开发中去,因为PCL库也是从原来的ROS中剥离出来的;大家都一致的认为点云结构是离散的N维信…...
clickhouse数据结构和常用数据操作
背景, 大数据中查询用mysql时间太长, 使用clickhouse 速度快, 数据写入mysql后同步到clickhouse中 测试1千万数据模糊搜索 mysql 需要30-40秒 clickhouse 约 100ms 一 数据结构和存储引擎 1 查看clickhouse所有数据类型 select * from system.data_type_families; 2 …...
upload-labs关卡9(基于win特性data流绕过)通关思路
文章目录 前言一、靶场需要了解的知识1::$data是什么 二、靶场第九关通关思路1、看源码2、bp抓包修改后缀名3、检查是否成功上传 总结 前言 此文章只用于学习和反思巩固文件上传漏洞知识,禁止用于做非法攻击。注意靶场是可以练习的平台,不能随意去尚未授…...
C++过河卒问题
#include <iostream> #include <cstring> using namespace std;int board[20][20]; // 棋盘 int dp[20][20][20][20]; // 动态规划数组int main() {int x0, y0, x1, y1;cin >> x0 >> y0 >> x1 >> y1; // 输入卒的起点和终点memset(board,…...
【机器学习12】集成学习
1 集成学习分类 1.1 Boosting 训练基分类器时采用串行的方式, 各个基分类器之间有依赖。每一层在训练的时候, 对前一层基分类器分错的样本, 给予更高的权重。 测试时, 根据各层分类器的结果的加权得到最终结果。 1.2 Bagging …...
nodeJs基础笔记
title: nodeJs基础笔记 date: 2023-11-18 22:33:54 tags: 1. Buffer 1. 概念 Buffer 是一个类似于数组的 对象 ,用于表示固定长度的字节序列。 Buffer 本质是一段内存空间,专门用来处理 二进制数据 。 2. 特点 Buffer 大小固定且无法调整Buffer 性能…...
Skywalking流程分析_9(JDK类库中增强流程)
前言 之前的文章详细介绍了关于非JDK类库的静态方法、构造方法、实例方法的增强拦截流程,本文会详细分析JDK类库中的类是如何被增强拦截的 回到最开始的SkyWalkingAgent#premain try {/** 里面有个重点逻辑 把一些类注入到Boostrap类加载器中 为了解决Bootstrap类…...
矩阵的QR分解
矩阵的QR分解 GramSchmidt 设存在 B { x 1 , x 2 , … , x n } \mathcal{B}\left\{\mathbf{x}_{1},\mathbf{x}_{2},\ldots,\mathbf{x}_{n}\right\} B{x1,x2,…,xn}在施密特正交化过程中 q 1 x 1 ∣ ∣ x 1 ∣ ∣ q_1\frac{x_1}{||x_1||} q1∣∣x1∣∣x1 q k …...
STL总结
STL vector 头文件<vector> 初始化,定义,定义长度,定义长度并且赋值,从数组中获取数据返回元素个数size()判断是否为空empty()返回第一个元素front()返回最后一个数back()删除最后一个数pop_back()插入push_back(x)清空clear()begin()end()使用s…...
资深测试总结,现在软件测试有未来吗?“你“的底气在哪里?
目录:导读 前言一、Python编程入门到精通二、接口自动化项目实战三、Web自动化项目实战四、App自动化项目实战五、一线大厂简历六、测试开发DevOps体系七、常用自动化测试工具八、JMeter性能测试九、总结(尾部小惊喜) 前言 1、为什么会有 “…...
Scalable Exact Inference in Multi-Output Gaussian Processes
Orthogonal Instantaneous Linear Mixing Model TY are m-dimensional summaries,ILMM means ‘Instantaneous Linear Mixing Model’,OILMM means ‘Orthogonal Instantaneous Linear Mixing Model’ 辅助信息 作者未提供代码...
树莓派超全系列教程文档--(62)使用rpicam-app通过网络流式传输视频
使用rpicam-app通过网络流式传输视频 使用 rpicam-app 通过网络流式传输视频UDPTCPRTSPlibavGStreamerRTPlibcamerasrc GStreamer 元素 文章来源: http://raspberry.dns8844.cn/documentation 原文网址 使用 rpicam-app 通过网络流式传输视频 本节介绍来自 rpica…...
JVM垃圾回收机制全解析
Java虚拟机(JVM)中的垃圾收集器(Garbage Collector,简称GC)是用于自动管理内存的机制。它负责识别和清除不再被程序使用的对象,从而释放内存空间,避免内存泄漏和内存溢出等问题。垃圾收集器在Ja…...
MVC 数据库
MVC 数据库 引言 在软件开发领域,Model-View-Controller(MVC)是一种流行的软件架构模式,它将应用程序分为三个核心组件:模型(Model)、视图(View)和控制器(Controller)。这种模式有助于提高代码的可维护性和可扩展性。本文将深入探讨MVC架构与数据库之间的关系,以…...
VTK如何让部分单位不可见
最近遇到一个需求,需要让一个vtkDataSet中的部分单元不可见,查阅了一些资料大概有以下几种方式 1.通过颜色映射表来进行,是最正规的做法 vtkNew<vtkLookupTable> lut; //值为0不显示,主要是最后一个参数,透明度…...
Python 包管理器 uv 介绍
Python 包管理器 uv 全面介绍 uv 是由 Astral(热门工具 Ruff 的开发者)推出的下一代高性能 Python 包管理器和构建工具,用 Rust 编写。它旨在解决传统工具(如 pip、virtualenv、pip-tools)的性能瓶颈,同时…...
网站指纹识别
网站指纹识别 网站的最基本组成:服务器(操作系统)、中间件(web容器)、脚本语言、数据厍 为什么要了解这些?举个例子:发现了一个文件读取漏洞,我们需要读/etc/passwd,如…...
多模态图像修复系统:基于深度学习的图片修复实现
多模态图像修复系统:基于深度学习的图片修复实现 1. 系统概述 本系统使用多模态大模型(Stable Diffusion Inpainting)实现图像修复功能,结合文本描述和图片输入,对指定区域进行内容修复。系统包含完整的数据处理、模型训练、推理部署流程。 import torch import numpy …...
从面试角度回答Android中ContentProvider启动原理
Android中ContentProvider原理的面试角度解析,分为已启动和未启动两种场景: 一、ContentProvider已启动的情况 1. 核心流程 触发条件:当其他组件(如Activity、Service)通过ContentR…...
学习一下用鸿蒙DevEco Studio HarmonyOS5实现百度地图
在鸿蒙(HarmonyOS5)中集成百度地图,可以通过以下步骤和技术方案实现。结合鸿蒙的分布式能力和百度地图的API,可以构建跨设备的定位、导航和地图展示功能。 1. 鸿蒙环境准备 开发工具:下载安装 De…...
全面解析数据库:从基础概念到前沿应用
在数字化时代,数据已成为企业和社会发展的核心资产,而数据库作为存储、管理和处理数据的关键工具,在各个领域发挥着举足轻重的作用。从电商平台的商品信息管理,到社交网络的用户数据存储,再到金融行业的交易记录处理&a…...
