Feb 20
奇妙的心电图数列
icon1 Matrix67 |icon2 Brain Storm | icon4 2010-02-20 0:06 | icon321 Comments »

    心电图数列 (EKG Sequence) 的定义简单而有趣:第一项为 1 ,第二项为 2 ,以后的每一项都是最小的和前一项不互质并且不曾出现过的数。换句话说,数列 a(1)=1 , a(2)=2 ,且当 n>2 时取 a(n) 为所有满足以下两个条件的数中最小的那一个:该数与 a(n-1) 有大于 1 的公约数,并且该数与前面 n-1 项都不相等。心电图数列的前面 20 项为

      1, 2, 4, 6, 3, 9, 12, 8, 10, 5, 15, 18, 14, 7, 21, 24, 16, 20, 22, 11 ...

    为什么它叫做心电图数列呢?原因很简单——因为把它描绘在图象上时,看上去像一张心电图。

 

查看更多 »

Feb 11

刚才在这里看到了如题所说的图像,立即想到用 Mathematica 验证一下。我选出了几个个人比较感兴趣的 k ,再用一句话便可输出所有对应 k 的图像:

kArray = {2, 3, 4, 6, 8, 10, 12, 14, 16, 18, 20, 36, 50};
For[i = 1, i <= Length[kArray], i++,
 Export["F:\\" <> ToString[kArray[[i]]] <> ".png",
  ArrayPlot[Table[Boole[Length[Divisors[x*y]] == kArray[[i]]], {x, 1, 400}, {y, 1, 400}],
   PixelConstrained -> {1, 1}, Frame -> False]]];

 
当 k=2 时,由于只有素数才有两个约数,因此所有点都是形如 (p, 1) 或者 (1, p) 的点,其中 p 为某个素数:

查看更多 »

Dec 9

    Dan Christensen发现,把所有次数不超过5的、系数在-4到4范围内的整系数多项式的所有根描绘在同一个复平面上,你会看到一个异常壮观的画面。图中的每个灰色点代表某个二次多项式的一个根,蓝色点代表三次多项式的根,红色代表四次多项式的根,黑色代表五次多项式的根。水平线代表实轴,0和±1的地方有很明显的空洞;竖直方向是虚轴,每个单位根处也都有明显可辨的空洞。

查看更多 »

Dec 6

    今天学到了一个新的名词,Runge现象。1901年,Carl David Tolmé Runge意外地发现,用差值插值多项式逼近函数f(x)=1/(1+25x^2)时出现了一些反常的现象。如图,灰色的粗线就是Runge函数在[-1,1]上的图象。蓝色虚线是过[-1,1]上的6个等距点所得到的5次多项式,红色虚线是过[-1,1]上的10个等距点所得到的9次多项式。可以看到,当次数变高时,插值多项式反而变得更不准确。

查看更多 »

Nov 14

    你相信吗?对于任意一个可数集,总能找出不可数个子集,使得从中任取两个集合,其中一个都是另一个的真子集。乍看之下,这似乎是不可能的。如果任两个集合之间都具有“其中一个是另一个的真子集”的关系,那它们就能构成一个“集合序列”(准确地说是全序关系),使得每个集合都是由它前面那个集合添加进若干元素得到;换句话说,我们能通过不断往一个空集中添加新的元素依次得到所有这些集合。但是如果这些集合中的元素就只有可数个,那这个“集合序列”中怎么会有不可数个集合呢?然而,涉及到无穷的问题总是那样违背直觉。下面我们只用三行字就能说明,这个命题的的确确是成立的。

查看更多 »

Nov 6

    重要通告:最近多次发现我的tom邮箱发出的邮件被识别成了垃圾邮件,是什么原因我还不是很清楚。最近向我的tom邮箱发过邮件但迟迟没有收到回复的朋友麻烦检查一下垃圾邮件箱,或者重新给我发一次邮件,我换一个邮箱回复您。

    数学学习真正悲哀的就是,记住了某个神奇而伟大的定理,看懂了其最严密的推导过程,但却始终没能直观地去理解它。虽然严密的推导是必要的,直观理解往往是不准确的,但如果能悟出一个让定理一瞬间变得很显然的解释,这不但是一件很酷的事,而且对定理更透彻的理解和更熟练的运用也很有帮助。我惊奇地发现,国内的每一本高数课本上都严格地讲解了微积分基本定理的证明,但几乎没有任何一个课本上讲过积分等于函数下方的图形面积究竟是为什么。事实上,这几乎是显然的,但还是有不少人学完微积分后仍然没有意识到。每当谈到这个问题时,我更愿意首先提出一个非常有启发性的事实——圆的周长是2·pi·r,圆的面积就是pi·r^2,后者的导数正好就是前者。这个现象是很容易理解的,因为圆的半径每增加一点,面积增加的就是周长那么一圈,换句话说面积的变化就等于周长。类似地,如果你能找到一个函数g(x),它的导数正好就是f(x),那么当x每增加一点,g(x)就增加了一条小竖线段,显然g(x)就应当是f(x)下方的面积。看清了这一点之后,我们才能欣赏到微积分基本定理真正牛B的地方。原先大家都是用分割求极限的办法来求函数下方的面积,但Leibniz却把面积看作一个可变的整体,用一种办法“一下子”就把它求了出来。有趣的是,这种现在看来如此自然的神奇办法,一千多年来居然没有任何人想到。

查看更多 »

Nov 1

    Banach-Tarski悖论指出,你可以把一个三维球体分割成有限多份,然后拼合成两个和原来一模一样的球体。这个构造是Stefan Banach和Alfred Tarski在1924年发表的论文中给出的,不过我还从来没有完整地瞻仰过这个牛B的构造过程。今天我看到了一个Banach-Tarski悖论的弱化版,但它的反直觉性绝对不亚于Banach-Tarski悖论。通过这个弱化的结论,你或许会对Banach和Tarski的构造方法有了更多的理解。
    下面我将会给出这样一个神奇的构造:取出[0,2]的一个子集S,把它分割为可数个不相交的点集,对每个点集各自进行适当的平移后,可以让它们的并集变为全体实数集。

查看更多 »

Oct 27

 

    有一个放听装可乐的货架,它的宽度要比四听可乐的直径稍微大一些。把10听可乐放进这个货架里,堆叠成一个三角形。虽然底下三层可乐罐歪歪斜斜有高有低,但最顶上的那听可乐一定位于货架的正中心,也就是说它到货架两壁的距离是相等的。这是为什么呢?

查看更多 »

« 更早的日志