Minkowski代数R1,1\mathbb R_{1,1}R1,1同构于2×22\times22×2实矩阵代数L2(R)L_2(R)L2(R)。一般线性群与特殊线性群与R1,1\mathbb R_{1,1}R1,1的乘法子群的同构关系是:
{G∈R1,1∣G∗G†≠0}≅GL2(R)\{G\in\mathbb R_{1,1}|G^*G^\dag\neq0\}\cong GL_2(R){G∈R1,1∣G∗G†=0}≅GL2(R){G∈R1,1∣G∗G†=1}≅SL2(R)\{G\in\mathbb R_{1,1}|G^*G^\dag=1\}\cong SL_2(R){G∈R1,1∣G∗G†=1}≅SL2(R)
SL2(R)SL_2(R)SL2(R)群是一个三参数群,结构如下:G=KαTβUφG=K_\alpha T_\beta U_\varphiG=KαTβUφ,其中Uφ=eφE/2,Kα=1+αe0=Kα†,Tβ=1+βe=Tβ†U_\varphi=e^{\varphi E/2},K_\alpha=1+\alpha e_0=K_\alpha^\dag,T_\beta=1+\beta e=T_\beta^\dagUφ=eφE/2,Kα=1+αe0=Kα†,Tβ=1+βe=Tβ†。此式的标量参数α,β,φ\alpha,\beta,\varphiα,β,φ的定义区间是[−∞,∞][-\infty,\infty][−∞,∞],我们对此式的兴趣源于它与后面描述的保角群的关系。通篇的其余部分,我们将处理Rn+1,1\mathbb R^{n+1,1}Rn+1,1,经常分解成直和的形式:Rn+1,1=Rn⊕R1,1\mathbb R^{n+1,1}=\mathbb R^{n}\oplus\mathbb R^{1,1}Rn+1,1=Rn⊕R1,1,这种分解称为共形分裂,因为它在本质上与Rn\mathbb R^{n}Rn上的保角群有关。用黑体字表示向量在Rn\mathbb R^{n}Rn的分量是很方便的,并且使用{eo,e}\{e_o, e\}{eo,e}表示R1,1\mathbb R^{1,1}R1,1的空基。Rn+1,1\mathbb R^{n+1,1}Rn+1,1中的任何向量aaa都可以分解为:a=a+αe0+βea=\textbf{a}+\alpha e_0+\beta ea=a+αe0+βe。
共形分裂由R1,1\mathbb R^{1,1}R1,1的单位伪标量EEE唯一确定。令III是Rn+1,1\mathbb R^{n+1,1}Rn+1,1的单位伪标量,那么E~=EI−1=−EI†\widetilde E=EI^{-1}=-EI^\dagE=EI−1=−EI†就是Rn\mathbb R^nRn的伪标量,并且有如下的分裂:a=PE(a)+PE⊥(a)a=P_E(a)+P_E^\perp(a)a=PE(a)+PE⊥(a),其中的射影算子由下式给出:PE(a)=(a⋅E)E=αe0+βe∈R1,1P_E(a)=(a\cdot E)E=\alpha e_0+\beta e\in\mathbb R^{1,1}PE(a)=(a⋅E)E=αe0+βe∈R1,1PE⊥(a)=(a⋅E~)E~†=(a∧E)E=a∈RnP_E^\perp(a)=(a\cdot \widetilde E)\widetilde E^\dag=(a\wedge E)E=\textbf{a}\in\mathbb R^{n}PE⊥(a)=(a⋅E)E†=(a∧E)E=a∈Rn
从上面的式子可以看出,Minkowski平面R1,1\mathbb R^{1,1}R1,1是由EEE唯一确定的,因此称之为EEE平面。投影PE⊥P_E^\perpPE⊥可以视为EEE平面的排斥(rejection)。Rn\mathbb R^nRn中的点a\textbf{a}a用仨向量(trivector)a∧Ea\wedge Ea∧E来表示,这两种表示方法各有优点,但是它们对Rn\mathbb R^nRn的表示是同构的,所以通常忽略它们之间的差别。
Rn\mathbb R^nRn中“点”的齐次坐标的基本思想是通过将Rn\mathbb R^nRn嵌入高维空间来消除原点的影响。为做到这一点,一个有效的方法是射影分裂。令eee是EEE平面上的向量,Rn+1,1\mathbb R^{n+1,1}Rn+1,1中任意不垂直于eee的向量aaa关于eee的射影分裂定义为:ae=a⋅e+a∧e=a⋅e(1+a∧ea⋅e)ae=a\cdot e+a\wedge e=a\cdot e\left(1+\dfrac{a\wedge e}{a\cdot e}\right)ae=a⋅e+a∧e=a⋅e(1+a⋅ea∧e)
这是用双向量a∧ea⋅e\dfrac{a\wedge e}{a\cdot e}a⋅ea∧e来表示向量aaa,这个表示方法与标量无关,因此可以采用一个特定的标量:a⋅e=e0⋅e=−1a\cdot e=e_0\cdot e=-1a⋅e=e0⋅e=−1,这对aaa在Rn\mathbb R^nRn中的部分没影响,因此我们把e∧a=−a∧ee\wedge a=-a\wedge ee∧a=−a∧e称为aaa的射影表示。齐次坐标的经典方法对应于关于非空向量的射影分裂。我们将看到,用空向量表示分裂有很多好处,结果是一种“广义”齐次坐标。
以向量n\textbf nn为法向量且经过aaa点的超平面记为Pn+1(n,a)P^{n+1}(\textbf n,a)Pn+1(n,a),里面的元素xxx满足:n⋅(x−a)=0(x∈Rn+1,1)\textbf n\cdot(x-a)=0(x\in\mathbb R^{n+1,1})n⋅(x−a)=0(x∈Rn+1,1),还可以表示为:n~∧(x∗a)=0(x∈Rn+1,1)\widetilde{\textbf n}\wedge(x*a)=0(x\in\mathbb R^{n+1,1})n∧(x∗a)=0(x∈Rn+1,1)其中n~=nI−1\widetilde{\textbf n}=\textbf nI^{-1}n=nI−1是对偶于n{\textbf n}n的n+1n+1n+1向量。超平面的定义是:Pn+1(e,e0)={x∈Rn+1,1∣e⋅(x−e0)=0}\mathbb P^{n+1}(e,e_0)=\{x\in\mathbb R^{n+1,1}|e\cdot(x-e_0)=0\}Pn+1(e,e0)={x∈Rn+1,1∣e⋅(x−e0)=0}
用e0e_0e0代替Rn\mathbb R^nRn的原点实现了齐次坐标。要使xxx在Rn\mathbb R^nRn对应唯一的x\textbf xx,还需要一个条件。
Homogeneous Model of Euclidean Space
Rn+1,1\mathbb R^{n+1,1}Rn+1,1中所有空向量的集合Nn+1\mathbb N^{n+1}Nn+1称为空锥(null cone)。把点视为空向量,并且使之位于Nn+1\mathbb N^{n+1}Nn+1和Pn+1(e,e0)\mathbb P^{n+1}(e,e_0)Pn+1(e,e0)的交集里面,我们完善了Rn\mathbb R^nRn中点的广义齐次坐标的定义。这里的交集是R2,1\mathbb R^{2,1}R2,1里面的抛物面:Nen=Nn+1⋂Pn+1(e,e0)={x∈Rn+1,1∣x2=0,x⋅e=−1}\mathbb N_e^n=\mathbb N^{n+1}\bigcap\mathbb P^{n+1}(e,e_0)=\{x\in\mathbb R^{n+1,1}|x^2=0,x\cdot e=-1\}Nen=Nn+1⋂Pn+1(e,e0)={x∈Rn+1,1∣x2=0,x⋅e=−1}
它在高维空间的双曲几何中推广为极限球面(horosphere),应用条件x2=0,x⋅e=−1x^2=0,x\cdot e=-1x2=0,x⋅e=−1,我们有:x=x+x2e2+e0x=\textbf x+\dfrac{\textbf x^2e}{2}+e_0x=x+2x2e+e0
这定义了一个由x∈Rn\textbf x\in\mathbb R^nx∈Rn到x∈Nenx\in\mathbb N_e^nx∈Nen的双射。于是,我们可以证明En,Nen,Rn\mathbb E^n,\mathbb N_e^n,\mathbb R^nEn,Nen,Rn是同构的,把Nen\mathbb N_e^nNen称为En,Rn\mathbb E^n,\mathbb R^nEn,Rn的齐次模型,把Nen\mathbb N_e^nNen的点称为齐次点。
我们用“齐次坐标”来描述我们的讨论,因为这是一个标准概念。然而,几何代数使我们能够将一个点表示为一个单独的向量,而不需要将一个向量分解成一组坐标。因此,更可取的说法是“齐次点”而不是“齐次坐标”。
令x=0\textbf x=0x=0,可以发现Rn\mathbb R^nRn的原点的齐次点是e0e_0e0;根据limx2→0x−x⋅e0=limx2→0e+2(x+e0x2)=e\lim\limits_{\textbf x^2\rightarrow0}\dfrac{x}{-x\cdot e_0}=\lim\limits_{\textbf x^2\rightarrow0}e+2\left(\dfrac{\textbf x+e_0}{\textbf x^2}\right)=ex2→0lim−x⋅e0x=x2→0lime+2(x2x+e0)=e可以发现eee是无限远点的齐次点。点的投影表示可以写为:e∧x=e∧x−e⋅x=ex+e∧e0e\wedge x=\dfrac{e\wedge x}{-e\cdot x}=e\textbf x+e\wedge e_0e∧x=−e⋅xe∧x=ex+e∧e0注意,当e⋅x=0e\cdot\textbf x=0e⋅x=0,有e∧x=ex=−xee\wedge\textbf x=e\textbf x=-\textbf xee∧x=ex=−xe。
直线、线段
在讨论几何对象的一般处理方法之前,我们先考虑齐次模型如何描述欧几里德几何中最简单的对象和关系。aaa和bbb的几何积可以展开为:ab=ab+(a−b)e0−(a2+b2)+(ba2−ab2)e+(b2−a2)E2ab=\textbf{ab}+(\textbf a-\textbf b)e_0-\dfrac{(\textbf a^2+\textbf b^2)+(\textbf b\textbf a^2-\textbf a\textbf b^2)e+(\textbf b^2-\textbf a^2)E}{2}ab=ab+(a−b)e0−2(a2+b2)+(ba2−ab2)e+(b2−a2)E根据双向量部分,我们得到:e∧a∧b=e∧(a+e0)∧(b+e0)=ea∧b+(b−a)Ee\wedge a\wedge b=e\wedge(\textbf a+e_0)\wedge(\textbf b+e_0)=e\textbf a\wedge \textbf b+(\textbf b-\textbf a)Ee∧a∧b=e∧(a+e0)∧(b+e0)=ea∧b+(b−a)E我们知道a∧b=a∧(b−a)\textbf a\wedge\textbf b=\textbf a\wedge(\textbf b-\textbf a)a∧b=a∧(b−a)是通过点a\textbf aa、以a−b\textbf a-\textbf ba−b为切线的直线的矩,因此e∧a∧be\wedge a\wedge be∧a∧b完全表征了这条直线。因此,我们将e∧a∧be\wedge a\wedge be∧a∧b解释为一条通过点aaa和bbb的直线,或者更具体地说,是一个一维单形。
标量部分a⋅b=−(a−b)2/2a\cdot b=-(\textbf a-\textbf b)^2/2a⋅b=−(a−b)2/2,因此,两个齐次点的内积可以直接由它们之间的欧氏距离的平方得到。当a2=b2=0a^2=b^2=0a2=b2=0的时候,有:(a−b)2=−2a⋅b=(a−b)2(a-b)^2=-2a\cdot b=(\textbf a-\textbf b)^2(a−b)2=−2a⋅b=(a−b)2。另外,这表明Rn\mathbb R^nRn在Nen\mathbb N_e^nNen中的嵌入是等距的。线段长度的平方是:(e∧a∧b)2=−(e∧a∧b)⋅(e∧a∧b)=[(b∧a)⋅e]⋅[e⋅(a∧b)]=[a−b]⋅[a−b]=(a−b)2\begin{aligned}{} (e\wedge a\wedge b)^2&=-(e\wedge a\wedge b)\cdot(e\wedge a\wedge b)\\ &=[(b\wedge a)\cdot e]\cdot[e\cdot(a\wedge b)]\\ &=[a-b]\cdot[a-b]=(a-b)^2 \end{aligned}(e∧a∧b)2=−(e∧a∧b)⋅(e∧a∧b)=[(b∧a)⋅e]⋅[e⋅(a∧b)]=[a−b]⋅[a−b]=(a−b)2这雨线段的欧几里德距离一致。我们还发现:e∧a∧b∧c=ea∧b∧c+E(b−a)∧(c−a)e\wedge a\wedge b\wedge c=e\textbf a\wedge\textbf b\wedge\textbf c+E(\textbf b-\textbf a)\wedge(\textbf c-\textbf a)e∧a∧b∧c=ea∧b∧c+E(b−a)∧(c−a)这告诉我们可以把a∧b∧c\textbf a\wedge\textbf b\wedge\textbf ca∧b∧c看作是带有切向量(b−a)∧(c−a)(\textbf b-\textbf a)\wedge(\textbf c-\textbf a)(b−a)∧(c−a)的平面的力矩。因此,e∧a∧b∧ce\wedge a\wedge b\wedge ce∧a∧b∧c表示通过点a,b,ca,b,ca,b,c的平面,或者更具体地说,是用这些点作为顶点的三角形(二维单形)。这个三角形的平方是(e∧a∧b∧c)2=[(b−a)∧(c−a)]2(e\wedge a\wedge b\wedge c)^2=[(\textbf b-\textbf a)\wedge(\textbf c-\textbf a)]^2(e∧a∧b∧c)2=[(b−a)∧(c−a)]2。
圆、球
Rn\mathbb R^nRn中以ρ\rhoρ为半径、p\textbf pp为圆心(球心)的圆球可以写为:(x−p)2=ρ2(\textbf x-\textbf p)^2=\rho^2(x−p)2=ρ2,等价的齐次方程式是:x⋅p=−ρ2/2x\cdot p=-\rho^2/2x⋅p=−ρ2/2。利用x⋅e=−1x\cdot e=-1x⋅e=−1,可以把上式化简为x⋅s=0x\cdot s=0x⋅s=0,其中向量sss有如下性质:{s=p−ρ2e/2=p+e0+p2−ρ22es2=ρ2>0s⋅s=−1\begin{cases} s=p-\rho^2e/2=\textbf p+e_0+\dfrac{\textbf p^2-\rho^2}{2}e\\ s^2=\rho^2>0\\ s\cdot s=-1 \end{cases}⎩⎪⎪⎨⎪⎪⎧s=p−ρ2e/2=p+e0+2p2−ρ2es2=ρ2>0s⋅s=−1因此Rn+1,1\mathbb R^{n+1,1}Rn+1,1中每一个向量sss都唯一确定了一个Rn\mathbb R^nRn中的圆球。从另一个角度来看,圆球可以用sss的对偶s~=sI−1\widetilde s=sI^{-1}s=sI−1来表示。因为I†=(−1)ϵI=−I−1,ϵ=(n+1)(n+2)/2I^\dag=(-1)^\epsilon I=-I^{-1},\epsilon=(n+1)(n+2)/2I†=(−1)ϵI=−I−1,ϵ=(n+1)(n+2)/2,我们可以改写sss的性质:{s2=−s~†s~=ρ2s⋅e=e⋅(s~I)=(e∧s~)I=−1\begin{cases} s^2=-\widetilde s^\dag\widetilde s=\rho^2\\ s\cdot e=e\cdot(\widetilde sI)=(e\wedge\widetilde s)I=-1 \end{cases}{s2=−s†s=ρ2s⋅e=e⋅(sI)=(e∧s)I=−1上述球面还可以写为对偶形式:x∧s~=0x\wedge\widetilde s=0x∧s=0。正如后面所提到的,s~\widetilde ss的优点是它可以直接从球面上的点计算出来。然后通过对偶得到sss,求出球的中心。这种球面的对偶性在计算和概念上都非常强大。
欧几里德几何
对于E2\mathbb E^2E2来说,齐次模型的优点在一个例子中体现得最为明显:
Simson定理\textbf{Simson定理}\quadSimson定理给定△ABC\triangle ABC△ABC,DDD是平面上的一点,从DDD作三角形三边的垂线,垂足分别是A1,B2,C1A_1,B_2,C_1A1,B2,C1。那么点A1,B2,C1A_1,B_2,C_1A1,B2,C1在一条直线上当且仅当DDD在△ABC\triangle ABC△ABC的圆周上。
暂时搁置用小写字母表示向量的约定,我们将上中的标签解释为E2\mathbb E^2E2中的齐次点。我们有几何代数来表述关系以便于分析。我们把△ABC\triangle ABC△ABC记为e∧A∧B∧Ce\wedge A\wedge B\wedge Ce∧A∧B∧C,它的边ABABAB记为e∧A∧Be\wedge A\wedge Be∧A∧B。这种记法这融合了综合几何的表达优势和几何代数的计算能力。在证明Simson定理之前,我们引用了欧几里德几何中的一些基本结果。首先,三角形e∧A∧B∧Ce\wedge A\wedge B\wedge Ce∧A∧B∧C的外接圆是:s~=A∧B∧C\widetilde s=A\wedge B\wedge Cs=A∧B∧C后文将给出一个一般性的方法,证明它确实是这三个点的外接圆。下面我们来计算半径:ρ2=s2(s⋅e)2=s~†s~(s∧e~)2=(C∧B∧A)⋅(A∧B∧C)(e∧A∧B∧C)⋅(e∧A∧B∧C)\rho^2=\dfrac{s^2}{(s\cdot e)^2}=\dfrac{\widetilde s^\dag\widetilde s}{(s\wedge\widetilde e)^2}=\dfrac{(C\wedge B\wedge A)\cdot(A\wedge B\wedge C)}{(e\wedge A\wedge B\wedge C)\cdot(e\wedge A\wedge B\wedge C)}ρ2=(s⋅e)2s2=(s∧e)2s†s=(e∧A∧B∧C)⋅(e∧A∧B∧C)(C∧B∧A)⋅(A∧B∧C)
上式右边是两个行列式的比值,当展开时,ρ2\rho^2ρ2就表示为三角形边的长度的关系式。分子可以展开为:
(A∧B∧C)2=−∣0A⋅BA⋅CB⋅A0B⋅CC⋅AC⋅B0∣=−2(A⋅B)(B⋅C)(C⋅A)=(A−B)2(B−C)2(C−A)24=(A−B)2(B−C)2(C−A)24\begin{aligned}{}
(A\wedge B\wedge C)^2&=-\begin{vmatrix}
0&A\cdot B&A\cdot C\\
B\cdot A&0&B\cdot C\\C\cdot A&C\cdot B&0
\end{vmatrix}=-2(A\cdot B)(B\cdot C)(C\cdot A)\\
&=\dfrac{(A-B)^2(B-C)^2(C-A)^2}{4}\\&=\dfrac{(\textbf A-\textbf B)^2(\textbf B-\textbf C)^2(\textbf C-\textbf A)^2}{4}
\end{aligned}(A∧B∧C)2=−∣∣∣∣∣∣0B⋅AC⋅AA⋅B0C⋅BA⋅CB⋅C0∣∣∣∣∣∣=−2(A⋅B)(B⋅C)(C⋅A)=4(A−B)2(B−C)2(C−A)2=4(A−B)2(B−C)2(C−A)2分母展开为:(e∧A∧B∧C)2=−4S△ABC2=[(B−A)⋅(C−A)]2−(B−A)2(C−A)2=[(B−A)⋅(C−A)]2−(B−A)2(C−A)2\begin{aligned}{}
(e\wedge A\wedge B\wedge C)^2&=-4S_{\triangle ABC}^2\\
&=[(\textbf B-\textbf A)\cdot(\textbf C-\textbf A)]^2-(\textbf B-\textbf A)^2(\textbf C-\textbf A)^2\\&=[(B-A)\cdot(C-A)]^2-(B-A)^2(C-A)^2
\end{aligned}(e∧A∧B∧C)2=−4S△ABC2=[(B−A)⋅(C−A)]2−(B−A)2(C−A)2=[(B−A)⋅(C−A)]2−(B−A)2(C−A)2归一化A∧B∧CA\wedge B\wedge CA∧B∧C并取其对偶,我们可以确定圆心PPP,因此:−A∧B∧C~e∧A∧B∧C~=P−ρ2e/2\dfrac{-\widetilde{A\wedge B\wedge C}}{\widetilde{e\wedge A\wedge B\wedge C}}=P-\rho^2e/2e∧A∧B∧C−A∧B∧C=P−ρ2e/2
为了把圆A∧B∧CA\wedge B\wedge CA∧B∧C与点DDD联系起来,我们使用∨\vee∨操作:(A∧B∧C)∨D=A∧B∧C~⋅D=−A∧B∧C∧D~{(A\wedge B\wedge C)\vee D=\widetilde{A\wedge B\wedge C}}\cdot D=-\widetilde{A\wedge B\wedge C\wedge D}(A∧B∧C)∨D=A∧B∧C⋅D=−A∧B∧C∧D进而有A∧B∧C∧D=ρ2−δ22e∧A∧B∧C(其中δ2=−2P⋅D)A\wedge B\wedge C\wedge D=\dfrac{\rho^2-\delta^2}{2}e\wedge A\wedge B\wedge C(\text{其中}\delta^2=-2P\cdot D)A∧B∧C∧D=2ρ2−δ2e∧A∧B∧C(其中δ2=−2P⋅D)
左边=0当且仅当DDD在圆A∧B∧CA\wedge B\wedge CA∧B∧C上,相应的,右边的ρ=δ\rho=\deltaρ=δ。
为了用代数方法构造Simson三角形,我们需要求出点DDD在e∧A∧Ce\wedge A\wedge Ce∧A∧C的垂足B1B_1B1。利用非齐次点,我们可以把垂足的条件写成:(B1−D)⋅(C−A)=0\bf(B_1-D)\cdot(C-A)=0(B1−D)⋅(C−A)=0,因此(B1−D)(C−A)=(B1−D)∧(C−A)=(A−D)∧(C−A)\bf(B_1-D)(C-A)=(B_1-D)\wedge(C-A)=(A-D)\wedge(C-A)(B1−D)(C−A)=(B1−D)∧(C−A)=(A−D)∧(C−A),两边除以C−A\bf C-AC−A:B1−D=[(A−D)∧(C−A)]⋅(C−A)−1=A−D−(A−D)⋅(C−A)−1(C−A)\begin{aligned}{} \bf B_1-D&=\bf[(A-D)\wedge(C-A)]\cdot(C-A)^{-1}\\&=\bf A-D-(A-D)\cdot(C-A)^{-1}(C-A) \end{aligned}B1−D=[(A−D)∧(C−A)]⋅(C−A)−1=A−D−(A−D)⋅(C−A)−1(C−A)所以B1=A+(D−A)⋅(C−A)(C−A)2(C−A)\bf B_1=A+\dfrac{(D-A)\cdot(C-A)}{(C-A)^2}(C-A)B1=A+(C−A)2(D−A)⋅(C−A)(C−A)我们可以很容易地把它转换成齐次点之间的关系。然而,我们只对Simson三角形e∧A1∧B1∧C1e\wedge A_1\wedge B_1\wedge C_1e∧A1∧B1∧C1感兴趣:e∧A1∧B1∧C1=E(B1−A1)∧(C1−A1)=E(A1∧B1+B1∧C1+C1∧A1)\begin{aligned}{} e\wedge A_1\wedge B_1\wedge C_1&=E\bf(B_1-A_1)\wedge(C_1-A_1)\\ &=E\bf(A_1\wedge B_1+B_1\wedge C_1+C_1\wedge A_1) \end{aligned}e∧A1∧B1∧C1=E(B1−A1)∧(C1−A1)=E(A1∧B1+B1∧C1+C1∧A1)不失一般性,把D\bf DD视为Rn\mathbb R^nRn的原点,计算得以简化。注意到δ2=−2P⋅D=p2\delta^2=-2P\cdot D=\textbf p^2δ2=−2P⋅D=p2,令D=0\textbf D=0D=0,有:e∧A1∧B1∧C1=(ρ2−δ24ρ2e∧A∧B∧C)e\wedge A_1\wedge B_1\wedge C_1=\left(\dfrac{\rho^2-\delta^2}{4\rho^2}e\wedge A\wedge B\wedge C\right)e∧A1∧B1∧C1=(4ρ2ρ2−δ2e∧A∧B∧C)计算中唯一棘手的部分是得到此式右侧的系数。利用上文关于δ,ρ\delta,\rhoδ,ρ的展开形式,以及A∧B∧C∧D=ρ2−δ22e∧A∧B∧CA\wedge B\wedge C\wedge D=\dfrac{\rho^2-\delta^2}{2}e\wedge A\wedge B\wedge CA∧B∧C∧D=2ρ2−δ2e∧A∧B∧C,我们可以得到下面的恒等式:e∧A1∧B1∧C1=A∧B∧C∧D2ρ2e\wedge A_1\wedge B_1\wedge C_1=\dfrac{A\wedge B\wedge C\wedge D}{2\rho^2}e∧A1∧B1∧C1=2ρ2A∧B∧C∧D
这证明了Simson定理,当且仅当D在圆上时,右边消去了,当且仅当三个点在同一条直线上时,左边消去了。