当前位置: 首页 > news >正文

酉矩阵(Unitary Matrix)和随机矩阵

    先讨论酉矩阵(Unitary Matrix)的性质。

1. 酉矩阵定义

酉矩阵(Unitary Matrix)是复数域上的方阵,满足以下条件:

        U^* U = U U^* = I
其中:

U^* 是 U 的共轭转置(即 Hermitian 转置,U^* = \overline{U}^T )。

I 是单位矩阵。

特殊情形(实数域):
如果 U 是实矩阵,则 U^* = U^T,此时酉矩阵退化为正交矩阵(Orthogonal Matrix),满足:

U^T U = U U^T = I


2. 酉矩阵性质

酉矩阵具有以下重要性质:

可逆性
U 可逆,且 U^{-1} = U^*

保内积
对任意向量 x, y \in \mathbb{C}^n,有:

                \langle Ux, Uy \rangle = \langle x, y \rangle
其中 \langle \cdot, \cdot \rangle 是标准 Hermitian 内积

保范数
\| Ux \| = \| x \|,即酉变换不改变向量的长度

特征值
酉矩阵的所有特征值 \lambda 满足 |\lambda| = 1,即位于复平面的单位圆上。

行列式
|\det(U)| = 1,即行列式的模为 1。

列(行)向量正交性
酉矩阵的列(或行)向量构成一组标准正交基。

乘积封闭性
两个酉矩阵的乘积仍是酉矩阵。

3. 重要定理


谱定理(Spectral Theorem)
任何正规矩阵(A^* A = A A^*)可被酉对角化:

                 A = U \Lambda U^*
其中 \Lambda 是对角矩阵,U 是酉矩阵。

QR 分解
任意矩阵 A 可分解为:

               A = QR
其中 Q 是酉矩阵,R 是上三角矩阵。

量子计算中的酉变换
量子门操作必须由酉矩阵表示,以保证概率守恒(即 \|\psi\|^2 = 1 始终成立)。

4. 计算示例

例1:验证酉矩阵

验证矩阵 U  是否为酉矩阵。

U = \frac{1}{\sqrt{2}} \begin{bmatrix} 1 & i \\ i & 1 \end{bmatrix}

解:

计算共轭转置 U^* :

U^* = \frac{1}{\sqrt{2}} \begin{bmatrix} 1 & -i \\ -i & 1 \end{bmatrix}

计算 U^{^*} U

U^* U = \frac{1}{2} \begin{bmatrix} 1 \cdot 1 + (-i) \cdot i & 1 \cdot i + (-i) \cdot 1 \\ -i \cdot 1 + 1 \cdot i & -i \cdot i + 1 \cdot 1 \end{bmatrix} = \frac{1}{2} \begin{bmatrix} 2 & 0 \\ 0 & 2 \end{bmatrix} = I

因此 U 是酉矩阵。

例2:酉矩阵的特征值


求矩阵 U = \begin{bmatrix} 0 & 1 \\ -1 & 0 \end{bmatrix} 的特征值。

解:

验证 U 是实正交矩阵(实数酉矩阵):

U^T U = \begin{bmatrix} 0 & -1 \\ 1 & 0 \end{bmatrix} \begin{bmatrix} 0 & 1 \\ -1 & 0 \end{bmatrix} = I

特征方程为 \det(U - \lambda I) = 0

\det \begin{bmatrix} -\lambda & 1 \\ -1 & -\lambda \end{bmatrix} = \lambda^2 + 1 = 0 \implies \lambda = \pm i

特征值为 i 和 -i,模均为 1,符合酉矩阵性质。

例3:量子计算中的酉矩阵

Hadamard 门是量子计算中的基本酉矩阵:

H = \frac{1}{\sqrt{2}} \begin{bmatrix} 1 & 1 \\ 1 & -1 \end{bmatrix}

验证 H^* H = I(留作练习)。

5. 应用方面

量子力学中,酉矩阵描述封闭量子系统的演化。

信号处理里,离散傅里叶变换(DFT)矩阵是酉矩阵。

数值线性代数里,用于稳定化算法(如 QR 迭代法求特征值)。

接下来先看随机矩阵,再看随机矩阵与 酉矩阵的关系。

6.随机矩阵(Stochastic Matrix)的定义

    随机矩阵(也称为概率矩阵或马尔可夫矩阵)是指满足以下两个条件的非负实矩阵 P \in \mathbb{R}^{n \times n} :

    行和为1(行随机矩阵):

            \sum_{j=1}^n P_{i,j} = 1 \quad \forall i \in \{1, \dots, n\}

元素非负:

           P_{i,j} \geq 0 \quad \forall i,j

变体:

        列随机矩阵:列和为 1(即 P^T 是行随机矩阵)。

        双随机矩阵:行和与列和均为 1。

7. 酉矩阵的元素模平方矩阵是否为随机矩阵

    设 U 是一个酉矩阵(U^* U = I),定义矩阵 S 为 U 的元素的模平方:

        S_{i,j} = |U_{i,j}|^2

问题:S 是否一定是随机矩阵

7.1. 行和的性质

    由于 U 是酉矩阵,其列向量是标准正交的,因此:

            \sum_{i=1}^n |U_{i,j}|^2 = 1 \quad \forall j

S 的列和为 1(而非行和)。因此:

    S 的列和满足随机矩阵的条件,但行和不一定。

    若 U 的行向量也是标准正交的(即 U 是双酉矩阵,如置换矩阵DFT 矩阵),则 S 的行和也为 1,此时 S 是双随机矩阵。

7.2. 验证

例子1:

考虑酉矩阵:

        U = \frac{1}{\sqrt{2}} \begin{bmatrix} 1 & 1 \\ 1 & -1 \end{bmatrix}

其模平方矩阵为:

        S = \begin{bmatrix} \frac{1}{2} & \frac{1}{2} \\ \frac{1}{2} & \frac{1}{2} \end{bmatrix}

行和与列和均为 1,因此 S 是双随机矩阵。

例子2:

考虑非双酉的酉矩阵:

        U = \begin{bmatrix} 1 & 0 \\ 0 & e^{i\theta} \end{bmatrix}

其模平方矩阵为:

        S = \begin{bmatrix} 1 & 0 \\ 0 & 1 \end{bmatrix}

行和与列和均为 1,但这是特殊情况(对角酉矩阵)。

例子3:

一般酉矩阵的行和可能不为 1。例如:

        U = \frac{1}{2} \begin{bmatrix} 1 & \sqrt{3} \\ -\sqrt{3} & 1 \end{bmatrix}

其模平方矩阵:

        S = \begin{bmatrix} \frac{1}{4} & \frac{3}{4} \\ \frac{3}{4} & \frac{1}{4} \end{bmatrix}

列和为 1,但行和也为 1(巧合)。更复杂的酉矩阵可能破坏行和条件。

7.3. 酉矩阵与随机矩阵的关系

    S 的列和恒为 1(因酉矩阵的列正交性),因此 S^T 是行随机矩阵。

    S 的行和不一定为 1,除非 U 的行向量也正交(即 U 是双酉矩阵)。

    若 U 是双酉矩阵(如置换矩阵、DFT 矩阵),则 S 是双随机矩阵。

7.4. 进一步讨论

量子力学中的概率解释

    在量子测量中,|U_{i,j}|^2  表示从状态 j 跃迁到状态 i 的概率,因此 S 的列和必须为 1(概率守恒)。

双随机矩阵与 Birkhoff-von Neumann 定理

     双随机矩阵可分解为置换矩阵的凸组合,而酉矩阵的模平方矩阵若为双随机,则对应量子通道的“均匀混合”性质。

非双酉矩阵的例外

    若 U 的行向量不正交(如随机生成的酉矩阵),S 的行和可能不为 1,此时 S 仅是列随机矩阵。

最终答案
不一定。酉矩阵 U 的模平方矩阵 S = [|U_{i,j}|^2] 总是列随机矩阵(列和为 1),但仅当 U 的行向量也正交时(即双酉矩阵),S 才是(行)随机矩阵。

典型双酉矩阵(如 Hadamard 矩阵DFT 矩阵)的 S 是双随机矩阵。

http://www.dtcms.com/a/318920.html

相关文章:

  • 在开发板上画出一个2048棋盘的矩阵
  • Jenkins全链路教程——Jenkins用户权限矩阵配置
  • 什么是键值缓存?让 LLM 闪电般快速
  • 面向远程智能终端的超低延迟RTSP|RTMP视频SDK架构与实践指南
  • 动手学深度学习(pytorch版):第一节——引言
  • web前端结合Microsoft Office Online 在线预览,vue实现(PPT、Word、Excel、PDF等)
  • 美食广场: 城市胃的便利店
  • JAVA,Maven继承
  • 开源大模型实战:GPT-OSS本地部署与全面测评
  • Postman接口测试详解
  • SpringBoot微头条实战项目
  • OpenCV入门:图像处理基础教程
  • 【题解】洛谷P3768 简单的数学题[杜教筛]+两种欧反公式解析
  • UDP网络编程chat
  • CompletableFuture的基础用法介绍
  • 技术优势铸就行业标杆:物联网边缘计算网关凭何引领智能变革?
  • 施耐德 Easy Altivar ATV310 变频器:高效电机控制的理想选择(含快速调试步骤及常见故障代码)
  • Flutter 局部刷新方案对比:ValueListenableBuilder vs. GetBuilder vs. Obx
  • 齐护机器人小智AI_MCP图形化编程控制Arduino_ESP32
  • 亚远景-ISO 42001:汽车AI安全的行业标准新趋势
  • 网站 博客遭遇DDoS,CC攻击应该怎样应对?
  • crew AI笔记[2] - 如何选型
  • MCU-TC397的UCB初识
  • 初识 MQ:从同步到异步,聊聊消息队列那些事
  • OpenCv对图片视频的简单操作
  • 深度学习(2):自动微分
  • 学深度学习,有什么好的建议或推荐的书籍?
  • MobileNetV3: 高效移动端深度学习的前沿实现
  • 从“炼金术”到“工程学”:深度学习十年范式变迁与未来十年路线图
  • 深度学习之opencv篇