Mandelbrot Set

人懒的时候容易被一些事情所烦,需要找些事情放松放松。看看图形是个不错的选择,比如经典的Mandelbrot集。虽然R画的没有绚烂,但是还能看。下面是搜集的稍作修改的一段代码和图形。

1.ggplot2的图形,循环迭代较多,运行确实很慢。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
library(ggplot2)
max_iter = 25
cl = colours()
step = seq(-2, 0.8, by = 0.005)
points = array(0, dim = c(length(step)^2, 3))
t = 0
for (a in step) {
    for (b in step + 0.6) {
        x = 0
        y = 0
        n = 0
        dist = 0
        while (n < max_iter & dist < 4) {
            n = n + 1
            newx = a + x^2 - y^2
            newy = b + 2 * x * y
            dist = newx^2 + newy^2
            x = newx
            y = newy
        }
 
        if (dist < 4) {
            color = 24
        }
        else {
            color = n * floor(length(cl)/max_iter)
        }
 
        t = t + 1
        points[t, ] = c(a, b, color)
    }
}
df = as.data.frame(points)
ggplot(data = df, aes(V1, V2, color = cl[V3])) + geom_point(size = 0.5) +
    opts(panel.background = theme_blank(), panel.grid.major = theme_blank(),
        panel.grid.minor = theme_blank(), axis.ticks = theme_blank(),
        axis.text.x = theme_blank(), axis.text.y = theme_blank(),
        axis.title.x = theme_blank(), axis.title.y = theme_blank(),
        legend.position = "none")

渲染的图形还是非常漂亮的,如下(点击图形看大图,否则看不出颜色渲染的效果

2. 动态看图形变化,此转于此处,代码简洁易理解。不过代码做出来的图与原网页展示的图形不同,颜色怎么调都没有他们的漂亮,主要是迭代次数还不够,也不够精细,直接看原链接的吧,搭配颜色是个细致活,伤脑细胞。

爱的声音

今天是母亲节。我是记不住这个节日的,可能本身这个节日在并不算做传统节日的缘故吧。长这么大,我一共陪母亲度过两个母亲节。一个是高三的时候,一个高四的时候。

那两年高考的最后几个月母亲过来陪读,很辛苦,很煎熬。但那是我第一次有了那么多时间与母亲交流,也是从那时候起,我才开始懂事,开始关心母亲,懂得一个小小的表达母亲都会很感动。每天的散步闲谈,让我深深体会到了父母对我的无条件爱,母亲说的有句话我一直刻在在心底“:只要你好,不管有没有出息,过的健康快乐,我和你爸就什么都好,都安心。”每次我一想起这句话,生活学习中的压力顿时烟消云散,母亲诚挚的话语,让我懂得了生活意义。虽然我无数次的质疑生活的意义,但是我从未放弃母亲告诉我的两样东西:爱与理想。

还有个场景我一回想就会热泪盈眶:高考前夕,烈日炎炎下,一群家长焦急等待自己孩子下课,手里都提着自己亲手做的可口饭菜,母亲因为身材娇小,而且怕热,每次都不得不站在人群后面,然后踮着脚,仰着头,远远的眺望,满脸的焦急与期盼。一想起这个场景,我心里都发酸,母亲为我倾注了那么多爱,而我却回报的如此至少,内心充满了愧疚与感动。

一晃,三年就过去了,如此之快。记得第一次我给母亲送了百合花,愿她健康快乐;第二次我送给她一只小乌龟,祝愿她长命百岁。我是一个不善于表达情感的人,满肚子的话,到了嘴边都会被憋回来。那两次送礼物的时候,我只是傻傻的笑,连句祝福不敢说,男孩子的倔强让我说不出那些肉麻的话。至今我也没对母亲说一句“妈,我爱你!”我知道,没有这些话母亲都很感动了,如果说了,怕我们两个都会哭出来,都是天蝎座的人,感情很丰富,很敏感,很脆弱,受不住这样的爱的表达。

三年的每个这个星期,虽然没法再送个我精心准备的礼物,但我都会或早或晚的打电话回去问候下母亲,今天也是一样。不过但是今天看到这个网页时候,我想到了一种表达爱的方式。我要一种比较特别方式送给母亲我一直不敢说的那几个字。不多言说,附图附代码。

1
2
3
4
5
6
7
8
9
10
11
12
library(seewave)
library(tuneR)
woaini = readMP3("我爱你.mp3")
woaini1=resamp(woaini, g = 2000)
muqinjie = readMP3("母亲节快乐.mp3")
muqinjie1 = resamp(muqinjie, g = 2000)
iloveyou = readMP3("i_love_you.mp3")
iloveyou1 = resamp(iloveyou, g = 2000)
op = par(mfrow = c(3,1),mar = c(4.5,4,2,2))
oscillo(woaini1, f = 2000, from = 3, to = 6,cexlab = 0.75,colwave = "red",title = "妈-我爱你")
oscillo(muqinjie1, f = 2000, from = 2, to = 5,cexlab = 0.75,colwave="gold",title = "妈-母亲节快乐")
oscillo(iloveyou1, f = 2000, from = 0, to = 4,cexlab = 0.75,colwave="lightblue",title = "i love you so much")

图解:

  • 最上是我用家乡话(信阳罗山话)说的:“妈,我爱你!”(罗山话中,“你”发音类似于“恩”,故振幅小)
  • 中间一幅依然用家乡话:“妈,母亲节快乐!”(录得不好,貌似有杂音)
  • 最下面的那副是英文:“I love you so much!”

其实从图的走势还是可以看出我羞涩的内心,爱的表达,对我来说,竟然真不是一件容易的事情,但是对许多男生来说,又何尝不是呢?男人有着内心的倔强,都懂得~