25 Dec

《新理解矩阵5》:体积=行列式

在文章《新理解矩阵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$维体。比如:平行四边形就是“平行二维体”,平行六面体就是“平行三维体”,高阶的只需要相应类比,不需要真正想象出高维空间的立体是什么样。

点击阅读全文...

26 Dec

体积与阿达马不等式

阿达马不等式
设有$n$阶实矩阵$\boldsymbol{A}=(a_{ij})_{n\times n}$,那么它的行列式满足阿达马(Hadamard)不等式
$$\begin{equation}
\left(\det \boldsymbol{A}\right)^2 \leq \prod\limits_{i=1}^{n}\left(a_{1i}^2+a_{2i}^2+\dots+a_{ni}^2\right)
\end{equation}$$

这是阿达马在1893年首先发表的。根据体积就是行列式的说法,上述不等式具有相当明显的几何意义。当$n=2$时,它就是说平行四边形的面积不大于两边长的乘积;当$n=3$时,它就是说平行六面体的体积不大于三条棱长的乘积;高维可以类比。这些结论在几何中几乎都是“显然成立”的东西。因此很难理解为什么这个不等式在1893年才被发现。当然,代数不会接受如此笼统的说法,它需要严格的证明。

点击阅读全文...

29 Dec

有质动力:倒立单摆的稳定性

前几天在“宇宙的心弦”浏览网页时,发现他更新了一篇很有趣的文章,叫《倒立单摆的稳定性与Ponderomotive Force》(果然,物理系的能接触到各种各样有趣的现象),里边谈到通过施加一个运动在单摆上面,倒立的单摆也可以是稳定的。这勾起了我的兴趣,遂也计算了一番。

点击阅读全文...

11 Mar

一维弹簧的运动(上)

我们通常用一个波动方程来描述弦的振动,但是,弦的振动是二维的,也就是说,它的“波”是在垂直方向的位移。让我们来考虑一根一端固定的一维理想弹簧,胡克系数为$k$,它的松弛状态是均匀的,线密度是$\rho$,长度是$l$,质量是$m$。

如何弹?
我们要分析这根弹簧的运动,即给定弹簧的初始状态,看弹簧的密度如何变化,这种情况类似于“横波”。但是,弹簧本身是连续介质,这是我们不熟悉的,但是我们可以将它离散化,将它看成无数个小质点的弹簧链。如下图

离散的弹簧

离散的弹簧

点击阅读全文...

4 Mar

平面曲线的曲率的复数表示

开学已经是第二周了,我的《微分几何》也上课两周了,进度比较慢,现在才讲到平面曲线的曲率。在平面曲线$\boldsymbol{t}(t)=(x(t),y(t))$某点上可以找出单位切向量。
$$\boldsymbol{t}=\left(\frac{dx}{ds},\frac{dy}{ds}\right)$$
其中$ds^2 =dx^2+dy^2$,将这个向量逆时针旋转90度之后,就可以定义相应的单位法向量$\boldsymbol{n}$,即$\boldsymbol{t}\cdot\boldsymbol{n}=0$。

常规写法

让我们用弧长$s$作为参数来描述曲线方程,$\boldsymbol{t}(s)=(x(s),y(s))$,函数上的一点表示对$s$求导。那么我们来考虑$\dot{\boldsymbol{t}}$,由于$\boldsymbol{t}^2=1$,对s求导得到
$$\boldsymbol{t}\cdot\dot{\boldsymbol{t}}=0$$

点击阅读全文...

13 Mar

一维弹簧的运动(下)

在上一篇文章中,我们得到了一维弹簧运动的方程
$$m\frac{\partial^2 X}{\partial t^2}=k\frac{\partial^2 X}{\partial \xi^2}$$
并且得到了通解
$$X=F(u)+H(v)=F(\xi+\beta t)+H(\xi-\beta t)$$
或者
$$X(\xi,t)=\frac{1}{2}\left[X_0(\xi+\beta t)+X_0(\xi-\beta t)\right]+\frac{1}{2\beta}\int_{\xi-\beta t}^{\xi+\beta t} X_1 (s)ds$$
在文章的末尾,提到过这个解是有些问题的。现在让我们来详细分析它。

点击阅读全文...

19 Jan

宇宙驿站服务器升级完毕

这一周科学空间时断时续的,原因是原来的服务器两个内存条坏了,内存不够用。

后来天文台决定给我们换一台服务器,这两天主要在转移数据,从而不能访问。

目前,基本上已经转移好了,服务器升级工作基本完成。新服务器的升级,CPU从原来的8核升级为48核,内存从16GB升级为64GB。再次感谢国家天文台宇宙驿站给予我们的服务^_^感谢各位技术人员的努力,让我们一起把中文科普事业做得更好~

14 Feb

情人节?元宵节!

元宵节快乐

元宵节快乐

今天是2014年2月14日,农历正月十五,传统的元宵佳节,祝大家元宵节快乐!

不过虽说是元宵佳节,但是我们这里的习俗却没有闹元宵的,好像在我们这里元宵节就像普通的初一十五一样,惯例地上个香,祭下神而已,唯一特别的地方就是早上妈妈放了个鞭炮,什么汤圆、灯笼、灯谜都没有呢。不过这并不妨碍我欣赏元宵节,印象里好像上学以来这是第一次在家过元宵节。幸好没有参加美国数学建模,不然又少了半个月的假期,少了一次难得的元宵,而多了得不偿失的劳动...

今天也是西方的情人节,但在这里我只强调元宵节。首要原因却不是我目前单身(当然这也是原因之一^_^),而是元宵节是中国传统节日。我这个人有个奇怪的“嗜好”,反正越潮流的东西我越不跟。于是乎,既然那么多人都庆祝着西方节日(什么万圣节、圣诞节、情人节),那么我就偏不凑这个热闹。我又想起了去年圣诞前夕有个师弟过来向我们宣传和推销圣诞的东西,被我批了一顿,我直言说“你为什么不等元旦再来?”。我想,如果哪一天,我也有机会庆祝情人节,我也只是庆祝中国的情人节,总感觉中国的情人节美多了:七夕,Qixi Festival,多美!不论是典故还是习俗都更美~

点击阅读全文...