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。最终我们得到了

点击阅读全文...

25 Dec

矩阵化简二次型(无穷小近似处理抛物型)

(阅读本文最好有一定的线性代数基础,至少对线性代数里边的基本概念有所了解。)

这学期已经接近尾声了,我们的《解析几何》已经讲到化简二次曲线了。可是,对于没有线性代数的其他同学们,直接用转轴和移轴这个计算公式来变换,那计算量会让我们很崩溃的;虽然那个“不变量”方法计算上有些简单,却总让人感到很诡异,总觉得不知从何而来,而且又要记一堆公式。事实上,如果有线性代数的基础,这些东西变得相当好理解的。我追求用统一的方法求解同一种问题,即用统一的方式处理所有的二次型,当然也希望计算量简单一点。

一般的模型

一般的二次型可以写成
$$x^T A x + 2 b^T x + c=0$$

其中$x,b$都是n维列向量(各元素为$x_i$和$b_i$),A是n阶方阵(各元素为$a_{ij}$),c是常数。在这里,我们只讨论n=2和n=3的情况。化简二次型的过程,可以归结为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$$

点击阅读全文...

26 Dec

费曼路径积分思想的发展(一)

注:这是郝刘祥前辈的一篇论文,98年的时候发表在《自然辩证法通讯》上,里边讲述了费曼以及路径积分的相关故事。我从网上下载下来,原文是很粗糙的pdf文件,我特意将它转化为网页文件,供大家欣赏。有些公式很模糊,所以我已经到图书馆查找了原文,但是由于作者非理论物理专业人员,还不确定部分公式是否正确,请读者慎读。原文较长,将分开几篇来发。如果涉及到版权问题,请作者告之(bojone@spaces.ac.cn),我将会尽快处理掉。

自然辩证法通讯(JOURNAL OF DIALECTICS OF NATURE)
第二十卷总115期,1998第3期

郝刘祥

摘要:该文首先阐述了 Richard Feynman为解决经典电动力学的发散问题而做的艰苦努力,进而论述了这种努力的副产品何以使他偏爱作用量表述,以及他是如何在Dirac文章的启发下得到非相对论量子力学的第三种形式--作用量量子化方案的。文章的第三部分叙述了费曼将其方案推广到相对论情形的尝试和费曼图的由来。最后,该文试图就路径积分方法在量子场论等领域中的广泛应用以及费曼对量子场论的重大疑惑作一简要的说明。

关键词:费曼,作用量,几率幅,路径积分

点击阅读全文...

20 Jan

评论功能修复了

一直都挺奇怪为什么这些天Blog都没有人来评论,难道我的文章质量下降了?可是统计显示访问人数比以往有所增加呀。一直郁闷但却无从找到答案。

今天收到了4anpan读者的邮件,表示评论解锁功能失效了。我自己测试才发现真的如此!!原来不一定是没有人来评论,而是评论插件出了问题。经排除,是最新的公式插件(ASCIIMathMLwFallback2)与滑动解锁的评论插件(IQapTcha)有冲突。所以只好把评论插件换回旧的验证码插件了。正在构思新一步的插件方案。

万分感谢4anpan读者的反馈,希望有更多的读者来帮忙完善科学空间。这里有你更精彩!!

期待你的声音^_^

1 Feb

纠缠的时空(一):洛仑兹变换的矩阵

我现在是越来越佩服爱因斯坦了,他的相对论是他天才的思想的充分体现。只有当相对论提出之后,宏观物理的大多数现象和规律才得到了统一的描述。狭义相对论中爱因斯坦对我们速度叠加常识的否定已经显示了他莫大的勇气,而一项头脑风暴性的工作——广义相对论则将他惊人的创造力体现得完美无瑕。我是被量子力学的数学吸引的,于相对论则是被相对论美妙的逻辑体系吸引。当然,其中也有相当美妙的数学。

狭义相对论中的核心内容之一就是被称为洛仑兹变换的东西,这在相对论发表之前已经由洛仑兹推导出来了,只不过他不承认他的物理意义,也就没有就此进行一次物理革命,革命的任务则由爱因斯坦完成。很久前我就已经看过洛仑兹变换的推导,那是直接设一种线性关系来求解的。但是我总感觉那样的推导不够清晰(也许是我的理解方式有问题吧),而且没有说明狭义相对论的两条原理如何体现出现。所以在研究过矩阵之后,我就尝试用矩阵来推导洛仑兹变换,发现效果挺好的,而且我觉得能够体现出相对论中的对称性。

两条原理

1、狭义相对性原理:在所有惯性系中,物理定律有相同的表达形式。这是力学相对性原理的推广,它适用于一切物理定律,其本质是所有惯性系平权。

2、光速不变原理:所有惯性系中,真空中的光速都等于c=299 792 458 m/s,与光源运动无关。迈克耳孙-莫雷实验是其有力证明。

点击阅读全文...

6 Feb

轻微的扰动——摄动法简介(2)

为了让大家更加熟悉摄动法的基本步骤,本文再讲一个用摄动法解代数方程的例子。这是从实际研究中出来的:
$$\begin{eqnarray*} x=\frac{k(1+k^2+k^4+l^2)}{2(1+k^2)^2} \\ k=\frac{dy}{dx}\end{eqnarray*} $$

这是一道微分方程。要求解这道方程,最好的方法当然是先从第一式解出$k=k(x)$的形式然后再积分。但是由于五次方程没有一般的显式解,所以迫使我们要考虑近似解。当然,一般来说熟悉mathematica的人都会直接数值计算了。我这里只考虑摄动法。

我们将原方程变为下面的形式:
$$x=\frac{k}{2}[1+\frac{l^2}{(1+k^2)^2}]$$

点击阅读全文...

7 Mar

高斯型积分的微扰展开(二)

为什么第二篇姗姗来迟?

其实要写这系列之前,我已经构思好了接下来几篇的内容,本来想要自信地介绍自己想到的一些积分展开的技巧;而且摄动法我本身就比较熟悉,所以正常来说不会这么迟才有第二篇。然而,在我写完第一篇,准备写第二篇的期间,我看到了知乎上的这篇回复:
http://www.zhihu.com/question/24735673

这篇文章大大地拓展了我对级数的认识。里边谈及到了积分的展开是一个渐近级数。这让我犹豫了,怀疑这系列有没有价值,因为渐近级数意味着不管怎样的展开技巧,得到的级数收敛半径都是0。

后来再想想,就算是渐近级数,也有改进的空间,有加速收敛的方法,所以我想我这几篇文章,应该还有一点点意义吧,还可以顺便介绍一下渐近级数和奇点的相关理论。嗯,就这么办吧。

点击阅读全文...