22
May
当Matlab遇上牛顿法
By 苏剑林 | 2013-05-22 | 61134位读者 | 引用牛顿法是求方程近似根的一个相当有用而且快捷的方法，我们最近科学计算软件课程（Matlab）的一个作业就是编写求方程近似解的程序，其中涉及到牛顿法。我们要实现的目标是，用户输入一道方程，脚本就自动求出根来。这看起来是一个挺简单的循环迭代程序，但是由于Matlab本身的特殊性，却产生了不少困难。
Matlab是为了数值计算（尤其是矩阵运算）而生的，因此它并不擅长处理符号计算。这就给我们编程带来了困难。在网上随便一搜，就可以发现，网上的Matlab牛顿法程序都是要求用户同时输入方程及其导函数，这显然是不方便的，因为Matlab本身就具备了求导功能。下面我们来分析一下困难在哪里。
我们要实现的最基本功能是定义一个函数，然后可以根据该函数求具体的函数值，并且自动求该函数的导数，接着求导数值。这些看起来很基本的功能在Matlab中却很难调和，因为Matlab的“函数”定义很广，一个具有特定功能的M文件叫“函数”，一个运算式$f(x)$也可能是一个函数，显然后者是可以求导的，前者却不行，所以Matlab一刀砍——不能对函数求导！！
欢聚兴隆，畅言科普
记信息时代的天文科普研讨会暨第三届宇宙驿站站长联谊会
在信息时代的今天，利用互联网相互交流以及查找各种资讯已经成为了许多天文爱好者的必经之道。同好们也许都浏览过牧夫天文论坛、星友空间站、空间天文网等天文科学网站，事实上，它们都源于一个共同的科普网站群体——宇宙驿站。正如她的名字所言，宇宙驿站是我们一大群天文爱好者在互联网上的“家”，她为我们这群热衷于网络科普的站长免费提供了稳定的网站空间。
宇宙驿站发起于2002年，是国家天文台LAMOST项目之一，迄今已经有近百位站长在上面“安家”。2013年6月28日到6月30日，我们这群站长齐聚兴隆，开展了一次别开生面的会议——“信息时代的天文科普研讨会暨第三届站长联谊会”。
农村的孩子免不了常做家务，当然我家也没有什么特别沉重的家务，通常都是扫地、做饭、洗菜这些简单的活儿。说到洗菜，洗完菜后总喜欢边放水边搅水，然后就在水面上形成一个颇为有趣的漩涡。现在我们从数学物理的角度来分析一下这个漩涡。
在讲洗手盆的漩涡之前，我们先来看一下一个比较类似的、更古老的问题——牛顿的旋转液面问题。牛顿假设有一个水桶（假设为圆柱形吧，但这不重要），水桶在绕自己的中轴线匀角速度旋转，直到桶内的水也随着匀角速度旋转（即水与水桶相对静止），此时水的液面形状是凹的，我们来看看该液面的形状。
牛顿的水桶
要分析形状，我们还要回顾之前提到过的流体静力学平衡：
http://kexue.fm/archives/1964/
31
Aug
暑假结束了，上学去~
By 苏剑林 | 2013-08-31 | 23003位读者 | 引用一个多月的暑假已经结束了，又回到了学校来。准确地说，昨天已经来到了学校，只是着搞卫生、社团等工作，无暇到blog上写点什么。早晨起来，一时无聊，就随便唠叨几句。
暑假就这样过去了，这也意味着大一完全过去了，我已经成为了师兄。曾不止一次感叹“光阴似箭，日月如梭”，而我越发地体味到这一点。不少人到了大学之后才明白高中生活的美好，而我有点不同，我在高中已经懂得大学并没有我们想象中的完美，所以我对大学和高中都抱有同样的眷恋和期待。大一过去了，从外边看来，我唯一的变化就是瘦了，沧桑了吧。还记得时隔一年的体检，我的体重居然少了十斤，以至于让我不得不怀疑那个秤的准确性；还记得多少次被小孩子喊做“叔叔”，被师兄称作“师兄”......
27
Sep
数学基本技艺之23、24（下）
By 苏剑林 | 2013-09-27 | 25399位读者 | 引用在上一篇文章中我们得到了第23题的解，本来想接着类似地求第24题，但是看着23题的答案，又好像发现了一些新的东西，故没有继续写下去。等到今天在课堂上花了一节课研究了一下之后，得到了关于这种拟齐次微分方程的一些新的结果，遂另开一篇新文章，与大家分享。
一、特殊拟齐次微分方程的通解
在上一篇文章中，我们求出了拟齐次微分方程$\frac{dy}{dx}=x+\frac{x^3}{y}$的解：
$$(2y+x^2)(x^2-y)^2=C$$
或者写成这样的形式：
$$(y+\frac{1}{2} x^2)(y-x^2)^2=C$$
5
Jan
不确定性原理的矩阵形式
By 苏剑林 | 2014-01-05 | 44055位读者 | 引用作为量子理论的一个重要定理，不确定性原理总是伴随着物理意义出现的，但是从数学的角度来讲，把不确定性原理的数学形式抽象出来，有助于我们发现更多领域的“不确定性原理”。
本文中，我们将谈及不确定性原理的n维矩阵形式。首先需要解释给大家的是，不确定性原理其实是关于“两个厄密算符与一个单位向量之间的一条不等式”。在量子力学中，厄密算符对应着无穷维的厄密矩阵；而所谓厄密矩阵，就是一个矩阵同时取共轭和转置之后，等于它自身。但是本文讨论一个更简单的情况，那就是n维实矩阵，n维实矩阵中的厄密矩阵就是我们所说的实对称矩阵了。
设$\boldsymbol{x}$是一个$n$维单位向量，即$|\boldsymbol{x}|=1$，而$\boldsymbol{A}$和$\boldsymbol{B}$是n阶实对称矩阵。在量子力学中，$\boldsymbol{x}$就是波函数，但是在这里，它只不过是一个单位实向量；并记$\boldsymbol{I}$是$n$阶单位阵。
考虑
$$\bar{A}=\boldsymbol{x}^{T}\boldsymbol{A}\boldsymbol{x},\bar{B}=\boldsymbol{x}^{T}\boldsymbol{B}\boldsymbol{x}$$
从这些记号可以看出，这些量对应着可观测量的期望值。当然，如果不懂量子力学，可以只看上面的矩阵形式。
在学车的时候，我堂大哥曾问我一道作圆的问题：
三圆的外切圆和内切圆 (1)
平面上给出三个两两相切的圆以及它们的圆心，求作一个圆与这三个圆都相切（尺规作图）。
如果从纯几何的途径入手，我们甚至很难判断这样的圆是否存在。但是我之前似乎已经看过类似的题目，于是很快想到一个名词：反演。反演可以将圆反演成直线（圆过反演点），也可以将圆反演成圆（圆不过反演点），而其他的相切、相交等关系保持不变。对反演后的图形进行相同的反演，就变回原来的图形。本题的难点在于圆太多，利用反演，我们可以将它变为两条直线和一个圆的问题。
假设读者已经有了反演的基本知识，如果没有，请到
http://zh.wikipedia.org/wiki/反演
阅读相关内容。
25
Mar
一本对称闯物理：相对论力学(二)
By 苏剑林 | 2014-03-25 | 18855位读者 | 引用从这个系列的第一篇文章到本文，已经隔了好多天。其实本文的内容是跟第一篇的内容同时完成的，为什么这么久才更新呢？原因有二，其一是随着春天的到来人也开始懒起来了，颓废呀～；其二，我在思考着规范变换的问题。按照朗道《场论》的逻辑，发展完质点力学理论后，下一步就是发展场论，诸如电磁场、引力场等。但是场论中有个让我比较困惑的东西，即场论存在着“规范不变性”。按照一般观点，我们是将规范不变性看作是电磁场方程的一个结果，即推导出电磁场的方程后，“发现”它具有规范不变性。但是如果用本文的方法，即假定场有这种对称性，然后就可以构建出场方程了。可是，为什么场存在着规范不变性，我还未能思考清楚。据我阅读到的资料来看，这个不变性似乎跟广义不变性有关（电磁场也是，这似乎说明即使在平直时空的电磁场理论中也暗示了广义不变性？）。还有，似乎这个不变性需要在量子场论中才能得到比较满意的解释，可是这样的话，就离我还很远了。
好吧，我们还是先回到相对论力学的推导中。
“无”中生有
上一篇文章我们已经构建了相对论力学的无穷小生成元，并进行了延拓。我已经说过，仅需要无穷小的变换形式，就可以构建出完成的相对论力学定律出来（当然这需要一些比较“显然”的假设）。这是个几乎从“无”到有的过程，也是本文标题的含义所在。另一方面，这种从局部到整体的可能性，也给我们带来一些启示：假如方法是普适的，那么可以由此构造出我们需要的物理定律来，包括电磁场、引力场方程等。（当然，我离这个目标还有点远。）
最近评论