Mar 17

    显然,过 Pizza 的圆心作四条直线,把一个周角平分成八等份,则整个 Pizza 饼也被分成了八等份。我们也很容易联想到,如果过圆心外的一点做出四条直线,并且同样满足每两条相邻直线夹 45 度角,那么这八块 Pizza 饼显然是不一样大的。考验你直觉的时候到了:你认为蓝色面积之和与红色面积之和相比,哪个大一些呢?

  

查看更多 »

Mar 8

    Mathematica 提供了一个看上去毫无用途的无厘头函数 Rasterize ,它可以以图片的格式输出运算结果。比如,下面这个句子可以打印出 (x+1)^n 的展开式的“倒影”:

   

    今天我突然想到,我们可以利用这个函数很方便地分析汉字在图象上的性质。函数 Binarize 可以把图象转换为单色单通道, ImageData 则可以把图象转换成数组的形式,以便我们定量分析。因此,下面这句话就可以把一个汉字转换成 12*12 的 01 矩阵:

   

查看更多 »

Mar 8

 

来源:MathOverflow
不得不说,确实很妙!

Feb 12

    这是一个经典智力问题,不知道大家见过没。下图是一辆自行车在泥地中驶过留下的痕迹,你能据此判断出这辆自行车是从左往右行驶的还是从右往左行驶的吗?

   

    提示:题目条件是充分的,根据这两道车轮印我们足以判定车行方向。这和图中的线条粗细、边缘锯齿没有关系,你完全可以把两道痕迹当作没有粗细之分的理想曲线;为了解决这个问题,必须仔细分析自行车驶过后两道车轮印一定会满足的几何性质。

查看更多 »

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次多项式。可以看到,当次数变高时,插值多项式反而变得更不准确。

查看更多 »

Oct 20

    连续函数f(x)满足f(0)=0且f(1)=0。证明,总能在[0,1]中找到两个数a和b满足b-a=1/2且f(a)=f(b)。换句话说,我们总能画出一条长为1/2的水平线段,它的两个端点都在函数f(x)上。
    这个证明再次用到了我们上次提及的零点定理。考虑f(1/2)的值,如果它也等于0,我们的问题就直接解决了。无妨设f(1/2)>0,那么考虑f(x+1/2)-f(x)的值:当x=0时,该值为一个正数;但当x=1/2时,这个值变成了一个负数。这表明,在x从0增长到1/2的过程中,一定有某一刻使得f(x+1/2)-f(x)恰好为0。

    我们接下来的问题是,除了长为1/2的横线段始终存在以外,还有哪些长度值具有相同的性质?下面我们证明,对任意一个正整数n,长为1/n的横线段也总是存在的。

查看更多 »

« 更早的日志