南箕北有斗,牵牛不负轭。

——汉·《明月皎夜光》


一、引言

  在没有深挖这个词之前,只是一直觉得它听起来既优雅又神秘。直到最近看一些Paper,里面涉及到了共轭梯度,再想起之前学过的共轭曲线、共轭复数、共轭转置,忽然觉得这些花狸狐哨的名词背后可能隐藏者一个不被我所知的理论体系,出于好奇和知识回顾的需要,就有了这篇文章。

1.1 轭及共轭

  这个词汉语里面是有的,但是共轭(conjugate/conjugation)这个词(概念)却是从英语里面翻译过来的。与熵(entropy)这类词所不同的是,这个词是有真真切切的实物所对应的,是有非常具象的含义而非译者自己凭空创造出来的。下面这个图就是一副轭:

图 1-1  轭(图片摘自网络)

  上图实际上是一副牛轭,分别将两只牛的牛头塞进那两个木套子,就能驭使它们犁地、耕田、拉货。如果只有一个套子的,就只能驭使一头牛,即为单轭。像上面这样有两个套子的,即为双轭(共轭)。

1.2 共轭关系

  有了图就很好解释和理解了。所谓的共轭关系是什么关系?共轭关系就是这幅双轭里面两头牛之间的关系:既相互制衡,相互约束,相互对立,又相互支撑,相互依存,相互统一。由于科学技术的长足发展,许多原本统一的概念在各个领域有了自己专属的含义。共轭也一样,在不同的科学领域有着不同的定义,共轭这个术语出现的学科包括但不限于以下这些:

  • 1. 数学;
  • 2. 化学;
  • 3. 物理;

  本文主要还是针对其数学领域而言,其他领域的请自行参考相关资料。

二、共轭(数学)

2.1 共轭的类型

  目前数学里面定义的共轭的类型包括但不限于以下内容:

  上述定义还不是数学里面全部的概念,零零总总加起来也已经有十几个了。这也正好验证了我们最开始的猜想:共轭的背后其实有一个庞大的家族,分布在各个研究方向中。
  要把这些概念都介绍一遍是非常耗时耗力的,也没有必要(其实是不了解),有许多高深的理论比如群论暂时我也用不上。因此只选择一些平时用的比较多的来讲,后面有需要的时候再来补充。

2.2 共轭复数

2.2.1 共轭复数的定义

  共轭复数(复共轭)是指实数部分相等,虚数部分互为相反数的一对复数,形如:

(2 - 1)z=a+bjconjz¯=abj

  或者用极坐标的形式表示为:

(2 - 2 - 1)z=reφjconjz¯=reφj

  上述两种形式通过欧拉公式相互连接起来,如下图所示:

图 1-2  共轭复数示意图(图片摘自wikipedia)

  注:对于任意复数z,其共轭复数通常记作z¯,也有的资料记作z

2.2.2 共轭复数的性质

  对任意给定的两复数zw,有:

(2 - 2 - 3){z±w¯=z¯±w¯zw¯=z¯w¯,w0zw¯=z¯w¯zn¯=(z¯)n|z¯|=|z|ez¯=ez¯log(z¯)=log(z)¯,z0

  此外,对任意给定的函数\varphi(x),若它为一个全纯函数(Holomorphic function),则恒有:

(2 - 2 - 4)φ(z¯)=φ(z)¯

2.3 共轭根式

2.3.1 共轭根式的定义

  共轭根式(Conjugate-squre-roots)可能是这里面最简单的了,它是指满足如下形式的一对含根号的式子:

(2 - 3 - 1)a+bdconjabd

  实际上共轭复数是共轭根式的一个特例(a=0, b=1, d=1)。

2.3.2 共轭根式的特性

  共轭根式的主要特性是:其和、积不包含平方根项,如下所示:

(2 - 3 - 2){(a+bd)(abd)=a2db2(a+bd)+(a+bd)=2a

  通常我们会用共轭根式的这个特性来消除分母中的平方根,如下示例:

(2 - 3- 3)1a+bd=abda2db2

  这个技巧相信大家平时几乎都用的炉火纯青了。这里需要提醒的是,对于d的取值,并没有特殊要求,它甚至可以是一个复数

2.4 共轭转置

2.4.1 共轭转置的定义

  共轭转置(也称埃尔米特共轭、埃尔米特转置)定义为满足如下等式的复数矩阵:

(2 - 4 - 1)(A)ij=Aji¯

  式中矩阵A即表示矩阵A的共轭转置矩阵。不同的学科中有不同的记法,但都表示同一个意思。线性代数中通常用A或者AH来表示。

  注:按照Wikipedia的说法,某些情况下A仅仅表示对原矩阵A的元素取复共轭(而不做转置)后形成的矩阵,具体属于哪种情况请结合上下文判断。

  上式等价的语言定义就是:共轭转置矩阵的每一个元素是原矩阵先做转置后对应元素的复共轭。共轭转置矩阵就是原矩阵先取转置矩阵,再将转置后的矩阵的每一个元素取复共轭后形成的矩阵。

  也即上述定义可以等价的表示为如下等式:

(2 - 4 - 2)A=(A¯)T=AT¯

  从上式也可以看出,无论是先对原矩阵做转置,再取复共轭。还是先对原矩阵取复共轭,再做转置。最终的结果都相同。

2.4.2共轭转置的示例

  共轭转置稍微复杂一些,因此我们举个例子来说明,给定如下复数矩阵:

(2 - 4 - 3)A=[1+j232j55j2+jj43+j]

  则其对应的共轭转置为:

(2 - 4 - 4)A=[1j5j25+j43+2j2j3j]

  由上可知,对称矩阵实际上是共轭转置矩阵的一种特殊情况(原矩阵所有元素均为实数)。

2.4.3 共轭转置的性质

  共轭转置有如下运算特性:

(2 - 4 - 5){(A+B)=A+B(rA)=r¯A,r(AB)=BA(A)=Adet(A)=(det A),Atr(A)=(tr A),Aeig(A)=(eig A)¯Ax,y=x,Ay

  上式中涉及矩阵乘法、矩阵与向量乘法的地方默认满足维度关系,不赘述。

2.4.4 共轭转置的延伸

  如果矩阵A是一个方阵,并且其共轭转置矩阵与其自身之间满足一些特定条件,我们就给A一些特殊的称谓。

2.4.4.1 埃尔米特矩阵

  如果矩阵A与其共轭转置满足如下关系:

(2 - 4 - 6)A=Aoraij=aji¯

  我们就称矩阵A为埃尔米特(Hermitian)矩阵(或自伴随矩阵(self-adjoint))。

2.4.4.2 斜埃尔米特矩阵

  如果矩阵A与其共轭转置满足如下关系:

(2 - 4 - 7)A=Aoraij=aji¯

  我们就称矩阵A为斜埃尔米特(skew Hermitian)矩阵(或反埃尔米特矩阵(antihermitian))。

2.4.4.3 正规矩阵

  如果矩阵A与其共轭转置满足如下关系:

(2 - 4 - 8)AA=AA

  我们就称矩阵A为正规(normal)矩阵。

2.4.4.4 酉矩阵

  如果矩阵A与其共轭转置满足如下关系:

(2 - 4 - 9)A=A1

  即矩阵的共轭转置为其逆矩阵,我们就称矩阵A为酉(unitary)矩阵。

2.5 共轭先验/分布

2.5.1 共轭先验/分布的定义

  根据Wikipedia的定义,在贝叶斯概率论中,如果后验概率(poeterior distributions)的概率分布p(θ|x)与其对应的先验概率(prior probability)的概率分布p(θ)属于相同类型的概率分布,则这对先验/后验概率分布被称为共轭分布

  此时,这个先验分布p(θ)就称为似然函数(likelihood function)p(x|θ)共轭先验。共轭先验一定是相对于似然函数而言的。

  这里再对先验概率、后验概率、似然函数做一个说明。

  • 1. p(θ|x):后验概率,是指在给定(已知)数据(样本)x的情况下,模型参数为θ的概率;
  • 2. p(x|θ):似然函数,刚好与后验概率相反,是指在给定(已知)模型参数θ的情况下,获得(观测到)数据x的概率;
  • 3. p(θ):先验概率,是指我们在获得(观测到)任何数据(样本)之前,预先对模型参数的一个估计(经验值),比如模型属于什么类型的概率分布。

  注:只有共轭分布、共轭先验的说法,并没有共轭后验这个术语。

2.5.2 常见的共轭分布

  按照似然函数分布的类型,我们有如下常见的共轭分布表:

2.5.2.1 离散分布

表 2 - 1 常见离散共轭分布函数表

似然函数 共轭先验 模型参数 先验超参数 参数注释
伯努力分布 Beta分布 p α, β αβ的详细解释见这里
二项分布 Beta分布 同上 同上 同上
几何分布 Beta分布 p0 同上 同上
超几何分布 Beta-二项分布 M N, α, β α, β同上
M——抽样数
N——样本总数
泊松分布 Gamma分布 λ k, θ
or
α, β
1θ(β)——时间间隔
k(α)——间隔内总次数
多项分布 Dirichlet分布 p——概率向量
k——概率向量维度
α αi——第i个元素发生次数
2.5.2.2 连续分布

表 2 - 2 常见连续共轭分布函数表

似然函数 共轭先验 模型参数 先验超参数 参数注释
正态分布
(已知τ)
正态分布 μ μ0, σ02 μ0——样本均值
σ02——样本方差
正态分布
(已知τ)
同上 同上 μ0, τ0(1σ02) μ0——同上
τ0——准确率
正态分布
(已知μ)
逆Gamma分布 σ02 α, β σ0=2β
正态分布
(已知μ)
Gamma分布 τ α, β 同上
多元正态分布
(已知)
多元正态分布 μ μ0, 0 μ0——样本均值向量
0——样本协方差阵
指数分布 Gamma分布 λ α, β αβ的详细解释见这里

  注:如果似然函数为指数函数,则其共轭先验必然存在,且通常也属于指数函数。

2.5.3 共轭分布示例

  举个实际的例子可能更容易理解。我们想了解某个大学的男女生比例,于是随机在校园里面进行了抽样,在遇到的100个学生里面,发现有66个是男生。问:该学校男生的比例?
  从频率学派的观点来看,有:

(2 - 5 - 1)θ(pboy)=66100=0.66

  从频率学派的角度来看,这个(概率)值存在唯一真值,这个值是不会变化的。频率派的谬误之处在于,假设该校总共有1000人,我抽了其中800人,恰好800人全是男生,于是该校的男生比例就是100%?
  从贝叶斯学派的角度出发,这个值θ本身属于某个概率分布,还需要结合我们的先验知识(比如往年的统计数据、学校是文科还是理科多等等)来进行修正,确定一个合理的数据。

  为了方便叙述,有如下定义:

  • m——抽样中男生的人数;
  • N——抽样的总人数(样本总数);
  • θ——在该高校中,男生比例;
  • p(θ)——在该高校中,男生比例θ的概率分布(先验概率);
  • X——在该高校中,随机抽取N个学生,其中有m个是男生的事件。

  根据贝叶斯公式:

(2 - 5 - 2)P(θ|X)=P(X|θ)P(θ)P(X)

  由上式可知,在获得了学生数据(观测样本)的情况下,我们要去推断该高校的男生的比例(即后验概率P(θ|X)),需要知道三个值:似然函数、先验概率(分布)、边缘概率。

  边缘概率可以视作归一化因子,那么后验概率最终就只取决于似然函数和先验概率的形式。这个例子的似然函数实际上是二项分布(每次实验只有男生、女生两种结果,)。

  由此有似然函数的表达式:

(2 - 5 - 3)P(X|θ)=CNmθm (1θ)Nm

  查离散表2-1可知,当似然函数为二项分布时,其对应的共轭先验是Beta分布。即如果我们将共轭先验选择为Beta分布,最终可以推导出后验概率也为Beta分布。

  因为θBeta(α, β),由此有似然函数表达式:

(2 - 5 - 4)P(θ)=P(θ;α, β)=Beta(α, β)=θα1(1θ)β101uα1(1u)β1du=θα1(1θ)β1B(α, β)

  式中:B(α, β)=01uα1(1u)β1du=Γ(α)Γ(β)Γ(α+β)Γ(x)伽玛函数

  由于先验分布为连续函数,因此有边缘概率的积分表达式(如果是离散函数则为累积求和表达式):

(2 - 5 - 5)P(X)=01p(X|θ)p(θ)dθ=01CNmθm (1θ)Nmθα1(1θ)β1B(α, β)=CNmB(α, β)01θα+m1(1θ)β+Nm1=CNmB(α, β)B(α+m, β+Nm)

  将上述式子2-5-3、2-5-4、2-5-5代入式2-5-2有:

(2 - 5 - 6)P(θ|X)=CNmθm (1θ)Nmθα1(1θ)β1B(α, β)CNmB(α, β)B(α+m, β+Nm)=CNmθm (1θ)Nmθα1(1θ)β1B(α, β)CNmB(α+m, β+Nm)B(α, β)=θα+m1(1θ)β+Nm1B(α+m, β+Nm)=Beta(α+m, β+Nm)

  我们算出来的后验概率P(θ|X)Beta(α+m, β+Nm),即后验概率仍然服从Beta分布!

  以这个题为例,假设我们抽样的高校是一个综合性大学,男女比例比较接近,由如下的Beta分布的概率密度函数图:

图 1-3  Beta分布概率密度函数图(图片摘自wikipedia)

  由上图推出α=2, β=2(图中紫色曲线)比较合理,即θBeta(2, 2),再结合N=100, m=66可知后验概率P(θ|X)Beta(68,46),如下图所示:

图 1-4  先验/后验概率密度对比函数图

  由上图可知,由于观测数据的加入,原本我们估计的男女比例均衡变成了男生比例偏高。这更符合我们的预期,模型参数会根据观测数据的变化而动态调整,同时结合了已有经验。

  表中剩下的共轭分布感兴趣的可以自己验证。共轭分布最大的优点就是,后验概率会保持与先验概率相同的分布形式。

参考文章

  1. https://blog.csdn.net/baimafujinji/article/details/51374202