Asymptotic Freedom

gnuplot所有的colornames及范例

2012-01-17T14:39:00.000-08:00

这里列出gnuplot 4.4所有colorname及其代表的颜色：

white black dark-grey red web-green web-blue dark-magenta dark-cyan dark-orange dark-yellow royalblue goldenrod dark-spring-green purple steelblue dark-red dark-chartreuse orchid aquamarine brown yellow turquoise grey0 grey10 grey20 grey30 grey40 grey50 grey60 grey70 grey grey80 grey90 grey100 light-red light-green light-blue light-magenta light-cyan light-goldenrod light-pink light-turquoise gold green dark-green spring-green forest-green sea-green blue dark-blue midnight-blue navy medium-blue skyblue cyan magenta dark-turquoise dark-pink coral light-coral orange-red salmon dark-salmon khaki dark-khaki dark-goldenrod beige olive orange violet dark-violet plum dark-plum dark-olivegreen orangered4 brown4 sienna4 orchid4 mediumpurple3 slateblue1 yellow4 sienna1 tan1 sandybrown light-salmon pink khaki1 lemonchiffon bisque honeydew slategrey seagreen antiquewhite chartreuse greenyellow gray light-gray light-grey dark-gray slategray gray0 gray10 gray20 gray30 gray40 gray50 gray60 gray70 gray80 gray90 gray100

数码单反（DSLR）相机raw格式图片的读取和处理

2011-05-27T14:53:00.000-07:00

随着消费用数码单反（DSLR）相机成本的下降和性能的提升，在很多要求不太高的科研场合DSLR相机可以用作昂贵的高档CCD数码相机的替代品。现在最新的数码相机可以提供14位A/D转换，而动态范围可以到达11档光圈以上【参考链接：1，2】，很多时候这已经可以满足要求了。一般摄影爱好者使用的JPEG图片格式只能存储8位的层次，而且JPEG图片压缩过程中会有失真，所以作为科研使用的DSLR相机必须使用raw图片格式。

一般各个消费用相机厂商都使用自己专有的raw格式，所以处理起来不像科研级相机那么直接。dcraw是一个Linux平台上的能将各主要品牌相机的raw格式图片转为Netpbm格式图片的小工具。Netpbm提供了一些工具软件，另外很多图像处理软件也包含对Netpbm格式的支持，所以可以在此基础上对图片进行进一步处理。

很多时候我们可能希望自己写代码来处理图片，因为：

很多图像处理软件只支持8位（或24位RGB）而不支持16位图片；
科研上用到的算法大都比较复杂，现有的软件不满足需求。

Netpbm可以把图片存成二进制文件或者纯文本（ascii）文件。无论哪种方式，开始的几行都是标识文件格式的纯文本信息，后面跟着是文件内容。

纯文本文件的好处是数值大小一目了然，但是文件比较大。二进制文件比较小，但是处理的时候需要注意：Netpbm存储文件时使用的是Big-endian方式，而一般x86机器都是使用Little-endian方式，所以在读取数据时要进行转换。glibc的endian.h提供这样的转换函数。二进制文件存储顺序是从上到下按行存储，每行内按从左到右顺序存储。

dcraw转换图片时默认使用8位方式，如果要使用16位方式，需要使用参数 -6。

《Maxima快速参考手册1.1》

2011-05-24T23:54:00.000-07:00

这是之前写的一份Maxima中文文档的更新版本。

Vim笔记（二）

2011-04-15T16:18:00.000-07:00

分割窗口
:split file2    把窗口分割成上下两半，在新开窗口打开文件file2
:new    上下分割窗口，在新开窗口打开新文件
CTRL-W w    切换窗口
CTRL-W +    增大窗口尺寸
CTRL-W -    缩小窗口尺寸
:vsplit file2    左右分割窗口，在新开窗口打开文件file2
CTRL-W h    切换到左边窗口，其他方向按j, k, l类推
CTRL-W H    移动窗口到最左边，其他方向按J, K, L类推
:qall    退出所有窗口，类似的有:wall
反复执行复杂命令
q{register}    开启寄存器，开始记录命令，register可以是a到z中任何一个
q    结束记录命令
@{register}    执行寄存器{register}内所有命令
查找替换一段字符
:[range]s/from/to/[flags] 将from替换为to，[range]控制替换范围，[flags]表示一些控制参数
常用[range]
- . 当前行
- $ 文件最后一行
- 1,5 第一行到第五行
- % 全部文件
- .+3 当前行之下第三行
常用[flags]
- g 替换所有匹配字符串，否则只替换第一个匹配
- c 执行每一个替换之前寻求确认
可视化块模式
CTRL-V    开启可视化块模式，用h,j,k,l选择块
I{string}<Esc>    在块的左边每一行插入一段字符串
A{string}<Esc>    在块的右边每一行加入一段字符串
c{string}<Esc>    替换一个块的文字
U    换为大写字母
u    换为小写字母
~    大小写互换
数行合并
J 数行合并，去除换行符

Vim笔记（一）

2011-04-14T22:54:00.000-07:00

常用的Vim命令小结，最基本的就不说了，例如如何安装，如何启动和退出，两种基本模式（正常模式和插入模式），用h,j,k,l移动光标，简单的复制粘贴等等。

移动光标

以字符为单位
f 右移至指定字符，例如fh移动至右边第一个h
F 左移至指定字符
以词为单位
w    右移至词首
b    左移至词首
e    右移到词尾
ge    左移到词尾
以行为单位
0    左移至行首
^    左移至行首第一个非空字符
$    右移至行尾
gg    移动至第1行
G    移动至最后1行
nG    移动至第n行
n%    移动至n%位置
以屏为单位
H    移动至屏幕顶端
M    移动至屏幕中间
L    移动至屏幕底端
CTRL-U    向上滚半屏
CTRL-D    向下滚半屏
CTRL-F    向下滚一整屏
CTRL-B    向上滚一整屏

简单搜索
/string    向后搜索字符串
?string    向前搜索字符串
n    下一个匹配字符串
N    上一个匹配字符串
/\<word\>    匹配整个单词，\<和\>分别匹配词首和词尾
%    寻找匹配的括号
书签
ma 标记此处为a，共可使用a到z二十六个书签
`a 回到标记为a处

几个关于Gnuplot小技巧的网页

2011-03-04T12:25:00.000-08:00

Gnuplot手册不能告诉你的小技巧。Gnuplot也可以画出很惊艳的图片。

not so Frequently Asked Questions
http://t16web.lanl.gov/Kawano/gnuplot/index-e.html
Gnuplot tricks blog
http://gnuplot-tricks.blogspot.com/
Impossible gnuplot graphs
http://www.phyast.pitt.edu/~zov1/gnuplot/html/intro.html
Gnuplot surprising
http://gnuplot-surprising.blogspot.com/
Gnuplotting
http://www.gnuplotting.org/

Octave笔记

2011-01-12T14:24:00.000-08:00

Octave可以看作是开放源代码版本的Matlab，它的命令几乎和Matlab完全一样，因此可以用Matlab的书籍和文档来学习Octave。另外，绝大多数Matlab的代码也可以不经修改在Octave中运行。对于初学者来说，下面这些笔记可以帮助快速上手Octave。更多的功能，可以参考Octave文档或者其他Matlab的资料。

脚本文件，以下面一行开始：
```
#!/usr/bin/octave -qf
```
注释，以%或者#开始
虚数单位：i, j, I, J都可以两种方法定义复数：3 + 4i或者complex(3, 4)

定义矩阵：

octave:1> a = [1, 2; 3, 4]
a =

1   2
3   4

随机数：rand()或者rand(m, n)

定义一段连续序列：

octave:1> 1:5
ans =

1   2   3   4   5

还可以指定步长：

octave:2> 1:2:10
ans =

1    3    5    7    9

字符串：单引号或者双引号

定义Data Structure：

octave:1> x.name = ”Yusuf”
octave:2> x.age = 33
octave:3> x.phone = “46103”
octave:4> x
x =
{
name = Yusuf
age = 33
phone = 46103
}

定义Cell Array：

octave1:> a = {”random matrix”, rand(2)}
a =

{
[1,1] = random matrix
[1,2] =

0.417140   0.460380
0.337105   0.078835

}

octave2:> a{2}(1,2)
ans =  0.460380

声明全局变量：
```
global c = 3e8
```

算术符号

`x + y`	加法
`x - y`	减法
`x * y`	矩阵乘法
`x .* y`	元素相乘
`x / y`	右除，相当于`(inverse (y’) * x’)’`
`x ./ y`	元素右除
`x \ y`	左除，相当于`inverse (x) * y`
`x .\ y`	元素左除，每个y元素被相应的x元素除
`x ^ y` 或 `x ** y`	乘方
`x .^ y` 或 `x .** y`	元素乘方
`x’`	转置复共轭
`x.’`	转置

比较符号

x < y 小于

x <= y 小于等于

x == y 等于

x >= y 大于等于

x > y 大于

x != y 或 x ~= y 不等于
if语句

形式一：
```
if (condition)
  then-body
endif
```

形式二：

if (condition)
  then-body
else
  else-body
endif

形式三：

if (condition)
  then-body
elseif (condition)
  elseif-body
else
  else-body
endif

switch语句

switch expression
  case label
    command_list
  case label
    command_list
  …
  otherwise
    command_list
endswitch

while语句
```
while (condition)
  body
endwhile
```
do-until语句
```
do
  body
until (condition)
```
for语句
```
for var = expression
  body
endfor
```
定义函数

形式一：无自变量 function name body endfunction
形式二：带自变量 function name (arg-list) body endfunction
形式三：单个返回值 function ret-var = name (arg-list) body endfunction
形式四：多个返回值 function [ret-list] = name (arg-list) body endfunction

可变自变量长度函数 function val = smallest (varargin) val = min ([varargin{:}]); endfunction
函数Handles（指向函数的指针） @function-name
匿名函数 @(argument-list) expression
终端输出 disp (x)
写读包含分隔符的数据文件 dlmwrite (file, a) dlmwrite (file, a, delim, r, c) dlmwrite (file, a, key, val ...) dlmwrite (file, a, "-append", ...) a是待写矩阵，delim是分隔符，默认为逗号。r, c分别为起始的空行和空列。key包含下列参数：

"append" 即上面的"-append"，"on"或者"off"
"delimiter" 即上面的delim
"newline" 换行符
"roffset" 即上面的r
"coffset" 即上面的c
"precision" 精确度，可以用类似fprintf的格式字符串，或者有效数字位数

data = dlmread (file) data = dlmread (file, sep, range)

range为四元素的矢量：[R0, C0, R1, C1]，为数据块左上角和右下角的标记，从0开始

Octave可以使用大多数类似C的I/O函数（fprintf，fopen等）

常用数学函数

exp (x)	指数函数	sinh (x)	双曲正弦
log (x)	自然对数	cosh (x)	双曲余弦
log10 (x)	以10为底的对数	tanh (x)	双曲正切
log2 (x)	以2为底的对数	coth (x)	双曲余切
sqrt (x)	平方根	sech (x)	双曲正割
abs (z)	绝对值	csch (x)	双曲余割
arg (z)	辐角	asinh (x)	反双曲正弦
conj (z)	复共轭	acosh (x)	反双曲余弦
real (z)	复数实部	atanh (x)	反双曲正切
imag (z)	复数虚部	acoth (x)	反双曲余切
sin (x)	正弦	asech (x)	反双曲正割
cos (x)	余弦	acsch (x)	反双曲余割
tan (x)	正切	sum (x)	求和
cot (x)	余切	prod (x)	求积u
sec (x)	正割	ceil (x)	不小于x的最小整数
csc (x)	余割	floor (x)	不大于x的最大整数
asin (x)	反正弦	factorial (n)	阶乘
acos (x)	反余弦	max (x)	最大值
atan (x)	反正切	min (x)	最小值
acot (x)	反余切	round (x)	四舍五入
asec (x)	反正割	sign (x)	符号函数
acsc (x)	反余割
atan2 (y, x)	计算atan (y / x)
sind (x),…	以度数为单位的三角函数

特殊函数

[a, ierr] = airy (k, z, opt)	Airy函数
[j, ierr] = besselj (alpha, x, opt)	第一类Bessel函数
[y, ierr] = bessely (alpha, x, opt)	第二类Bessel函数
[i, ierr] = besseli (alpha, x, opt)	变形第一类Bessel函数
[k, ierr] = besselk (alpha, x, opt)	变形第二类Bessel函数
[h, ierr] = besselh (alpha, k, x, opt)	第一类（k=1）和第二类（k=2）Hankel函数
beta (a, b)	Beta函数
gamma (z)	Gamma函数
bincoeff (n, k)	二项式系数
erf (z)	误差函数
legendre (n, x, normalization)	n阶Legendre函数，m = 0…n

线性代数
det (a) 行列式
eig (a) 本征值
dot (x, y) 点乘
inv (a) 逆矩阵
解非线性方程 [x, fval, info] = fsolve (fcn, x0, options) fcn是包含待解方程组的函数，x0是自变量初始值，x为返回的数值解，fval为返回的函数值，如果结果收敛，info为1，否则为其他值。例子： function y = f (x) y(1) = -2*x(1)^2 + 3*x(1)*x(2) + 4*sin(x(2)) - 6; y(2) = 3*x(1)^2 - 2*x(1)*x(2)^2 + 3*cos(x(1)) + 4; endfunction [x, fval, info] = fsolve (@f, [1; 2]) 结果： x = 0.57983 2.54621 fval = 6.1872e-08 -3.2708e-07 info = 1
数值积分（各函数对应不同算法） [v, ier, nfun, err] = quad (f, a, b, tol, sing) quadl (f, a, b, tol) quadgk (f, a, b, prop, val, ...) quadl (f, a, b, tol)

统计笔记（一）

2010-01-21T12:30:00.000-08:00

协方差（covariance）定义：

$Cov (Y_{1}, Y_{2}) = E [(Y_{1} - μ_{1}) (Y_{2} - μ_{2})]$

其中 $E ()$ 表示期望值， $μ_{1} = E (Y_{1})$ ， $μ_{2} = E (Y_{2})$ 。
相关系数（coefficient of correlation）定义：

$ρ = \frac{Cov (Y_{1}, Y_{2})}{σ_{1} σ_{2}}$

其中 $σ_{1}$ 和 $σ_{2}$ 分别是 $Y_{1}$ 和 $Y_{2}$ 的标准差。
任意分布样本的平均值
$Y_{1}$ ， $Y_{2}$ ，...， $Y_{n}$ 是 $E (Y_{i}) = μ$ ， $V (Y_{i}) = σ^{2}$ 的独立随机变量，则平均值

$\bar{Y} = \frac{1}{n} \sum_{i = 1}^{n} Y_{i}$

的期望值 $E (\bar{Y}) = μ$ ，方差 $V (\bar{Y}) = σ^{2} / n$ 。
正态分布样本的平均值
$Y_{1}$ ， $Y_{2}$ ，...， $Y_{n}$ 是来自正态分布（期望值为 $μ$ ，方差为 $σ^{2}$ ）的样本，其平均值

$\bar{Y} = \frac{1}{n} \sum_{i = 1}^{n} Y_{i}$

也呈正态分布，期望值 $μ_{\bar{Y}} = μ$ ，方差 $σ_{\bar{Y}}^{2} = σ^{2} / n$ 。
$Y_{1}$ ， $Y_{2}$ ，...， $Y_{n}$ 同上，那么

$\sum_{i = 1}^{n} Z_{i}^{2} = \sum_{i = 1}^{n} {(\frac{Y_{i} - μ}{σ})}^{2}$

呈自由度为 $n$ 的 $χ^{2}$ 分布。
$Y_{1}$ ， $Y_{2}$ ，...， $Y_{n}$ 同上，那么

$\frac{(n - 1) S^{2}}{σ^{2}} = \frac{1}{σ^{2}} \sum_{i = 1}^{n} (Y_{i} - \bar{Y})^{2}$

呈自由度为 $n - 1$ 的 $χ^{2}$ 分布。其中，

$S^{2} = \frac{1}{n - 1} \sum_{i = 1}^{n} (Y_{i} - \bar{Y})^{2}$

是样本方差。 $S^{2}$ 和 $\bar{Y}$ 为独立随机变量。
学生 $t$ 分布(简称 $t$ 分布)定义：
如果 $Z$ 是正态分布随机变量， $W$ 是自由度为 $ν$ 的 $χ^{2}$ 分布随机变量，如果 $Z$ 和 $W$ 互相独立，则

$T = \frac{Z}{\sqrt{W / ν}}$

呈自由度为 $ν$ 的 $t$ 分布。
在上面的定义中，如果让 $Z = \sqrt{n} (\bar{Y} - μ) / σ$ ， $W = (n - 1) S^{2} / σ^{2}$ ，则由上面第4和第6条结论可以得出

$\sqrt{n} (\frac{\bar{Y} - μ}{S})$

呈自由度为 $n - 1$ 的 $t$ 分布。
$F$ 分布定义：
如果 $W_{1}$ 和 $W_{2}$ 为独立 $χ^{2}$ 分布变量，自由度分别为 $ν_{1}$ 和 $ν_{2}$ ，则

$F = \frac{W_{1} / ν_{1}}{W_{2} / ν_{2}}$

被称作具有 $ν_{1}$ 分子自由度和 $ν_{2}$ 分母自由度的 $F$ 分布。
中心极限定理（Central Limit Theorem）
$Y_{1}$ ， $Y_{2}$ ，...， $Y_{n}$ 是独立的具有同样分布的随机变量，期望值 $E (Y_{i}) = μ$ ，方差 $V (Y_{i}) = σ^{2} < \infty$ 。定义

$U_{n} = \sqrt{n} (\frac{\bar{Y} - μ}{σ})$

其中

$\bar{Y} = \frac{1}{n} \sum_{i = 1}^{n} Y_{i}$

那么，当 $n \to \infty$ 时， $U_{n}$ 趋近于正态分布。

写了一份Maxima的中文文档

2010-01-18T16:56:00.000-08:00

以前我为计算机代数系统Maxima写了一些简介，这次抽时间写了一个比较完整的中文文档，希望能为开放源代码软件的推广做点贡献。

《Maxima快速参考手册》：

写的仓促，欢迎反馈。

在eps图片上添加LaTeX公式

2009-05-11T01:16:00.000-07:00

利用MetaPost的exteps宏包，可以在metapost文件中加入eps图片，同时利用metapost的label命令在图片任何地方加入LaTeX公式。下面是一个使用exteps的例子：

input exteps
verbatimtex
\documentclass{article}
\begin{document}
etex
prologues:=2;
beginfig(1);
begineps "figure.eps";
endeps;
label.rt(btex $f(x) = \sqrt x$ etex, (2cm,3cm));
endfig;
end

为了使用LaTeX命令，需要设置环境变量

TEX=latex

GnuPlot在PostScript终端下使用任意字体

2009-04-02T15:00:00.000-07:00

有时候我们在GnuPlot中想使用一些其他字体，例如LaTeX字体或者一些艺术字体，这时候我们可以用fontfile选项把自己的字体内嵌进生成的PostScript文件：

set term postscript fontfile "myfont.ext" "FontName"

这里的字体可以是Type 1字体，也可以是TrueType字体。命令必须包含完整字体文件名，包括扩展名。字体文件可以在当前目录，也可以用GNUPLOT_FONTPATH环境变量控制字体路径。FontName（字体名）不是字体文件名，通常不那么容易直接得到，尤其是TrueType字体。最简单的办法是在GnuPlot里执行

set term postscript fontfile "myfont.ext"

GnuPlot会自动给出字体名。

这样，由于可以使用LaTeX字体，配合PostScript终端的enhanced选项，我们几乎可以在GnuPlot中输入任意特殊字符。另外，为使用TrueType字体，需要安装ttf2pt1工具进行字体自动转化。

sed编辑器的使用

2009-03-16T12:06:00.000-07:00

sed被称为流编辑器（stream editor）。它和一般写字板之类的文本编辑器不同，它通过逐行读入文本文件，按照用户指定的方式对文本进行过滤。用户通过脚本命令告诉sed该如何处理文本文件，所以sed可以看作一种程序语言。这看起来没有使用鼠标和键盘编辑文件的编辑器方便，但是当我们需要自动处理大量文本文件时，这会非常方便。例如，我需要把很多数据文件的某些行注释掉，也就是在所有数据文件的某几行开始处加上字符“# ”。如果一个个文件编辑，工作量就太大了。下面这个命令可以很方便的实现这个功能：

sed -i '1,5 s/^/# /' *.txt

这个命令把每一个.txt作为后缀的文件前5行开始处加上“# ”。下面把常用的一些sed命令总结一下。

替换一段文字

sed

sed 's/a/A/' <old >new

old

new

old

sed 's/a/A/g' <old >new

/

sed

执行多个命令

-e

sed -e 's/a/A/' -e 's/b/B/' <old >new

指定对某些行操作

sed '1,10 s/A/a/'

另外，这里有一个很有用的sed常用单行命令参考。

GnuPlot使用png terminal时找不到字体的问题

2009-01-09T12:19:00.000-08:00

GnuPlot不同的终端（terminal）使用不同的处理字体的方法。png/jpeg/gif使用libgd，所以字体目录由GDFONTPATH这个环境变量控制，可以直接使用TrueType字体。例如把需要用到的字体放在$HOME/fonts目录下，在.bashrc中添加：

export GDFONTPATH=$HOME/fonts

当路径设置好后，可以用GNUPLOT_DEFAULT_GDFONT这个环境变量设置默认字体。

在USB drive（U盘）上安装Linux

2009-01-02T13:17:00.000-08:00

用UNetbootin可以很方便的把Linux装到U盘上，只要主板支持，就可以从U盘启动Linux。由于U盘可读写，比Live CD更方便。目前UNetbootin支持Fedora，Debian，Ubuntu等各种主流Linux版本。如果你的U盘比较小，可以试试Puppy Linux，挺小巧的一个Linux版本，只有100M左右大小。

Mie散射计算器

2008-11-21T11:59:00.000-08:00

Mie散射理论为球形颗粒对于电磁波的散射提供了一套严格解，其数值计算十分复杂。Bohren和Huffman的书里提供了用于数值计算的Fortran代码。另外，Oregon Medical Laser Center的网页提供了一个在线的Mie散射计算器，可以快速的计算出一些常用参数，并且给出一些图线。

GNU screen的几个常用命令

2008-11-19T21:46:00.000-08:00

GNU screen是命令行界面下的窗口管理器，允许用户在命令行界面下同时在几个虚拟窗口运行不同的程序。GNU screen有一个很有用的功能：如果用户登录到远程机器上利用screen运行程序，该程序在用户切断连接之后可以继续运行，并且下次用户可以登录继续上次的操作。下面是一些常用的命令。

首先在终端执行screen命令进入程序。执行之后虽然看不到任何变化，其实screen已经启动，我们看到的就是一个虚拟终端。

建立一个新的窗口（虚拟终端）：C-a c
也就是按Ctrl+a然后再按c
切换到之后的窗口：C-a n
切换到之前的窗口：C-a p
切换的第N个窗口：C-a N
显示现有窗口列表，用上下箭头（或者类似于vi中用j和k）移动光标，按Enter选择：C-a "
改变当前窗口的名字：C-a A
关闭目前的screen但是保留其中运行的程序（detach）：C-a d
如果直接关闭窗口也可以有同样作用。
运行screen继续上次在screen里运行的程序（reattach）：screen -r
进入拷贝模式：C-a [
在screen里不能向上翻屏，进入拷贝模式可以使用光标上下移动翻屏。用空格键可以选择要拷贝的文字，第一次按空格选择开始的字符，第二次选择结束的字符。
粘贴选择的文字：C-a ]

Linux下用sopcast看电视

2008-11-16T03:39:00.000-08:00

Sopcast是一个基于P2P技术的网络电视软件，在Linux和Windows上都可以运行。Sopcast网站提供的Linux版本倒挺简单，直接下载解压就可以用了。但是这只是一个命令行的软件。有几个图形界面可以用，觉得gsopcast不错，但是源代码在Fedora和Debian上编译都有问题，header文件需要改动几处。可能因为大家用的编译器不一样？编译好的文件在这里，是我在Debian下用gcc4.3编译的，gsopcast版本是0.4.0，可以直接下载下来用，应该不依赖特定的发行版。

如果自己在桌面建立一个gsopcast快捷方式图标的话，似乎当从图标启动的时候，$PATH变量不起作用，所以gsopcast可能找不到sopcast的执行文件。另外，新的sopcast执行文件是sp-sc-auth，而不是以前的sp-sc。要保证正常使用，最好在/usr/local/bin目录建立一个符号链接，取名作sp-sc。

光圈档数（aperture stop）的计算方法

2008-11-10T17:31:00.001-08:00

摄影的时候计算曝光量经常用到几“档”曝光这个术语。这个曝光的“档数”（stop）也经常被称作EV（exposure value）。这个“档数”到底怎么算出来的？曝光量是由光圈和快门共同决定的。曝光量增加一倍，就叫做曝光增加一档。当光圈一定时，曝光量和快门时间（也就是曝光时间）成正比，例如1/30秒就比1/60秒曝光增加一倍，也就是增加一档。这比较容易理解。比较容易迷惑的是光圈。光圈一般用F数（F number或者F/#）来表示，其定义为：

f / # = \frac{f}{D}

其中f是镜头焦距，D是镜头入射孔径。如果拿起一个镜头，从前面看进去，所看到的那个圆形小孔的直径，就是镜头的入射孔径了。对于一般相机镜头，这个孔径大小是可调的，由f/#的设置来控制。有些人可能会想，既然计算曝光量，为什么不直接拿孔径大小D来算，还要定义一个什么f/#。这是因为，感光器件（胶片或者CCD）的曝光量不仅和孔径大小有关，还和焦距有关。用f/#能更准确的表示曝光量。但是这有一个前提，就是物体离镜头比较远，如果物体离镜头很近（例如显微镜），就不能用f/#来算曝光量了。这时候一个更常用的量是数值孔径（numerical aperture或NA）。因为孔径D在分母上，所以f/#数值越小，光圈越大。例如光圈f/2比f/4要大。但是f/2并非比f/4大一档光圈，而是二档。因为镜头通光量和孔径面积成正比，也就是和孔径大小的平方成正比。所以，比f/4大一档的光圈是

\frac{4}{\sqrt{2}} ≅ 2.8

实际上还会经常碰到分数表示，例如半档，1/3档。对于普遍情况下，可以用下面的公式计算两个光圈究竟差几档：

N = 2 \log_{2} \frac{f / # 2}{f / # 1}

例如，f/1.4和f/1.8相差

N = 2 \log_{2} \frac{f / 1.8}{f / 1.4} ≅ \frac{2}{3}

LaTeX CJK生成PDF文件的中文搜索以及复制粘贴

2008-11-04T10:53:00.003-08:00

以前的blog介绍过如何在Fedora下安装配置LaTeX CJK，其他的Linux发行版也可以用类似的方法。但是如果直接用ps2pdf或者dvipdfm生成PDF文件，这样的PDF文件没办法搜索中文，也不能复制粘贴中文内容。可以用dvipdfmx将生成的dvi文件转换为PDF，这样生成的PDF可以完成上述功能。另外，我把自己的LaTeX CJK的字体和配置文件放到了网上，地址是http://www.mediafire.com/file/dwutwxjjmdn/texmf.tar.gz，有需要的朋友可以去下载。这个文件包含CJK包和三个字体，分别是Fedora自带的文鼎宋体，楷体，还有文泉驿的黑体，应该满足大多数应用。下载后解压到$HOME目录，执行texhash $HOME/texmf。不需要进行任何其他配置。下面的文件可以用来测试安装结果：

\documentclass{article}
\usepackage{CJKutf8}
\begin{document}
\begin{CJK}{UTF8}{song}
CJK 测试。

\emph{CJK 测试。}

\textbf{CJK 测试。}
\end{CJK}
\end{document}

用WebPlayer在网页里播放mp3

2008-10-31T22:30:00.003-07:00

WebPlayer利用Flash技术，允许你在你的网页里直接播放mp3文件，而不需要任何其他软件。WebPlayer使用非常简单，只要下载一个很小的文件放在你的网站上，把要播放的mp3文件也放上去，再在你的网页的HTML文件里加几行代码，就可以了。下面是我用口琴演奏的《天鹅》。

祝老婆纪念日快乐！

光束质量因子M2（M squared）的测量

2008-10-29T16:53:00.027-07:00

各种文献经常会用M²（M squared）因子来评价光束的质量。到底什么是M²？一般文献会告诉你，M²就是一个光束的腰半径和发散角的乘积跟TEM₀₀高斯光束相比的比值。M²≥1，因为TEM₀₀高斯光束的该乘积具有最小值。要想从实验得到M²的值，不可避免的要涉及光束半径测量的问题。

光束半径

f (x) = a e^{- \frac{x^{2}}{2 σ^{2}}}

I (x) = e^{- \frac{2 x^{2}}{w_{x}^{2}}}

₁

₂

D = | x_{2} - x_{1} |

₀₀

TEM₀₀高斯光束

误差函数

\frac{\int_{- \frac{D}{2}}^{\frac{D}{2}} e^{- \frac{2 x^{2}}{w^{2}}} d x}{\int_{- \infty}^{\infty} e^{- \frac{2 x^{2}}{w^{2}}} d x} = \erf (\frac{\sqrt{2} D}{2 w})

D = 1.2816 w = 2.5631 σ

₀₀

M²因子

₀₀

w^{2} (z) = w_{0}^{2} + {(\frac{λ}{π w_{0}})}^{2} (z - z_{0})^{2}

₀

₀₀

w^{2} (z) = w_{0}^{2} + M^{4} \times {(\frac{λ}{π w_{0}})}^{2} (z - z_{0})^{2}

刀锋法测量

₀₀

Siegman的一篇论文

₀₀

让blogger支持MathML公式输出

2008-10-28T15:58:00.004-07:00

我想在我的blog里显示公式怎么办？我可以把每个公式变成图片文件，然后把公式作为图片显示。我不需要截屏软件，OpenOffice可以很方便的把包含公式的文档转换为包含图片的HTML文件。但是这么做太笨了。更好的办法是用MathML在网页上显示公式。 MathML是XML的一个应用，是W3C推荐使用的网页数学公式显示方案。遗憾的是，目前不是所有的浏览器都支持MathML。Firefox本身包含对MathML的支持，但是你可能需要安装显示数学公式需要的字体。如果你使用Linux，你的发行版很可能包含（甚至已经安装）了该字体，请以mathml为关键词搜索以确定是否已经安装。如果你不幸使用Windows，请到这里按照说明自行安装字体。如果你更不幸使用IE浏览器，那只好自生自灭了。（开个玩笑，请安装MathPlayer插件）。现在你可以到这里测试一下你的浏览器是不是可以正确显示MathML公式了。能显示，很好。可是，blogger并不能自动支持MathML的显示，因为blogger的网页扩展名都是html，浏览器不会自动按照xml来显示。这个问题可以这么解决：到这里下载一个小javascript文件，把这个文件直接放到你的blogger template文件的head里面，或者放到互联网上某个地方，在template的head里面加入下面一行：

<script type="text/javascript" src="javascript文件的链接"></script>

这样blogger就可以显示你的MathML公式了。 MathML语法非常复杂，要想直接用MathML写出一个公式并不容易。TtM是个把LaTeX公式转换成MathML代码的小工具。这下就方便多了。把转换好的MathML代码复制到你的blog就可以了。下面是一个MathML显示公式的例子：

U (r) = A_{0} \frac{W_{0}}{W (z)} \exp [- \frac{ρ^{2}}{W^{2} (z)}] \exp [- jkz - jk \frac{ρ^{2}}{2 R (z)} + j ζ (z)]

不用重启切换runlevel

2008-10-21T23:15:00.003-07:00

Linux下telinit命令允许不重启的情况下切换runlevel。

魔术组合键SysRq

2008-10-19T00:23:00.014-07:00

很多用户不熟悉SysRq键到底有什么用，甚至很多用户根本没注意过电脑键盘上还有这么一个键存在。在Linux系统出现问题时，只要内核没有完全崩溃，SysRq就可以为我们提供一个诊断和控制内核的途径。SysRq键的使用方法是同时按下：

Alt + SysRq + CommandKey

有些键盘无法同时处理那么多键，例如一些笔记本的键盘，很可能还需要按下额外的Fn键。这时候可以按照下面的顺序：

按住Alt不放 => 按下并放开SysRq => 按下CommandKey => 放开所有键

Fedora默认是不开启SysRq功能的。如果想开启这项功能，可以修改/etc/sysctl.conf文件，把kernel.sysrq设置为1：

kernel.sysrq = 1

常用的CommandKey有如下这些：

CommandKey	功能
`m`	目前内存信息
`t`	目前任务信息
`s`	同步所有挂载文件系统
`u`	重新挂载所有文件系统为只读
`b`	立即重启（不sync或者unmount）
`r`	键盘从raw模式切换到XLATE模式
`k`	Secure Access Key杀死现在终端上所有程序
`c`	故意让系统崩溃
`e`	向除init外所有进程发出SIGTERM信号
`i`	向除init外所有进程发出SIGKILL信号

如果syslogd仍然起作用，SysRq触发结果将会保存在/var/log/messages。

如何使用Maxima（四）

2008-09-30T14:30:00.004-07:00

矩阵运算

输入矩阵

entermatrix(m,n)

m=n

(%i1) m: entermatrix(2,2);

Is the matrix  1. Diagonal  2. Symmetric  3. Antisymmetric  4. General
Answer 1, 2, 3 or 4 :
4;
Row 1 Column 1:
a;
Row 1 Column 2:
1;
Row 2 Column 1:
0;
Row 2 Column 2:
b;

Matrix entered.
                                  [ a  1 ]
(%o1)                             [      ]
                                  [ 0  b ]

matrix

(%i2) matrix([a,b,c],[1,2,3],[d,e,f]);
                                  [ a  b  c ]
                                  [         ]
(%o2)                             [ 1  2  3 ]
                                  [         ]
                                  [ d  e  f ]

矩阵的转置，逆和行列式

(%i3) transpose(m);
                                   [ a  0 ]
(%o3)                              [      ]
                                   [ 1  b ]
(%i4) invert(m);
                                 [ 1     1  ]
                                 [ -  - --- ]
                                 [ a    a b ]
(%o4)                            [          ]
                                 [      1   ]
                                 [ 0    -   ]
                                 [      b   ]
(%i5) determinant(m);
(%o5)                                 a b

矩阵的代数运算

(%i6) m . %o4;
                                   [ 1  0 ]
(%o6)                              [      ]
                                   [ 0  1 ]
(%i7) %o6 + %o4;
                               [ 1         1  ]
                               [ - + 1  - --- ]
                               [ a        a b ]
(%o7)                          [              ]
                               [        1     ]
                               [   0    - + 1 ]
                               [        b     ]

提取矩阵的某一行或某一列

(%i8) row(m,1);
(%o8)                              [ a  1 ]
(%i9) col(%o7,2);
                                   [    1  ]
                                   [ - --- ]
                                   [   a b ]
(%o9)                              [       ]
                                   [ 1     ]
                                   [ - + 1 ]
                                   [ b     ]

为矩阵增加一行（列）或几行（列）

(%i10) addcol(m,[0,1],[c,1]);
                                 [ a  1  0  c ]
(%o10)                           [            ]
                                 [ 0  b  1  1 ]
(%i11) addrow(%o10,[1,0,0,1]);
                                [ a  1  0  c ]
                                [            ]
(%o11)                          [ 0  b  1  1 ]
                                [            ]
                                [ 1  0  0  1 ]

另外，+, -, *, /, ^, exp(x)等运算都可以应用于矩阵，这些运算符分别作用于每一个元素。

(%i12) m/c;
                                   [ a  1 ]
                                   [ -  - ]
                                   [ c  c ]
(%o12)                             [      ]
                                   [    b ]
                                   [ 0  - ]
                                   [    c ]

Vim输入特殊字符（例如希腊字母）

2008-09-25T16:43:00.004-07:00

Digraphs可以用来在Vim中输入特殊字符。Digraphs使用两个一般字符来表示一个特殊字符，输入之前先按Ctrl+K，然后紧接着键入两个一般字符。例如要键入拼音字母ü，可以先按Ctrl+K，然后键入u:。如果要键入希腊字母θ，可以先按Ctrl+K，然后键入h*。这个双字符组合比较有规律，例如第二个字符

*表示希腊字母
：表示字母上面加两点
R表示罗马数字

等等。默认的组合列表可以在Vim里用:dig命令查看。另外，还可以手动添加字符。详细说明可以参考Vim Digraphs文档。

改变Linux终端分辨率

2008-09-22T14:35:00.004-07:00

在grub.conf文件里kernel启动选项加入参数

vga=xxx

其中xxx是代表显示模式的数字，具体多少可以从这里查出来。下面是一些常用的数值：

	640x480	800x600	1024x768
256 color palette	769	771	773
16-bit	785	788	791
24-bit	786	789	792

在Linux下合法使用微软字体

2008-09-20T22:17:00.002-07:00

corefonts提供了绕开Windows授权使用微软字体的方法，可以很方便的得到rpm安装包。这里是我在Fedora 9上编译好的rpm文件，以方便以后使用。

在Xfig或者Gnuplot生成的图画中加入公式或label

2008-08-20T16:53:00.003-07:00

Xfig和Gnuplot都可以把图画输出为MetaPost文件。用MetaPost的label命令就可以方便的插入公式或者标签。修改后的MetaPost文件可以借助mpost和mps2eps生成eps文件，具体方法可以参考这里。

新的Linux Flash Player

2008-04-23T13:44:00.002-07:00

仅仅安装Adobe的flash-plugin不能保证新的Linux Flash Player工作，还需要安装libflashsupport。最简单的办法是到Adobe的下载网站使用YUM for Linux。

新安装Fedora 9 Preview

2008-04-20T14:42:00.002-07:00

KDE4一开始不太适应，过了一段时间觉得还不错。不知道什么原因，每次log out之后，键盘就不管用了，用鼠标restart X就没有问题了，大概是kdm的bug。可以设置kdmrc让kdm每次log out后重启X。在kdmrc中

[X-*-Core]

后面增加一行：

TerminateServer=true

矢量绘图语言Asymptote（一）

2007-12-07T16:02:00.000-08:00

Asymptote是从MetaPost扩展而来的矢量数学绘图语言。它使用浮点数运算，支持三维图形，使用类似C++语法。Asymptote的很多语法和MetaPost类似，但是也有不同。下面简单介绍它的使用方法。假设Asymptote命令存在文件test.asy中，则执行下面的命令可以生成eps文件：

asy -V test

其中-V参数会自动打开一个PostScript浏览窗口，这样可以立即浏览生成的图形。

简单直线图形

Asymptote代码	生成的图像
draw((0,0)--(100,0)--(100,100) --(0,100)--cycle);
size(0,3cm); draw(unitsquare);
unitsize(2cm,3cm); draw(unitsquare); label("$A$",(0,0),SW); label("$B$",(1,0),SE); label("$C$",(1,1),NE); label("$D$",(0,1),NW);

Asymptote代码

生成的图像

draw((0,0)--(100,0)--(100,100)
--(0,100)--cycle);

size(0,3cm);
draw(unitsquare);

unitsize(2cm,3cm);
draw(unitsquare);
label("$A$",(0,0),SW);
label("$B$",(1,0),SE);
label("$C$",(1,1),NE);
label("$D$",(0,1),NW);

曲线图形

Asymptote代码	生成的图像
draw((0,0){up}..(100,25){right}..(200,0){down});
draw((0,0){up}::(100,25){right}::(200,0){down});
size(0,100); path unitcircle=E..N..W..S..cycle; path g=scale(2)*unitcircle; filldraw(unitcircle^^reverse(g),yellow,black);
size(200); real mexican(real x) {return (1-8x^2)exp(-(4x^2));} int n=30; real a=1.5; real width=2a/n; guide hat; path solved; for (int i=0; i<n; ++i) { real t=-a+iwidth; pair z=(t,mexican(t)); hat=hat..z; solved=solved..z; } draw(hat); dot(hat,red); draw(solved,dashed);

graph module

Asymptote代码	生成的图像
import graph; size(150,0); real f(real x) {return exp(x);} pair F(real x) {return (x,f(x));} xaxis("$x$"); yaxis("$y$",0); draw(graph(f,-4,2,operator ..),red); labely(1,E); label("$e^x$",F(1),SE);

计算折射率和其他光学变换的网站

2007-11-13T14:37:00.000-08:00

计算折射率和其他光学变换的网站

http://www.luxpop.com/

在Fedora 7下安装LaTeX CJK

2007-09-03T14:48:00.000-07:00

LaTeX CJK的安装一向都很麻烦，主要是字体安装和配置不容易。这里主要根据Edward Lee的Blog（http://blog.bs2.to/post/EdwardLee/7672）里的方法进行配置。

下载安装CJK包

http://cjk.ffii.org/

$HOME/texmf/tex/latex/CJK/

生成Type1字体和配置文件

FontForge

subfonts.pe

$HOME/song/

uming.ttf

utils

subfonts.pe

http://delloye.free.fr/Unicode.sfd

Unicode.sfd

fontforge -script subfonts.pe uming.ttf song Unicode.sfd

for i in *.tfm
do
cat >> song.map << EOF
${i%.tfm} ${i%.tfm} < ${i%.tfm}.pfb
EOF
done

c70song.fd

$HOME/texmf/tex/latex/CJK/texinput/UTF8/

% This is the file c70song.fd of the CJK package
%   for using Asian logographs (Chinese/Japanese/Korean) with LaTeX2e
%
% created by Werner Lemberg 
%
% Version 4.7.0 (17-Oct-2006)

\def\fileversion{4.7.0}
\def\filedate{2006/10/17}
\ProvidesFile{c70song.fd}[\filedate\space\fileversion]


% character set: Unicode U+0080 - U+FFFD
% font encoding: Unicode

\DeclareFontFamily{C70}{song}{\hyphenchar \font\m@ne}

\DeclareFontShape{C70}{song}{m}{n}{<-> CJK * cyberb}{}
\DeclareFontShape{C70}{song}{bx}{n}{<-> CJKb * cyberb}{\CJKbold}

\endinput

cyberb

song

$HOME/texmf/fonts/map/dvips/CJK/song.map
$HOME/texmf/fonts/tfm/CJK/song/*.tfm
$HOME/texmf/fonts/type1/CJK/song/*.pfb

生成ls-R数据库

texhash

更新字体map

updmap --enable Map song.map

测试

\documentclass{article}
\usepackage{CJKutf8}
\begin{document}
\begin{CJK}{UTF8}{song}
\LaTeX CJK 测试。
\end{CJK}
\end{document}

如何使用Maxima（三）

2007-03-30T16:07:00.001-07:00

复数运算

(%i1) a:2-3*%i$
(%i2) b:5+4*%i$
(%i3) a+b;
(%o3)                               %i + 7
(%i4) a*b;
(%o4)                        (2 - 3 %i) (4 %i + 5)

%i1

%i2

$

;

(%i5) expand(%);
(%o5)                              22 - 7 %i
(%i6) realpart(%);
(%o6)                                 22
(%i7) imagpart(%o5);
(%o7)                                 - 7

(%i8) abs(%o5);
(%o8)                              sqrt(533)
(%i9) carg(%o5);
                                        7
(%o9)                             - atan(--)
                                        22

(%i10) polarform(%o5);
                                     - %i atan(7/22)
(%o10)                    sqrt(533) %e
(%i11) rectform(%);
                               7                            7
(%o11)       sqrt(533) cos(atan(--)) - sqrt(533) %i sin(atan(--))
                               22                           22
(%i12) %,numer;
(%o12)                           22.0 - 7.0 %i

方程数值解

f(x)=0

[a, b]

find_root(f(x), x, a, b);

在线收听古典音乐频道

2006-11-28T13:50:00.000-08:00

WCPE古典音乐频道，位于Rolesville，NC，24小时播音。这里是在线收听地址： http://audio-ogg.ibiblio.org:8000/wcpe.ogg.m3u

Linux下复制Video DVD

2006-10-26T15:51:00.000-07:00

首先制作iso镜像文件：

mkisofs -dvd-video -o dvd.iso /media/dvd

刻录DVD：

growisofs -dvd-compat -Z /dev/hdc=dvd.iso

Gnuplot的Palette Map作图

2006-09-11T19:31:00.000-07:00

例子（TEM22高斯型光束强度分布存在gaussian_profile文件中）：

gnuplot> set view map
gnuplot> unset surface
gnuplot> set style data pm3d
gnuplot> set style function pm3d
gnuplot> set pm3d at b
gnuplot> set palette gray
gnuplot> set size square
gnuplot> splot "gaussian_profile"

生成的图像：

从YouTube，Google Video等站点下载视频

2006-08-17T15:39:00.000-07:00

YouTube，Google Video等网站利用新的FLV（Flash Video）流媒体视频技术允许用户在互联网上传播录像，但有时候我们希望把录像保存到自己电脑上，而不是仅仅在网上浏览。Firefox的扩展工具VideoDownloader能够让用户把这些录像下载到本地。下载之后的文件是FLV格式，如果想在一般媒体播放器里看，可以在Linux下用ffmpeg把FLV文件转化成MPG或其他格式：

 ffmpeg -i video.flv -ab 128 -ar 44100 -b 200 -s 320x240 video.mpg

上面命令中，-b表示视频比特率，-ab表示音频比特率，-ar表示音频取样率，-s表示视频画面大小（宽x高）。

学用绘图工具MetaPost（二）

2006-08-13T20:53:00.000-07:00

线性方程

z1

z2

z3

z3=1/3[z1,z2]

z1

z2

whatever[z1,z2]

z1

(x1,y1)

MetaPost代码	生成的图像
beginfig(7); z1=(-1.5cm,0); z2=(-3.2cm,-2cm); z3=(2cm,0); z4=(3cm,-2cm); z5=whatever[z1,z2]=whatever[z3,z4]; draw z2--z1--z5; draw z4--z3--z5; pickup pencircle scaled 4pt; draw z1; draw z2; draw z3; draw z4; draw z5; endfig;

MetaPost代码

生成的图像

beginfig(7);
z1=(-1.5cm,0);
z2=(-3.2cm,-2cm);
z3=(2cm,0);
z4=(3cm,-2cm);
z5=whatever[z1,z2]=whatever[z3,z4];
draw z2--z1--z5;
draw z4--z3--z5;
pickup pencircle scaled 4pt;
draw z1; draw z2; draw z3;
draw z4; draw z5;
endfig;

图像里插入文字和符号

label

lft|rt|top|bot|ulft|urt|llft|lrt

btex...etex

beginfig(8);
a=.8in; b=.6in;
z0=(0,0);
z1=-z3=(a,0);
z2=-z4=(0,b);
draw z1..z2..z3..z4..cycle;
draw z1--z0--z2;
label.top("a", .5[z0,z1]);
label.lft("b", .5[z0,z2]);
dotlabel.bot("(0,0)", z0);
endfig;

beginfig(9);
numeric u;
u=1cm;
draw (0,2u)--(0,0)--(4u,0);
pickup pencircle scaled 1pt;
draw (0,0){up}
   for i=1 upto 8: ..(i/2,sqrt(i/2))*u endfor;
label.lrt(btex $\sqrt x$ etex, (3,sqrt 3)*u);
label.bot(btex $x$ etex, (2u,0));
label.lft(btex $y$ etex, (0,u));
endfig;

用gnuplot和xfig输出MetaPost文件

set terminal mp
set output "filename.mp"

学用绘图工具MetaPost（一）

2006-08-11T21:22:00.000-07:00

MetaPost是贝尔实验室的John Hobby基于Knuth的METAFONT设计的一种绘图语言。它适于绘制精确的数学图形，可以方便的生成EPS文件，供插入到LaTeX文档中。

Linux下可以用mpost命令处理编辑好的MetaPost文件，生成EPS文件。生成的EPS文件可以作为图片插入到LaTeX文件中，但是没办法直接用gv等PostScript浏览器直接观看，这里有关于这个问题的讨论和解决办法。另外，可以用mptopdf命令代替mpost，可以解决这个问题，而且可以直接生成PDF文件。

下面是MetaPost的简介。beginfig...endfig之间是绘图命令。实际使用中，好几个beginfig...endfig可以放在一个文件里，文件最后需要有一个end命令。

画直线和曲线

MetaPost代码	生成的图像
beginfig(1); z0=(0,0); z1=(3cm,2cm); z2=(2cm,4.5cm); z3=(0.5cm,3.5cm); z4=(1.5cm,2.5cm); draw z0--z1--z2--z3--z4; endfig;
beginfig(2); u=1cm; z0=(0,0); z1=(3u,2u); z2=(2u,4.5u); z3=(0.5u,3.5u); z4=(1.5u,2.5u); draw z0..z1..z2..z3..z4; endfig;
beginfig(3); z0=(0,0); z1=(3cm,2cm); z2=(2cm,4.5cm); z3=(0.5cm,3.5cm); z4=(1.5cm,2.5cm); draw z0..z1..z2..z3--z4--cycle; endfig;

MetaPost代码

生成的图像

beginfig(1);
z0=(0,0);
z1=(3cm,2cm);
z2=(2cm,4.5cm);
z3=(0.5cm,3.5cm);
z4=(1.5cm,2.5cm);
draw z0--z1--z2--z3--z4;
endfig;

beginfig(2);
u=1cm;
z0=(0,0);
z1=(3u,2u);
z2=(2u,4.5u);
z3=(0.5u,3.5u);
z4=(1.5u,2.5u);
draw z0..z1..z2..z3..z4;
endfig;

beginfig(3);
z0=(0,0);
z1=(3cm,2cm);
z2=(2cm,4.5cm);
z3=(0.5cm,3.5cm);
z4=(1.5cm,2.5cm);
draw z0..z1..z2..z3--z4--cycle;
endfig;

表示方向的dir和用于循环的for

dir a

dir 0

right

dir 90

up

left

down

beginfig(4); draw ( for i=0 upto 7: 80dir (135i)-- endfor cycle); endfig;
beginfig(5); for a=0 upto 7: draw (0,0){dir 45}..{dir 10a}(6cm,0); endfor endfig;

改变画笔粗细

beginfig(6);
u=1cm;
draw (-u,0)--(u,0);
draw (0,-u)--(0,u);
pickup pencircle scaled 3;
draw (-u,u)--(u,u)--(u,-u)--(-u,-u)--cycle;
endfig;

pdf文件处理工具：pdftk

2006-05-19T15:44:00.000-07:00

pdftk是一个处理pdf文件的小工具，可以进行合并pdf文档，提取其中几页成为新pdf文件，加密解密等操作。

合并多个pdf文件

pdftk 1.pdf 2.pdf 3.pdf cat output 123.pdf

提取文件某几页合成为新文件

pdftk A=one.pdf B=two.pdf cat A1-7 B1-5 A8 output combined.pdf

128位加密文件

pdftk mydoc.pdf output mydoc.128.pdf owner_pw foo user_pw baz

加密但是允许打印

pdftk mydoc.pdf output mydoc.128.pdf owner_pw foo user_pw baz allow printing

解密文件

pdftk secured.pdf input_pw foopass output unsecured.pdf

如何使用Maxima（二）

2006-05-05T16:10:00.001-07:00

三角函数变换

常用的三角函数变换有下面几种：

trigexpand: 利用和差化积公式展开
trigreduce: 利用积化和差等公式变成sin或cos的和
trigsimp: 利用sin²(x) + cos² = 1等公式简化
trigrat: 简化分数形式

下面是一些例子：

(%i1) sin(2*x)/cos(x)+cos(2*x);
                              sin(2 x)
(%o1)                         -------- + cos(2 x)
                               cos(x)
(%i2) trigexpand(%o1);
                             2                    2
(%o2)                   - sin (x) + 2 sin(x) + cos (x)
(%i3) trigreduce(%o1);
(%o3)                     sec(x) sin(2 x) + cos(2 x)
(%i4) trigsimp(%o1);
                          sin(2 x) + cos(x) cos(2 x)
(%o4)                     --------------------------
                                    cos(x)
(%i5) trigrat(%);
(%o5)                         cos(2 x) + 2 sin(x)

微积分

微分：

(%i1) f: sqrt(x-1)*exp(2*x^2);
                                               2
                                            2 x
(%o1)                         sqrt(x - 1) %e
(%i2) diff(f, x);
                                                     2
                                         2        2 x
                                      2 x       %e
(%o2)               4 sqrt(x - 1) x %e     + -------------
                                             2 sqrt(x - 1)
(%i3) diff(sin(x)*exp(x^2), x, 2);
                       2             2                 2
                  2   x             x                 x
(%o3)          4 x  %e   sin(x) + %e   sin(x) + 4 x %e   cos(x)

不定积分：

(%i4) integrate(%o2, x);
                                       2   log(x - 1)
                                    2 x  - ----------
                                               2
(%o4)                     (x - 1) %e

定积分：从这个例子可以看出，maxima是一个很聪明的代数系统，它不会像其他数学程序那样想当然的默认一个结果，而是主动问你必要的条件，然后给出正确的结果。

(%i5) integrate(exp(-a*x^2), x, 0, inf);
Is  a  positive, negative, or zero?

positive;
                                   sqrt(%pi)
(%o5)                              ---------
                                   2 sqrt(a)

泰勒展开：这个例子给出函数g的前三阶泰勒展开。

(%i6) g: sin(2*x)*exp(a*x);
                                  a x
(%o6)                           %e    sin(2 x)
(%i7) taylor (g, x, 0, 4);
                                2       3     3         4
                        2   (3 a  - 4) x    (a  - 4 a) x
(%o7)/T/     2 x + 2 a x  + ------------- + ------------- + . . .
                                  3               3

刘再复这样概括巴金

2006-05-04T21:52:00.000-07:00

‘文革’结束了，伸冤的伸冤，平反的平反，乌纱帽掉了又带上了，而那些在牛棚里受难的作家与革命家们个个都在讲述自己的苦难过去和向社会索取报酬，甚至讴歌‘二次解放’。只有一个受折磨的灵魂，在寓所的墙角下对过去自己的著作与行为感到深深不安，感到揪心的惭愧，灵魂和双手都在战栗，这就是《真话集》(即《随想录》)的伟大作者……惟一,惟一，Only One！中国的大地这么辽阔，大地上所发生的历史事件这么惨烈，卷入历史运动的人群数亿之众，但我们在大喧哗之后听到的良心呻吟与自我鞭挞却只有一个！

如何使用Maxima（一）

2006-05-03T20:22:00.001-07:00

Maxima是一个优秀的开放源代码代数运算系统，符号运算能力可以和著名商业数学软件Mathematica相媲美。这里是它的使用方法的一个简介。

运行和退出
运行：
```
maxima
```
进入maxima环境后，会看到开始处有(%i1)，这是行标。%i表示输入行，%o表示输出行，后面的数字表示第几个输入或输出。退出：
```
quit();
```
可以看出，maxima语法和C语言类似，语句结束需要加分号。maxima里面是要区分大小写的。
算术运算
符号+, -, *, /, ^分别表示加、减、乘、除、乘方。maxima是一个精确计算的代数系统，例如分数，开方等运算会保持符号形式，不会被自动转化为小数。
```
(%i1) 2/5+3/4;
                                      23
(%o1)                                 --
                                      20
```
如果要输出小数，可以用numer命令：
```
(%i2) %,numer;
(%o2)                                1.15
```
百分号%代表最后一次输出的结果。另外，可以用bfloat输出任意精确度的小数，有效数字位数由变量fpprec控制，默认是16位。
```
(%i3) fpprec:100;
(%o3)                                 100
(%i4) bfloat(%pi);
(%o4) 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117068b0
```
(%i3)中的冒号表示赋值。

代数运算

多项式展开：

(%i5) (x^2+2*x*y+3*y^2)^3;
                                 2            2 3
(%o5)                        (3 y  + 2 x y + x )
(%i6) expand(%);
            6         5       2  4       3  3       4  2      5      6
(%o6)   27 y  + 54 x y  + 63 x  y  + 44 x  y  + 21 x  y  + 6 x  y + x

等量代换：

(%i7) %o6, x=2/z;
                 5        4        3        2
            108 y    252 y    352 y    336 y    192 y   64       6
(%o7)       ------ + ------ + ------ + ------ + ----- + -- + 27 y
              z         2        3        4       5      6
                       z        z        z       z      z

提取公分母：

(%i8) ratsimp(%);
          6  6        5  5        4  4        3  3        2  2
      27 y  z  + 108 y  z  + 252 y  z  + 352 y  z  + 336 y  z  + 192 y z + 64
(%o8) -----------------------------------------------------------------------
                                         6
                                        z

因式分解：

(%i9) factor(%);
                                2  2             3
                            (3 y  z  + 4 y z + 4)
(%o9)                       ----------------------
                                       6
                                      z

解非线性方程组：

(%i10) a+b+2*c=0;
(%o10)                          2 c + b + a = 0
(%i11) 2*a-b=1;
(%o11)                            2 a - b = 1
(%i12) a-2*b*c=5;
(%o12)                           a - 2 b c = 5
(%i13) solve([%o10,%o11,%o12],[a,b,c]);
             sqrt(7) + 1      2 sqrt(7) - 1      sqrt(7) - 14
(%o13) [[a = -----------, b = -------------, c = -------------],
                  3                 3            4 sqrt(7) - 2
                          sqrt(7) - 1        2 sqrt(7) + 1      sqrt(7) + 14
                   [a = - -----------, b = - -------------, c = -------------]]
                               3                   3            4 sqrt(7) + 2

反驳对Linux的两种无知的指责

2006-04-24T16:55:00.000-07:00

经常听到对Linux有若干指责，其实完全是无知的表现。第一，说“用Linux有法律隐患”。其实，目前针对Linux的法律起诉非常少，而且据我所知，Linux还没有一起败诉的。相反，Microsoft招致的起诉远远多于Linux，而且很多以败诉告终。以今天的软件复杂度，有知识产权纠纷是非常正常的，可以说任何一个软件企业都有陷入法律纠纷的可能，单单指责Linux是没有道理的，更何况Linux根本不是一个企业。Windows和MS Office都有败诉的例子，而人们还在使用他们，这些人凭什么去指责Linux有法律问题。第二，说“Linux版本众多，相互不兼容”。Linux下绝大部分软件是以源代码形式发布的，安装时根据不同的平台进行编译，兼容众多的硬件平台，绝对比只运行在PC上Windows兼容性强。为了大家方便，也有很多编译后的二进制代码以软件包的形式发布。虽然Linux有很多版本，但是软件包的管理方式主要有两种：一种是Red Hat的rpm，一种是Debian的dpkg。通过alien软件，这两种包可以互相转换，所以不存在不兼容的问题。一些硬件厂商以Linux版本多，支持起来不容易为由拒绝为Linux提供驱动，就更是可笑了。他们只要提供内核级的驱动支持就可以了，Linux发行版再多跟他们也没关系。有些人抱怨Linux支持的硬件少，这是因果颠倒，他们应该去抱怨这些硬件厂商不支持Linux。

莫扎特作品下载

2006-04-21T18:17:00.000-07:00

这里有好多莫扎特作品mp3下载。虽然音质不是特别好，但是作为莫扎特作品入门完全可以了。

Monte Carlo方法

2006-04-18T15:45:00.001-07:00

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

谈谈Super64吧

2006-04-14T13:23:00.001-07:00

生日礼物耶！我不太舍得买，但是很仰慕这把琴（或者说比较仰慕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

2006-04-12T13:30:00.000-07:00

这里假设是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图形中插入文字

2006-04-11T18:30:00.000-07:00

% 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文件

2006-04-10T17:48:00.000-07:00

有时候需要把多个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。各个文件中间部分是描述图形的，不需要管。只要改动一下每个文件的开头和结尾部分就可以了。

每个eps文件开头都有类似于这样一行：
```
%!PS-Adobe-2.0 EPSF-2.0
```
只要保留第一个文件中的就行了，其余文件开头的这一行都删去。
每个eps文件结尾都有这样的命令：
```
showpage
```
只需要保留最后一个就可以了，其余文件结尾处的该命令都删去。
每个文件开头都有类似于这样一行：
```
%%BoundingBox: 0 0 100 120
```
这表示图形的大小，四个数字分别是左下角和右上角的坐标。我们需要根据自己的要求把数字改一下。例如如果每个eps大小都是100x120，我们希望得到2x2的图形列表，那么合成之后的图应该是200x240的。那么我们把第一个eps文件开头的BoundingBox改为
```
%%BoundingBox: 0 0 200 240
```
，然后把其余各eps文件开头的BoundingBox行删掉。
在第一个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

2006-04-06T21:04:00.000-07:00

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

设置Linux router

2006-04-01T17:07:00.000-08:00

如果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图

2006-03-30T20:33:00.000-08:00

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

生日快乐！

2006-03-24T16:42:00.001-08:00

祝自己生日快乐！

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

Fedora Core 5安装flash插件

2006-03-23T14:45:00.000-08:00

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

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

2006-03-22T13:04:00.000-08:00

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

用gnuplot拟合Gaussian曲线

2006-02-22T14:03:00.000-08:00

在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下载流媒体文件

2006-01-30T16:18:00.000-08:00

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

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

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

古典音乐作品编号标记

2006-01-24T14:03:00.000-08:00

很多时候我们见到一些古典音乐作品的名字后面跟着由字母和数字组成的编号,例如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文件多页合并到一页的小工具

2005-10-23T19:01:00.000-07:00

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

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

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

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

2005-10-20T14:39:00.001-07:00

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

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

2005-10-18T20:25:00.000-07:00

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

在Linux下把wma文件转化为mp3

2005-09-21T17:43:00.000-07:00

需要用到mplayer和lame

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

transcode编码DVD的bug

2005-09-09T19:25:00.000-07:00

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里作拼写检查

2005-08-18T19:14:00.000-07:00

:w
:!aspell -c %

一个介绍激光基础知识和各种激光器的网页

2005-06-29T13:53:00.000-07:00

Sam's Laser FAQ

Fedora使用kdm显示管理器

2005-06-24T01:05:00.000-07:00

Fedora默认显示管理器是gdm，对于习惯使用KDE桌面环境的用户，无法在退出桌面时直接选择关机或重启。在/etc/sysconfig/desktop文件里做如下更改可以使用kdm显示管理器：

DISPLAYMANAGER="KDE"

注意：

DESKTOP="GNOME"

是默认桌面环境，和显示管理器是两个概念。也就是说，完全可以在kdm下使用GNOME或KDE任何一种桌面。

Linux下刻录CD失败的解决方案

2005-06-19T17:48:00.000-07:00

最近安装Fedora Core 4,刻录的CD总是无法通过安装时的mediacheck检测。Google发现这是一个很普遍的问题，刻录之后的CD经常和原始的iso文件 md5sum不一样。凭我个人的经验，K3B刻出的所有CD都有这个问题，X-CD-Roast时好时坏，其中原因我目前还不清除。其实这样的CD并不意味着文件有问题不能使用，但是不能通过md5sum总让人无法放心。有人提出使用-pad参数可以避免这样的问题，可是我使用X-CD-Roast刻盘时使用了这个参数的，一样会出问题。这里有一个方案，提出使用padsize=63s。例如：

cdrecord dev=/dev/cdrom speed=16 padsize=63s -pad -v -eject /home/myuid/myiso.iso

尝试了一下这个方法，刻的CD成功通过了md5sum测试。其效果有待进一步测试。

Adobe Reader 7.0在Linux下打印无法选择打印机的问题

2005-06-18T18:51:00.000-07:00

虽然Linux下有不少很好的PDF文件浏览器，例如xpdf，但是某些功能例如form还是只有Adobe Reader才能实现。Adobe Reader 7.0的Linux版本功能很好，但是如果有多个打印机，Adobe Reader不允许用户选择打印机，只能用默认的打印机。LinuxFocus提供了一个用kprinter的解决方案：

#!/bin/sh
# save this as /usr/local/bin/printerScript
# and make it executable with:
# chmod 755 /usr/local/bin/printerScript
## create temporary file
tmp="/tmp/postscript-print-$$.ps"
[ -f "$tmp" ] && exit 1
## save the file to a new name which will not be deleted
## by acroread. the first argument is the file name of
## the postscript file
cp  "$1"  $tmp
## start the KDE printer frontend on this file
kprinter  $tmp
## clean up afterwards
rm -f     $tmp
##

将上面的文件存为/usr/local/bin/printerScript，然后把Adobe Reader的printer Command一栏由/usr/bin/lp改为/usr/local/bin/printerScript，这样在打印时会自动启动kprinter，允许用户选择使用的打印机。

Fedora Firestarter

2005-06-14T20:04:00.000-07:00

Fedora Core 4发布。Fedora Extra被默认加入yum镜像列表。Firestarter是Fedora Extra的一个防火墙配置工具，比Fedora Core默认的system-config-securitylevel略为灵活一些，然而，生成的iptables规则有些过于复杂。另外，卸载之后，使用DHCP动态IP的系统重启时会提示firestarter文件找不到，这是因为/etc/dhclient-exit-hooks文件有关firestarter的内容并没有随着卸载而删除掉。

Guarddog

2005-06-10T23:37:00.000-07:00

Guarddog是一个KDE环境下的iptables防火墙配置工具。它的图形界面非常友好，使用者不需要太多网络方面的知识。 Guarddog默认分为两个网络区域，分别是Internet（外面的世界）和Local（你的计算机）。软件默认首先阻挡所有网络连接，然后由用户指定哪个端口需要开启，这样最大可能的保证了系统安全性。指定端口的方式很简单。如果用户的计算机不运行任何服务（一般用户的电脑都是这样的），那么用户的计算机（Local区域）的所有端口都是关闭状态。如果用户要访问互联网网页，只要把Internet区域中的http开启。如果要访问ftp，就把Internet中的ftp开启。常用的端口包括http, https, ftp, telnet, ssh, POP3, SMTP, DNS等等，另外还有MSN, Yahoo等即时通讯软件和Quake等游戏服务器端口，基本上满足了绝大多数用户的要求。 Guarddog目前功能上的局限也很明显，例如无法配置NAT，处理多个网络界面时缺乏区别对待的灵活性，iptables的很多参数还无法实现等等。但是对于一般的非技术用户，应该足够了。

Debian 3.1 发布

2005-06-08T02:26:00.000-07:00

三年多了，Debian终于发布了新版本。安装程序比以前有了大改观，但是易用性比Fedora等仍然差一些，目前还是只能定位于高级用户。

发现一个光衍射干涉分析软件DIFFRACT

2005-05-31T14:59:00.000-07:00

MM Research Inc. website: http://www.mmresearch.com/

决定自己学修车了

2005-05-25T16:45:00.000-07:00

又要花几百美元修车了。以前觉得费事，不如多花些钱让车铺去修，但是这次觉得实在不值，很多问题可以自己解决的，并没有那么麻烦。也好，让我终于下决心要尝试自己修车了。

今天的感想

2005-05-20T23:55:00.000-07:00

以后一定要google搜索之后再开始自己的计算推导，可能会省很多精力。

Fedora Core 4 Test　3　发布

2005-05-11T18:26:00.000-07:00

Fedora Core 4 Test　3　昨天发布。今天安装试用。界面改进很明显，或许是托新版KDE和Gnome的福。难得的是除了bluecurve之外其他的主题改进也很明显。现在FC已经成了一个可用性非常强的系统。当然bug还是有的。例如Gimp，还有中文支持，但是应该比较容易改正，总体讲比FC2的测试版的bug少了许多。OpenOfffice.org 2.0　Beta版本增加了数据库软件Base，因此和MS Office的竞争力又增强了一些。另外，主菜单里多了一些新面孔，例如新的电影播放器Totem，新的pdf和postscript文档浏览器 Evince。其他软件正在试用中。

标题解读

2005-05-05T16:56:00.000-07:00

瑞典皇家科学院将2004年诺贝尔物理学奖授予美国科学家戴维·格罗斯(David Gross)、戴维·波利策(David Politzer)和弗兰克·维尔切克(Frank Wilczek)，以表彰他们发现了粒子物理的强相互作用理论中的“渐近自由”(asymptotic freedom)现象。 3位科学家在上个世纪70年代提出的理论认为，强作用力会随着夸克彼此间距离的增加而增大，因此没有夸克可以从原子核中向外迁移，获得真正的自由。通俗地说，这一现象有点像拉一根具有弹性的橡皮筋：橡皮筋拉得越长，其产生的力量越大，人拉起来也更为费劲。同样根据“渐近自由”理论，强作用力会随着夸克间距离的变小而减弱，这意味着，约束在质子等内部的夸克在彼此距离足够小时将近乎自由地进行运动。三位科学家于1973年完成的工作所采用的数学形式是后来描述强相互作用理论——量子色动力学(QCD)的基础。

`x < y`	小于
`x <= y`	小于等于
`x == y`	等于
`x >= y`	大于等于
`x > y`	大于
`x != y` 或 `x ~= y`	不等于

det (a)	行列式
eig (a)	本征值
dot (x, y)	点乘
inv (a)	逆矩阵

Asymptotic Freedom

gnuplot所有的colornames及范例

数码单反（DSLR）相机raw格式图片的读取和处理

《Maxima快速参考手册1.1》

Vim笔记（二）

Vim笔记（一）

几个关于Gnuplot小技巧的网页

Octave笔记

统计笔记（一）

写了一份Maxima的中文文档

在eps图片上添加LaTeX公式

GnuPlot在PostScript终端下使用任意字体

sed编辑器的使用

GnuPlot使用png terminal时找不到字体的问题

在USB drive（U盘）上安装Linux

Mie散射计算器

GNU screen的几个常用命令

Linux下用sopcast看电视

光圈档数（aperture stop）的计算方法

LaTeX CJK生成PDF文件的中文搜索以及复制粘贴

用WebPlayer在网页里播放mp3

光束质量因子M2（M squared）的测量

让blogger支持MathML公式输出

不用重启切换runlevel

魔术组合键SysRq

如何使用Maxima（四）

矩阵运算

Vim输入特殊字符（例如希腊字母）

改变Linux终端分辨率

在Linux下合法使用微软字体

在Xfig或者Gnuplot生成的图画中加入公式或label

新的Linux Flash Player

新安装Fedora 9 Preview

矢量绘图语言Asymptote（一）

计算折射率和其他光学变换的网站

在Fedora 7下安装LaTeX CJK

如何使用Maxima（三）

在线收听古典音乐频道

Linux下复制Video DVD

Gnuplot的Palette Map作图

从YouTube，Google Video等站点下载视频

学用绘图工具MetaPost（二）

学用绘图工具MetaPost（一）

pdf文件处理工具：pdftk

合并多个pdf文件

提取文件某几页合成为新文件

128位加密文件

加密但是允许打印

解密文件

如何使用Maxima（二）

三角函数变换

微积分

刘再复这样概括巴金

如何使用Maxima（一）

运行和退出

算术运算

代数运算

反驳对Linux的两种无知的指责

莫扎特作品下载

Monte Carlo方法

谈谈Super64吧

体积：

吹嘴：

盖板：

按键：

簧片：

音色：

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

在eps图形中插入文字

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

谢旭东的Blog

设置Linux router

用blender做的第一个3D图

生日快乐！

祝自己生日快乐！

Fedora Core 5安装flash插件

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

用gnuplot拟合Gaussian曲线

用Mplayer下载流媒体文件

古典音乐作品编号标记

Fedora Core 4 Test　3　发布