余姓

今天我来考古一下。以历史资料开始，从我们村结束。

根据2005年12月公布的《新百家姓》，余姓人口580万，排在第40位大姓。
在潮汕为第20位大姓，占潮汕人口1.01%，约13.7万人。
饶平县余姓有7万多人。当年高考分数出来后，饶平有个宗亲组织，还给我发了个奖状(只要上了大学都给奖状)。一直被我扔在墙角边。

余姓起于春秋秦之由余，由余做为余姓的开始，已成定论。
由余又从那来？由余是罕的第37世孙。轩辕黄帝之五世孙夏禹的三子少康氏，又名罕，封于无余，为无余王。没想到挖到黄帝，囧rz，再挖就变神话故事了。

由余，原是晋国人，后来避乱，跑到西戎，西戎是比较野蛮的民族，周封秦国就是让秦国在西边挡西戎的。后来由余奉命出使秦国，为秦穆公所谋用，拜为上卿。秦穆公很厉害的一个人，秦国能强大起来，全靠他，打下了西戎，吞并了晋国，有个词叫秦晋之好，有时候别人结婚会用一下，祝人家永结秦晋之好，不过我觉得很奇怪，虽然是通婚了，但后来，女婿打丈人，丈人打女婿，后世怎么会出来这个词，觉得好奇怪啊？！扯远了，囧，回到讲余姓。秦穆公对由余委以重任，由余帮助他谋划征讨西戎，奠定秦国西方霸主的地位（春秋战国，打来打去，其实就是四个国家，北方晋国，西方秦国，南方楚国，东方齐国，其它的小打小闹，没啥作为，中原地带那时候很无能，尽出亡命之徒，搞谋杀）。

后代以由余之名为姓，始有余氏，后来秦大夫百里奚、蹇叔，颂扬由余功绩，为余氏撰写世数辈序。以由余为余姓第一世，辈序诗以一字一世计，延续使用至第52世。 这是由来，且对后世影响很深，全国各地的余姓，几乎都排辈序。不过不同地方的余姓，后来都各自发展，没有统一的辈序。

唐代高宗总章二年（公元669年），河南光州固始人（或说山西人）余元华入闽南粤东地区平“獠乱”而落籍于福建漳州，成为入南方历史较早并且有史可考的余氏先祖之一（根据清代《漳州余氏族谱》记载）。

到了北宋，宋名臣蔡襄为余姓辈序续撰14字:
渊源积厚自流芳，远近攸分尽炽昌。

到了57世“自”字辈广东余姓出了个牛人余靖，广东后世余姓多以余靖为太祖，以余靖为一世，排到我这辈就是35世。近1千年。

下面这个余靖的生平简介来自于百度百科，http://baike.baidu.com/view/220929.htm
余靖（公元1000–1064）本名希古，字安道，号武溪。北宋韶州曲江（今属广东韶关）人。生于宋真宗咸平三年，卒于英宗治平元年，年六十五岁。出身于仕宦之家。少时随舅父师从进士黄正读书，性聪慧，过目不忘，后师林和靖先生，学业大进。天圣二年（公元1024）进士及第，起家为赣县尉。官至朝散大夫，守工部尚书，集贤院学士，知广州军州事兼广南路兵马都铃辖经略安抚使，柱国，始兴郡开国公。赐紫金鱼袋，赠刑部尚书。北宋英宗治平元年（1064年）余靖回京述职途次南京，偶染风寒，病重不起，卒于江宁府秦淮亭，享年六十五岁。英宗闻讯恻然，辍朝一日追赠刑部尚书，谥曰“襄”后人尊称忠襄公。余靖生三子，六女。
　　北宋大文学家欧阳修撰襄公余靖道碑谓：“公为人资重刚劲，而言语恂恂，不见喜怒。自少博学强记，至于历代史记，杂家小说，阴阳律历，外鐕浮屠，老子之书，无所不通。”
　　数上书，建言班固汉书舛谬。命与五洙并校司马迁、范晔二史，﹝按:五洙疑为尹洙之讹﹞奏擢集贤校理。景祐三年（1036），因向皇帝上疏为被贬的礼部员外郎范仲淹辩护，与尹洙、欧阳修同被贬，降职为监筠州酒税。由是益知名。庆历三年（1043），复起任擢升为谏院右正言，专司向皇帝进谏奏事，余正直敢谏，曾多次为建言“轻徭薄赋”整顿户政，去除贪残之吏，抚疲困之民事而向皇帝抗声力争，以致唾液飞溅至皇帝的“龙颜”上仍意犹未尽。他的建议大多为仁宗所接受，因而与欧阳修、王素、蔡襄，同被誉为朝廷敢于进谏的“四谏”，时论重之。庆历四年（1044），受命为出使契丹的使者。他巧妙地运用外交手段折服了雄据一方的辽主，从而适时的在复杂的宋辽夏三角关系中维护了宋朝的利益。
　　三使契丹，习外国语。侬智高叛，经制南事；为帅十年，不载南海一物。广州有八贤堂，靖为其中之一。官至工部尚书。卒，谥襄。靖著有《武溪集》二十卷，见《四库总目》。
————–百度百科转载完毕———————–

余靖的著作也很多，主持校勘《史记》、《汉书》、《后汉书)，着有《三史勘误》四十卷，诗文集《武溪集》等。《粤东诗海》评其诗“骨格清苍，吾粤宋诗无出其右。”《岭南历史诗选》说他的五言诗“有骚雅遗音”。

余靖被尊称为忠襄公，潮汕余姓系此一脉，所以有“无余不忠襄”的说法。余靖的第三子余叔英到潮州当官，余姓开始进入潮汕地区。

从第8世以后，因为以前排的辈序用完了，从此各个余姓分支开始另尊开基始祖，并自创辈序。

传到第14世余克轩，跑到我们饶平海山岛，成为小岛上余姓的始祖。

我们村叫美宅村，全村都姓余。我们村是从第22世开始自己排辈序的。排到第53世:
元衍思文，绍美忠功，斯昌维远，克光其宗，修齐治平，传启和同，丰盈亿万，可继乃中。

我的光字，排在第35位。

讲方言的软件

前几天聊起了方言，然后就上网找到了这个小软件。

功能不多，就是给它一段文字(只要是文本的就行，网页上选定文字，也OK的)，

这个软件能够把它念出来，其实就是带了一个语言库。

我把它放在dropbox上，http://dl.getdropbox.com/u/774072/software/zgfy_setup.zip

听歌也是一个不错的学习途径。

推荐个学习普通话和粤语的网站: http://www.putonghuaonline.com/

use Maxima之初等代数篇

准备写个maxima的小教程。maxima是一个计算机代数系统(Computer Algebra System)，用于公式推导、符号计算、数值计算等，Maxima的前身是Macsyma，Macsyma由MIT在60年代末所开发，用LISP实现。Maxima应该是开源CAS里最广为人知的，不过它的功能不是最强的。最强大的是IBM的Axiom, 此外还有SINGULAR和Yacas等其它几个CAS。

功能最强大不一定是最好的，适合自己的才是最好的。比如你要做一个同学的通讯录，你需要的是excel而不是sql。

maxima可以从sourceforge下载。有Linux、MacOS X和Windows多种版本。本身是文本的界面，有多种图形界面的实现，最初的是xmaxima，我觉得很差，一直在emacs里用maxima，现在有一个界面wxMaxima，感觉做得不错，推荐使用。

首先，它可以做为一个计算器使用。
(%i1) 1+2*3+4/5+6^7;
(%o1)
maxima的命令行i代表输入，o代表输出，语句以;结束，回车后，自动输出。用wxmaxima回车不输出，shift+Enter才输出。
结果可重用。如下，将上面的输出，换成浮点数。
(%i2) float(%o1);
(%o2) 279943.8

命名，使用label，名字和表达式之间用:分隔。
(%i3) radius: 11;
(%o3) 11
(%i4) area: %pi * radius^2;
(%o4)

(%i5) y: 3*x^2+5*x+7;
(%o5) 3*x^2+5*x+7
(%i6) solve(y=1,x);
(%o6)
不单变量可以命名，表达式也可以，(%i5)是命名，不是函数定义。
函数定义的语法是:=
(%i7) y(x) := x^2;
(%o7) y(x):=x²
(%i8) y(4);
(%o8) 16
解方程使用solve命令。
(%i9) solve(y(x)=3, x);
(%o9)

另外一个重要的命令是ev，ev可以强制使用数值输出，还可以定义局部变量。
(%i10) exp(5/3);
(%o10) e^5/3

(%i11) ev(exp(5/3), float);
(%o11) 5.29449005047003
ev是evaluate的缩写，float参数表示使用浮点数输出，还有另一个参数numer，float是numer的一种类型，使用float不会对numerval的变量求值。
(%i12) numerval(b,5);
(%o12) [b]
(%i13) a:3;
(%o13) 3
(%i14) ev(a*b, float);
(%o14) 3*b
(%i15) ev(a*b, numer);
(%o15) 15

ev还可以用于定义局部变量。下面的例子a,b,c,d都是局部变量。
(%i16) ev(solve(a*x^2+b*x+c=d,x),a=3,b=4,c=5,d=6);
(%o16)

数列求和:
(%i17) sum(k^3,k,1,10);
(%o17) 3025
在表达示前加上’将会列出相应的表达式子，而不会进行求值。
(%i17) ’sum(k^3,k,1,10);
(%o17)

先写到这，下次写微积分。
另外一个小小的tip，输出我用latex写的，blog上看到的是图片，但是选中图片，复制，然后到文本编辑器里粘贴，会看到latex的代码。

Fences-桌面图标分类工具

使用不同的操作系统会有不同的操作习惯，用linux我喜欢用终端，用windows我喜欢把东西都往桌面上丢。时间长了，桌面看起来一团糟。

如果你也跟我有一样的习惯，显然我们需要分类和组织桌面图标。

Fences便是这样一个小工具，它可以将不同的图标放到不同的半透明“容器”当中，你还可以移动和拉伸这个“容器”。双击桌面，隐藏所有图标；再次双击桌面，显示所有图标。没其它功能了。

Fences可以从官方站点上下载，免费的。

http://www.stardock.com/products/fences/downloads.asp

截屏一张：

用Taylor公式估计函数值

最近在看微积分，刚看到Taylor theorem，本科用的是阉割版的农科教材，竟然没讲这么个超级经典的方法。

 efv <- function(f, value,variable="x", a=0, eps=0.001) {
     #estimate function value using Taylor theorem
     assign(eval(variable), a)
     fv.old <- eval(f)
     k <- 1     
     repeat {
         df <- D(f, variable)
         if (df == 0)
             break
         fv.new <- fv.old + eval(df)*(value-a)^k/factorial(k)
         if (fv.new - fv.old < esp)
             break
         fv.old <- fv.new
         f <- df
         k <- k+1
     }
     return (fv.new)
 }

#估算一下e值，跟系统函数算的差不多。
> efv(expression(exp(x)),1)
[1] 2.718254
> exp(1)
[1] 2.718282