图形技 Archives - 境由心生

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

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

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

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

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

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

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

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

^?View Code RUBY

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！”

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

境由心生

统计边，生活处

Category Archives: 图形技

Mandelbrot Set

爱的声音