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

酉矩阵(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.xdnf.cn/news/17343.html

相关文章:

  • 拥抱云原生:从传统架构到云原生架构的演进与实践
  • 慢SQL优化实战:从一例线上慢SQL探究执行引擎工作过程
  • 如何快速掌握Excel公式?14天轻松通关
  • 大疆前端笔试题目详解
  • Kafka数据生产和发送
  • 如何将Dubbo从Zookeeper平滑地迁移到Nacos?
  • OpenCV图像处理入门实战指南
  • 【办公自动化】使用Python来自动化处理Excel表格有哪些方法?
  • 急危重症专科智能体”构建新一代急诊、手术与重症中心的AI医疗方向探析
  • 图像认知与OpenCV——图像预处理4
  • Java并发与数据库锁机制:悲观锁、乐观锁、隐式锁与显式锁
  • NFS 服务器与iSCSI 服务器
  • 安卓系统属性之androidboot.xxx转换成ro.boot.xxx
  • PostgreSQL 强制索引:当重复数据让优化器“失明”时的解决方案
  • 对话访谈|盘古信息×冠捷科技:全球制造标杆的智能化密码
  • 2025 年华数杯全国大学生数学建模竞赛C 题 可调控生物节律的 LED 光源研究--完整成品、思路、模型、代码、结果分享
  • 本地组策略编辑器无法打开(gpedit.msc命令异常)
  • Windows中安装rustup-init.exe以及cargo build报错443
  • Mac 电脑安装 ADB 环境完整指南
  • 【自动驾驶】《Sparse4Dv3》代码学习笔记
  • golang的面向对象编程,struct的使用
  • DLedger(自动选举)
  • 【多重BFS】Monsters
  • 人工智能——自动微分
  • Spring Boot + ONNXRuntime CPU推理加速终极优化
  • 02电气设计-安全继电器电路设计(让电路等级达到P4的安全等级)
  • PostgreSQL面试题及详细答案120道(61-80)
  • 仁懋电子MOT11N45——音响电路的卓越选择
  • 亚马逊广告运营:有什么好用的辅助工具
  • 接口自动化-pytest