当不会说谎的人去说谎

明天见分晓吧, 看我是否会更新这篇BLOG. 09/01/2010

拿去吧, 大叔大妈们的水手服!

今天下了动新84期的MV, 发现里面有武道馆演唱会版的Lucky Star OP 拿去吧, 水手服
好欢快滴说…
于是我很手贱的把音轨替换了成了 JAMがもってった! セーラーふく
违和啊…
其实我主要是想试一下CUDA编码H.264, 但是貌似失败了, 速度和CPU速度不应该没差别吧, 囧

附原版

PS: 的确, 论高清的效果, Vimeo比Youtube好很多, 2个视频都是2800K码率的, 但是V的颜色更加饱满, 下面的字幕也更加清晰

密码保护：嗯 酒 黑历史

这是一篇受密码保护的文章。您需要提供访问密码：

密码：

LILIUM

今天看了 空之音 第2话, 由于第1话是没有OP的, 所以当看到OP动画的时候, 我猛地一惊, 以第1话来看, 这应该是一部有卖萌+治愈的动画吧, 但是这个画面, 怎么如此的诡异 — 仿佛暗示着死亡, 完全不搭界的样子, 而这种风格的OP, 似乎勾起了我某个回忆

第一部把我带入”宅”圈的动画, 是妖精的旋律, 15禁, 至于我看的年纪, 大概也是15岁吧, 这部作品所给我带来的影响我现在也没发说清楚, 但是可以说我心理的黑暗面, 有很可观的一部分是这部作品带来的

妖精旋律的OP是LILIUM, 画面灵感来自Gustav Klimt的作品(Wiki):

Gustav Klimt画作特色在于特殊的象征式装饰花纹, 并在画作中大量使用性爱主题. 喜欢他画作的人称赞有其”沉闷美感”与大胆象征寓意.
Klimt画作另一特色为画中主角大部分都是女人, 主题则为”爱”, “性”, “生”与”死”的轮回宿命.

而空之音的OP, 画面依然是Klimt的作品的风格, 这自然不是巧合, 因为这2部作品的监督都是神戸守, 我其实还看过他的一部”鬼公子炎魔”, 也是黑暗系的, 他最近的一部作品除了空之音外, 是去年9月的新番信蜂, 考拉评价这部作品很高, 这样我也想看一下了, 等完结吧

附: 空之音OP动画 + LILIUM音乐

正多面体共有5个 的 2种证明

今天BBC的一个节目提到了这一点, 觉得蛮有趣的, 于是找了一下如何证明这一点, 记录下来

1. 一个初等证明(不利用欧拉公式, 出自一則關於正多面體祇有五種的簡單證明)

先想以下问题: 现有一纸张ABCDEF, 我们希望能将它折成一个立体图形的一部分, 使G成为该部分的顶点. 要做到这点, 必须将这张纸的一部分(即ΔDEG)切去, 然后才可把这张纸折成所需的立体, 而且没有重叠的地方.

如果将这个问题倒过来说, 一个立体图形的任何一个顶点的四周的面, 当它们一块块的互相拼接在一起(就好像将这个顶点压平一样), 都会成为一个有缺口的平面.

故此, 这些面在这顶点的角的和必少于360°

现在谈谈有关正多面体, 我们先假设某个多面体每一顶点是由P个正n边形连接着.

正n边形内角 = (180°n -360°)/n = 180° – 360°/n

由于每个顶点由P个正n边形连接, 将它们一块块地拼在一起, 由上述的讨论中得出:
P(180° – 360°/n) < 360°
P(1 – 2/n) < 2
nP − 2P < 2n
nP − 2P − 2n + 4 < 4
(n − 2)(P − 2) < 4

(P, n)都是大于2的整数, 上式子只有5组解 — (3,3),(3,4),(3,5),(4,3),(5,3), 分别对应5个正多面体 — 正4,6,12,8,20面体

2. 利用欧拉公式的证明

欧拉公式: 对于任意多面体(即各面都是平面多边形并且没有洞的立体), F, E 和V分别表示面, 边 ,顶点的个数, 那么F-E+V=2.

首先证一下吧, 证法蛮多的, 比如这里就提供了多达19种的证明, 我最早看到这个定理貌似是在 从一到无穷大 的书上, 印象中它里面提到的应该是柯西的证法, 挖掉一个面, 展开, 然后分割为三角形…
具体如下:

1. 把多面体(图①)看成表面是薄橡皮的中空立体.
2. 去掉多面体的一个面, 就可以完全拉开铺在平面上而得到一个平面中的直线形, 像图中②的样子. 假设F′, E′和V′分别表示这个平面图形的(简单)多边形、边和顶点的个数, 我们只须证明F′-E′+V′=1.
3. 对于这个平面图形, 进行三角形分割, 也就是说, 对于还不是三角形的多边形陆续引进对角线, 一直到成为一些三角形为止, 像图中③的样子. 每引进一条对角线, F′和E′各增加1, 而V′却不变, 所以F′-E′+V′不变. 因此当完全分割成三角形的时候, F′-E′+V′的值仍然没有变. 有些三角形有一边或两边在平面图形的边界上.
4. 如果某一个三角形有一边在边界上, 例如图④中的△ABC, 去掉这个三角形的不属于其他三角形的边, 即AC, 这样也就去掉了△ABC. 这样F′和E′各减去1而V′不变, 所以F′-E′+V′也没有变.
5. 如果某一个三角形有二边在边界上, 例如图⑤中的△DEF, 去掉这个三角形的不属于其他三角形的边, 即DF和EF, 这样就去掉△DEF. 这样F′减去1, E′减去2, V′减去1, 因此F′-E′+V′仍没有变.
6. 这样继续进行, 直到只剩下一个三角形为止, 像图中⑥的样子. 这时F′=1, E′=3, V′=3, 因此F′-E′+V′=1-3+3=1.
7. 因为原来图形是连在一起的, 中间引进的各种变化也不破坏这事实, 因此最后图形还是连在一起的, 所以最后不会是分散在向外的几个三角形, 像图中⑦那样.
8. 如果最后是像图中⑧的样子, 我们可以去掉其中的一个三角形, 也就是去掉1个三角形, 3个边和2个顶点. 因此F′-E′+V′仍然没有变.

即F′-E′+V′=1 成立, 于是欧拉公式: F-E+V=2得证.

对于正多面体, 假设它的各面都是正n边形, 而且每一个顶点处有r个边相遇。这样就有:
nF=2E  — (1)
rV=2E  — (2)
(1)的右边系数2是因为每边出现在2面中, (2)的右边系数2是因为每边通过2个顶点。把(1)和(2)代入欧拉公式中, 就得到：

即

显然n≥3, r≥3, 因为多边形至少有三边, 而在每顶点处也至少有三边. 但n＞3, 且r＞3又是不可能的, 因为那样就要有

可是E＞0. 所以r和n中至少有一个等于3.
设n=3,

因此r=3,4,5,由是E=6, 12, 30, 而F=4, 8, 20, 这就给出了正四面体, 正八面体和正二十面体.
设r=3,

因此n=3,4,5, 由是E=6, 12, 30, 而F=4, 6, 12, 这就给出了正四面体, 正六面体和正十二面体.

星之所在

星の在り処
星之所在

歌：う～み
コーラス：伊藤和子
作词：浜田英明
编曲：和田耕平

君 の 影 星 の ように 朝 に 溶 け て 消 え て い く
你的影子像星星一样 随着黎明的到来而消逝
kimi no kage hoshi no you ni asa ni to ke te ki e te i ku

行 き 先 を 失 く し た ま ま 想い は 溢 れ て く る
失去目的地 思念溢满心头
iki sa ki wo na ku shi ta ma ma omoi wa afu re te ku ru

强さにも弱さにも この心は向き合えた
强大也好软弱也好 迎向这颗心
tsuyo sa ni mo yowa sa ni mo ko no kokoro wa mu ki a e ta

君と な ら ど ん な 明 日 が 来 て も 怖 く な いの に
和你在一起的话 无论明天会怎样都无所畏惧
kimi to na ra don na ashita ga ki te mo kuwa ku na i no ni

二 人 歩 い た 时 を 信 じて い て 欲 し い
二人同行之时 希望你能相信我
futa ri a ru i ta toki wo shinjite i te ho shi i

真 実 も 嘘 も な く 夜 が 明 け て 朝 が 来 る
不是真实也不是谎言 夜晚过去早晨到来
shin jitsu mo uso mo na ku yuru ga a ke te asa ga ku ru

星 空 が 朝 に 溶 け て も 君 の 辉 き は わ か る よ
星空溶入朝霞 知晓了你的光辉
hoshisora ga asa ni to ke te mo kimi no kagaya ki wa wa ka ru yo

思 い出 を 羽 ば た か せ 君 の 空 へ 舞 い 上 が る
张开思念的羽翼 飞上你的天空
omoi de wo ha ba ta ka se kimi no sora e ma i a ga ru

星 空 が 朝 に 溶 け て も 君 の 辉 き は わ かる よ
星空溶入朝霞 知晓了你的光辉
hoshi sora ga asa ni to ke te mo kimi no kagaya ki wa wa ka ru yo

爱してる た だ そ れ だ け で 二人 は い つか ま た 会 え る
爱你 只要爱你 二人就一定会再次相见
aishiteru ta da so re da ge de futari wa i tsuka ma ta a e ru

—

流水账似的回忆并非我的性格呢, 罗列一本本书, 一场场电影之类之类, 有些无聊呢
所以呢, 作为2009年的结束, 只是贴一首玩过的游戏的主题曲吧
愿自己能以美好的心境活着

跟客人接触的时候 那里就是决斗的场所 绝对不能动怒 不能轻敌 看穿对手的心 调动四肢 全身贯注 使出全部的技巧 还有不能缺少笑容和气度 不管发生什么 别忘了扑克脸

—黑羽盗一 柯南TV 219话

猛然记起这么一段话, 在昨晚玩21点的时候想到的, Poker face, 哈~

结束了哦…

猛然滴发现, 自己MSN签名还是3个月前, 刚开学时的”开始了哦…” 但是这个学期就这样结束了 没有一点征兆的 当我还在等待数值方法什么时候涉及到微分方程 或者统计分布什么时候讲极限理论 或者回归模型什么时候讲非线性模型的时候 居然发现这个学期已经结束了 随之而来圣诞的假期 1月10号开始的考试 然后便是下个学期了
很欣喜地, 我在这个假期里可以”享受”一下专心读书的乐趣了, 因为新房子在1月9号之前不会有互联网
同样欣喜滴, 达达和夫人会过来玩
虽然前段时间关于房子有点忧郁, 过去了都…

B’z 衝動

日文

长い夜に目覚めて 青白い部屋の中
不意に大事な何かを 伤つけたくなる
暗暗のfreedom
扉の前で立ち尽くす
ちっぽけな背中に気づいて欲しい
仆にも谁かを爱せると その手を重ねて知らせて
あなたの温もりがくれる 冲动
真実ばかり追いかけて 揺れ动くeveryday
（きりがない）
强がる自分の中身も 丸っきり谜だらけ
燃え上がるfreedom
钢のような绊だけが へこんでく心を踏ん张らせる
仆にも何かを変えられる さりげない言叶で嗫いて
あなたの声が明日への 冲动
谁もが无限の可能性を
抱きしめて生まれて来たんでしょうね
仆にも谁かを爱せると その手を重ねて知らせて
希望とは目の前にある道
何処かに行けると信じよう
あなたの全てが仆の 冲动
go for it go for it 爱情こそが冲动

罗马拼音

Nagai yoru ni mezamete aojiroi heya no naka
Fui ni daiji na nanika wo kizutsuketaku naru
(Kurayami no freedom)

Tobira no mae de tachi tsukusu
Chippoke na senaka ni kiduite hoshii

Boku ni mo dareka wo aiseru to
Sono te wo kasanete shirasete
Anata no nukumori ga kureru shoudou

Shinjitsu bakari oikakete
Yure ugoku everyday
Tsuyogaru jibun no nakami mo
Marukkiri nazo darake
(Moeagaru freedom)

Hagane no you na kizuna dake ga
Hekondeku kokoro wo funbaraseru

Boku nimo nanika wo kaerareru
Sarigenai kotoba de sasayaite
Anata no koe ga ashita he no shoudou

Daremo ga mugen no kanousei wo
Dakishimete umarete kitan deshou? Nee

Boku nimo dareka wo aiseru to
Sono te wo kasanete shirasete
Kibou to wa menomae ni aru michi
Dokoka ni ikeru to shinji you
Anata no subete ga boku no shoudou

go for it go for it
Aijou koso ga shoudou

中文

漫漫长夜中醒来 苍白空漠的房间里
突然想要破坏重要的什么
黑暗的freedom
久久伫立在门前 渴望被留意到这渺小的背影
请握住我的手让我明了 我也可以爱恋上谁
你的温暖给我去爱恋的冲动
一味追逐着真实 动荡摇摆的everyday
永无止境
故作坚强的心底 其实充满困惑
熊熊燃烧的freedom
钢铁般的情思 令我萎靡屈服的心振作起来
请平静地向我悄声诉说 我也可以改变什么
你的声音给我奔向明天的冲动
任谁都是紧拥着无限的可能降生的 对吧
请握住我的手让我明了 我也可以爱恋上谁
所谓希望既是呈现眼前的远途 要坚信可以通向何方吧
你的全部给我前进的冲动
go for it go for it
正是爱情给我冲动
go for it go for it
正是爱情给我冲动

我不是一个情绪化的人, 但是我要活得真实, 笑~

rpy2 – 结合Python与R

Python被称作胶水语言, 她和N多语言都能结合的很好, 在Python中引入R使用的模块是rpy2, 可以从http://rpy.sourceforge.net/下载
安装后尝试
import rpy2.robjects as robjects
如果提示找不到R.exe, 需要在Windows环境变量中添加R_HOME, 值为R的安装目录, 我的是C:\Program Files (x86)\R\R-2.9.2
如果提示没有win32api模块, 需要安装pywin32
大致上, 有2种调用R的方法
其一是直接使用R代码, 比如随机选取10个标准正态分布的变量
R Code
rnorm(10)
Python

1. import rpy2.robjects as robjects
2. r = robjects.r
3. y = r(‘rnorm(10)’)

注意这里得到的结果总会是一个RVector(尽管有时候只有一个值), 有点类似于python的list, 她当然也可以转换为list
其二是调用R的函数, 依然是上面的例子
Python

1. import rpy2.robjects as robjects
2. r = robjects.r
3. y = r.rnorm(10)

下面的程序, 以标准正态分布为例, 演示了用小样本方差去估计总体方差时, 虽然使用Bessel’s correction后的无偏估计的确在均值方面更加接近真实值, 但是用通常的距方法得到的有偏估计, 却可以提供更小的MSE, Wiki中有更加精确的数学表达

1. import rpy2.robjects as robjects
2. r = robjects.r
3. ss = r.rnorm
4.  
5. def test(n):
6.     y = ss(n)
7.     avg = sum(y)/n
8.     s=0
9.     for i in y: s=s+(i-avg)**2
10.     v1=s/n
11.     v2=s/(n-1)
12.     return v1,v2
13.  
14. def test0(m,n):
15.     v1,v2=[],[]
16.     for i in range(m):
17.         v1.append(test(n)[0])
18.         v2.append(test(n)[1])
19.     avg1=sum(v1)/m
20.     avg2=sum(v2)/m
21.     s=0
22.     for i in v1: s=s+(i-avg1)**2
23.     s1=s/m
24.     s=0
25.     for i in v2: s=s+(i-avg2)**2
26.     s2=s/m
27.     return avg1,avg2,s1,s2
28.  
29. t1,t2=0,0
30. for k in range(10):
31.     v=test0(100,5)
32.     if abs(1-v[0])<abs(1-v[1]): t1=t1+1
33.     if v[2]<v[3]: t2=t2+1
34.  
35. print "有%d次,有偏估计的期望更接近真实值" % t1
36. print "有%d次,有偏估计的方差更小" % t2

Google API for LaTeX to Image

Here is a Google Chart API which is used to render your LaTeX code to an image.
http://chart.apis.google.com/chart?cht=tx&chs=1×0&chf=bg,s,FFFFFF00&chco=000000&chl=$$$e^{i\pi}%2B1%3D0$$$
Use the HTML IMG tag wraping above url and it looks

For someone who doesn’t use WordPress XD

不相关且不独立的2组随机变量的构造 连续&离散

今天在课本上瞄到这样一句话, 随手写了2组均为(0,1)分布的随机变量, 发现如果不相关的话一定就是相互独立的, 于是我就纠结了, 后来才知道, 对于只有2个取值的2组随机变量, 不相关和独立是等价的, 但是通常情况下: 独立=>不相关 但是 不相关 =/=> 独立, 下面是一些例子

连续型:
设X为均值为0的服从正态分布的随机变量
设Y=X^2
由Cov(X,Y)=E(XY)-E(X)E(Y)=E(X^3)-0=E(X^3)=0 可知 X, Y不相关
但是
f(y|x) =/= f(y)
所以X, Y也不独立
类似的例子WIKI上也有: Normally distributed and uncorrelated does not imply independent

离散型:
这个就有些复杂了, 刚才说过, 如果随机变量必须有3种或者以上的取值才能构造出来的, 下面是一个例子

E(X)=3.8
E(Y)=0.6
E(X)E(Y)=2.28
E(XY)=2.28
所以 不相关
P(X=5|Y=2)=0.0845/0.4=0.2124 =/= P(X=5)=0.2
所以 不独立
有一篇文章专门介绍这种离散型随机变量的构造, 文章在这里不相关且不独立二维离散型随机向量的构造方法
文中附带了一个生成程序, 但是是一个试卷生成软件用的代码, 也不复杂, 下面是我改写的Python版

1. from numpy import array,dot,min,mat
2. import random
3.  
4. def gcd(a,b):
5.     while b > 0: a,b = b, a%b
6.     return a
7.    
8. x2=random.randint(1,3)
9. x1=random.randint(3,5)
10. x3=random.randint(3,5)
11. q=gcd(gcd(x3-x2,x1-x3),gcd(x1-x3,x2-x1))
12. y1=random.randint(1,3)
13. y2=random.randint(-1,1)
14. y3=random.randint(-1,1)
15. r=gcd(gcd(y3-y2,y1-y3),gcd(y1-y3,y2-y1))
16. t=random.randint(1,2)
17. s=random.randint(1,2)
18. a1=random.randint(1,2)
19. a2=random.randint(1,2)
20. a3=5-a1-a2
21. xi=0.2*array([[a1,a2,a3]])
22. xi=mat(xi)
23. b1=random.randint(1 ,2)
24. b2=random.randint(1 ,2)
25. b3=5-b1-b2
26. eta=0.2*array([[b1],[b2],[b3]])
27. eta=mat(eta)
28. X=array([x1,x2,x3])
29. Y=array([y1,y2,y3])
30. H=array([[(y3-y2)/r,(y1-y3)/r,(y2-y1)/r],\
31.         [t,0.0,-t]])
32. G=array([[s,(x3-x2)/q],\
33.         [-s,(x1-x3)/q],\
34.         [0.0,(x2-x1)/q]])
35. GH=dot(G,H)
36. etaxi=dot(eta,xi)
37. c=min(GH)
38. d=min(etaxi)
39. b=-d/c
40. A=b*GH+etaxi
41. print A
42. print x1,x2,x3
43. print y1,y2,y3
44. raw_input("\nXD")
45.