科学空间|Scientific Spaces

昨天清晨，台风“韦森特”正式来袭我们新兴，话说凌晨三点我已经被风声吵醒了。大概7点钟起床，刚好是台风最抢镜的时候，猛烈地刮呀刮，声音有点像卡车启动的声音......

昨天一整天断电，上午还断了固定电话（农村地区是这样的啦，断电是整体的，台风刮倒了电线杆；断电话是我自己家的问题），中午的时候，固话却自动连上了，郁闷中。下午风雨都基本停下来了，妈妈和我们就出去收拾“残局”，被风刮倒的东西可真不少，尤其是我家门口的两个小棚，惨不忍睹；还有门前的一些盆栽、菜、树等，都倒下了不少。三个人爬上爬下，慢慢维修、收拾。

晚上还是没来电，也好，很久没有尝试过烛光晚餐了。九点多钟的时候，电来了，但是又是一番故障——其他人家中的电都很正常，就我这里灯泡很暗、日光灯启动不了，明显是电压不够的问题。没办法，只好硬着头皮抢修了，排除了很多原因，最后甚至从隔壁家搭电过来，发现我们家的灯还是那么暗（电压不足的问题没有解决）——这说明只有一个可能了，外部电路都没有故障，是我家的内部电路出了问题，猜想某个地方串联了一个用电器分去了电压。但是电线都镶进墙里了，这么黑根本维修不了，没办法，先睡觉了。

点击阅读全文...

在天文爱好者眼中，黑洞是一个球体，其半径为$\frac{2GM}{c^2}$；这是广义相对论的施瓦兹黑洞的结果，也从经典力学推导推导出来，虽然用经典力学是错误的，但是对于多数的天文爱好者（包括笔者）来说，这是目前唯一的一种可行的理解方法（广义相对论那些复杂推导会让我们很崩溃的）。当然，事实上，黑洞不是一个球体，它只是一个密度很大的点。至于密度有多大，目前公认的说法是无穷大，但是严格的物理是不接受这个说法的，或者说，物理是不会接受任何无穷大的说法，所以现在积极发展量子引力理论来统一相对论和量子力学，不过这是另话了。$\frac{2GM}{c^2}$只不过是黑洞的视界，视界之内，我们就什么也不知道了。本文主要就从经典力学的角度探讨一下两个黑洞的合并过程中其视界的变化。读者将会发现，这些视界的形状相当有趣。

经典力学中的黑洞是这样定义的：天体表面的逃逸速度超过了光速，于是连光都无法逃脱，所以这个“洞”就很黑。也就是说，光子的总能量（引力势能与动能之和，经典力学意义下的）要为负，负数表示受到束缚。用数学公式来讲，就是：

$$\frac{1}{2}mc^2 - \frac{GM_1 m}{r_1}-\frac{GM_2 m}{r_2}-...-\frac{GM_n m}{r_n} \leq 0$$

点击阅读全文...

马里乌斯·索菲斯·李

在这篇日志发表之前，科学空间在整个十月就只是在国庆期间发了一篇小感想，这是比较少见的。一个小原因是这学期社团（广播台）方面的活动有点多，当然这不是主要的，其实这个月我大多数课余时间放到了两件事情上：一是无线电路的入门，二就是本文所要讲的《求解微分方程的李对称方法》。

李对称方法主要是通过发现微分方程的对称性来求解微分方程。我首次接触到这个方法是在一本叫《微分方程与数学物理问题》的书上边，书中写得很清晰易懂，后来我还买了类似的《微分方程的对称与积分方法》，后者相对抽象一些，讨论也深入一些。在我目前发现的中文书籍中，这是唯一的两本以李对称方法求解微分方程为主题的书。这两本书还有一个共同特点，就是它们都是外国教材的翻译版。

点击阅读全文...

在文章《新理解矩阵3》：行列式的点滴中，笔者首次谈及到了行列式的几何意义，它代表了n维的“平行多面体”的“体积”。然而，这篇文章写于我初学矩阵之时，有些论述并不严谨，甚至有些错误。最近笔者在写期末论文的时候，研究了超复数的相关内容，而行列式的几何意义在我的超复数研究中具有重要作用，因此把行列式的几何意义重新研究了一翻，修正了部分错误，故发此文，与大家分享。

一个$n$阶矩阵$A$可以看成是$n$个$n$维列向量$\boldsymbol{x}_1,\boldsymbol{x}_2,...,\boldsymbol{x}_n$的集合
$$A=(\boldsymbol{x}_1,\boldsymbol{x}_2,\dots,\boldsymbol{x}_n)$$
从代数的角度来看，这构成了一个矩阵；从几何的角度来看，这$n$个向量可以建立一个平行$n$维体。比如：平行四边形就是“平行二维体”，平行六面体就是“平行三维体”，高阶的只需要相应类比，不需要真正想象出高维空间的立体是什么样。

点击阅读全文...

向量之间的运算有点积和叉积（Cross Product，向量积、外积），其中点积是比较简单的，而且很容易推广到高维；但是叉积不同，一般来说它只不过是三维空间中的东西。叉积的难以推广在于它的多重含义性，如果将向量及其叉积放到张量里边来看（这属于微分形式的内容），那么三维以上的向量叉积是不存在的；但是如果只是把叉积看成是“由两个向量生成第三个与其正交的向量”的工具的话，那么叉积也是可以高维推广的，而且推广的技巧非常巧妙，与三维空间的叉积也非常相似。

回顾三维空间

为了推广三维空间的叉积，首先回顾三维空间的叉积来源是有益的。叉积起源于四元数乘法，但是从目的性来讲，我们希望构造一个向量$\boldsymbol{w}=(w_1,w_2,w_3)$，使得它与已知的两个不共线的向量$\boldsymbol{u}=(u_1,u_2,u_3),\boldsymbol{v}=(v_1,v_2,v_3)$垂直（正交）。从普适性的角度来讲，我们还希望构造出来的向量没有任何“奇点”，为此，我们只用乘法构造。至于叉积的几何意义，则是后话，毕竟，先达到基本的目的再说。

点击阅读全文...

这也是我的期末论文之一...全文共17页，包括了四元数的构造方法，初等应用等。附录包括行列式与体积、三维旋转的描述等。使用LaTex进行写作（LaTex会让你爱上数学写作的）

几何的数与数的几何
――超复数的浅探究

摘要
今天，不论是数学还是物理的高维问题，都采用向量分析为基本工具，数学物理中难觅四元数的影子。然而在历史上，四元数的发展有着重要的意义。四元数（Quaternion）运算实际上是向量分析的“鼻祖”，向量点积和叉积的概念也首先出现在四元数的运算中，四元数的诞生还标记着非交换代数的开端。即使是现在，四元数还是计算机描述三维空间旋转问题最简单的工具。另外，作为复数的推广，四元数还为某些复数问题的一般化提供了思路。

本文把矩阵与几何适当地结合起来，利用矩阵行列式$\det (AB) =(\det A)(\det B)$这一性质得出了四元数以及更高维的超复数的生成规律，并讨论了它的一些性质以及它在描述旋转方面的应用。部分证明细节和不完善的思想放到了附录之中。

点击阅读全文...

进入期末了，很多同学都开始复习了，这学期我选的几门课到现在还不是很熟悉，本想也在趁着这段时间好好看看。偏生五天前我在浏览数学研发论坛的编程擂台时看到了这样的一道题目：

设对于给定的$L$，方程
$$\frac{1}{x}+\frac{1}{y}=\frac{1}{n}$$
满足$0 < x < y \leq L$的正整数解共有$f(L)$种情况。比如$f(6)=1,f(12)=3,f(1000)=1069$，求$f(10^{12})$。

这道题目的来源是Project Euler的第454题：Diophantine reciprocals III（丢潘图倒数方程），题目简短易懂，但又不失深度，正符合我对理想题目的定义。而且最近在学习Python学习得不亦乐乎，看到这道题目就跃跃欲试。于是乎，我的五天时间就没有了，而且过程中几乎耗尽了我现在懂的所有编程技巧。由于不断地测试运行，我的电脑发热量比平时大了几倍，真是辛苦了我的电脑。最后的代码，自我感觉已经是我目前写的最精彩的代码了。在此与大家共享和共勉~

上述表达式是分式，不利于编程，由于$n=\frac{xy}{x+y}$，于是上述题目也等价于求$(x+y)|xy$（意思是$x+y$整除$xy$）的整数解。

点击阅读全文...

含有$n$个一阶常微分方程的一阶常微分方程组
$$\dot{\boldsymbol{x}}=\boldsymbol{A}\boldsymbol{x}$$
其中$\boldsymbol{x}=(x_1(t),\dots,x_n(t))^{T}$为待求函数，而$\boldsymbol{A}=(a_{ij}(t))_{n\times n}$为已知的函数矩阵。现在已知该方程组的$n-1$个线性无关的特解$\boldsymbol{x}_1,\boldsymbol{x}_2,\dots,\boldsymbol{x}_{n-1}$（解的列向量），求方程的通解。

这是我的一位同学在6月5号问我的一道题目，我当时看了一下，感觉可以通过李对称的方法很容易把解构造出来，当晚就简单分析了一下，发现根据李对称的思想，由上面已知的信息确实足以把通解构造出来。但是我尝试了好几天，尝试了几何、代数等思想，都没有很好地构造出相应的正则变量出来，从而也没有写出它的显式解，于是就搁置下来了。今天再分析这道题目时，竟在无意之间构造出了让我比较满意的解来~

点击阅读全文...