3.2 埃尔米特转置

news/2024/11/24 9:53:53/

定义

  对于复矩阵,转置又不一样,常见的操作是共轭转置,也叫埃尔米特转置Hermitian transpose。埃尔米特转置就是对矩阵先共轭,再转置,一般来说用三种符号表示埃尔米特转置:

  1. 第一种符号是AHA^HAH,这是国内教材通用的做法,H是埃尔米特名字首字母;
  2. 第二种符号是A∗A^*A,这是国外教材喜欢用,这个符号在国内教材表示伴随矩阵,如以下文档:
    在这里插入图片描述
  3. 第三种符号是匕首符号A†A^{\dagger}A,但是有时候也用来表示矩阵的加号逆。

  求埃尔米特转置的代码比较简单,python就一行代码:

    # 埃尔米特转置def hermitian_transpose(self):return Matrix([[e.conjugate() for e in v] for v in self.__vectors]).transpose_matrix()

  测试了一个矩阵:
(1−i6−i2−8i2+i5+i4−i)H=(1+i2−i6+i5−i2+8i4+i)\begin{pmatrix}1-i & 6-i & 2-8i\\ 2+i & 5+i & 4-i\\ \end{pmatrix}^H= \begin{pmatrix}1+i & 2-i\\ 6+i & 5-i\\ 2+8i & 4+i\\ \end{pmatrix} (1i2+i6i5+i28i4i)H=1+i6+i2+8i2i5i4+i

埃尔米特阵

  如果一个矩阵,埃尔米特转置后还是它自己,这样的矩阵就是埃尔米特阵。毫无疑问,矩阵必须得是一个方阵。所以它的判断方式也很简单,首先判断是否为方阵,再以对角线为对称轴判断就完事了,但是要注意数据类型,把复数和其他类型区分开来,所以代码会稍微长一点:

# 是否埃尔米特阵def is_hermitian(self):m = len(self.__vectors[0])n = len(self.__vectors)if m != n:return False# 遍历每一行对角线以上的元素for i in range(m):for j in range(i+1, n):e = self.__vectors[j][i]f = self.__vectors[i][j]if isinstance(e, complex):if e != f.conjugate():return Falseelse:if e != f:return Falsereturn True

  比如以下矩阵就是一个埃尔米特阵:
(1−i2−i3+i2+i5+i3+i3−i3−i3−i)\begin{pmatrix}1-i & 2-i & 3+i\\ 2+i & 5+i & 3+i\\ 3-i & 3-i & 3-i\\ \end{pmatrix} 1i2+i3i2i5+i3i3+i3+i3i

酉矩阵

  一个方阵的逆矩阵恰好是自己的埃尔米特转置,这样的矩阵被称为酉矩阵unitary matrix,也就是:
AAH=AHA=IAA^H=A^HA=I AAH=AHA=I
  AAH=AHAAA^H=A^HAAAH=AHA这个定义就限制了必须为方阵。所以它的判断也比较简单:

   # 是否为酉矩阵def is_unitary(self):m = len(self.__vectors[0])n = len(self.__vectors)if m != n:return Falsex = self * self.hermitian_transpose()return x.is_identity()# 是否为单位矩阵def is_identity(self):m = len(self.__vectors[0])n = len(self.__vectors)if m != n:return Falsefor i in range(n):for j in range(n):if i == j:if self.__vectors[i][j] != 1:return Falseelif self.__vectors[i][j] != 0:return Falsereturn True

  比如以下两个矩阵就是一个酉矩阵:
(0.5−0.5i−0.5+0.5i0.5i0.50.5+0.5i0.5+0.5i−0.5+0.5i0)(0.5−0.5i0.5−0.5i0.5i0.5−0.5−0.5i−0.5−0.5i0.5−0.5i0)=(100010001)\begin{pmatrix}0.5 & -0.5i & -0.5+0.5i\\ 0.5i & 0.5 & 0.5+0.5i\\ 0.5+0.5i & -0.5+0.5i & 0\\ \end{pmatrix}\begin{pmatrix}0.5 & -0.5i & 0.5-0.5i\\ 0.5i & 0.5 & -0.5-0.5i\\ -0.5-0.5i & 0.5-0.5i & 0\\ \end{pmatrix}\\=\begin{pmatrix}1 & 0 & 0\\ 0 & 1 & 0\\ 0 & 0 & 1\\ \end{pmatrix} 0.50.5i0.5+0.5i0.5i0.50.5+0.5i0.5+0.5i0.5+0.5i00.50.5i0.50.5i0.5i0.50.50.5i0.50.5i0.50.5i0=100010001


http://www.ppmy.cn/news/23604.html

相关文章

java 自定义json解析注解 复杂json解析

java 自定义json解析注解 复杂json解析 工具类 目录java 自定义json解析注解 复杂json解析 工具类1.背景2、需求-各式各样的json一、一星难度json【json对象中不分层】二、二星难度json【json对象中出现层级】三、三星难度json【json对象中存在数组】四、四星难度json【json对象…

麻省理工出版 | 2023年最新深度学习综述手册

UCL Simon Prince的新书:《Understanding Deep Learning》 ,在2023年2月6日由MIT Press出版。他之前写过很受欢迎的《Computer Vision: Models, Learning, and Inference》。 关于这本最新的深度学习手册,作者这样介绍它: 正如书…

target31对传感器采样率限制

target31谷歌做了一些列的适配,其中对看似无关紧要的传感器采用率做了限制,最近几年很火的摇一摇开屏广告就是基于传感器去实现的。这篇文章,我们就一起去看一下target31对传感器采样频率的限制。 一、背景 先说一下背景:近期&am…

「Vue源码学习」常见的 Vue 源码面试题,看完可以说 “精通Vue” 了吗?

Vue源码面试题一、行时(Runtime) 编译器(Compiler) vs. 只包含运行时(Runtime-only)webpackRollupBrowserify二、Vue 的初始化过程(面试关问:new Vue(options) 发生了什么&#xff1…

python学习笔记-查看数据结构、均值、中位数、分位数、众数、离中趋势(标准差、方差、求和、偏态系数、风险系数)正态分布pdf、cdf、ppf

①引入pandas包,命名为pd。 import pandas as pd ②读入HR.csv数据 dfpd.read_csv(“./data/HR.csv”) ③查看是什么结构 type(df) ④查看单个类别satisfaction_level的数据结构 type(df[“satisfaction_level”]) ⑤查看均值的数据结构 type(df.mean()) …

Vue3.0 setup的使用及作用

目录开篇:1.什么是setup2.setup怎么使用3.setup中包含的生命周期函数3.setup相关参数4.setup特性总结总结开篇: 从vue2升级 vue3,vue3是可以兼容vue2。所以v3可以采用v2的选项式api,但是v2不能使用v3的组合式api,由于…

【项目精选】基于Java的银行排号系统的设计与实现

银行排号系统是为解决一些服务业营业大厅排队问题而设计的,它能够有效地提高工作人员的工作效率,也能够使顾客合理的安排等待时间,让顾客感到服务的公平公正。论文首先讨论了排号系统的背景、意义、应用现状以及研究与开发现状。本文在对C/S架…

一篇文章读懂阿里云企业级数据库最佳实践

今天阿里数据库不再是简单的电商业务,而是涵盖了视频娱乐、IM、地图、在线零售、新零售、物流、在线旅游、音乐、IoT等纵多领域。2017年双十一交易额达1682亿,数据库交易峰值也以数十倍的速度在增长。超大规模的业务压力,在阿里巴巴内部淬炼出…