2008-May 15Thursday

地动山摇

2008年5月15日11:10　 书虫之语, 学而时习　4条评论 || 4 Comments

　　承蒙信息时代传递之速，尽管我并没有感觉到地震，但是在四川大地震主震发生之后不到半小时，我就知道了这一可怕的事件。昨天尽自己微薄之力，捐款400大洋，总算觉得有所作为，心安定了一些。媒体传来传去的事实，比如毁坏情况，天国人数，以及可歌可泣的事件等，死难者的惨状令人悲哀，幸存者的坚毅令人感动，但这留给媒体和社会讨论，我想也无必要反复叙述。本人天文可能知一点，地理的话大体是不知的，在这里由于科学好奇心发作，尽力避开大家耳熟能详的事实，也斗胆谈一下吧。
　　地震的定义很简单：由地震波所造成岩石圈的震动。地震波可由自然现象和人为原因造成，因此想体验地震很容易：去工地或者铁路旁边，你就能感受得到。我倒是感受过一次真的地震，那是2004年日本新泻大地震之后几天，在东京迪士尼乐园的高架火车站上感受到的一场里氏6.0级的余震（美国地质调查局国家地震信息中心/USGS测定结果）。幸好只是猛震了一次，而且那座高架车站很坚固，不然今天我恐怕就没法给大家在这里吹牛皮啦。

2008年5月14日上午10:36，在大学城拍到疑似地震云。注：地震云的出现不代表地震会在该处发生

　　我主要想说说的是关于地震的预测，虽然我也完全是外行，只是刚刚去学习了一些相关文献。因为这几天看到很多人对地震学家发起猛烈的攻击，更有不少人和社会各界开恶意玩笑，宣称“此地、某时”会有强烈余震。2006年底的台湾海峡地震之后，便有人宣称“凌晨2时30分大学城将有地震，请大家起来避险”，结果有许多人被骗得半夜三更爬起来胆战心惊。
　　从各种资料来看，地震预测的难度非常高，以至于“地震学家一般不接受地震预测”。与气象预测非常不同，地底的状态人类所知甚少。就拿台风来说，台风初期的热带扰动远在陆地几千公里以外，我们也可以通过卫星图像来准确测定其位置，并估计其中心的湍流情况；但目前人类最深的钻井还不到25公里（壳牌公司，2004年），我国最深的钻井仅为5公里（http://zhidao.baidu.com/question/22268227.html），而本次大地震的震源深度就有19公里（系主震震源，根据USGS更新后的测定结果）。所以，即使是提地震预测，也仅仅是统计学意义上的地震预测，比如有科学家预测2003-2032年间旧金山湾区发生高危地震的概率为62%。
　　在国际范围内认可的、准确的地震预测记录，历史上似乎只有两次，均在中国，一次是1975年2月的海城7.3级地震（Richard G., “Earthquake prediction: Haicheng, China - 1975″, Earth Science），当时政府疏散了约百万人，从而使伤亡人数降到未疏散估计值的1/5；另一次是1999年11月的海城5.4级地震。对于本次四川大地震，似也有个人成功预报者（例如http://www.meteoquake.org/cin.html，提到一位中国科学家在一个月前并预测了地震，误差为8天；另有网络上流传的提前预测），但由于无法查证更多的文献资料，故暂时不算在内。亦有以前国外科学家的多篇文献进行某些地震预测，但大多失败。另有一位华裔科学家利用地震云来进行地震预测（http://quake.exit.com），准确性似乎颇高，但由于地震云的认定暂时还不具有科学意义上的普适性，科学界并不接受利用地震云来进行地震预测的方法。

中大学子踊跃提供助灾物资，一个下午便收集到如此多袋。回宿舍的路上还看到十来个宿舍倾巢而出，提着大包小包向募捐地点走去

一个下午加一个晚上，楼下的休息室已经放了不少助灾物资

　　有很多人指出四川大地震前出现气象以及生物学上的诸多异常，然而，科学上的、准确的说法是“尚无法证明动物能预知地震”。假使动物能感觉地震前的纵波，那么地震仪对此更敏感。通过包括气象和生物学上的异常现象来预测地震，“技术上的可靠性无法确立，也无法反演”，我们应当了解，地震预测很大程度上“难以满足实践上和科学上的可靠性”。因此，我们必须接受一个事实：资料的严重不足以及技术方面的极其匮乏，使得地震预测的难度非常大。尤其提醒各位，不要把气象预测的难易度和地震预测的难易度混为一谈。
　　最后想谈的一个问题是，大地震会在一个地方短时间内连续发生吗？因为看来有许多朋友在这次大地震发生以后，夜晚睡眠质量大降，担心很快会有另一次大地震发生。本人查证了USGS统计的自1973年以来的538次里氏7.0级以上的地震，结论是：这种可能性很低。一次大震之后虽然余震不断，但威力相当的大余震似乎在短时间内发生的概率很小。四川大地震以后最强的余震仅为6.0级，发生在主震之后15分钟。
　　好吧，就先谈这么多，本人才学疏浅，不正之处，还请各位指正。地震无情人有情，各位有钱出钱，有力出力，顺为受难者祈祷。

参考文献（按引用顺序）：
“世界上最深的地址钻井有多深，在哪里”，http://zhidao.baidu.com/question/22268227.html
EARTHQUAKE CLOUDS AND SHORT TERM PREDICTION ，http://quake.exit.com/
Latest Earthquakes M5.0+ in the World - Past 7 days，http://earthquake.usgs.gov/eqcenter/recenteqsww/Quakes/quakes_big.php
Richard G., “Earthquake prediction: Haicheng, China - 1975″, Earth Science
Message from Prof. Song Qi, http://www.meteoquake.org/cin.html
USGS/NEIC (PED) Earthquakes 1973 - Present, USGS/NEIC
Earthquake prediction, http://en.wikipedia.org/wiki/Earthquake_prediction
Earthquake cloud, http://en.wikipedia.org/wiki/Earthquake_cloud
Abercrombie R., Seismology: The start of something big?, Nature V. 438 N. 7065, pp. 171-172

2008-Apr 22Tuesday

《齐格弗里德》听后记

2008年4月22日15:26　 学而时习, 椰子笔谈　评论/Comment.

　　齐格弗里德和布伦希尔德是《指环》全局的中心人物。与后者不同，齐格弗里德直到本剧中才出现，尽管在序剧和第一剧中我们可以模糊的感到这一人物的存在。第一幕，前奏曲沉稳平缓但略嫌阴暗的旋律表现了“沉思的动机”，暗示着一个阴谋正在进行之中。弦乐奏起了“尼伯龙根的动机”，把它和前者联系起来，我们可以大概猜到后面将要表现什么；音乐慢慢生动起来，接着是管乐的“指环的动机”，事件的脉络慢慢清晰起来，到达高潮的时候小号吹出的“剑的动机”，渐缓。
　　铁砧发出叮叮当当的声音：在森林深处一个洞穴里，尼伯龙根族侏儒、阿尔伯里希的弟弟米梅正在铸剑，可是不管他用什么方法，他每铸好一把剑就被齐格弗里德弄断。米梅希望齐格弗里德能够杀死变成大蛇看护指环的巨人法弗纳，为他取得指环，因此他急于早日造出新剑。这时，小号吹响了齐格弗里德的动机，披着兽皮的齐格弗里德进来了，见他牵着捕获的一头熊，米梅吓得滚倒在地。看到米梅如此惶恐，齐格弗里德大笑起来，直到米梅表示新剑已铸好，齐格弗里德才将大熊放回森林，准备一试新剑的力量。
　　然而米梅担心的事情又发生了：新剑又被齐格弗里德击碎，乐队奏起了充满年青人气力的乐音，齐格弗里德怒气冲冲地责难米梅无用。米梅的动机忧愁地响起，伤心的米梅哭着讲起他是如何养育齐格弗里德的，但代表齐格弗里德的铜管乐器不耐烦地打断了米梅忧愁的叙述：“你这套故事我早就听腻了！我一点都不喜欢你。你面貌是如此丑陋，我长得和你一点都不一样！我看到林中的野兽都成双成对的出现，你肯定隐瞒了什么。我现在要你把有关我父母的事情统统告诉我！”米梅无奈，终于缓缓地说出实情。多年前他救了一个在森林中昏倒的女人，她在山洞中生下齐格弗里德后不久便死了，那女人的名字叫齐格林德，至于他父亲则不知是被谁杀害的。“没有证据，你要我怎么相信你呢？”齐格弗里德说。这时，乐队奏起了我们在《女武神》中已经耳熟能详的一段音乐—— 诺顿宝剑的旋律。米梅找出齐格蒙德（齐格弗里德的父亲—— 我加个注释，因为那俩兄妹的名字实在很容易搞混）留下的断剑交给齐格弗里德对证。齐格弗里德仔细的审视起断剑，命令米梅必须将剑铸造成新的，说完吹着号角出门去了。

米梅营救齐格林德，A. Rackham作

　　齐格弗里德的号角渐渐远去。乐队奏起了“绝望的动机”，我们看到无法铸好剑的米梅沮丧地坐在地上。这时，扮作流浪者的众神之王沃坦出现了（微弱的“瓦尔哈拉神殿”的动机提醒着我们这个流浪者的身份），他请米梅问他三个问题，答不出的话愿以自己的头颅作为代价。本剧开始时“沉思的动机”再次浮现，米梅略加思考便问：“哪个种族住在黑暗的地下？”“是尼贝尔海姆的尼伯龙根族的侏儒。他们的首领阿尔伯里希抢走了莱茵的黄金，制成魔力无穷的指环。”“很对。”米梅说，又问：“地上住着哪一族类？”“巨人。”流浪者回答说，——乐队奏出那让人印象深刻的“气势汹汹但只有两个音组成的简单旋律”—— 巨人的旋律，应合着沃坦的答案—— “受到阿尔伯里希的诅咒，弟弟法弗纳为了指环杀了哥哥法索尔特，现在变成一条巨蛇看守尼伯龙根的财宝。”“没错。那么，天上为谁的领地？”“沃坦，他与众神住在瓦哈拉统治世界。”沃坦面不改色地回答，背景是“瓦尔哈拉动机”的再现。流浪者答对了米梅的三个问题，现在轮到他问米梅：“这世上哪一族的人血亲间互相残杀？”米梅回答道：“维尔塞族。齐格蒙德是沃坦的儿子，但亲手葬送了他的却是沃坦。现在维尔塞族中最勇敢的人是齐格蒙德的儿子齐格弗里德。”“那么，现在有一个尼伯龙根人将齐格弗里德养育成人，想靠他杀死巨人取得指环，齐格弗里德要用什么武器才能打赢大蛇？”“齐格蒙德留下的一把断剑。”米梅不假思索地回答（这个沃坦，明知故问），乐队轻柔地奏出“剑的动机”。“谁能够重新铸好这把剑？”沃坦追问道。乐队尝试给米梅提醒，奏出了“齐格弗里德的动机”，不过米梅当然是听不到的—— 于是乐队再度奏出“绝望的动机”。“我…… 我不知道。”“那么，听我的一个预言吧。铸好剑的会是一个英雄，但你的命将会失在他手里。”说完，流浪者就走了，剩得米梅一个人坐在地上发呆，心怀恐惧。
　　过了一会儿，齐格弗里德哼着歌回来了，他看到米梅坐在地上不铸剑，大为生气，责问米梅为什么不铸剑。神志不清的米梅却只一味说着恐惧。“什么叫‘恐惧’？”齐格弗里德问道。这时米梅清醒过来，他赶紧讲起火焰与变作大蛇的法弗纳来掩饰自己的惶恐。“那些东西我都不怕，”齐格弗里德恼火的说，“快点把剑弄好，这样我就可以出去闯荡了，也可以学习一下什么叫‘恐惧’。”然而米梅却仍然无意铸剑。急不可耐的齐格弗里德就自己动手，他一边生火一边唱起了“铸剑之歌”，这是号角和诺顿神剑旋律的混合。接着，以铁砧为主角的打击乐也登场了，在当当当的响声和“铸剑之歌”的应和之下，齐格弗里德奋力地挥动铁锤，每一下敲打都激得火星四溅。乐队不可抑制地显得越来越宏伟和明亮，弦乐器短促的泛音描述着剑浸入水里发出嗤嗤的声音，火光照亮了齐格弗里德充满喜悦的脸。在一旁的米梅看到剑将铸成，联系起流浪者的预言，他心怀深一步的恐惧地开始盘算，想等齐格弗里德杀死大蛇后就骗他喝下迷药，自己再将他杀死，于是在一旁开始煮一锅加入迷药的汤。剑终于铸成了！齐格弗里德高举起剑欢呼，扩大版的铜管阵容奏出能令人产生非理性的兴奋的巨响，兴奋的齐格弗里德一剑朝铁砧劈去，铁砧竟然被劈成两半！
　　第二幕讲述了齐格弗里德杀死大蛇法弗纳并夺取指环的故事。前奏曲，大提琴与定音鼓奏出“巨人的动机”的变型，以此代表法弗纳，接着低音号加入进来，奏出预示危险的“蛇的动机”，说明法弗纳已经变为巨蛇。之后是木管“指环的动机”与铜管“诅咒的动机”，“憎恨的动机”将旋律推向高潮。随后“憎恨”逐渐减弱，代以弦乐和木管奏出黎明前森林的幽幽声响。在法弗纳洞穴前的阴影中，阿尔伯里希心怀仇恨地窥探被夺去的宝藏，一心想要再夺回来。风吹开了树枝，一缕光照射过来，另一个幽幽的旋律加进来—— 瓦尔哈拉旋律，沃坦穿着流浪者的衣服出现。阿尔伯里希认出了沃坦，怨恨的咒骂他，沃坦告诉阿尔伯里希关于米梅的阴谋后离开。过了一会儿，齐格弗里德与米梅一起走来，阿尔伯里希立刻躲起来。米梅醉翁之意不在酒，对齐格弗里德大谈大蛇的可怕之处。“够了，你给我滚开。”齐格弗里德厌烦透顶，“让你看看我是怎么干掉他的。”米梅退到不远处，心里暗暗希望齐格弗里德与法弗纳同归于尽。齐格弗里德坐在一棵树下等待大蛇醒来，一边忧愁地想起自己的身世。这时，黎明森林的幽幽声响打断了他的思绪，长笛和小提琴领衔乐队，轮番奏出清晨万物苏醒的天籁之音。齐格弗里德深为此着迷，屏住气欣赏着。少顷，单把小提琴奏出“小鸟的旋律”，原来小鸟聚在齐格弗里德身边，只是唱的歌他都听不懂。为了表示愿意和小鸟做朋友，齐格弗里德吹起了用芦苇作的笛子（听不出乐队用的是什么乐器—— 好像是简易喇叭？），却难听得很，无法模仿鸟鸣声，于是他吹响了自己的号角—— 单支圆号快速地交替奏起了“号角的动机”和“齐格弗里德的动机”。但圆号所演奏的很快变成了阴暗而恐怖的旋律，之后由数支长号来演奏—— 不，现在所演奏的是“蛇的动机”，原来嘹亮的号角声惊醒了大蛇，它爬出洞穴。定音鼓捶出了“巨人的动机”。齐格弗里德毫不畏惧，走上前去问大蛇法弗纳：“你知道什么是恐惧吗？”定音鼓再捶出“巨人的动机”：“你竟然见到我都不知道什么叫恐惧！”大蛇张开嘴露出獠牙威胁他。齐格弗里德拔剑刺向大蛇，一阵扭打之后，齐格弗里德把剑插入大蛇的心窝，乐队忙不迭地交替演奏“齐格弗里德”、“剑”和“憎恨”，后者描述的是大蛇痛苦地卷曲，蛇血溅在齐格弗里德的手上，他把手指举到唇边舔，一瞬间，他听懂了小鸟的语言。弦乐演奏起久别的“指环”和“隐身盔”的旋律，原来小鸟告诉齐格弗里德关于尼伯龙根宝藏的事情，齐格弗里德谢过它们，走进洞中寻找。
　　齐格弗里德找到宝藏，走出山洞，这时阿尔伯里希和米梅正在争吵。听到“指环”和“莱茵黄金”的旋律，阿尔伯里希心虚地溜走了。米梅端着迷药走过来，小鸟叮咛齐格弗里德要当心米梅。米梅强迫齐格弗里德喝汤，看穿他的齐格弗里德一剑将他刺死，之后疲倦的靠着树坐下，他感到十分孤独—— “号角的动机”略显疲乏。坐了一会儿，齐格弗里德抬头问树上的小鸟在哪里能找到一个朋友，于是小鸟唱起关于沉睡在被烈火包围的岩石山上的布仑希尔德的故事，齐格弗里德听得心花怒放。木管乐曲快速地演奏愉悦的旋律（这在瓦格纳的音乐中可不多见），告诉我们齐格弗里德迫不及待地要出发去找布伦希尔德。太阳越升越高，森林弥漫着生气，于是小鸟带齐格弗里德出发了。
　　第三幕的前奏曲将第二幕结尾原本愉快的旋律又消隐得无影无踪，“契约”和“命运”交替出现，暗示着《指环》的大结局。阴暗的旋律继续延续，定音鼓和大镲诠释着此时的背景—— 暴风雨。这是岩石山荒凉的山脚处，黑夜。扮作流浪者的沃坦唤醒在地底沉睡的智慧女神埃尔达，询问众神的未来。智慧女神睁开眼睛，只说齐格弗里德与布仑希尔德能够摆脱诅咒，乐队用演奏“爱情的动机”诠释着—— 这是《女武神》中用于表达齐格蒙德和齐格林德的旋律，只不过稍微做了一点转换，—— 接着是略为阴暗的“瓦尔哈拉动机”，原来智慧女神对于众神的命运避不作答。沃坦说他已经不怕众神的灭亡，因为这是“他的意志”，他的女儿布伦希尔德将会拯救世界。埃尔达对此并不感兴趣，无奈的沃坦只得将埃尔达送回梦乡，他已厌倦了与命运抗争，隐隐的感到有一种新的力量将改变现状，神的时代正在衰败。

齐格弗里德邂逅布伦希尔德，A. Rackham作

　　过了一会儿，“小鸟的动机”响了起来，原来小鸟引着齐格弗里德来到岩石山。沃坦横枪加以阻拦，不让齐格弗里德上山；齐格弗里德并未认出他的祖父，又急于见到布仑希尔德，因此对沃坦毫不客气，这时乐队又奏起“剑的动机”：哐！齐格弗里德一剑劈断了象征众神之王沃坦权利与力量的长枪。沃坦平静地看着齐格弗里德离去，拾起长枪的碎片离开了。山顶上一圈火焰熊熊燃烧，红色的火光照亮了黑夜，但乐队奏起了坚毅的旋律：齐格弗里德吹响了他的号角，无畏的踏进烈火中，火焰在他的面前消退，最后云烟飘散。天空变得晴朗，齐格弗里德发现了在盾牌下沉睡的布仑希尔德，乐队轻柔的奏起“少女的动机”和“安眠的动机”，他将盾牌搬开，终于看见美丽的布仑希尔德，心里又惊又喜—— 这里，瓦格纳巧妙的把前两剧中爱情之神福瑞娅和婚姻之神弗里卡的旋律，和代表维尔塞族和喜悦的旋律让乐队同时演奏出来，预示着情节的发展。布仑希尔德苏醒过来，唱起了赞美阳光的歌。凝固的美丽慢慢生动起来，布伦希尔德眼神透露出高贵，金色的头发随风飞扬，齐格弗里德不由心动。布仑希尔德问将她唤醒的英雄是谁；对于我们听众来说，这个问题是乐队代齐格弗里德回答了—— 先演奏出坚毅的“齐格弗里德的动机”，然后又演奏“命运的动机”。接下来乐队不可自拔地陷入充满热情的旋律，描述着俩人相互凝视，打心里欢喜，齐格弗里德的心已经陷入了爱情，但这种陌生的感情让他感到为难；而布仑希尔德也在为被永远的驱逐出女武神行列的命运感到悲伤，但强烈的爱情驱散了他们心中的阴影，人性中美好的感情如一道暖流注入彼此的心中，两个人立下了永远相爱的誓言。至此，乐队奏出洪亮的“幸福的动机”和“爱情的动机”，《齐格弗里德》剧终。

2008-Mar 14Friday

[zz]庆祝圆周率日

2008年3月14日10:24　 学而时习, 鸡毛蒜皮　2条评论 || 2 Comments

圆周率的历史

　　古希腊欧几里得《几何原本》（约公元前3世纪初）中提到圆周率是常数，中国古算书《周髀算经》（ 约公元前2世纪）中有“径一而周三”的记载，也认为圆周率是常数。历史上曾采用过圆周率的多种近似值，早期大都是通过实验而得到的结果，如古埃及纸草书（约公元前1700）中取π=（4/3）^4≒3.1604 。第一个用科学方法寻求圆周率数值的人是阿基米德，他在《圆的度量》（公元前3世纪）中用圆内接和外切正多边形的周长确定圆周长的上下界，从正六边形开始，逐次加倍计算到正96边形，得到(3+(10/71))<π<(3+(1/7)) ，开创了圆周率计算的几何方法（亦称古典方法，或阿基米德方法），得出精确到小数点后两位的π值。
　　中国数学家刘徽在注释《九章算术》（263年）时只用圆内接正多边形就求得π的近似值，也得出精确到两位小数的π值，他的方法被后人称为割圆术。他用割圆术一直算到圆内接正192边形。
　　南北朝时代数学家祖冲之进一步得出精确到小数点后7位的π值（约5世纪下半叶），给出不足近似值3.1415926和过剩近似值3.1415927，还得到两个近似分数值，密率355/113和约率22／7。其中的密率在西方直到1573才由德国人奥托得到，1625年发表于荷兰工程师安托尼斯的著作中，欧洲称之为安托尼斯率。
　　阿拉伯数学家卡西在15世纪初求得圆周率17位精确小数值，打破祖冲之保持近千年的纪录。
　　德国数学家柯伦于1596年将π值算到20位小数值，后投入毕生精力，于1610年算到小数后35位数，该数值被用他的名字称为鲁道夫数。
　　1579年法国数学家韦达给出π的第一个解析表达式。
　　此后，无穷乘积式、无穷连分数、无穷级数等各种π值表达式纷纷出现，π值计算精度也迅速增加。1706年英国数学家梅钦计算π值突破100位小数大关。1873 年另一位英国数学家尚可斯将π值计算到小数点后707位，可惜他的结果从528位起是错的。到1948年英国的弗格森和美国的伦奇共同发表了π的808位小数值，成为人工计算圆周率值的最高纪录。
　　电子计算机的出现使π值计算有了突飞猛进的发展。1949年美国马里兰州阿伯丁的军队弹道研究实验室首次用计算机（ENIAC）计算π值，一下子就算到2037位小数，突破了千位数。1989年美国哥伦比亚大学研究人员用克雷－2型和IBM－VF型巨型电子计算机计算出π值小数点后4.8亿位数，后又继续算到小数点后10.1亿位数，创下新的纪录。至今，最新纪录是小数点后12411亿位。
　　除π的数值计算外，它的性质探讨也吸引了众多数学家。1761年瑞士数学家兰伯特第一个证明π是无理数。1794年法国数学家勒让德又证明了π^2也是无理数。到1882年德国数学家林德曼首次证明了π是超越数，由此否定了困惑人们两千多年的“化圆为方”尺规作图问题。还有人对π的特征及与其它数字的联系进行研究。如1929年苏联数学家格尔丰德证明了eπ 是超越数等等。

背诵口诀

3  .   1 4  1 5 9      2 6     5  3 5     8  9 7     9  3 2     3  8 4     6  2  6
山巅一寺一壶酒，尔乐。苦煞吾，把酒吃，酒杀尔，杀不死，乐尔乐。
4  3  3    8  3  2 7      9 5  0  2 8      8 4  1 9  7      1  6 9  3 9      9 3  7
死珊珊，霸占二妻。 救我灵儿吧！ 不只要救妻， 一路救三舅， 救三妻。
5  1 0  5 8      2 0  9     7 4 9  4 4  5  9 2     3  0 7
我一拎我爸，二拎舅（其实就是撕我舅耳）三拎妻。
8  1  6    4  0  6    2  8  6    2  0 8     9  9  8  6
不要溜！司令溜，儿不溜！儿拎爸，久久不溜！
2  8 0     3  4  8     2 5  3 4  2 1     1  7 0  6     7  9  8
饿不拎，闪死爸，而我真是饿矣！要吃人肉？吃酒吧！
（作者华罗庚）

圆周率前2008位

3.1415926535 8979323846 2643383279 5028841971 6939937510 5820974944 5923078164 0628620899 8628034825 3421170679 8214808651 3282306647 0938446095 5058223172 5359408128 4811174502 8410270193 8521105559 6446229489 5493038196 4428810975 6659334461 2847564823 3786783165 2712019091 4564856692 3460348610 4543266482 1339360726 0249141273 7245870066 0631558817 4881520920 9628292540 9171536436 7892590360 0113305305 4882046652 1384146951 9415116094 3305727036 5759591953 0921861173 8193261179 3105118548 0744623799 6274956735 1885752724 8912279381 8301194912 9833673362 4406566430 8602139494 6395224737 1907021798 6094370277 0539217176 2931767523 8467481846 7669405132 0005681271 4526356082 7785771342 7577896091 7363717872 1468440901 2249534301 4654958537 1050792279 6892589235 4201995611 2129021960 8640344181 5981362977 4771309960 5187072113 4999999837 2978049951 0597317328 1609631859 5024459455 3469083026 4252230825 3344685035 2619311881 7101000313 7838752886 5875332083 8142061717 7669147303 5982534904 2875546873 1159562863 8823537875 9375195778 1857780532 1712268066 1300192787 6611195909 2164201989
—– [1000] —–
3809525720 1065485863 2788659361 5338182796 8230301952 0353018529 6899577362 2599413891 2497217752 8347913151 5574857242 4541506959 5082953311 6861727855 8890750983 8175463746 4939319255 0604009277 0167113900 9848824012 8583616035 6370766010 4710181942 9555961989 4676783744 9448255379 7747268471 0404753464 6208046684 2590694912 9331367702 8989152104 7521620569 6602405803 8150193511 2533824300 3558764024 7496473263 9141992726 0426992279 6782354781 6360093417 2164121992 4586315030 2861829745 5570674983 8505494588 5869269956 9092721079 7509302955 3211653449 8720275596 0236480665 4991198818 3479775356 6369807426 5425278625 5181841757 4672890977 7727938000 8164706001 6145249192 1732172147 7235014144 1973568548 1613611573 5255213347 5741849468 4385233239 0739414333
4547762416 8625189835 6948556209 9219222184 2725502542 5688767179 0494601653 4668049886 2723279178 6085784383 8279679766 8145410095 3883786360 9506800642 2512520511 7392984896 0841284886 2694560424 1965285022 2106611863 0674427862 2039194945 0471237137 8696095636 4371917287 4677646575 7396241389 0865832645 9958133904 7802759009
  —– [2000] —–
94657640

按：传说今天是“白色情人节”，叹~ 人们太善于编造让两头动物有理有据地呆在一起的氛围了，难道一个Valentine’s Day还不够吗！？看来上帝设计了相当多的诡计，以便使（大部分）人们在一些闭合的圆圈上耗费时间，不能去探寻他那王国的奥秘。正感慨上帝的狡诈，忽然发现今天原来是“圆周率日”！哈哈，对于这个可爱数字的节日，怎可以不庆祝庆祝呢？各位，友情提示：今天是“圆周率日”！以上内容转自http://baike.baidu.com/view/3287.htm

2007-Jun 20Wednesday

MPD计划和它的新进展

2007年6月20日23:22　 学而时习, 鸡毛蒜皮　评论/Comment.

　　MPD计划是我在进行的一个雄心勃勃的、具有国际先进水平的大计划，将动用我的全部编程功底来做一个大程序。呵呵—— 那具体是什么哈？总之涉及天体轨道计算一类，具体内容暂时保密……

　　今天利用午饭时间仔细思考了一下，花了大约一堂课时间导出了指定条件下Ω和ω值的算式，使得那个有9个变量的微分方程变成了只有7个变量…… 其中好像有5个是指定值的？那明天的任务就是考完英语以后导出剩下两个变量了。

　　为了完成这第一个步骤，还要研究Dr.Jedicke等的论文，到时候再来涂两把……

2007-Apr 21Saturday

Cygwin下架设php

2007年4月21日23:34　 学而时习, 鸡毛蒜皮　评论/Comment.

　　今天一早起来很悲哀地发现MPC的另一个小bug，这个小bug浪费了另一个六小时。什么bug呢？具体来说就是忽略了23小时59分和1天之间的区别。

　　千里之堤，溃于蚁穴呀……

　　随后更郁闷的发现，专门给我做登陆跳板的服务器挂了，管理员说星期一才会把它搞好。于是我不得不决定：自己架一个linux的php环境，把MPC移到本地来进行！

　　最佳选择自然是cygwin。安装cygwin还是挺简单的，不过要在这个linux+windows的杂合体中安php却是件难事。搞不懂为何php这等流行语言都进不了cygwin的大门，还要我们来麻烦地编译？

　　搞到最后终于搞定了，就把足迹记下，方便后来人，因为我也是好麻烦才搞定的。虽然cygwin论坛上有帖子介绍（http://www.cygwin.com/ml/cygwin/2004-01/msg00383.html），不过这个方法太麻烦了。我的方法虽然可能不是那么正宗，但现在好像还没出什么大问题（貌似有一点小问题，也许后来人可以解决一下）。

　　首先先要把gcc和flex这两个插件装上，cygwin里面都有，不过默认是不装的。否则对php做configure的时候会出错。而且php不能用5版本的，我用的是4.4.6版。做configure，大约要三五分钟，之后将make中的libphp4.so更为libphp4.dll，make，再到libs文件夹执行如下命令：

gcc -shared -o libphp4.dll -Wl,–out-implib=libphp4.dll.a -Wl,–export-all-symbols -Wl,–whole-archive=libphp4.a /usr/local/apache/libexec/libhttpd.dll -Wl,–no-whole-archive -lcrypt

　　执行完毕以后，回到根目录，再做make install。现在就好啦，测试一下吧？php run-test.php，运行得还可以吧？

　　只是我发现无论运行什么都会出现X-powered by php之类的信息，难道是哪里出错了？谁有兴趣可以试试看。

　　好了，之后可以在这个临时搭建的cygwin+php下运行我的MPC了，出不少错误而且速度很慢（哪里比得上国台的4芯机子……）。上面的过程让我从晚5点忙到晚8点，终于搞好了，开始运行MPC。结果晚9点的时候管理员通知我服务器恢复，我兴冲冲的关了cygwin准备再次享受4芯的畅快，结果发现：服务器还没有好！

　　>_< 懒得重开cygwin，只好再次艰难的等待！

2007-Feb 21Wednesday

傻傻的电脑

2007年2月21日23:39　 学而时习, 泉志咨文　评论/Comment.

　　为了第三代7Timer!的自动运行，我花了两小时啃入了BASH的门。我的程序中某一个地方要循环体，需要定义一个循环相加的变量。于是我定义
TEST = 9
　　下一个循环中它要是12。根据学C、jsp、php之类的经验，我毫不犹豫地写
TEST = $TEST+3
echo $TEST
　　各位猜怎么着？电脑居然显示
9+3
　　我狂晕！觉得又好气又好笑。当然，这等问题本来不应该拿到桌面上说，不然肯定被老手笑话你这个菜鸟不见世面，可我的确狂笑了1分钟：电脑啊电脑，你究竟还是只能听人指令的傻东西！哈哈哈哈哈哈哈哈！
　　后来查了一下，正确的写法应该是
let “TEST = $TEST+3″
echo $TEST
　　电脑就能够显示“12”了。

2007-Feb 5Monday

第四代7Timer!—— 5公里精度？

2007年2月5日15:28　 学而时习, 鸡毛蒜皮　评论/Comment.

　　野心太大，有时真是一件坏事！嘿嘿！

　　昨天在7Timer!-INFO N.63里兴致勃勃地发布了GRIB2数据成功解码的消息，又顺利地把首页的图弄成了Mollweide投影，已经觉得是“重要成就”了，一个晚上都乐得屁颠屁颠的。可今天太阳还没落，GRIB2 decode技术好像已经是恐龙时代的事物一样。为什么呀？今天仔细一看才知道，NCEP打算今年8月1日开始全面转用GRIB2格式的数据，GRIB1到了那时候就不再发布。也就是说，7Timer!最迟在今年8月1日就得开始采用0.5deg（今天换算了一下，0.5deg应该是35km。55km是对应赤道的）数据了。今天才知道可以用NCEP自己发布的cnvgrib(http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2/)把GRIB2转为GRIB1，再用GrADS来处理。至此，7Timer!采用35km数据理论上已经扫清了最后一步障碍了。采用GRIB2数据的7Timer!将是第三代7Timer!，其精度可以大大提高。如果第三代7Timer!应用的话，那我的地点预测就可以恢复了。这项业务将命名为7Timer!-Panel2+，而准备投入民用的7Timer!-Weatherman（暂定名）也在开发中。

　　最不妙的就是很快又发现，大名鼎鼎的WRF模型原来是可以用GFS/ETA的数据来算出精度更高的预报。其实么，本来是想再看看MM5的可行性。最后发现MM5是可以做到的，因为现在我对Linux/Unix下的操作也算基本入了门，而MM5原来用的是NCEP Reanalysis的数据（还可以用国内的T213？不管，反正不好找到）。用了MM5，精度可以再上一个档次，就可以和北美CSC叫板了，正是我梦寐以求的啊。然而MM5毕竟要把现在的平台全部推倒重来，代价似乎太高了。用WRF的话，可以顺用GFS数据，只是用WRF内核来处理，然后再用WRF2GrADS转到GrADS下处理，现在的平台还可以用上约一半（不知道我理解有没有错误）。精度呢？哇哇哇，WRF可是中小尺度预测的，目标是1-10km。想象一下，5km精度的7Timer!，让人口水都要流出来了！

　　但精度做到这么高，代价肯定很高—— 服务器绝对吃不消。我的设想是，先对一些小地区（比如，北京上海广州）进行中小尺度的预测，如果效果不错，说不定可以向崔博士那边申请一台专门的机子……？哈哈，别想那么远了。现在第一步就是要上手WRFSI…… 万里长征第一步还是要迈出去才行啊。

　　至于原来的计划么？新版本ArP？唉，后延后延…… 又一个星期没出去了，还是要找天出去溜溜啊。各位，罗马不是一天建成的，但高精度的第四代7Timer!已经在地平线上出现了……

2007-Feb 4Sunday

认真研究Kreutz

2007年2月4日23:37　 学而时习, 鸡毛蒜皮　一条评论哦 || 1 Comment

　　最近要做的事情真多，不过也犯不着和别人说“忙”，因为这些事情么，都是召之即来的，要是不喜欢可以拖延，十分自由。

　　昨天拿到天爱《增刊》，发现里面介绍有位MM弄Kreutz轨道的。呵呵，牛啊，不看都不知道北方牛人原来还不少。自我批评一下，应该向dreamh学习，d调。

　　话说回来到Kreutz，一直很低调地研究克族的性质，现在也可公开一些了。地面研究虽然不是没有可能，但看来难度还是很高的，决定先准确掌握Kreutz性质再说。今天翻到Sekanina的一篇重要论文：http://www.journals.uchicago.edu/ApJ/journal/issues/ApJ/v607n1/59666/59666.html，因为实在太长又太复杂，决定慢慢啃，顺便放在这里看看谁有兴趣也来研究。晚上，朱头把孙同学的论文发了来，可谓雪中送炭哈！原来她把轨道根数已经算了个平均值出来。这些天我还偷懒不想算呢…… 这下有人送上门来了。孙同学还列出了若干重要数据，明天慢慢研究！现在先睡觉。

2007-Jan 21Sunday

梯瑟兰参数

2007年1月21日11:00　 学而时习　评论/Comment.

　　今天在讨论疑彗小行星2007 BJ的时候遇到一个新名词：Tisserand parameter，似乎是讨论天体似彗程度的一个指标。学习了一下，具体内容如下：

如果Tj小于3，则目标的轨道呈彗星状。

参考文献：

http://farside.ph.utexas.edu/teaching/336k/lectures/node122.html

http://tech.groups.yahoo.com/group/mpml/message/18685

2006-Dec 6Wednesday

致泰勒公式的抗议书

2006年12月6日22:27　 中东物语, 学而时习　评论/Comment.

致泰勒公式代表处：
　　我谨代表所有深受N项泰勒公式(N>1)展开复杂性之问题伤害之受害者，对泰勒公式表示严正抗议。泰勒公式以其虚假宣传、欺骗受害者，引诱受害者上钩后又用以二次等差数列式增长的混合求微商运算对受害者进行虐待，导致受害者使用后身心受到严重伤害，同时又有拉格朗日余项、皮亚诺余项等狼狈为奸之辈，为非作歹，加重了受害者的痛苦，而许多情况下又达不到其所宣称的无限迫近所求值之所谓“神奇”效果。以上几点严重伤害受害者感情、破坏人们为征服高等数学作出努力的行为，我们要求泰勒公式立即向所有受害者郑重道歉，并给予高数满分的赔偿。
签名：xxx,xxx,xxx,xxx,xxx,xxx,xxx,…,xxx,o(xxx)
S年o月m日