30 Nov

算子与线性常微分方程(上)

简介

最近在学习量子力学的时候,无意中涉及到了许多矩阵(线性代数)、群论等知识,并且发现其中有不少相同的思想,其中主要是用算子来表示其对函数的作用和反作用。比如我们可以记$D=\frac{d}{dx}$,那么函数$f(x)$的导数就可以看作是算子D对它的一次作用后的结果,二阶导数则是作用了两次,等等。而反过来,$D^{-1}$就表示这个算子的反作用,它把作用后的函数(像)还原为原来的函数(原像),当然,这不是将求导算子做简单的除法,而是积分运算。用这种思想来解答线性微分方程,有着统一和简洁的美。

线性微分方程是求解一切微分方程的基础,一般来说它形式比较简单,多数情况下我们都可以求出它的通解。在非相对论性量子力学的薛定谔方程中,本质上就是在求解一道二阶偏线性微分方程。另一方面,在许多我们无法求解的非线性系统中,线性解作为一级近似,对于定性分析是极其重要的。

一阶线性常微分方程

这是以下所有微分方程求积的一个基础形式,即$\frac{dy}{dx}+g(x)y=f(x)$的求解。这是通过常数变易法来解答的,其思想跟天体力学中的“摄动法”是一致的,首先在无法求解原微分方程的时候,先忽略掉其中的一些小项,求得一个近似解。即我们先求解
$$\frac{dy}{dx}+g(x)y=0$$

点击阅读全文...

30 Nov

算子与线性常微分方程(下)

不可交换

很自然会想到把这种方法延伸到变系数微分方程的求解,也许有读者回去自己摆弄了一下却总得不到合适的解而感到困惑。在这里群的非Abel性就体现出来了,首先用一个例子来说明一下,我们考虑算子的复合
$$(D-x)(D+x)=D^2-x^2+(Dx-xD)$$

我们要谨慎使用交换律,我们记$[P,Q]=PQ-QP$

其中P和Q是两个算子,此即量子力学中的“对易式”,用来衡量算子P和算子Q的可交换程度,当然,它本身也是一个算子。我们先来求出$[D,x]$给出了什么(要是它是0的话,那就表明运算可以交换了)。究竟它等于什么呢?直接看是看不出的,我们把它作用于一个函数:
$$[D,x]y=(Dx-xD)y=D(xy)-xDy=yDx+xDy-xDy=y$$

由于“近水楼台先得月”,所以$Dxy$表示x先作用于y,然后D再作用于(xy);而$xDy$表示D先作用于y,然后x再作用于Dy。最终我们得到了

点击阅读全文...

23 Dec

“末日”的快乐!

传说中的“末日”,正好是中国传统的冬至节。它的到来并没有引起我们的恐慌,反而,让我有一颗更加平静的心去享受“过冬”的乐趣。

冬至在我们农村是很重要的一个节日,一般来说全家人都要聚在一起吃顿饭,还会包粽子等等。来到广州之后,回去就自然不大方便了,不过我还是想找找那种节日的感觉。于是,下午我就跑到华师西门那里,看看有没有粽子卖。不过发现西门那里基本上都是一些风味小吃,没有那久违的粽子香。不过,忽然想起小飞说她们那里冬至是吃汤圆的,好吧,入乡随俗,我也买了两包汤圆回宿舍煮啦。

点击阅读全文...

14 Dec

关于“微分”的理解

我从来不想在教科书上的定义上纠结太多,因为我知道,真正对定义的理解,需要在长期的实践应用中慢慢感悟的,所以我们唯一需要做的是继续我们的研究。

但是前些天有些朋友问到我关于微分的理解,比如“dx是不是一定很小”等等,所以决定在此写写我的理解。

与微分联系很紧密的,也是我们很熟悉的东西,当然是“增量 ”啦,比如$\Delta y$、$\Delta x$等等,增量显然是可以任意大的(只要自变量还在定义域内)。那么考虑一个函数$y=f(x)$,函数的微分是怎么出现的呢?那是因为我们直接研究函数的增量是比较麻烦的,所以就引入了微分dy,当$\Delta x$很小时,它代表增量的主项:$\Delta y=dy+o(\Delta x)=A \Delta x+o(\Delta x)$,A是一个常数。

点击阅读全文...

18 Dec

黑洞融合的简单模拟

在天文爱好者眼中,黑洞是一个球体,其半径为$\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$$

点击阅读全文...

5 Jan

又折腾数学公式插件了

从2013.11.15开始,使用MathJax插件。主要原因是MathJax在兼容性方面比ASCIIMath Image Fallback Scripts做得好很多。而且从长远考虑,用MathJax也是应该的。
官方网站:http://www.mathjax.org/

复制数学公式:http://www.mathjax.org/demos/copy-and-paste/

-------以下内容已经过时(写于2013.01.05)--------

原来一直是使用“数学研发论坛”完善的数学公式插件来显示数学公式的,使用很简单,载入速度很快,这样一下子就用了三年了。

不过进入大学后,学习的东西越来越多,数学符号也越来越多,郭大哥的插件的不足也暴露出来了。最要命的是它居然无法显示$\hbar$,这叫我这个学习量子力学的孩子情何以堪...(不过郭大哥新版的插件已经加入了这个符号)。还有另外一个不足的地方,就是郭大哥的插件进行了大量的化简,使得数学公式的输入简单了不少,但是反而对标准的Latex代码支持不足了。久而久之,会带来一个弊端,就是迁移性不强。万一哪天这个插件无法使用了,就难以找到替代品了。考虑到这些,我写latex代码的时候总是用标准的语法而不用简化语法,后来$\hbar$的问题出来后,一下子用上了MathJax这个强大的插件(考虑过JsMath,但是发现它的行内公式显示效果不大好)。

点击阅读全文...

7 Jan

角的疑惑——为什么使用弧度?

也许当我们从小学数学进入中学数学的过程中,让我们最郁闷的事情就是课本上把用的好好的角度制改为弧度制了,那个好好的360°的周角无端端变成了一个无理数$2\pi$,为此还多了一堆转换公式,那时这可把我折腾了好一阵子。为什么一个完美的360°不用,反而转向一个无理数$2\pi$?这里边涉及到了相当多的原因,在这些原因中,重新体现了数学体系的一致与简约。当然,文章里的观点只是我自己的看法,仅供大家参考。

弧度制:简约的要求

如果读者已经学过了极限理论,那么我就可以直接说,引入弧度制,是为了在这样的一种角的度量体制下,满足:
$$\lim_{x\to 0} \frac{\sin x}{x}=1$$

点击阅读全文...

23 Jan

校外通过VPN通道访问华师资源

在学校使用校园网时,是很容易访问到华师内部的各个网站的,比如教务系统、图书馆电子资源等,但是如果使用校外网或者是在家时就不那么容易了。另外一种情况是,期末同学回家了,很想早点知道成绩,一般我们会上http://jwc.scnu.edu.cn查询,这个网站在校内校外都可以登陆,但是通常来说为了成绩的录入,会把成绩查询功能关闭掉一段时间,事实上,大部分的成绩都已经录入了。在校园网时,心急的朋友可以访问http://222.201.93.5:211来查询,可是这个网址在校外是不能用的。VPN通道就是为这些校外需求而开通的。

使用方法很简单,打开
https://121.8.171.37/

输入校园卡的账号密码登陆就行了。登录后就会出现校内网的各种链接,包括图书馆资源、教务系统等等。

点击阅读全文...