Monte Carlo方法

The Review of Particle Physics有一篇关于Monte Carlo方法的Review，在Mathematical Tools里面，介绍了Monte Carlo方法和任意分布随机数的生成，地址是 http://pdg.lbl.gov/2005/reviews/monterpp.pdf。 另外还有一些关于概率和统计的Review，对进行计算机模拟实验很有帮助。

谈谈Super64吧

生日礼物耶！我不太舍得买，但是很仰慕这把琴（或者说比较仰慕Franz Chmel），经常念叨，结果老婆就作为生日礼物送给我，很惊喜啊。现在把玩了有一段时间了，谈谈感受吧。

• 体积：

  的确比12孔口琴大了不少，也重了一些，但是在乐器中仍然属于非常小巧的，所以没什么抱怨了。前些时候跟一个邻居介绍口琴有三到四个八度的音域，他还很吃惊。他虽然不精通音乐，但是也玩过很多乐器，所以对乐器有些了解，知道一般管乐器没有这么大的音域。这一点小小口琴的确挺让人吃惊的。

• 吹嘴：

  非常顺滑。一开始感觉比270吹嘴厚，但是习惯之后觉得比270更容易含多孔演奏八度和声，大概是因为形状更符合人体工程学吧。吹快速的曲子时比较容易感觉嘴唇粘连，似乎吹嘴有吸水效果似的，大概和镀银有关。但是如果稍微有点水，就会变得非常滑溜。至于圆形吹孔，我并没有感觉和方形的有太大不同，因为我很少真的把舌头贴在吹孔上的，除了极少数吹和声的时候。孔距比270略大一些，虽然差别不大，但是我需要花一段时间去适应。

• 盖板：

  以前觉得Super64外观不如270好看，但是真的拿到手，觉得还是很不错的。尤其是新琴，外壳很亮，造型很有现代感，这点和新的270Deluxe比较相似，和270的古典美不一样。拿给一位学小提琴的朋友看，他一个劲的夸漂亮。外壳似乎比较薄，捏起来像日本汽车外壳的感觉，比较软。难道是我手劲太大了？

• 按键：

  比270行程长，比255行程短。弹簧没有270那么硬，花了一段时间去适应。新琴按键声音比较大，拆开后在按键推杆边缘抹了一点点半固态的润唇膏，声音小了许多。这个用量一定要小，用手指头沿着边缘抹一点点就好，多了就容易沾灰了。

• 簧片：

  明显比270容易驱动，同样的力度变化，在Super64上显的更容易一些。270上发音略微滞后的现象没有了，用老婆的话说，就是听我吹曲子“不赶了”。高音比270省力，低音比270容易发声，多出的低音四孔需要点功力驱动。有少数几个簧片发音有问题。例如多出的低音四孔里有个别簧片有颤音，但是吹了一段时间后好了些。高音有些音发音不通畅，需要重新调整一下簧片和面板间距。有一个c音至今仍然有问题，吹起来似乎有关风薄膜被粘住那种“扑扑”声，但是明显不是薄膜的问题。快速或大力的吹没问题，但是温柔一点的吹问题就来了。我估计是簧片平直度的问题，但是因为在内侧，不太好调整，我也不希望把新琴弄坏，反正同样音高的c有好几个呢，就这样吧。

• 音色：

  和270比较，270明亮中略带沙哑，Super64深沉饱满，类似于单簧管。感觉270比较有激情，Super64婉约一些。这一点上不同的人可能有不同的喜好，我老婆喜欢Super64的音色。不管他人如何看法，我觉得琴格材料对音色影响还是满大的。

Fedora Core 5英文环境下运行中文Firefox的Script

这里假设是KDE环境，Gnome没有试过，不知道行不行。

#!/bin/bash

export LANG="zh_CN.UTF-8"
export GDM_LANG="zh_CN.UTF-8"
export QT_IM_MODULE=scim
export XMODIFIERS=@im=SCIM
export GTK_IM_MODULE=scim
firefox &

在eps图形中插入文字

% Sample of printing text

/Times-Roman findfont   % Get the basic font
20 scalefont            % Scale the font to 20 points
setfont                 % Make it the current font

newpath                 % Start a new path
72 72 moveto            % Lower left corner of text at (72, 72)
(Hello, world!) show    % Typeset "Hello, world!"

showpage

参考链接： First Guide to PostScript

把多个eps文件合为一个eps文件

有时候需要把多个eps图形文件合为一个文件。在网上搜索了一下，似乎没有什么特别合适的工具。这里有个解决办法。 在linux下eps文件其实就是一个文本文件，可以用vim等文本编辑器手动编辑。例如要把01.eps, 02.eps, 03.eps, 04.eps四个图形文件合为一个2x2的图形文件，可以首先执行下面命令把它们合为一个文件

cat 01.eps 02.eps 03.eps 04.eps > outfile.eps

然后用vim打开outfile.eps。各个文件中间部分是描述图形的，不需要管。只要改动一下每个文件的开头和结尾部分就可以了。

1. 每个eps文件开头都有类似于这样一行：

   %!PS-Adobe-2.0 EPSF-2.0

   只要保留第一个文件中的就行了，其余文件开头的这一行都删去。
2. 每个eps文件结尾都有这样的命令：

   showpage

   只需要保留最后一个就可以了，其余文件结尾处的该命令都删去。
3. 每个文件开头都有类似于这样一行：

   %%BoundingBox: 0 0 100 120

   这表示图形的大小，四个数字分别是左下角和右上角的坐标。我们需要根据自己的要求把数字改一下。例如如果每个eps大小都是100x120，我们希望得到2x2的图形列表，那么合成之后的图应该是200x240的。那么我们把第一个eps文件开头的BoundingBox改为

   %%BoundingBox: 0 0 200 240

   ，然后把其余各eps文件开头的BoundingBox行删掉。
4. 在第一个eps文件开始处（BoundingBox之下，其余postscript命令之上）加入

   save gsave

   ，在其余eps文件开始处加入

   grestore restore
   save gsave
   xxx yyy translate

   xxx和yyy对应于每个eps图形平移的坐标。在所有eps文件最后的showpage之前加入

   grestore restore

然后就可以用postscript浏览器看效果了。 参考链接： http://www.debian.org/doc/manuals/user/ch11.html

谢旭东的Blog

今天和好朋友谢旭东联系上了。好久没见面了。这是他的Blog: 谢旭东的Blog

设置Linux router

如果Linux机器上有两块网卡，eth0（外部网）和eth1（内部局域网），可以如下设置让内部局域网的机器通过这台电脑连接外部网络：

iptables -t nat --flush
iptables -t nat --delete-chain
iptables -t nat -A POSTROUTING -o eth0 -j MASQUERADE
iptables --insert FORWARD -i eth1 -j ACCEPT
echo 1 > /proc/sys/net/ipv4/ip_forward

用blender做的第一个3D图

就是为了好玩，挺简单的一个3D图。

生日快乐！

祝自己生日快乐！

收到了老婆送的礼物，Hohner Super64口琴，非常非常高兴！一直梦寐以求的琴啊！确实比270的音色好，簧片反应快，但是孔距和按键的行程和270都有些不一样，需要一些时间去适应。

Fedora Core 5安装flash插件

从Macromedia直接下载安装flash插件在Fedora Core 5里行不通，firefox总提示找不到插件。可以从下面这个网址直接下载安装flash-plugin的rpm包， http://macromedia.mplug.org/ 或者干脆按照网页所示的方法把Repos加入/etc/yum.repos.d/用yum安装。

把文本文件中的逗号分隔符换成TAB

sed 's/,/\t/g' input-file > output-file

用gnuplot拟合Gaussian曲线

在gnuplot里定义Gaussian函数

f(x)=a*exp(-(x-b)*(x-b)/(c*c))

这个Gaussian函数的宽度 FWHM = 2*sqrt(ln2)*c = 1.665c 建立一个参数列表文件gaussian.par：

#three parameters
a = a0
b = b0
c = c0

其中a0, b0, c0是人为设置的初始值。在gnuplot里运行

fit [xrange] f(x) 'datafile' via 'gaussian.par'

上面[xrange]替换成相应的自变量区间，datafile是数据文件名。

用Mplayer下载流媒体文件

互联网上有些音像文件是以流媒体形式存在的，例如rtsp,mms协议。使用mplayer可以把这些文件下载到本地机器上。例如：

mplayer -dumpstream mms://hostname/directory/filename

下载之后的文件名是stream.dump，只要改回相应文件名就可以了。

古典音乐作品编号标记

很多时候我们见到一些古典音乐作品的名字后面跟着由字母和数字组成的编号,例如op37, kv525, RV4, D112等.究竟这些标号表示什么意思,来源于哪里,有什么作用？最常见的记号是op或opus,拉丁语的意思是"作品".它大概起源于17世纪初叶,出版音乐乐谱的出版商为了区别同一作曲家的不同作品而标上记号,这就是op记号. op记号通常是按同一作曲家出版作品的次序由小到大编号,但这又并不是绝对的.不过, op标记也有问题,例如很多作曲家的作品由不同出版商出版,那时不同出版商有不同的op序号.因此会出现某作品在一个出版商标号为op3,在另一出版商标号为op8,甚至有第三个出版商那里编号为op6,这样显然会造成混乱. 　　到了近两百年,越来越多学者研究前人的音乐作品.这些学者为研究这些作曲家的作品,列出他们所作的所有作品进行研究甚为重要.由于前面所讲的 op记号有时会造成混乱,或是研究者在研究过程中发现新的作品,又或者是一些作品未正式出版而没有op记号,这时对所有作品进行重新编号便成为这些学者的又一项工作.由此就产生如BWV, D等编号,这些编号往往是提出新作品编号的学者名字的代表字母.下面列出一些作曲家作品常用的记号: 　　作品记号/适用的作曲家 　　AV Richard Strauss 理查-斯特劳斯 　　B Dvorak 德沃扎克 　　BB Bartok 巴托克 　　BWV J.S. Bach 巴赫 　　D Schubert 舒伯特 　　F Vivaldi 韦瓦尔第 　　H Holst 侯斯特 　　HWV Handel 汉度 　　K Mozart 莫扎特 　　KV Mozart 莫扎特 　　L Debussy 德彪斯 　　M Vivaldi 韦瓦尔第 　　Op 多数作曲家 　　P Vivaldi 韦瓦尔第 　　R Vivaldi 韦瓦尔第 　　RV Vivaldi 韦瓦尔第 　　S J.S. Bach 巴赫 (与BWV相同) 　　S Liszt 李斯特 　　Sz Bartok 巴托克 　　WWV Wagner 韦格勒 　　Z Purcell 帕萧

pdfnup: 让pdf文件多页合并到一页的小工具

有时候做handout的时候，想把几页slides打印到一张纸上。pdfnup可以很方便的把一个pdf文件的几页纸合为一页。例如，如下命令

pdfnup --nup 2x3 --frame true --delta "1cm 1cm"
        --scale 0.9 handout.pdf

自动生成一个每页纸包含2x3个slides的pdf文件。

半音阶口琴笔记——持琴法

半音阶口琴基本的持琴方法是：左手拇指和食指捏紧琴身，右手食指抵住按键，双手手掌自然合拢，形成一个共鸣腔。这样出来的音色饱满甜润，方便使用手震音。还有一种半开放的持琴法，即双手不合紧，下面打开一定空间，这样出来的音色明亮活泼。我以前使用第一种方法多一些，但是后来发现封闭的持琴法音量毕竟小一些，不太适合演出时使用麦克风的需要，尤其麦克风不是很好时，所以现在用第二种方法多一些。这样需要尽量使用气震音。然而，即使是半开放的持琴法，手震音还是有效的。 按琴键的食指可以分别使用上中下三个关节。食指上关节比较软，一般抒情乐曲可以满足要求，但是演奏快速活泼的乐曲时，觉得按键速度不够快，这时候似乎用中下关节比较好。如果要演奏快速的颤音，还需要配合手腕的抖动。

沉痛悼念世纪的良心——巴金老人

今天上网惊闻巴金老人去世了。虽然早知道巴金老人年纪很大了，一直患病，但仍不免心情沉重。 巴金是我最喜爱的作家。他早期的作品充满了青年火一样的热情，对真理的渴望，以及对于劳动人民深切的同情。我正是在那火一样的年纪里接触到了巴金的作品，感觉到他说出了我一直以来想说却没能说出口的话，感受到了从绝望中寻找希望的那种强烈共鸣，从此以后一直深深喜欢巴金的作品。从他那里我学到了什么是善良，什么是青春。他是一个把自己的命运和社会的前途融入一起的人，他的散文作品有强烈的个人意识，几乎总是以第一人称出现，表达对社会和人生的思考。小说《家》虽然不算自传小说，但是有巴金本人深深的影子。巴金把自己的思想通过主人公觉慧淋漓尽致的表现了出来，表达了对封建家庭的厌恶和对新生活的渴望。 巴金老人一直主张讲真话，对自己在文革中未能坚持真理内咎了一辈子。五集《随想录》便是巴金对于自己这段人生的思考和解剖。《随想录》标志了巴金散文的一个新的高度。和青年时期的作品不同，《随想录》中的文章读起来宛如一位善良的老人坐在你的面前娓娓诉说他所经历的历史，以及他对于这段历史的思考。作品中反复强调说真话，解剖自己在这段历史中犯下的错误，思考人们如何避免这段悲剧历史再现。巴金老人一生有两个重要的愿望：一个是建立中国现代文学馆，这个愿望在巴金老人的不懈努力下已经实现，中国现代文学馆数年前已经正式开馆，文学馆的门把手就是依据巴金老人手的模型制造的，每位参观者仿佛都会被巴金老人亲手领进这个文学宝库；另一个愿望便是建立文革博物馆，思考悲剧产生的原因，警戒后人不要犯同样的错误，这个愿望在巴金老人的有生之年未能实现。具有讽刺意义的是，巴金老人仅仅因为自己在文革中没有讲真话而自责了一辈子，而真正制造悲剧的人却都活的好好的，至今没有任何人出来为文革的发生承担责任。 巴金老人给称为世纪的良心。如今，他离我们而去了，还有谁能够在官话套话横行的世界里举起讲真话的旗帜，凭自己的良心讲话呢？

在Linux下把wma文件转化为mp3

需要用到mplayer和lame

mplayer "filename.wma" -ao pcm:file="filename.wav"
lame -h "filename.wav" "filename.mp3"

transcode编码DVD的bug

dvd::rip在Linux下把DVD压缩编码为avi，当transcode时无法进行，提示“Job 'Transcoding video - title #1, pass 1' failed”。根据google的结果，是transcode在AMD处理器上的bug。transcode时手动加上transcode options "--accel 3dnow"可以解决问题。

用aspell在vim里作拼写检查

:w
:!aspell -c %