一些可能很有用的矩阵知识

一些可有可无的矩阵知识

  • 酉矩阵

酉矩阵

  • 一个服从正态分布的向量乘以一个酉矩阵,得到的向量仍然服从正态分布

酉矩阵是一个复数矩阵,满足其转置的共轭等于其逆矩阵。当一个向量通过一个酉矩阵进行线性变换时,它的模长保持不变,只是发生了旋转和缩放。这意味着如果原始向量服从正态分布,变换后的向量仍将服从相同的正态分布。

p r o o f : proof: proof:

当一个向量服从正态分布时,其概率密度函数(PDF)可以表示为: f ( x ) = 1 2 π σ e − ( x − μ ) 2 2 σ 2 f(x) = \frac{1}{\sqrt{2\pi}\sigma}e^{-\frac{(x-\mu)^2}{2\sigma^2}} f(x)=2π σ1e−2σ2(x−μ)2其中, μ μ μ 是均值, σ σ σ 是标准差。现在,我们有一个酉矩阵U,将向量 x x x 乘以U 得到 y : y = U x y: y=Ux y:y=Ux。

对于 y 的概率密度函数,首先,计算y 的均值 μ y μ_y μy: μ y = E ( y ) = E ( U x ) = U E ( x ) \mu_y = E(y) = E(Ux) = UE(x) μy=E(y)=E(Ux)=UE(x)由于 x x x 服从正态分布且期望是 μ μ μ,则 μ y = U μ x = U μ μ_y=Uμ_x=Uμ μy=Uμx=Uμ,然后,计算 y y y的协方差矩阵 Σ y Σ_y Σy: Σ y = E [ ( y − μ y ) ( y − μ y ) T ] = E [ ( U x − U μ ) ( U x − U μ ) T ] = U E [ ( x − μ ) ( x − μ ) T ] U T \Sigma_y = E[(y - \mu_y)(y - \mu_y)^T] = E[(Ux - U\mu)(Ux - U\mu)^T] = UE[(x - \mu)(x - \mu)^T]U^T Σy=E[(y−μy)(y−μy)T]=E[(Ux−Uμ)(Ux−Uμ)T]=UE[(x−μ)(x−μ)T]UT由于 x x x 服从正态分布且协方差矩阵是 Σ Σ Σ,则 Σ y = U Σ U T Σ_y=UΣU^T Σy=UΣUT,现在,我们可以得到 y y y的概率密度函数 f y ( y ) f_y(y) fy(y):
f y ( y ) = 1 ( 2 π ) n ∣ Σ y ∣ e − 1 2 ( y − μ y ) T Σ y − 1 ( y − μ y ) f_y(y) = \frac{1}{\sqrt{(2\pi)^n|\Sigma_y|}}e^{-\frac{1}{2}(y-\mu_y)^T\Sigma_y^{-1}(y-\mu_y)} fy(y)=(2π)n∣Σy∣ 1e−21(y−μy)TΣy−1(y−μy)将 μ y \mu_y μy和 Σ y \Sigma_y Σy带入可得: f y ( y ) = 1 ( 2 π ) n ∣ Σ ∣ e − 1 2 ( y − U μ ) T ( U Σ U T ) − 1 ( y − U μ ) f_y(y) = \frac{1}{\sqrt{(2\pi)^n|\Sigma|}}e^{-\frac{1}{2}(y-U\mu)^T(U\Sigma U^T)^{-1}(y-U\mu)} fy(y)=(2π)n∣Σ∣ 1e−21(y−Uμ)T(UΣUT)−1(y−Uμ)由于酉矩阵 U 具有单位行列式( ∣ U ∣ = 1 ∣U∣=1 ∣U∣=1)和单位逆矩阵( U − 1 = U T U^{−1}=U^T U−1=UT),上式可简化为: f y ( y ) = 1 ( 2 π ) n ∣ Σ ∣ e − 1 2 ( y − U μ ) T ( U μ ) − 1 ( y − U μ ) f_y(y) = \frac{1}{\sqrt{(2\pi)^n|\Sigma|}}e^{-\frac{1}{2}(y-U\mu)^T(U\mu)^{-1}(y-U\mu)} fy(y)=(2π)n∣Σ∣ 1e−21(y−Uμ)T(Uμ)−1(y−Uμ)这与正态分布的概率密度函数形式相同,只是参数变为 Σ Σ Σ和 U μ U_μ Uμ。因此, y y y 也服从正态分布,其均值为 U μ U_μ Uμ,协方差矩阵为 Σ Σ Σ。

相关推荐
Elastic 中国社区官方博客3 分钟前
Observability:适用于 PHP 的 OpenTelemetry:EDOT PHP 加入 OpenTelemetry 项目
大数据·开发语言·人工智能·elasticsearch·搜索引擎·全文检索·php
ytttr87329 分钟前
Landweber迭代算法用于一维、二维图像重建
人工智能·算法·机器学习
feifeigo12337 分钟前
Matlab编写压缩感知重建算法集
人工智能·算法·matlab
紫小米1 小时前
提示词(Prompt)工程与推理优化
人工智能·ai·prompt·ai agent
子非鱼9211 小时前
1 NLP导论及环境准备
人工智能·自然语言处理
狠活科技1 小时前
Claude Code 重大更新:支持一键原生安装,彻底别了 Node.js
人工智能·aigc·ai编程·claude·claude code
mwq301232 小时前
解密“混合专家模型” (MoE) 的全部魔法
人工智能·llm
能来帮帮蒟蒻吗3 小时前
深度学习(2)—— 神经网络与训练
人工智能·深度学习·神经网络
新加坡内哥谈技术3 小时前
从文字到世界:空间智能是人工智能的下一个前沿
人工智能
oil欧哟3 小时前
文心 5.0 来了,百度大模型的破局之战
前端·人工智能·百度·prompt