liuxqsmile‘s Blog

ActPiano:The Android App for learning piano and sheet music by MIDI

2012-02-15T10:39:00.002+08:00

The Piano App turns your MIDI music to Sheet music and Play it on your Android devices! The app can generate a sheet music from MIDI file. You can playing the music by yourself by just
following the automatically scrolling keyboard and score. This App can also be a music learning game!

Top Features:
1. Generate sheet music from MIDI file;
2. Play the music by the app;
3. Learn to play the music by yourself. You can chose to play the
Treble clef or the bass clef;
4. Customizable keyboard and score view.

ActPiano on Android Market: https://market.android.com/details?id=actpiano.liuxqsmile.android

MARIAGE D'AMOUR

2011-03-28T23:09:00.001+08:00

Mariage D'Amour

20100312 三多村

2011-03-24T21:55:00.002+08:00

65Km。

2011.01.01阳宗海 94KM

2011-01-28T12:51:00.001+08:00

2010.12.19长虫山 44km

2010-12-23T09:21:00.001+08:00

如图

2010.11.20棋盘山

2010-11-22T09:56:00.001+08:00

44KM; 16KM/h; 2:45

Milestone三

2010-11-10T19:23:00.001+08:00

The sound of Silence

就这个水平，将就着能听。

Milestone二

2010-11-09T19:05:00.001+08:00

这个只是花钱换来的

Milestone一

2010-11-08T17:14:00.002+08:00

From liuxqsmile‘s Blog

2010.10.16，9:30~19:00, 145KM，19.5KM/h，Giant ATX750.

PDF免费工具

2010-11-03T21:24:00.001+08:00

PDF Xchange Viewer

免费的PDF阅读器，可以加入注释、图章等元素。

PDFCreator

开源PDF打印机。不支持保存为中文名。

Word to PDF: Conversion made easy with PDF-T-Maker

实际上是一个Word插件，需要配合PDFCreator等PDF生成工具，可以从Word的标题等生成PDF书签。Word 2010里另存为PDF时可以加入书签，就不需要这个插件了。

PDFill

其中包含了很多免费的工具，合并、水印、格式转换等。

SPKtool加入Sourceforge

2008-10-07T20:38:00.001+08:00

作为我这几年的重要收获之一，我希望能够将这个工具长期维护、发展下去，即使今后不在从事相关的工作。

因此，我将SPKtool放在了Sourceforge上，地址在这里： http://spktool.sourceforge.net/。

如果您也从事相关的工作，有这方面的实用程序，欢迎拿出来共享，共同创造一个完善的电生理信号分析工具。

本博客和Google Group上放的版本仍然保留，但不会再更新，建议到新地址下载。

最近把版本号升到了1.0.1，包含了比较完整的中、英文使用手册。改进了Perievent Raster绘图。

这个工具包还有很多不足，最大的问题是各个MATLAB函数的内置说明不完善。我还想增加多电极分析、更好的自动聚类算法和更多的统计绘图功能，甚至移植到开源的软件平台。但现在时间太紧张，希望一年之内可以再开始更新吧。

SourceForge建立项目

2008-09-24T18:12:00.002+08:00

前几天，在Sourceforge 上建了一个项目。平时网络的“花边”新闻倒是看了不少，一操作起来可就立马撞墙了。SF的帮助在它的Wiki 上，help里面反倒没什么有用的信息。简单总结一下，供外行们参考。

Sourceforge支持OpenID登录，但必须要申请一个SF的ID才可以建立项目。
项目的名称是以后可以改的。但unix name不能改，必须是唯一的。
项目必须经过站方的审批。所以关于项目的描述应该是英文的（这是我猜的，SF的审核人员应该不懂中文吧）。审核需要一个工作日，通过后需要为自己的项目建立分类信息等等，这些也都可以在以后重新修改。

审核通过后就可以发布源代码等文件了。在Admin中选择File Release，建立一个新的Package，新的Release。仔细看一下里面的文字描述，一般有不少链接，指向更详细的说明。
可以通过很多方式将本地文件上传到SF服务器，Web直接上传或者使用一些工具。只有 rsync over ssh是支持断点续传的，推荐国内用户使用。Windows用户可以安装Cygwin ，它相当于Windows上的Linux虚拟机。Cygwin是在线安装的，选择一个镜像地址，在安装包里选择以rsync和openSSH开头的几个包。完成后运行，会在安装目录里生成用户文件夹，例如当前Windows用户是USER，安装目录中会有homeUSER文件夹。将需要上传的文件拷贝到这个目录，然后按照SourceForge.net Wiki上的描述，在Cygwin中执行相应的命令就可以了。
通过Rsync over SSH上传文件：

[jsmith@linux ~]# rsync -avP -e ssh FILE AcountName@frs.sourceforge.net:uploads/
AcountName@frs.sourceforge.net's password:
building file list ...
1 file to consider
FILE
15000000 100% 34.13kB/s 0:07:08 (xfer#1, to-check=0/1)

sent 15001925 bytes received 42 bytes 29560.53 bytes/sec
total size is 15000000 speedup is 1.00

等到最后一行出现，才表示上传完成，有时候需要点耐心。但总比网页上传一次次报错好多了。然后到文件发布页面，选择修改Release，里面应该已经有了上传的文件，附加上去就可以了。

还可以在SourceForge上为项目建立网站，SourceForge提供了100MB的空间。对静态网页，在本地将网站建好，上传上去就OK了。我只建了一个index.html，只需要上传到服务器上的htdocs目录下就可以了。

用Cygwin中的sftp上传：

 sftp AcountName,ProjectName@web.sourceforge.net
cd htdocs
put index.html
bye

然后在浏览器中输入http://ProjectName.sourceforge.net就可以访问项目主页了。

SPKtool使用手册

2008-09-04T19:18:00.003+08:00

SPKtool
A tool for neurophysiological data processing
锋电位分析工具

版本1.0 2008-8

1 简介及安装

2 数据文件

3 锋电位检测

5 聚类分析

6 绘图功能

7 帮助及反馈

1 简介及安装

1.1简介

SPKtool是一个在MATLAB平台上编写的神经元锋电位分析工具软件，提供了从锋电位检测、特征提取、聚类分析到放电序列分析的各种功能：

锋电位检测：

显示记录信号中锋电位峰值的分布曲线；
通过阈值检测锋电位；
通过原始信号、能量、非线性能量进行检测；
三种锋电位波形对齐方式：阈值、波谷、波峰。
GUI交互式界面。

特征提取：

多种常用特征：PCA，峰值、能量等；
小波方法；
特征分别图，特征密度分布图等。

聚类分析：

K-means；
模板匹配；
基于高斯混合模型的EM算法；
Valley seeking；
通过绘制轮廓线，手动聚类；
对聚类结果进行合并，删除等。

动作电位串分析：

相邻动作电位间隔分布图(ISI)；
Poincare map；
相关函数图；
锋电位分布光栅图；
放电频次直方图；
刺激前后放电直方图（PSTH）；
刺激前后放电光栅图。

1.2安装及界面

本程序在MATLAB 7.1 SP3上开发完成，可能需要该版本或更新的版本的MATLAB才能运行。

将程序压缩包拷贝到MATLAB中的work目录下，解压缩，产生SPKtool目录。在MATLAB中，选择菜单File->Set Path … ，在“Set Path”对话框中选择“Add with Subfolders…”，把该目录及子目录添加到MATLAB的路径中去，选择“Save”按钮保存。

在命令窗口输入：spktool，打开GUI界面如下图。

程序界面中包含三个部分，菜单、信息表格和底部的状态栏。软件的所有功能均通过菜单命令来操作。信息表格列出了分析数据的主要信息，如采样频率，记录长度，文件名称、锋电位数目、特征名称、各个分类包含的锋电位数目等。对于多通道数据，通过信息表格右上角的Channel下拉菜单选择当前通道。状态栏给出程序运行状态的简明信息。

2 数据文件

SPKtool没有和其它软件进行数据交换的接口，不能直接读入记录系统产生的数据。用户需要首先将数据（连续记录信号、锋电位波形矩阵等）导入MATLAB，并按照特定的变量名保存到.mat文件（MATLAB中默认的数据文件格式），才能在SPKtool中打开，进行分析。

在菜单中选择File->Load…，打开文件选择对话框。加载数据后，根据文件中的变量，自动设置各个菜单的状态。

SPKtool通过变量名区分数据，因此，对数据文件(.mat)中的变量名有严格的要求。

2.1原始记录数据

原始记录数据文件中至少应包含两个变量：

dataflow：连续记录数据，二维数组，行对应通道，列对应采样点；

fs：采样频率，数值。

默认情况下，SPKtool对读入的数据进行低通滤波，使用6阶Butterworth滤波器，截至频率4500Hz。可以打开MLoadRaw.m，找到butter函数对应的位置，修改滤波参数。

2.2锋电位波形数据

从原始信号中检测出的锋电位波形数据，至少应该包含以下三个变量：

posind：锋电位波形在连续数据中的采样位置标记。行向量。第n个波形发生的时刻为posind(n)/fs。

waveforms：锋电位波形数据，二维数组。每行对应一个锋电位的采样数据。例如，1000个波形，每个波形含40个数据点，则waveforms尺度为1000×40；

fs：采样频率。

数据文件中还可以包含如下一些变量：

featurespace：锋电位特征空间，二维数组，每行包含一个波形的特征。如1000个波形的PC1和PC2特征，则featurespace为1000×2。特征的数目不限。

featnames：特征名称，包含字符串的单元数组，每个字符串对应一个特征的名称。例如: 两个特征的名称分别为PC1和PC2，则 featnames 为{‘PC1’,’PC2’};

unitcell：包含分类信息的单元数组。每个单元对应一个类别中的锋电位序号列表。例如，unitcell{1}对应没有归类的锋电位编号。 unitcell{3} 是第二个类别中的锋电位编号数组。

对于多通道数据，上面这些变量需要放在一个名为ResultSet的结构体数组内，每个结构体对应一个通道。

刺激信息

stimmark: 刺激的开始和结束时间。行向量。例如[20 30 60 70]，表示两次刺激，第一次刺激从第20s开始，到30s结束，第二次刺激60s开始，70s结束。

上述变量可以共存于一个文件中。在程序中，如果找不到dataflow这个变量，将关闭锋电位检测功能（不能激活Detect菜单），找不到waveforms，将关闭特征提取功能（不能激活Features菜单）。

另外，程序中用到了两个全局变量，采样频率fs和定义各个类别色彩的矩阵colorary，所以在程序运行过程中，请勿执行clear all等能够清除全局变量的命令。

2.3数据保存和导出

选择菜单File->Save …，保存分析结果。文件格式同样是.mat文件，变量名称与读入数据相同。（注意，连续记录数据不包含在该数据文件内。）可以用SPKtool重新打开保存的文件，也可以直接用读入MATLAB，将特定的数据导出。

选择菜单File->Save Current Channel…，仅保存当前通道的分析数据。

菜单File-> Experimental Info …，设置试验信息，包含记录位置、时间、坐标、数据文件名称等信息。

菜单File-> Export to TXT file …，导出数据到文本文件。这将产生两个文本格式的文件，后缀名为.inf的文件包含上面的试验信息，以及数据信息（通道数、单元数等）。后缀名为.txt的文件是一个数据表格，第一列对应锋电位时间（即变量posind），第二列为类别编号，后面的各列为波形数据。

3 锋电位检测

如果读入的数据中包含原始连续记录信号（即dataflow变量），菜单Detect可用，其中包含以下两个命令：

3.1峰值密度分布曲线

菜单Detect->Plot Peak Histogram，绘出原始记录数据的峰值密度分布曲线。幅度为对数坐标。

3.2锋电位检测界面

菜单Detect->Load Detector…，加载锋电位检测，界面如下：

参数设置：

Detect On：该下拉菜单中包含了三种波形检测方式，即分别在原始信号（Raw Signal），原始信号的能量(Energy)和非线性能量(Nonlinear Energy)^¹曲线上进行检测；

Threshold：阈值。输入阈值后左图中的红色水平线也相应调整。信号中与阈值线有交点的波形将根据后面的参数进行截取；

Spike length：每个锋电位波形的长度，单位为ms。如果采样频率为40KHz，1ms长度的锋电位波形包含40个数据点；

Prethreshold ratio：对齐点以前数据长度比例。对齐点前保留的数据长度与波形数据总长度的比值。例如包含40个点的锋电位波形，对齐方式为valley时，在波谷位置前保留10个采样点，则该参数需要设置为10/40=0.25；

Dead time：不应期长度，单位为ms。例如，如果采样频率为40KHz，该参数设置为2时，程序在检测到一个波形后，将跳过该波形后紧邻的80个点，然后才开始寻找下一个波形。

Alignment：波形对齐方式，阈值位置（gate）、波谷 (valley)、和波峰 (peak) 对齐三种形式。阈值和波谷的对比如下：

a) gate b) valley

使用方法：

中间的坐标轴根据Detect On的设置，显示原始信号或其（非线性）能量曲线。底部坐标轴在（非线性）能量检测方式下显示原始信号波形。左侧坐标轴显示中间坐标轴中数据的峰值分布曲线（横坐标为对数坐标）。

可以利用顶部的工具栏，对中间坐标轴进行缩放和平移，左侧坐标轴的纵坐标和底部坐标轴的横坐标也将进行相应的调整。

参数设置完毕后，选择Detect进行波形检测，检测完成后，将打开一个单独的绘图窗口显示锋电位波形。如果不符合要求，可以修改参数，重新检测。最后，选择Done完成，并关闭该界面。

需要说明的是，由于原始记录信号中的数据量通常很大，为了加快运行速度，坐标轴中只绘出前0.4e6个采样点。当然，锋电位波形检测还是针对所有的采样数据。

右侧的参数中除Detect On外，具有自动保存功能，再次运行Detector时将自动加载前一次设置的参数，即使你前一次选择Cancel退出。

在切换不同的检测方式时，阈值将根据中坐标轴中的数据自动调整为：均值+/-3倍标准差。

4 锋电位特征提取

通过阈值检测出锋电位波形后，或者读入的数据中保护锋电位数据，这时，Feature菜单将可用。在该菜单中，可以进行锋电位特征的选择、计算、修改、绘图等功能。通常需要定义特征后，才能打开Sort菜单进行聚类操作。

4.1特征定义

菜单 Features->Feature Define…，打开如下的特征选择对话框：

左侧All Features 列表中列出了所有可用的特征。中间有两个特征列表：Features Vs Features 中的特征用于绘制特征比较图（菜单 Features-> Features Vs Features）；Features For Cluster 才是用于聚类分析的特征。本说明中，除非特别指明，特征均指用于聚类的特征。在右侧可以修改一些特征相关参数。

可用特征：

peak：锋电位波峰处的幅值；

peak.pos：峰值对应的位置；

valley：锋电位波谷处的幅值；

valley.pos：波谷对应的位置；

peak-valley：峰峰值，波峰到波谷的振幅。

width：波形宽度，即波峰位置到波谷位置的距离：peak.pos-valley.pos；

energy：波形的能量，即采样点的平方和。

nonlinear.energy：所有点的非线性能量之和，幅度平方与频率平方的乘积，其计算方法为：NEO(t) = v(t)²- v(t-1)*v(t+1)，（v(t)为记录信号）；

Slice：锋电位波形中某个时刻的幅值，对应参数Slice Parameters。例如每个锋电位含32个采样点，Slice Parameters为0.5，则Slice特征为锋电位中第32×0.5=16个采样点的数值。Slice Parameters可以包含多个值，提取多个采样点，例如 0.5，0.75 则第16和第24个采样点的数值都作为特征。

pca：主成份分析方法，可以在右侧PCA Parameters中设置要提取的主成份序号；

ridge.pca：基于小波脊的PCA特征提取方法，首先对波形逐个进行小波变换（使用mexh小波），然后提取出小波脊，再对小波脊上的小波系数进行PCA，提取特征。需要三个参数：小波分析的频率范围Frequency limits，小波脊的数目Number of Ridges，以及需要提取的主成份序号PCA Parameters。

完成特征选取后，选择Done，开始提取用于聚类分析的特征，并给出特征分布散点图。

4.2特征修改

特征比例缩放

Features->Scale Features，打开如下图的对话框，在其中，可以调整特征的相对范围，可能有利于改善K-Means的聚类结果。

删除特征

Features->Remove Features，删除部分使用中的特征。一般特征的计算不需要花很长时间，重新选择并计算即可。在同时使用Ridge.pca和PCA作为特征时，由于共用PCA参数，可以通过该工具重新选择。

4.3特征分布散点图

Features->Plot Features 绘制特征分布散点图，如果定义了两个特征，散点图是二维的，如果有三个或更多的特征，则只使用前三个，以三维形式显示，并同时画出特征向三个平面上的投影分布，如下图。

4.4特征密度图

Features-> Features Density：特征分布密度图。如下图，该图只能绘出前两个特征的分布密度。

4.5特征对比图

Features-> Features Vs Features，绘制特征分布的对比图，注意该功能中使用的特征与聚类分析的特征是独立的。并且，每次绘图都需要重新计算特征。

如果已经划分类别，不同类别的特征将以不同颜色区别显示。

4.6主成份

Features-> Principal Components：给出特征提取中主成份曲线。例如特征定义中使用第一个和第二个主成份，则给出PC1和PC2的形状曲线。

5 聚类分析

聚类功能集中于菜单Sort中。在选择具体的聚类方法之前，应该先设置聚类的参数。

5.1参数设置

菜单Sort->Options…, 打开参数设置对话框：

其中共有四个参数：

Unit Number：划分类别的数目，K-Means，期望最大化(EM)方法，以及模板匹配方法使用该参数；

K-Means Radius：K-Means方法的聚类半径，值越小，聚类的半径越小。模板匹配方法也受该参数影响；

Parzen Mult：寻谷聚类算法的初始半径大小，初始半径越小，获得的类别数目越多。

Sort On：设置聚类的锋电位范围，选择All Spikes时对所有锋电位的特征进行聚类；选择Unsorted Spikes时仅对未分类的锋电位进行聚类。

5.2手动聚类

菜单Sort->Contour…，对应手动聚类，选择该菜单后，弹出如下图的手动聚类对话框。

对话框中间给出了未分类的特征的分布图。底部有四个控制按钮。首先单击New Contour按钮，光标将变成图中所示的十字形，在特征分布图中单击左键，将出现一个红色的十字形标记点，这样依次点击，就可以确定一个类别的大致轮廓，之后，鼠标右击，这些点将自动连接成一条封闭的曲线。然后选择底部的Update按钮，曲线内的点就被划分到一个类别。再选择New Contour绘制下一个类别的边界，依此操作，直到聚类完成。

也可以连续绘制多条边界线，然后再选择Update，显示聚类结果。最后选择Done将显示的聚类结果传回主程序，选择Cancel放弃。

手工聚类的类别数目不受前面聚类参数对话框中的类别数目的影响。

5.3自动聚类方法

SPKtool中包含了几种常见的聚类算法，简单说明如下：

K-Means

根据特征之间的距离来进行聚类。用户指定聚类的数目（聚类参数对话框中），该算法首先从样本（特征）中随机抽取确定聚类中心，然后计算每个样本到各个聚类中心的距离，并将其归到距离最近的类中，第一次聚类完成后，计算各个样本的均值作为新的聚类中心，重新计算每个样本到新的聚类中心的距离，并将其划分到聚类最近的类别中，如此循环，直到聚类中心的位置稳定^[^²^]。

Template 模板匹配法

在SPKtool中，模板方法的与K-Means基本上是一样的，不同之处仅在于直接用锋电位的波形数据作为特征，而K-Means使用锋电位降维后的特征数据。

Gaussian EM 基于高斯混合模型的期望最大化算法

该算法首先用K-Means方法确定初始聚类，然后将锋电位特征的分布看成符合混合高斯分布，用期望最大化算法确定高斯混合模型的均值、协方差等参数，然后计算每个特征在各个高斯分布中的概率，并划分到概率最大的一个类别（高斯分布）中^[2]。

Valley Seeking 寻谷聚类算法

在特征空间中寻找特征分布密度的极小位置，作为聚类的边界。算法的流程可以简述如下^³:

1. 确定初始半径，可以取为样本标准差的某个倍数；初始半径的设置影响聚类的数目；

2. 在该半径范围内，确定每个特征的近邻数目，以及近邻的集合。如果特征的数目为N，需要初始化一个N×N的数组来存放每个特征的近邻编号，计算完成后，可以根据最大近邻集合中元素的数目，删去多余的元素。

3. 如果一个特征的近邻数目，超过它每一个近邻各自的近邻数目，则该特征确定了一个特征密度峰。如此，设置一个最小近邻标准，可以找出全部峰。

4. 五个峰确定5个类别，给每个峰的近邻特征设置类别标记；

5. 将所有特征按照近邻数目从大到小排序；

6. 从近邻数目最大的特征开始，看其近邻中属于哪一类别的最多，就将该特征也归入这一类；如此，完成全部特征的聚类。

5.4聚类结果调整

移除未分类的锋电位

菜单Sort->Edit units …中，选择复选框Remove unsorted spikes。未分类的锋电位将从数据中删去。

合并分类

同样在菜单Sort->Edit units …中，在Merge Units框中输入欲合并的单元编号；例如，输入“2 3”将单元2和单元3合并为单元2，输入“0 3”删除单元3（其中的波形转为未归类状态）。

移除特征离散度很大的锋电位

可以首先通过轮廓线聚类方法，将要保留的特征划分到单元1中，然后按照上面的方法，先移除没有归类的边缘特征，再将单元1与0（未分类）合并。

6 绘图功能

6.1锋电位波形

菜单：Plot->Spikes：绘出所有锋电位波形；

菜单：Plot-> Sorted spikes：分布绘出各个类别的波形。每个波形图左下的数字标明了该分类中波形的数目。

6.2特征分布图

菜单：Plot-> Features：绘出特征分布散点图；

菜单：Plot->Sorted Features：画出分类后的特征散点图；如下图。

6.3时间间隔(ISI)分布

菜单：Plot->ISI …，打开如下的ISI参数设置窗口：

Min. Interval: 最小间隔时间；

Max. Interval: 最大间隔时间；

Number of Bins: 阶梯图的分段数目，数值越大，ISI曲线越平滑；

Units: 选择单元，只画出指定单元的ISI分布图。

6.4 间隔-时间分布图

菜单：Plot->ISI versus time…，绘制相邻放电间隔的时间分布图。参数设置对话框：

X Min. :横坐标（时间）最小值；

X Max.：横坐标最大值；

Units：选择需要进行绘图的单元编号。

6.5庞加莱图

菜单：Plot-> Poincare maps：该图以与前一个放电的时间间隔为纵坐标，与后一个放电的时间间隔为横坐标，标出每一个放电的位置。坐标轴的范围（即最大时间间隔）由ISI中的Min. Interval和Max Interval 指定。如下图。

6.6放电频次阶梯图

菜单：Plot->Fire Rate…：绘制放电的频次直方图。参数设置对话框：

X Min. :横坐标（时间）最小值；

X Max.：横坐标最大值；

Time per Bin：每个阶梯表示的时间长度；

Units：选择需要进行绘图的单元编号。

Add stimulation markers：选中该项后，可以在下面的Markers中输入刺激标记（对应数据文件中的simmark变量）。例如，刺激开始的时刻为49.89885s，结束的时刻为262.4599s，在Markers中填写如下：

放电频次分布直方图如下，刺激时间将以一个黑框标示出来。

6.7放电光栅图

菜单：Plot->Rasters …：绘制放电的光栅分布图。参数设置对话框：

X Min. :横坐标（时间）最小值；

X Max.：横坐标最大值；

Units：选择需要进行绘图的单元编号。

Add stimulation markers：与频次直方图中相同。

6.8相关函数图

菜单：Plot->CrossCorrelograms…：绘制单元间的互相关和自相关函数。参数设置对话框：

Pre-ref. Time：横坐标（时间）最小值，计算相关函数时信号的偏移范围；

Post-ref. Time：横坐标最大值；

Number of Bins：相关函数分段数目；

Units：参与计算的单元编号；

Ref Units：作为参考的单元编号。

6.9刺激前后放电直方图PSTH

菜单：Plot->Peristimulus Histogram … ：打开如下的参数设置对话框：

T. pre-Event：刺激事件前保留的时间长度；

T. post-Event：刺激事件后保留的时间长度；

Time/Bin：直方图的每段包含的时间长度；

Units：选择需要进行绘图的单元编号。

Stimulus Phase：选择刺激事件的相位，“on”表示刺激开始；“off”表示刺激结束；

Stimulus No.：当有多个刺激周期时，可以选择其中的一个刺激周期。设置为0时，将所有刺激周期的的PSTH叠加。

6.10刺激前后放电光栅图

菜单：Plot->Peristimulus Raster … ：打开如下的参数设置对话框：

T. pre-Event：刺激事件前时间长度；

T. post-Event：刺激事件后时间长度；

Units：选择需要进行绘图的单元编号。

Stimulus Phase：选择刺激事件的相位，“on”表示刺激开始；“off”表示刺激结束；

Stimulus No.：当有多个刺激周期时，可以选择其中的一个刺激周期。设置为0时，将所有刺激前后放电都画出。

在该组信号中，共有三个刺激事件，Stimulus No.设置为0时，三个刺激前后的放电都绘出，如下图。纵坐标名称中，E*表示事件的编号，U*表示单元编号：

部分参考文献：

^{1Mukhopadhyay S, Ray GC. A new interpretation of nonlinear energy operator and its efficacy in spike detection. IEEE Transactions on Bio-Medical Engineering. 1998 Feb ;45(2):180-7.}

2 （美）迪达等著，李宏东等译. 模式分类（原书第2版）[M]. 北京: 机械工业出版社，2003.

3 Zhang C, Zhang X, Zhang MQ, Li Y. Neighbor number, valley seeking and clustering. Pattern Recogn. Lett. 2007 ;28(2):173-180.

其它文献：

梁培基，陈爱华. 神经元活动的多电极同步记录及神经信息处理 [M]. 北京：北京工业出版社, 2003.

Lewicki MS. A review of methods for spike sorting: the detection and classification of neural action potential [J]. Net Comput Neural Syst, 1998, 9(4): 53-78.

Answer.com: Expectation-maximization algorithm

锋电位分析工具SPKtool

2008-08-06T18:06:00.002+08:00

SPKtool是一个在MATLAB平台上编写的神经元锋电位分析工具软件，具有锋电位检测、聚类以及放电序列统计分析等功能，可以满足神经信号分析的大部分要求。并且带有批处理功能，可以加快大量数据的分析速度。SPKtool集成了锋电位分析中常用的成熟算法，并且代码是开放的，方便使用者进行二次开发，根据需要扩展其功能。软件具有直观的图形化界面和详细的说明文档，使用方便。

其功能包括：

锋电位检测：

Ø 显示记录信号中锋电位峰值的分布曲线；

Ø 通过阈值检测锋电位；

Ø 通过原始信号、能量、非线性能量进行检测；

Ø 三种锋电位波形对齐方式：阈值、波谷、波峰。

Ø GUI交互式界面。

特征提取：

Ø 多种常用特征：PCA，峰值、能量等；

Ø 小波方法；

Ø 特征分别图，特征密度分布图等。

聚类分析：

Ø K-means；

Ø 模板匹配；

Ø 基于高斯混合模型的EM算法；

Ø Valley seeking；

Ø 通过绘制轮廓线，手动聚类；

Ø 对聚类结果进行合并，删除等。

动作电位串分析：

Ø 相邻动作电位间隔分布图(ISI)；

Ø Poincare map；

Ø 相关函数图；

Ø 锋电位分布光栅图；

Ø 放电频次直方图；

Ø 刺激前后放电直方图（PSTH）；

Ø 刺激前后放电光栅图。

软件及使用手册下载地址：http://groups.google.com/group/liuxq-works/web

界面截图：

如何将曲线以不连续的形式画出来

2008-06-18T20:35:00.001+08:00

最近有些“忙”， blog也没有什么内容好更新。弄点简单的东西，凑个篇幅。
MATLAB中的数据，用plot画出来一般是一条连续的曲线，怎么用将它画成离散的形式呢？
例如下面的连续正弦曲线，怎么画成后面的离散形式呢？

这个方法很简单，这个上图中的连续曲线：
x = 1:100;
y = sin(x/20);
figure, plot(x,y);
不连续曲线：
y1 = y(1:2:100);
y2 = y(2:2:100);
x1 = x(1:2:100);
x2 = x(2:2:100);
figure, plot([x1;x2], [y1;y2]);
简单点来说，就是将每段曲线的起点和终点的横纵坐标分别放在一个2×N的数组中。
在这个不连续曲线中，每个线段都是一个绘图对象，所以具有不同的颜色。这种方法的限制是每条连续线段必须包含相同的数据点数。否则的话，只能通过循环来逐段绘制。
实际上，如果数据中有NaN（即：Not a Number），该点是不会绘出来的，这样也可以形成不连续效果。
仍然是上面的曲线：
yn = y;
yn(3:3:100) = NaN;
figure, plot(x, yn)
得到下图。这种方法就没有了各段尺度相同的限制。简单点来说，就是在需要断开的数据点（包括横纵坐标）前面插入NaN。这样的不连续线仍然是作为一个整体的。

十六国魏晋南北朝更迭

2008-05-30T20:37:00.003+08:00

最近读夏曾佑的《中国古代史》，制作了一幅从西晋到隋即十六国和南北朝时期的更迭示意图。可能会有一些错误。

参考资料：
http://tieba.baidu.com/f?kz=45314607 ；
http://www5c.biglobe.ne.jp/~paruwees/Histoly_Goko.html
这里有一幅十六国的地理位置示意图：
http://ntpicc.51.net/cgx/images/5-16.jpg 。
SVG源文件：http://liuxqsmile.blogbus.com/files/12121506310.svg

Inkscape中文化工作

2008-05-21T14:35:00.001+08:00

Chinese translation of Inkscape interface

Inkscape 0.46的界面中共有3500项字符串，翻译的工作量很大。Chris Jia 在2006、2007年做了很多工作。由于Inkscape仍然在不断修改完善，到2008年初，我开始参与这项工作时，只有60%左右的完成度。经过一段时间的努力，在0.46版发布之前，我将这个比例提高到了89%。
我并不从事图形图像方面的工作，专业知识不足，因此汉化后的界面中肯定存在各种问题。大家在使用Inkscape过程中，如果对一些名词或描述的译法不满意，欢迎来信指明。
如果你觉得某个中文译法很别扭，想看一看相应的英文，应该怎么做呢？在中英文界面间切换需要修改locale(Windows需要修改“区域和语言设置”)，太麻烦。
简单一些的方法是查看Inkscape的界面翻译文件。到下面的地址下载中文界面文件zh_CN.po：
http://inkscape.svn.sourceforge.net/viewvc/inkscape/inkscape/trunk/po/
Linux上的编辑器一般都支持这种格式的文件，Windows上可以用Poedit 。打开这个文件，用搜索功能查找相应的中文字符，就可以同时找到相应的英文。实际上，界面汉化就是在这个文件上进行的，有些时候很难确定一些词句在界面上的具体位置，错漏也就在所难免了。

另外，我还翻译了一些Inkscape教程，已经可以在Inkscape官方网站上看到：
http://www.inkscape.org/doc/index.php?lang=en
另外，我的博客中也有相应的内容。

锋电位分析软件MClust简介

2008-05-02T19:11:00.001+08:00

MClust 是一个MATLAB上的锋电位分类工具包，用于四电极胞外记录信号的分类，也可以用于单电极的情况。程序可以直接读入Neuralynx、Axona等记录软件的数据，也可以自己编写程序，读入特定的数据。程序本身只有分类功能，因此读入的数据必须是提取的锋电位，不能是连续记录信号。可以设置PCA、幅度等多种特征提取方法。分类方法有三种，除手动外，还有BubbleClust和KlustaKwik 两种自动聚类方法。另外这个软件的文档也很详细。
打开该软件后，首先设置好特征，以及通道数目，选择数据加载方法，然后选中“Creat/Load FD files”，在弹出文件框中选择数据文件。如果没有错误的话，就可以进行聚类了。如果要重新开始分析，或者设置不同的特征，需要先选择“Clear workspace”，然后再从头来过。
手动聚类界面：

两种自动聚类方法需要借助该软件中的批处理功能，参考说明文档中的“Batch Processing”部分。然后再用MClust对聚类的结果进行二次处理。

寻谷聚类算法

2008-04-20T21:29:00.002+08:00

valley seeking clustering method
所谓的valley seeking聚类方法，就是在特征空间中寻找特征分布密度较小的位置，作为聚类的边界。如下图中的2维特征空间，形成5个分布密度较高的峰，如果能确定峰-峰间的谷的位置，也就等于确定了各个类别的边界，从而完成聚类。算法的目标是非常直观的。

这种聚类方法的流程可以简单描述如下:
1. 确定初始半径，可以取为样本标准差的某个倍数；初始半径的设置影响聚类的数目；
2. 在该半径范围内，确定每个特征的近邻数目，以及近邻的集合。如果特征的数目为N，需要初始化一个N×N的数组来存放每个特征的近邻编号，计算完成后，可以根据最大近邻集合中元素的数目，删去多余的元素。
3. 如果一个特征的近邻数目，超过它每一个近邻各自的近邻数目，则该特征确定了上图中一个特征密度峰。如此，设置一个最小近邻标准，可以找出全部五个峰。
4. 五个峰确定5个类别，给每个峰的近邻特征设置类别标记；
5. 将所有特征按照近邻数目从大到小排序；
6. 从近邻数目最大的特征开始，看其近邻中属于哪一类别的最多，就将该特征也归入这一类；如此，完成全部特征的聚类。

在上面的聚类结果中，右上角有几个较密集的特征没有归类，反而归类了旁边更分散的几个特征。出现这种问题的原因是这几个特征的形成一个小的局部峰值，这些特征的聚类次序在周围特征之前，而这时，周围的这些特征还没有获得聚类信息。

一辆伪山地车的悲惨遭遇

2008-04-19T17:08:00.001+08:00

A Bicycle's miserable life

2005年的时候，我决定丰富自己的运动生活，跑步之外，再加点别的东西，为了不至于还是在原地转圈圈，就选择了自行车，去体验人造风的猎猎吹拂，转一个更大的圈圈。
在自行车论坛上泡了几周，又去商场转了几圈。苦于钱包不鼓，就买了现在这辆，捷马TXT 5500，新车当然是很漂亮的，可惜当初没有合影留念，网上google了一下，又百度了一下，连个图片都找不到。全钢的车身很重，锻炼身体嘛，重才有效果！四连杆的车架，多了个起副作用的后避震，但看上去总还有一点“运动”气质。

拆掉前后挡泥板，拔高车座，“运动”气质又调高了几分。那时候实验室在一楼，车子放在房间里，每天看它几十次，喜不自胜，俨然一自行车运动高手。学校离紫金山近，往山脚下转悠几次，踩踩点，然后就开始赶着车轱辘往山上滚了。上天文台，头陀岭，再从青年路下来。不定期地在这条路上颠来倒去，一个更大的圈形成了。
说起来，这辆伪山地车应该像它的伪自行车爱好者主人一样骄傲了，出身贫贱，但也可以经常呆立在头陀岭索道边的平台上，俯瞰南京灰蒙蒙的城市景观。这通常是身价数倍于它的同类才能享受的待遇。从这个意义上说，它反倒是很幸运的，尽管上山时链条和齿盘咬牙切齿，下山时车架和前叉咣当作响。

当然，不幸还是主流，它的两个轮子不动的时间就跟主人的臀部不动的时间一样长。更糟的是，它被请出了房间，只能在学校这个大工地的车棚下面呆着，忍受灰尘和雨水，以及某些人试图让它重回正常生命轨迹的努力。为了抵抗这种努力，我给它配了两把锁。
莫名其妙地前胎破了，趁机去买了胶水和微型打气筒，自己将就着补好了。过了几个月，后胎又破了，这次锁上有了碰撞的痕迹。后来，不骑的时候，我都会把轮胎的气放掉，这下总扎不穿了吧。去年，自行车也开始主动配合，前胎自动放气，打次气可以支持个五六小时，正好够绕上一圈。
长期放在外面，生锈是难免的，座垫杆和曲柄上没有上漆，早早地就生了锈，只能经常用WD40喷一下，让锈别堆得太过分。轮圈和辐条也失去了光泽，“运动”气质也只能靠高昂的车座撑撑面子了。

一次取车的时候，发现后轮锁上套了一个塑料袋，里面有一些浅绿色的黏糊糊的东西，我有些莫名其妙，什么人吃了东西垃圾乱扔，还故意扔到我锁上？等到发现锁口上锈迹斑斑，钥匙也拧不动时，我才明白过来，撬不开拿酸液来腐蚀？这也想得出！
幸好我有WD40。
前几天再去看时，后轮上的锁怎么又变了？再一看，原来是上面的塑料套没了。又遭黑手了，还好两把锁给面子。不过，丢车的危险身已经像那些闪闪发光的部件一样越来越少，坚持到我离开南京应该不是问题了。

Windows共享Internet连接中的RPC错误

2008-04-15T18:10:00.002+08:00

在Window XP中，通过Internet连接共享可以与其它计算机共享网络连接，一拖多共同上网。在网络连接中启用该功能时，可能会弹出错误信息：RPC Server unavaliable（RPC 服务器不可用）。
解决的方法是在控制面板->管理工具->“服务”中启动“DHCP Client”，并将启动方式改为“自动”。
在我的机子上，未启动该服务时，并不总是弹出错误提示，但一段时间后，网络连接自动断开，连主机也不能上网。在网络连接中先“禁用”外网网卡，再重新“启用”后又恢复正常。启用该服务可以解决这个问题。

参考链接:http://www.legend.ws/blog/tips-tricks/rpc-server-unavailable-ics/

本来想用Windows自带的ICS与另一机子共享上网的，弄了半天还是不行，一个很奇怪的现象是，主机可以ping通客户机，而客户机只有在主机 ping客户机后很短的时间内才能ping通主机。试了网上的各种解决方法，打开各种服务，关闭防火墙等等，还是不通。最后装了Sygate就马上解决问题了。

MATLAB中的正态分布检验

2008-04-13T13:50:00.001+08:00

要对一组样本进行正态性检验，在MATLAB中，一种方法是用normplot画出样本，如果都分布在一条直线上，则表明样本来自正态分布，否则是非正态分布。
MATLAB中也提供了几种更正式的检验方法：
kstest Kolmogorov-Smirnov 正态性检验，将样本与标准正态分布（均值为0，方差为1）进行对比，不符合正态分布返回1，否则返回0；该函数也可以用于其它分布类型的检验；
lillietest Lilliefors test。与kstest不同，检验目标不是标准正态，而是具有与样本相同均值和方差的正态分布。
jbtest Jarque-Bera test。与 Lilliefors test 类似，但不适用于小样本的情况。

MATLAB中的一些小技巧(2)

2008-03-31T19:57:00.001+08:00

MATLAB中的一些小技巧(2)

1. Ctrl+D打开子程序
在MATLAB的Editor中，将输入光标放到一个子程序名称中间，然后按Ctrl+D可以打开该子函数的m文件。当然这个子程序要在路径列表中（或在当前工作路径中）。实际上该快捷键对应右键菜单中“Open Selection”命令。

2. 列出依赖函数
在Editor菜单中，选择Tools->Show Dependency Report，可以列出当前m文件调用的子程序。如果需要将程序移动到别的电脑上，使用这个工具可以确保不会漏掉一些东西。这个工具也适用于文件夹。

3. 显示鼠标数据内容
调试程序过程中，程序运行到断点后，将鼠标移到一个变量名上，将会弹出一个方框，显示该变量的内容。不在调试状态下，怎么使用这个功能呢？在File菜单中打开设置选项（Preferences），选择 Editor/Debugger -> Display，在 "General Display Options" 中选择 "Enable datatips in Edit Mode "。

4. 单元数组的基本操作
清除cell中的元素

cellname(n1:n2) = [];

合并如下的两个cell数组

a = {'str1', 'str2', 'str3'};
b = [1:3];

c = [a(1:2), b, a(3)];

得到：c = { 'str1' 'str2' [1x3 double] 'str3'}

而 c = {a(1:2), b, a(3)}

则得到 c = {1x2 cell} [1x3 double] {1x1 cell}

cell array本质上是数组，前一种操作方式与普通数组的操作方式相同，而后一种方式则是将逗号分隔的各个元素分别一个cell，整体作为一个cell数组。

cell的两种引用方式
a(1) 得到 {'str1'}, 因为a是一个cell数组，所以其中的每一个元素都是一个cell，
而a{1} 得到 'str1'。

5. MATLAB 中的字体平滑
在File菜单中打开设置选项（Preferences），在 Fonts 中，选择底部的 “Using antialiasing to smooth desktop fonts” ，MATLAB 中默认字号比较小，可以调大一些（可以在下拉菜单中输入数字，更精细地控制文字大小）。有衬线字体打开反锯齿后可能会看起来“发虚”。

6. 同步缩放坐标轴linkaxes
linkaxes命令，可以将一个figure中的各个axes指定相同的坐标范围，实现各个绘图同步放大、缩小。在缩放一个绘图时，其它也会同步缩放。可以同时关联x、y轴，也可以只关联一个坐标轴。

7. 鼠标绘图
MATLAB Help中的Interactive Plotting一文提供了一种在figure中用鼠标绘图的方法，主要使用ginput函数。

8. GUI程序关闭按钮的callback
为GUI窗口右上角的关闭按钮添加回调函数：在GUIDE中打开GUI，在没有控件的地方单击，选中整个figure，右键菜单中选择 View Callbacks -> CloseRequestFcn，m文件中将自动添加该回调函数，输入代码就可以了。

开源矢量绘图软件Inkscape 0.46 发布！

2008-03-25T20:08:00.001+08:00

Inkscape社区3月25日发布了该矢量绘图软件的最新版本。Inkscape 0.46 最重要的特性是引入了原生的PDF支持，为PDF文档的编辑提供了一种简单的开源解决方案。

新版本的执行效率显著提升，同时引入了众多新特性。其中包括：

对话框可以停靠在编辑窗口中；
渐变的修改可以完全在画布上执行；
新的油漆桶工具，可以为封闭区域填色；
新的立方体工具，有助于使绘图更符合透视规律；
新的扭曲工具，为路径和对象提供了更直观的编辑方式；
新的实时路径效果，可以为路径提供“笔刷”等多种灵活的路径效果；
改进的色彩管理工具，除sRGB之外，还提供对其它色彩模式的支持；
实现了多数 SVG 滤镜，并提供了强大的滤镜堆栈编辑界面。

== 下载Inkscape 0.46 ==

Inkscape 0.46 已经包含在Ubuntu Hardy (8.04)中，可以正常安装；
Ubuntu Gutsy (7.10) 可以通过如下方式安装：
在Admin : Software Sources : Third-Party Software: 中添加
deb http://ppa.launchpad.net/inkscape.testers/ubuntu gutsy main
deb-src http://ppa.launchpad.net/inkscape.testers/ubuntu gutsy main

Macintosh OS X 用户可以下载 Leopard 通用安装包：
http://sourceforge.net/project/showfiles.php?group_id=93438

Fedora 、Debian 以及 Windows 等其它平台的安装包也将于近期发布。

Inkscape网站
http://www.inkscape.org/

完整的Release Notes见这里：
http://wiki.inkscape.org/wiki/index.php/ReleaseNotes046

新功能截图及介绍：
http://inkscape.org/screenshots/

MATLAB 中GUI子程序的参数传递

2008-03-25T13:40:00.002+08:00

输入参数传递：
比如子GUI的名称为subGUI, 设想的参数输入输出为：[out1, out2] = subGUI(in1, in2)
在subGUI的m文件中（由GUIDE自动产生）：
1.第一行的形式为：function varargout = subGUI(varargin)
该行不用做任何修改；varargin 和 varargout 分别是一个可变长度的cell数组(MATLAB帮助文件中有说明)。输入参数in1和in2保存在varargin中，输出参数out1，out2包含在varargout中；
2.在subGUI的OpeningFcn中，读入参数，并用guidata保存，即：
handles.in1 = varargin{1};
handles.in2 = varargin{2};
guidata(hObject, handles);

返回参数的设置：
1. 在GUI子程序的OpeningFcn函数的结尾加上uiwait(handles.figure1); figure1是subGUI的Tag；
2. subGUI中控制程序结束（如"OK”和"Cancel"按钮）的callback末尾加上uiresume(handles.figure1)，不要将delete命令放在这些callback中；
3. 在子GUI的OutputFcn中设置要传递出去的参数，如 varargout{1} = handles.out1；varargout{2} = handles.out2;末尾添加 delete(handles.figure1); 结束程序。

在GUI的OpenFcn中，如果不加uiwait，程序会直接运行到下面，执行OutputFcn。也就是说程序一运行，返回值就确定了，再在其它部分对handles.output作更改也没有效果了。

加上uiwait后，只有执行了uiresume后，才会继续执行到OutputFcn，在此之前用户有充分的时间设置返回值。

通过以上设置以后，就可以通过 [out1, out2] = subGUI(in1, in2) 的形式调用该子程序。

在一个GUI中调用另一个GUI时，主GUI不需要特别的设置，同调用普通的函数一样。在打开子GUI界面的同时，主程序还可以响应其它的控件。不需要担心子GUI的返回值被传错了地方。

磁共振成像MRI的简单介绍

2008-03-22T16:57:00.001+08:00

磁共振成像MRI的简单介绍

如果你看美剧《豪斯医生》，对MRI这个词一定感到很熟悉，在豪斯的诊断中几乎是必不可少的环节。那么什么是MRI，它又是怎样工作的呢？

因为工作的关系，我曾经接触过MRI一段时间，简单查阅过一些文献。在这里，我想从一个外行的角度，对MRI做一些形象的描述。

MRI的全称是Magnetic Resonance Imaging，即磁共振成像，它在医学上的应用有些类似于CT和X光。MRI设备对人体的某些部位（如大脑等），进行扫描，得到氢原子核H在组织中的密度分布。不同的组织，水的含量也是不同的，病变的组织的水含量又与正常的组织不同，通过对比，来了解组织的状况，以作为诊断的重要依据。

MRI图像中，灰度反映了组织中水的含量http://www.umshp.org/rs/mrict.htm

那么MRI设备是怎么得到人体内部的氢原子的含量呢，当然就像它的名字那样，通过核磁共振，本来这种成像方式应该成为核磁共振成像，为了消除病人的恐惧，去掉了里面的“核”字。
氢原子核实际上就是一个带电的质子，原子核绕其中心旋转，我们可以认为在原子核上产生了电流，这个环形的电流就使原子核产生了磁场。可以将H原子核看成类似于地球，具有南北磁极，同时也在不停地绕自身旋转（自旋）。

将这个小地球放入一个更强的磁场时，由于自转方向与磁场力方向不同，原子核在自转的同时，将绕着磁场方向公转，这种现象叫进动。磁场强度不同，进动的情况也不同。实际上由于组织内部分子的热运动，只有极少部分的氢原子核能够发生这种进度，但正是这极少数的原子核才是磁共振成像成为可能。当然磁场的强度越高，发生进度的原子核越多，对成像也越有利。

我们知道，电磁波也是有能量的，当用电磁波作用于发生进动的原子核时，原子核吸收电磁波的能量，进动的形式会发生变化。不过，由于原子核的量子特性，原子核的能量并不是平滑变化，而是会在特定的能级间跳跃，所以只有特定频率的电磁波才会被原子核吸收，发生能量的跃迁，也就是核磁共振。用来激发核磁共振的电磁波的频率很高，属于射频范围。磁场强度不同，氢原子核的进动就不同，吸收的电磁波的频率也就不同。

处于较高能级的原子核是不稳定的，在激发电磁波取消后，它会重新回到低能量状态，同时将能量仍以电磁波的形式释放出来，发射电磁波的频率同激发的频率是一致的。释放能量的原子核越多，产生的电磁波的强度也越高。通过检测设备，我们就可以将这个释放的电磁波检测出来，根据这个电磁波的性质，我们就可以判断出氢原子核的相对密度。

到这里我们应该很清楚了，要测出组织内的氢原子核的密度，需要一个很强的磁场，一个射频发射装置，一个射频检测装置，当然还要有相应的控制装置，计算机软件处理等。但这还是不够的，在这种情况下如果将一个组织放入其中扫描，由于磁场是均匀的，组织内每一处的氢原则核都发生进动，在射频下都共振，检测出来的信号也是全体组织的，将整个脑袋放进去也就是一个大亮点而已。
解决的方法是，在主磁场的基础上，再叠加一个空间上均匀变化的磁场（称为梯度场），使组织内不同位置的氢原子核的进动情况都不同，然后用不同频率的射频分别激励，各个位置的原子核依次响应，放射出不同频率的电磁波，空间上的差异就检测出来了。

通常在医院看到的MRI设备是下面这样的，一个大圆筒和一个病人躺的工作台。你肯定能猜到这个大圆筒是干什么的了。没错，它就是一个超级大磁铁。可以产生磁场的方法无外乎两种，永久性的磁铁和电磁效应产生磁场。永久性的磁铁磁场弱，而且又大又沉，一般用在低端设备上，场强低于0.3特斯拉。而更高场强的设备只能用电来驱动。但较高的场强，不仅需要很大的功率，线圈的散热也很成问题。
现在最多的是通过超导体产生的永磁场。比如图中这台GE的MRI，场强为3特斯拉，而地球磁场的强度仅仅为0.000005特斯拉。

图片来自：http://www.inst.msstate.edu/3tmri

超导需要极低的温度，因此MRI的励磁线圈浸泡在液态氦中，温度为4.2K，在线圈中通电后，切断电源，电流就在线圈中永无休止，永不衰减地流动，并且产生强大的磁场。人们早就发现了在更高温度下具有超导性的材料，为什么还要还用这么低的温度呢？原来，大部分超导材料在较强的磁场作用下都会失去超导性，直到1970年代，人们才发现了可以在强磁场中工作的超导体。
超导线圈要放在液氦中，而盛放液氦的容器与环境总是存在一些热交换，会导致液氦不断蒸发。为了减少液氦的消耗，盛放液氦的容器又放在液氮中，液氮的温度77K比液氦高，但比室温又低多了，液氦与环境的温度差降低，蒸发量也就下降了。当然，液氮也会不断蒸发，但它的价格比较便宜。

前面提到，需要一个空间上变化的磁场，来确定成像的位置信息。这个就是梯度磁场，梯度磁场是用普通线圈通过电磁效应产生的，强度要远小于主磁场。在快速成像中，梯度磁场实际上是变化的，这引起金属器件受力变化，产生变形，这是噪声的主要来源。梯度磁场的线圈位于主磁场的内侧，更靠近圆筒内壁的地方。

图片来自MRI: A Guided Tour

射频线圈的形式多种多样。用于大脑成像的射频线圈一般是圆筒形，扫描时套在头部。由于射频的发射和接受在时间上是分开的，所以可以用一个线圈同时充当发射线圈和接受线圈，在两种模式下交替工作。也可以用两个独立的线圈。射频信号被线圈接受后，经过复杂的计算机处理，生产影像。
MRI的主要缺点是，病人被限制在一个狭小的空间里，扫描的部位要一动不动，病人体验不够好。也有一些是开放式的结构，不过一般场强较低，能力有限。

参考资料：关于低温的概念、技术和应用

檢查脊髓的利器：核磁共振攝影 ( MRI )

黄继英梁星原：磁共振成像原理

Nonlinear energy:非线性能量

2008-03-18T13:10:00.002+08:00

Nonlinear energy:非线性能量
对神经信号进行锋电位提取时，一种方法是直接对采集信号设置阈值，提取波形；另一种方法是对信号的能量分布曲线设置阈值。常用的能量计算有两种，线性能量（Linear energy）和非线性能量(Nonlinear energy)。
线性能量相当于直接对数据点求平方，相当于即时功率。而非线性能量是指，幅度平方与频率平方的乘积，即(Amplitude^2) * (Frequency^2)，由于频率的引入，这种能量计算方法可以放大突变信号（频率较高）的能量，在锋电位提取中可能会获得更好的效果。
非线性能量的计算一般简化为：

NEO = v(t)²-v(t-1)*v(t+1)。v(t)表示 t 时刻的信号幅度，NEO 即非线性能量运算符 Nonlinear energy operator。

参考：http://www.nbb.cornell.edu/neurobio/land/PROJECTS/spikeSort/index.html
一个典型的锋电位波形的非线性能量：

不过，非线性能量对信号中的毛刺非常敏感，如下图：

Wallpaper：一张桌面背景

2008-03-13T21:06:00.001+08:00

Wallpaper：一张桌面背景
我自己做的一个桌面墙纸，内容取自《美女与野兽》（Beauty And The Beast）的招贴画：http://www.plos.asn.au/Productions/Beast/Beast.html

局部放大：

Rubicon:从苏拉到奥古斯都

2008-03-06T20:52:00.001+08:00

Rubicon:从苏拉到奥古斯都

图片来源：http://tspl.xzcn.com/book/25960

公元前1世纪，罗马将领们继续率领军团在地中海沿岸所向披靡。共和国的公民们或许还没与意识到，正是这些伟人们断送了他们的共和、自由和民主。每一个强盛的国家都伴随着对荣誉的追逐，但荣誉和共和国并不总能够相处无间。长期以来，罗马公民的资格仅局限于罗马周围的局部地区，而军事领袖的征服范围扩大了，个人的力量有了与共和国一争高下的机会。

在扫荡了亚平宁山区的叛乱后，两位将领马略和苏拉开始争夺下一界执政官的席位，胜者还将会获得在亚洲对本都战争的指挥权，这也意味着更大的荣耀和财富。年迈的马略尽管已经6次出任执政官，拥有了无数的战功和荣誉，却还是不满足，但这次年轻的苏拉赢了。马略心有不甘，与保民官策划了一场暴动，撤了苏拉的执政官搭档，苏拉自己勉强保住了位置，心存一些希望，能够在东方的战争中东山再起，挽回颜面。然而，对手不愿给他这样的机会，接替他帅印的人拍马赶到，今年的指挥权都没了，还想明年？
苏拉作出了选择，打破了共和国传统，率领5个军团回到了罗马，军官中只有一个响应，士兵中却没有人反对。二十年前，马略取消了军团士兵的财产限制，由国家发给军饷，现在他们只忠诚于发军饷的首领，而不是共和国。
苏拉虽然打倒了自己的对手，却没能扶植起自己的亲信，他还没有无耻到操纵选举的地步。选举产生的两个执政官没有一个是他的盟友，他也只能让他们发誓不要做得太过份，就登上了东去的舰队。
两个执政官很快发生内讧，其中一个秦纳很快和马略勾结起来，率领军队回到罗马，将同事的脑袋摆在了桌上。马略杀戮一番后，第七次当选了执政官，不久死于任上。秦纳一头独大，连续三年坐在执政官位子上，等苏拉的返回。
尽管还没有端掉本都王米特拉达提的老窝，苏拉也不能再对罗马坐视不理。秦纳不是久经战阵的苏拉的对手，很快被干掉了，在消灭秦纳和马略余党的过程中，富人克拉苏和年轻的“大将军”庞培也粉墨登场，站在苏拉一边。前两次军队入罗马只能算逼宫，而公元前83年这次，战火最终烧到了罗马的城门下，随后是对数千名俘虏的屠杀。
大权在握的苏拉少不了一番清洗，好给喽罗们中饱私囊的机会。作为秦纳女婿的凯撒也牵涉其中，苏拉让凯撒休妻，后者不从，宁愿逃跑。他当然也忘不了保民官所带来的麻烦，让他们不再具有推翻元老院决议的权力。尽管苏拉按照自己的意志重复分配共和国权力，但他却认为自己是为了挽救共和国，不让马略之流有兴风作浪的机会。五年后，苏拉去世，但他的手下继续把持着罗马。

苏拉在罗马敲敲打打，年轻的庞培正转战非洲，对马略余党赶尽杀绝。胜利回国后，又自愿去西班牙剿匪，元老们当然高兴让这个危险的25岁毛头小子走得越远越好。克拉苏也巴不得他在那个山险水恶的地方摔个大跟头。但是，西班牙渐渐安静下来，庞培的声势日渐壮大起来。克拉苏坐不住了，虽然手里钱很多，但胜利的荣誉是买不到的。
公元前73年，斯巴达克思起义了。一年前，卢卡拉斯率领罗马军队去了亚洲，罗马空虚。克拉苏顺利地从两位执政官手中拿到了指挥权，最终的结果大家都很清楚，斯巴达克思被克拉苏赌在了亚平宁半岛的脚后跟上。突围北上再次被围，全军覆没。庞培正从西班牙班师，赶上灭掉了最后一支义军，然后写信给元老院，吹嘘是自己救了罗马，卡拉苏狂怒，为了宣传自己的功绩，俘虏全部钉上了阿庇安大道两侧的十字架，绵延100多英里。
然而，卡拉苏的战争只是给罗马遮丑，庞培在西班牙的征服却是罗马的荣耀，元老院当然更愿意鼓吹后者。随后，克拉苏和庞培共同当选执政官。苏拉派的势力被削弱，保民官的否决权又回来了。

前74年，苏拉派的卢卡拉斯领军东征，希望能彻底吃掉东山再起的米特拉达提。卢卡拉斯占据了本都，米特拉达提逃到亚美尼亚避难，卢卡拉斯又将亚美尼亚打了个落花流水，但米特拉达提仍然在逃，卢卡拉斯被牵着团团转，五年后，他的士兵受不了了，反苏拉的人开始在家里抄他老底。他只好转攻为守，坐看老米又冒出头来。又两年后，他被庞培取代。
在此之前，庞培刚以三个月的时间扫荡了肆虐地中海数十年之久的海盗。本来，地中海的统治者是迦太基和希腊城邦，但罗马打残了这些国家，消灭了他们的舰队，而自己又没有建立起永久性的海上力量，海盗们才应运而生。卢卡拉斯只好眼睁睁地看着数年辛苦拱手让人，颜面扫地地回到罗马，从此意志消沉，沉溺于吃喝玩乐。
庞培将米特拉达提重新赶到山沟里，将亚美尼亚和叙利亚转为罗马的附属国。前64年，又南下攻入耶路撒冷，在前往佩特拉途中，传来了米特拉达提的死讯，这位屡败屡战的强人被自己的儿子囚禁，服毒未果后让卫兵刺死了自己，尸体成了儿子向庞培投降的见面礼。在米特拉达提的遗物中发现了一件亚历山大的红色斗篷，穿在身上，庞培觉得很合身。东方平定了，庞培准备回罗马，向公民们展示自己的斗篷。

前63年，凯撒靠借贷疯狂贿赂选民，获得了大祭司的职位。声名狼藉的喀提林依法炮制，竞选执政官，却败在了律师出身的西塞罗手上。眼看下一年依然不举，这下他到了破产的边缘。在社会的下层，同样有许多人负债累累，无以度日，相同的处境拉近了他们之间的距离。很快，西塞罗就察觉了他的动向，为他安排了罪名：喀提林打算来一场暴力革命，推翻共和国。喀提林逃出了罗马，拼凑起一支杂牌军，但在庞培回国前几个月就被消灭，他的同谋也被处死。西塞罗则成了共和国的拯救者。

在这个传统观念分崩离析的时代，既有像卢卡拉斯那样一蹶不振的反面人物，当然也有恪守共和国光辉理念的卫道士，加图就是这样的代表人物。通过对传统价值信条的固守和宣扬，他成了罗马正统的准绳，他的敌人也就是共和国的敌人。甚至连庞培这样的伟人，也渴望能获得他的尊重。
回到罗马后，庞培获得了第三个凯旋式，荣耀达到了顶峰。庞培离了婚，公开希望能和加图的家族扯上姻亲关系，然而得到的却是公然的拒绝和无情的嘲讽。庞培用从东方夺来的财富开始修建罗马第一座石质的大剧院，渴望自己的威望永存。共和国怎能容忍伟人？加图和卡拉苏团结起来，极力让他从公众的视线中消失。庞培当然不会坐以待毙。他需要一位重量级的政治盟友，而后者正远在西班牙。
一年前，凯撒担任了西班牙的行政长官，这位曾经站在罗马时尚最前沿的新新人类出乎预料地在任上干得相当不错，挣了不少钱，还得到一次凯旋式。现在，他知道自己的机会来了，前60年，为了参加执政官的选举，他放弃了凯旋式。在庞培的鼎立支持下，凯撒有大量的钱可以派给选民们，加图没有办法，也只好容忍己方的候选人依葫芦画瓢。双方的钱都没有白花，各占据一个执政官席位。
但凯撒很快就将克拉苏也拉到自己一方来，前三头联盟悄悄地成立了。尽管克拉苏对庞培很不待见，但凯撒安排他和庞培负责老兵的安置，显然大有油水可捞。而站在加图后面又可以得到什么呢？凯撒从执政官位置下来后，加图本来只希望他在意大利境内陪土匪玩玩，而现在他取得到了一个长达五年的总督职务，负责意大利北部高卢的大片领地，边境外还有无数的荣誉等他去摘取。
加图的优势在于控制了传媒，给当权者的一切都蒙上一层非法和卑鄙的外衣。凯撒将自己的女儿嫁给了庞培，加图当然会居高临下地评价一番。凯撒对别人的谩骂不屑一顾，庞培可就没有那么坦然了。传统的荣誉观念在庞培身上根深蒂固，这其中当然不能缺少共和国精英们的尊重和肯定。

在凯撒的成全下，贵族出身的大混混克洛迪乌斯被贬为平民，从而当选为保民官。他以底层民众的代言人自居，煽动暴力事件是他的拿手好戏，曾经指控过他的西塞罗被迫流亡。克拉苏和他的关系还不赖，但庞培却被整得忍气吞声，甚至整天缩在家里，元老院则在一旁看热闹。好在克洛迪乌斯的一年任期很快就结束了，庞培支持的保民官米洛走马上任，组织起人马和克洛迪乌斯分庭抗争。在这场争斗中，庞培和克拉苏间的矛盾开始浮上水面。反三头的元老以为机会来了，开始蠢蠢欲动。庞培成了暴力事件的替罪羊，再次陷入低谷。看着郁闷的庞培，元老们很开心，准备再拿凯撒开刀，拿掉他的高卢总督职位。
然而，三头们很快再次团结起来，甚至克洛迪乌斯也被拉到了他们的阵营中来，庞培和克拉苏再次和解，联手竞选前55年的执政官，有克洛迪乌斯在旁出手，连加图都可以被打回家去，还会有什么能挡住他们？
两位执政官的回报也是巨大的，他们分别获得叙利亚和西班牙的总督，任期同样为5年，而凯撒的任期结束后将再获得5年。两位总督的权力甚至包括了征召军队，开战，媾和等。以前这些权力只属于元老院，现在最杰出的公民成了共和国最大的敌人。
克拉苏任期一结束，对东方的掠夺就开始了，庞培参观过的耶路撒冷圣殿被洗劫一空，随后又度过幼发拉底河，进入帕提亚王国。事实证明，他的军事天赋远赶不上敛财天赋，前53年，他的部队被围，儿子战死，自己也被骗去谈判，混乱中被杀。
克拉苏死后，在罗马庞培一头独大，为了平衡他的势力，加图支持被庞培抛弃的前保民官米洛竞选执政官。米洛的死对头克洛迪乌斯岂能容忍，冲突中，克洛迪乌斯被米洛杀死，他的支持者为他举行的火葬，连元老院也被烧得精光。加图没有办法，只能求助于庞培，将军队开进罗马，两帮乌合之众作鸟兽散。庞培再次把握了罗马的权柄，这次连加图也不得不表示出一些敬意。

几乎整个公元前50年代，凯撒转战西欧大陆，不断将财富和荣誉送回罗马，同时也提醒健忘的公民们，不要忘了自己的存在。前52年，经过阿莱西娅城的决战，高卢领袖韦辛格托利克斯投降，凯撒最终平定了高卢，准备回罗马，然而迎接他的会是什么呢？
如果他放弃军队的指挥权，以平民的身份回去，很可能像庞培当年一样，被加图等踩在脚下。而且，庞培之前可是没有什么把柄的，而凯撒的违法勾当却太多了。
庞培尽管曾是凯撒的女婿，但比他年长几岁，一直将凯撒当成一个晚辈看待。现在情况不同了，手握重兵的凯撒对自己地位是一个巨大的威胁，他开始和加图站在一起，要求凯撒交出军队的指挥权。凯撒的对策是扶植起支持的自己的保民官，元老院一作出不利于自己的决议，保民官立即否决掉。这样反复几次后，元老们失去了耐心，庞培的军队再次进入罗马。凯撒的保民官乔装溜出了罗马，北上拉韦纳，凯撒和他的13军团正驻扎在那里，注视着罗马的一举一动。
接下来发生的事情就更广为人知了，凯撒迅速南下，庞培带领元老院退出罗马，后又被凯撒追到希腊。前48年，庞培本来有希望将凯撒困死，但在元老院的煽动下，他发起了决战，被凯撒久经沙场的军团打败，逃亡埃及时死在托勒密国王手里。凯撒随后赶到，扶植起国王的姐姐克里奥佩特拉。随后在亚历山大被困了一段时间，解围后，很快平定了亚洲的叛乱。前46年，击败了在非洲的共和派，加图绝望自杀。
凯撒得到了一个为期十年的独裁官任期，以及任命所有的行政官员的权力，在共和国之后，还从来没有一个人获得如此集中的权力，前44年，又被元老院任命为终身独裁官。这下，共和国的遗老们彻底坐不住了，一个月后，他被刺死在元老院内。

争权夺利的场面又出现了，但已经没有了共和派的位置。西塞罗曾一度屈服于凯撒的强权，但现在他成了共和派在罗马的领军人物。利用他的辩才和声望，他希望能够扶植起可以和安东尼对抗的人物。凯撒的两位老将被选为执政官，凯撒年轻的继承人屋大维也受到吹捧。
不幸的是，在和安东尼的战斗中，两位执政官同时战死。战败的安东尼随后和同为凯撒副手的李必达同盟。小凯撒没有希望战胜安东尼，于是调转方向，率军进入罗马，西塞罗的希望破灭了。随后三头政治正式建立起来，将罗马帝国分而治之，元老院和公民们只能靠边站。为了养活三头数目众多的军团，对罗马的掠夺又开始了，富人一个个成了国家公敌，人头落地，财产充公。三个人甚至还纳了投名状，各将一个自己人列上名单，安东尼的一个叔叔，李必达的一个兄弟，以及被屋大维曾成为“父亲”的西塞罗。
共和国最后的一批追随者，刺杀凯撒的布鲁图斯等人倒在东方的战场上，在那里，安东尼扮演起传统的总督角色。而在罗马，屋大维完全承担起破坏共和国的使命，为大量老兵们安置土地的计划使大批农民失去土地，意大利盗匪横行，在加上庞培的儿子率领的海盗控制了海上航线，罗马饥荒严重。经过几年的努力，屋大维打败了海盗，意大利的情况开始好转，共和国时的种种做法似乎又开始上演了，小凯撒俨然成了共和主义者。
而安东尼似乎完全被东方的传统征服了，开始梦想像埃及女王那样加入神的行列，一个高贵的罗马人被东方人同化了！罗马人最后的一点自尊心被屋大维激励起来，前31年，经过亚克兴角的海战后，安东尼和托勒密王国一起走到了历史的尽头。
小凯撒再无对手，前27年，他宣布恢复共和国，但在奥古斯都的称号下，仍然把权力牢牢地抓在手里，直到公元14年死去。在他之后，一个人的性格决定国家命运的东方模式在罗马正式登场了。

《卢比孔河:罗马共和国的胜利和悲剧》, [英] 汤姆.霍兰著，杨军译。
这本身还是不错的，语言比较流畅，趣味性较强。如果对这段时期的历史不是很熟悉，还是很值得一读的。

图片来源：http://ch.sysu.edu.cn/hope/sites/greeceroman/ppt/map.html

VibraCapture: 简单的振动数据采集分析软件

2008-03-01T17:44:00.002+08:00

VibraCapture: 简单的振动数据采集分析软件

我自己用LabVIEW做的一个简单的振动数据采集和分析程序，几个月前就已经放弃，不再继续开发。但毕竟也花了不少的精力，不舍得完全扔掉，在这里把代码共享出来，供有需要的朋友参考。程序的功能很不完整，也没有经过详细的测试。开发工具是 LabVIEW 7 Express。

具有的功能：
1. 触发式与连续数据采集，设计目标是针对PCI4472采集卡；采集参数控制：采样频率，采样点数，平均次数，保存参数和数据等。窗函数只使用了LabVIEW内置的几种。这一部分只在仿真卡上部分测试，可能问题较多。
2. 多窗口数据显示：可同时显示1、2、4或者8个窗口，在每个窗口可以独立配置其显示内容，包括：数据通道，时间历程，功率谱，频率响应函数等。在一个或两个graph时，可以调整光标。
3. 连续小波变换识别模态频率和阻尼比。

程序的使用：
程序的目的是进行一个模型的完整的振动数据采集，因此在此之前，应该有一个模型的测量方向表（.dirc文件）（程序中附有一组数据，可供参考）。在Parameters 面板中将 Project Path 指向测量表的位置。右侧的 Node Control 中将自动读入该测量表。
选择 Apply 后，采集参数保存在.cfg文件中，ICP的参数不保存（没有测试）。
进入 DAQ 面板，选择 Node Control 中的项目，如果该节点已有采集数据，则在graph中显示。选择START DAQ开始采集数据。数据自动保存在测量方向表所在的文件夹，文件名与节点名相适应。注意在底部方框中为graph设置合适的参数。当前graph的背景为紫色。
在 CWT Analysis 面板中，连续小波变换的程序是用 MATLAB 编写的，计算机中需要安装有MATLAB才能调用。选择节点，自动开始计算小波系数（可能需要等待几秒钟），在左上角的Graph中将光标移动到某一个峰值频率处，点击 Estimate，估计频率和阻尼。可以调节Morlet小波的参数，以及参数估计中曲线拟合的范围。

所有保存文件，包括参数文件、数据文件实际都是以文本格式存盘，可以用记事本打开查看。
一个简单的流程图

下载地址： http://liuxq-works.googlegroups.com/ 或 http://download.csdn.net/user/liuxqsmile
屏幕截图

惠普HP LaserJet 5L 打印机出纸口卡纸的维修

2008-02-23T21:03:00.001+08:00

惠普HP LaserJet 5L 打印机出纸口卡纸的维修

这种机子虽然已经很老了，可能还有不少人在用。出纸口卡纸的时候，一般都可以拽出来，但要是卡的太紧，把纸拽碎了，留在里面的就很难取出。要送修就太浪费时间和金钱了。下面介绍一下这种情况时的处理方法，也不是很复杂，需要一些一字和十字的螺丝刀。
原文地址，http://www.elliotfan.com/temp/HP.html ，所有图片也来自该网站，我加了一些中文注释。此文章中的方法可能会对你或你的设备造成的危害和损失！原文作者和翻译者对此概不负责;-)

1. 切断电源，最好等机器冷却下来后再继续操作。
2. 取出硒鼓，避光保存。注意不要碰到硒鼓下面的滚筒，皮肤表面的油性物质可能对其造成永久的损伤。

2. 将前盖拉开，按下左侧与底座连接的地方，将前盖取下。

3. 取下背板上的三个螺丝，背板和其它部分还通过四个卡子勾在一起，两个在顶部，用螺丝刀伸进去顶一下就开了。另两个比较隐蔽，可以将螺丝刀从底部的空隙伸进去，将卡扣撬开。将背板上撑纸的板拉起，取下背板。

4. 拧下打印机顶部的两个螺丝，取下左下方的扩展内存盖，剩下的塑料机体盖通过前部的两个暗扣与机体相连，用螺丝刀从底部空隙伸进去将其撬开，取下塑料机体盖。

5. 将出纸滚轮两端的塑料卡子向内弯，旋转90度，可以很容易地将该滚轮取下。

6. 按下右侧碳粉盒导轨上的卡子，将导轨取下。

7. 拧下里面的螺丝，取下加热元件的盖子，到这里一般都可以把碎纸取出来了，如果不行，继续拆。

8. 取下加热组件上的金属盖。金属盖下面是用弹簧撑着的，拧下螺钉时最好能记下螺钉旋进的深度（我觉得这里不需要拧紧，可能有一个很小的间隙约0.2mm左右，最好记一下深度，待会恢复时也拧同样的深度）。

9. 清理完毕后，按照相反的顺序将各部分再装好。

MATLAB中的一些小技巧

2008-02-19T18:43:00.001+08:00

MATLAB中的一些小技巧

1. Ctrl+C 中断正在执行的操作
如果程序不小心进入死循环，或者计算时间太长，可以在命令窗口中使用Ctrl+c来中断。MATLAB这时可能正疲于应付，响应会有些滞后。

2. figure命令新建一个绘图窗口
figure 可以打开一个空的绘图窗口，接下的绘图命令可以将图画在它里面，而不会覆盖以前的绘图窗口。当有多个figure窗口时，在命令窗口中执行如Plot等命令将覆盖当前figure窗口中的对象。所谓的当前figure窗口，也就是最后一次查看的窗口(可以用命令gcf得到)。
figure(N)，将编号为N的窗口置为当前figure，出现在所有窗口的最前面，如果该窗口不存在，则新建一个编号为N的空白figure。

3. 注释掉一段程序：%{、%}。
经典方法是用 if 0，但缺点是不够直观，注释掉的内容仍然保持代码的颜色。现在可以用 %和{ 的组合。输入%{后，后面的代码都会变绿，在注释结束的地方再加上 %} 。
也可以选中要注释内容，在右键菜单中选择Comment (Uncomment去掉注释标记)，或使用快捷键Ctrl+R。
将光标放在需要注释的行中，按Ctrl+R，将该行变为注释。取消注释也是一样的，快捷键为Ctrl+T。

4. doc 命令名，打开命令的帮助文档
与help命令不同，帮助文档MATLAB Help中对命令的描述更详细，往往还有一些例子，更便于理解。

5. clc 清屏
清除命令窗口中的所有输入和输出信息，不影响命令的历史记录。

6. clear 和clear all
clear 变量名：可以清除workspace中的无用的变量，尤其是一些特别大的矩阵，不用时及时清理，可以减少内存占用。
clear all 清除所有的变量，使workspace一无所有，当重新开始一次算法验证时，最好执行一次，让workspace中的变量一目了然。

7. close all
如果开了多个绘图窗口，用这个命令将它们一起关掉。

8. 上下光标键↑↓
在命令窗口中，上下光标键可以将历史记录中的命令复制到输入位置，便于快速重新执行。如果输入命令的前几个字母如 [row, col] =，再使用光标键，则只会选择以这些字母开始的命令。

9. Tab补全
对名称记得不太清楚的命令，可以输入开头的几个字母，然后按Tab键，当只有一个以这些字母开头的命令时，将自动补全命令名，否则显示一个命令名列表，方便从中选择。当然，只在命令窗口中有效。

10. cell模式
在一个长长的脚本m文件中，可能需要对其中的一段反复修改，查看执行效果，这时，cell模式就非常有用了。cell模式相当于将其中的代码拷贝到命令窗口中运行。两个%后接一个空格(%% )开始一个cell。将输入光标放到一个cell中时，背景将变为浅黄色，Ctrl+Enter执行cell中的代码。
执行cell中代码时不需要保存m文件，该m文件可以不在路径列表中。
cell模式中，断点不起作用，当然，调用的子程序中的断点还是正常的。

11. 获取文件列表，批处理
MATLAB内置了一些文件操作命令，例如cd(切换工作目录)，dir(同ls，显示目录内文件列表)等。dir命令可以返回目录中的文件和文件夹列表，存在一个结构体数组中。如果需要对一些数据文件进行批处理，而文件名又没有一定的规律，你可能需要借助于这个命令。

Inkscape tutorial advanced：高级教程

2008-02-13T18:18:00.001+08:00

Inkscape tutorial advanced：高级教程

本文翻译自矢量绘图软件Inkscape帮助文件中的高级教程，图片来自于http://www.inkscape.org/ 。由于原文档的格式是SVG，可以直接在帮助文件中操作绘图对象，所以文中的叙述方式可能有些特别。
翻译：liuxqsmile (liuxqsmile@gmail.com)

高级

作者：bulia byak, buliabyak@users.sf.net ； josh andler, scislac@users.sf.net

本教程包括：复制/粘贴、节点编辑、手绘和Bezier曲线、路径操作、布尔操作、偏移、简化、以及文本工具。

通过Ctrl+arrows, 滚轮, 或者中键拖动将绘图页面向下卷动。绘图对象的创建、选择、变换等基本操作，请参考帮助Help > 教程Tutorials中的基础教程。

粘贴操作

当用Ctrl+C复制对象或Ctrl+X剪切对象后，通常的粘贴Paste命令(Ctrl+V)将复制的对象粘贴到鼠标光标处，如果光标在绘图窗口外，则粘贴到文档窗口的中心。实际上，剪贴板中的对象仍然记着它的原始位置，你可以用原位粘贴Paste in Place将它粘回原始位置(Ctrl+Alt+V)。

另一个粘贴命令，粘贴样式Paste Style(Shift+Ctrl+V)，将复制对象的样式应用到所选对象。样式包括：填充、轮廓、以及字体设置，但不包括形状、大小、以及与该形状相关的参数，如星形的角数等。

命令粘贴大小Paste Size，将复制对象的大小应用到所选对象上。该命令包括：粘贴大小、宽度、高度，以及分别粘贴大小、宽度、高度。

粘贴大小Paste Size将全部选择的总大小缩放到剪贴板中对象的总大小。粘贴宽度Paste Width/粘贴高度Paste Height则仅影响水平和竖直方向上的尺寸。这些命令依据复制对象的长宽比是否锁定（选择工具控制栏，W和H的中间），如果复制对象的长宽比锁定，目标对象的另外一个方向上的尺寸将根据该比例自动缩放；否则，另一个方向的尺寸将不改变。带有“分别Separately”的相应命令也是类似的，不同之处在于将每个选择对象都分别缩放以适应复制的对象。

注意，Inkscape使用自己的内部剪贴板，除了在文本工具中复制/粘贴文本外，不使用操作系统的剪贴板。

手绘和规则路径

创建任意形状的最简单的方法是使用铅笔(手绘)工具(F6):

对于更规则一些的形状，可以用钢笔(Bezier)工具(Shift+F6):

在钢笔工具中，每次单击创建一个没有曲线控制柄的尖锐点，所以，一系列的单击产生一串线段，点击然后拖动产生一个光滑的Bezier节点，两边各有一个共线的控制柄。拖动一个控制柄时，按住Shift可以保持另一个不动。同样，Ctrl限制当前线段或Bezier控制柄的方向为15度的整数倍。Enter结束绘制，Esc取消。如果只取消上一段，使用Backspace。

在手绘和bezier工具模式下，选中路径的两端都会显示一个方形的锚点anchors，在这些锚点上可以继续绘图，从而延长路径，或使其封闭（从一个锚点画到另一个锚点），而不产生新的路径。

编辑路径

形状工具创建的是形状，而钢笔和铅笔工具创建的是路径。路径由直线和Bezier曲线构成，像其他对象一样，路径也可以设置任意类型的填充和轮廓属性。但与形状不同的是，修改路径时可以随意调整节点和(直线或曲线)段，而不是预先设置好的控制柄。切换到节点工具(F2)，然后选择下面的路径：

你会看到上面有一些灰色的方形节点。通过点击、Shift+点击、和拖出弹性选框，来选择这些节点。也可以单击路径中的一段来选择相邻的节点。选中的节点将高亮显示，并显示节点控制柄：一个或两个与该节点相连的小圆圈。!在当前子路径范围内反选节点(子路径上至少选择一个节点)。Alt+!在整个路径范围内反选节点。

路径的编辑通过拖动节点、节点控制柄、或路径段来进行。请在上面的路径上练习一下。Ctrl仍然有限制移动和旋转的作用。光标arrow 键, Tab, [, ], <, >的作用与选择工具中一样，但作用于节点，而不是整个绘图对象。在路径上的任意位置双击或Ctrl+Alt+Click可以添加节点。

选中节点后按Del ，或者Ctrl+Alt+单击一个节点，可以将其删去。删除节点时，将尽可能保持路径的形状，如果要保持相邻的节点不变（形状将发生变化），用Ctrl+Del删除。另外，Shift+D可以再制选中的节点，Shift+B可以将路径在选中的节点处打开，如果选中的是路径的两个终点，则可以将其结合在一起。

Shift+C可以使节点尖锐Cusp ，它的两个控制柄独立，可以各自调整角度；Shift+S可以使节点平滑smooth ，两个控制柄共线；Shift+Y可以让节点对称symmetric ，两个控制柄共线并且等长。改变节点的类型时，将鼠标悬停在一个控制柄上，可以保持该控制柄不变，仅另一个控制柄相应地改变。

通过Ctrl+click控制柄，可以将节点的控制柄收回（到节点上），如果相邻两个节点的控制柄都被收回，它们中间将变为直线。在节点上Shift+drag可以将控制柄重新拉出。

子路径和结合

一个路径可以包含数个子路径subpath。每个子路径中的节点互相连接，子路径与子路径之间则是断开的。左下图，三个子路径组合为一个路径，右下图中则互相独立，各自为一个路径：

要注意的是，复合路径并不等同于群组，它是一个单独的对象。如果你选中左上的对象，然后切换到节点工具，将会看到，三个子路径上的节点都显现出来，而在右侧，每次只能选中一个路径进行节点编辑。

通过对几个路径进行结合Combine可以形成一个复合路径(Ctrl+K)，也可以将一个复合路径分解为几个独立的路径 (Shift+Ctrl+K)。在上图中练习一下。由于一个对象只能有一种填充和轮廓样式，结合后的复合路径将继承第一个对象（处于叠放次序的底层）的属性。

在合并有填充的路径时，如果路径之间有重叠区域，合并后，重叠部分的填充将消失：

这是创建内部有孔的形状的最简单的方法。路径工具的高级操作请参考下面的“布尔操作”。

转换为路径

任何的形状和文本都可以转为路径 (Shift+Ctrl+C)。这个步骤不改变对象的外观，但对象的原本所具有的编辑方式（例如矩形倒圆，改变文本内容等）都将不复存在，而变为用节点工具进行编辑。这里有两个星形，左边的一个是形状，而右边的一个已经转为路径，切换到节点工具模式，选择这两个对象，看看他们的区别：

而且，任何对象的轮廓stroke都可以转换为路径(“outline”)。下图中第一个是原始路径（无填充，黑色轮廓），第二个是执行轮廓转为路径Stroke to Path后（黑色填充，无轮廓）：

布尔操作

路径Path菜单中命令可以将多个路径以布尔操作boolean operations的方式结合到一起：

布尔操作对应的快捷键也与相应的运算相适应(合并union对应加号，相减difference对应减号，等)。命令相减Difference和排除Exclusion 只针对两个路径，其它操作可以应用于任意数量的对象。操作后的对象总是保留参与操作的底层对象的样式。

排除Exclusion 与结合Combine 操作有些类似，只不过，排除Exclusion 在原始对象相交的地方添加节点。分割Division 和剪切路径Cut Path命令的区别在于前者用顶层路径将底层路径完全剪切，而后者只剪切轮廓，填充则完全完全删除(适用于将不用填充的轮廓分为数段)。

嵌入与扩展

Inscape不仅可以通过缩放，也可以通过偏移offsetting来扩展和收缩形状，即将路径上的点沿法线方向移动。相应的命令为：嵌入Inset (Ctrl+() 和扩展Outset (Ctrl+))。下图中给出了原始路径(红色)以及通过嵌入和扩展产生的新路径：

正常情况下，嵌入Inset 和扩展Outset命令生成的对象是路径（如果原始对象不是路径，将先转为路径）。通常，更方便的命令是动态偏移Dynamic Offset (Ctrl+J)，通过一个拖动控制柄（同形状的控制柄类似）来控制偏移量。选中下面的对象，切换到节点工具，拖动控制柄到一个合适的位置：

这种动态偏移对象dynamic offset object会记录原始位置，多次调整偏移时不会产生退化(degrade)。如果不需要再调整，可以将偏移对象转为路径。

也许，更有效的是关联偏移linked offset，与动态偏移类似，但原始对象仍然保留，并且可以编辑。一个原始对象可以有多个关联偏移。下图中，原始对象是红色的，其中一个关联偏移轮廓是黑色的，没有填充，另一个有黑色填充，但没有轮廓。

选择红色的对象，编辑其节点，观察关联偏移对象的变化。选择关联对象，拖动控制柄，调节偏移量。你会注意到，移动和改变原始对象影响到关联偏移对象，而偏移对象的移动和变换是独立的，同时保持和源对象的链接关系。

简化

简化Simplify (Ctrl+L)命令在尽量保持形状的情况下减少路径上的节点。铅笔工具创建的对象，节点数目往往过多，需要这个工具来简化。下图中，左侧的形状是通过手绘工具创建的，右侧是简化后的。原始对象有28个节点，简化后只有17个(节点工具编辑时更容易一些)，而且更平滑。

简化的程度（称为阈值threshold）取决于选区的大小。所以，如果选择路径的同时也选择了较大对象，简化的程度将更大。并且，简化的速度将加快。也就是所，如果连着按几次Ctrl+L（间隔不超过0.5秒），每次简化的阈值将递增。（如果等一会再执行，阈值又会还原原始大小。）通过这种方法可以比较精确地控制简化的程度。

除了简化手绘对象，简化命令还可以产生许多创造性的效果。显得和尖锐和呆板的对象经过简化经常后产生更柔和的效果。锐角变得平滑，引入更自然的变形效果，显得更有趣，更美观。下面是一个的对象经过简化后的效果：

创建文本

Inkscape可以创建复杂的文本。也可以很方便地绘制简短的文字对象，例如标题，标识，标语，流程图等中的文字。本节介绍Inkscape中文本工具的基本功能。

切换到文本工具(F8)，在页面上的任意位置点击，然后输入文字。打开文本和字体对话框Text and Font dialog(Shift+Ctrl+T)，可以修改文字的字体，样式，大小和对齐方式。这个对话框里也有一个文字输入框，可以修改选中的文本的内容。在这个对话框里输入文本可能比在画布上更方便（而且支持拼写检查）。

像其他的工具一样，文本工具模式下可以选择其自身类型的对象——文本对象——点击选择，将输入人光标放到文本中的任意位置（比如这个段落）。

文本编辑中常用的一个操作是调整文字间距和行间距，Inkscape中同样有对应的键盘操作方式。当编辑文本时，Alt+< 和Alt+>改变当前行的字间距letter spacing，该行的长度在当前缩放级别上每次改变一个像素（选择工具中，同样用这些键实现像素级别的缩放）。通常，如果字体比默认的大，字间距紧凑一些看起来更协调。例如：

紧凑一些的作为标题看起来更好一些，但仍然不是很完美：字间距并不一致，例如，“a"和"t"的间隔比"t"和"i"的间距大。在一些质量比较差的字体中，（尤其是字体比较大的情况下）这种不均衡的紧排更明显；但是，不管任何字体，总会存在这种文本组合，需要手工调整松紧。

在Inkscape中调整起来是很方便的，将光标放到需要调整的两个字符的中间，Alt+arrows键移动光标右侧的文字。与上面相同的文字，手动调整字符间距后：

除了Alt+Left 和Alt+Right将文字左右移动，Alt+Up 和 Alt+Down也可以将文字上下移动：

当然也可以将文字转为路径(Shift+Ctrl+C)，像路径一样移动。但是让其保持为文字无疑是更好的选择，仍然可以编辑，并且改变字体时不会丢失间距，同时文件的体积也小一些。保留文本的一个缺点是，当你将该SVG文件拿到别的计算机上打开时，这个机子上必须安装有相应的字体。

与字间距类似，在多行文本中也可以调整行间距line spacing。在本教程的任意段落中，Ctrl+Alt+< 和Ctrl+Alt+>来增大和缩小行间距，每次调整，整个文本的高度在当前缩放级别上改变一个像素。与选择工具类似，配合Shift键，行间距和字间距的调整量扩大十倍。

XML编辑器

Inkscape中的终极工具是XML编辑器(Shift+Ctrl+X)，可以实时显示整个文档的XML树形图。修改绘图时，你可以注意一下XML树形图中的变化。也可以在XML编辑器中修改文本、元素或者节点属性，然后在画图上查看效果。这是一个非常形象化的学习SVG格式的交互式工具。并且可以实现一些通常的编辑工具无法完成的功能。

小结

这个教程只展示了Inkscape功能的一小部分，我们希望你能喜欢。欢迎探索它的功能，展示你的灵感。更多信息，最新版本，以及寻求用户社区的帮助，请登录www.inkscape.org。

翻译：柳小勤Liu Xiaoqin, 2008.

16000米

2008-02-12T17:03:00.000+08:00

16000米

2008年2月6日，农历2007年腊月三十
这两天从运动场旁边经过时，看到跑道上的积雪已经给工人除去了，堆在跑道的两边。想想这个无所事事的假期，而我的长跑极限还停留在05年12月的10000米，这个除夕，还不如做点有意义的事，来一次两万米长跑吧。

五点多吃了学校的年夜饭。七点半回宿舍换了衣服，运动袜脚掌部分已经被磨出洞了，李宁的慢跑鞋鞋帮上的护垫海绵也已经破了，这次功成之后，他们也许可以安息了。
到操场上是八点，活动一下，8：05开始，跑了半圈才发现，跑道上的雪只除了一半，另半段上只有一条道上没有雪，但上面也已经结冰，走上去都要小心翼翼，更不用说跑了。犹豫了一下，既然来了，那就在这半段来回跑吧，跑上100次，不也就是20000米吗。
8:15 正式开始，从东南方向的弯道的开始，沿第一道跑到西北方向的直道的结束，然后掉过头，沿第六道跑回去，每一圈的实际长度超过400米。气温虽然在零度以下，但是没有风，也不感觉冷。跑道上只有我一个人，周围已经尽是连绵不绝的烟火了。
9:15 跑完25圈，10000米，我停下来，拿起手机看时间，感觉腿有些不灵敏了。跑步时膝盖下侧和脚腕上感觉有些沉重外，心肺并没有感到很大的压力。于是就挺乐观，再跑10000米应该没有太大的问题吧。
不过，形势很快急转直下，右脚掌用力时感觉有些疼，不知道是不是冻疮的问题。手感觉越来越冷，我戴的是半指手套，手指冰凉。速度明显下降了，转弯后速度很难再提上去。腿更沉重，心肺的压力虽然不像快跑时感觉那样剧烈，但更深重，似乎是一种从内到外的精疲力竭。只有大脑还比较清醒，体力透支过度会要人命的，比如马拉松。我觉得这种危险性似乎再增长。
9:55，15圈，我停下来，想稍微休息一下，然后接着跑完最后的10圈。腿疲劳得快不能响应神经的指挥了，慢走了两圈，重新再跑，情况仍然很严重，更糟糕的是，刚才还只是手冷，现在整个上半身都冷，只跑了一圈就放弃了。没有完成20000米，好歹也有16000，也算将个人记录向前推进了一大步。

回去后腿就开始疼了，右脚掌上磨出了一个水泡，还好没破，休息了三天也就全好了。我的体能这两年都没有提高，毕竟运动量都维持在原始的水平上：跑步还是隔一天一次，差不多3000米，近几个月自行车也骑的少了。以后的节奏也许该变一下了。

Inkscape tutorial: shape --Inkscape教程：形状工具

2008-02-02T12:26:00.001+08:00

Inkscape tutorial shape：Inkscape形状工具教程
Inkscape是一款开源的跨平台（Windows，Linux，）矢量绘图软件，类似于Coreldraw，Adobe Illustrator等。不仅自由免费，而且功能强大，官方网站及下载：http://www.inkscape.org/。
本文翻译自该软件附带的教程文档。

形状

作者：bulia byak, buliabyak@users.sf.net

本教程涵盖四种形状工具：矩形、椭圆、星形和螺旋线。我们将向你展示Inkscape在形状绘制上的能力，并通过实例演示这些工具的用法和用途。

通过Ctrl+Arrows, 滚轮, or 中键拖动将绘图页面向下卷动。绘图对象的创建、选择、变换等基本操作，请参考帮助Help > 教程Tutorials中的基础教程。

Inkscape有四种通用的形状工具，可以创建和编辑相应的形状。每种形状类型都具有共同的特征，可以通过拖动式的控制柄和一些数值参数来调节，使其具有不同的外观。

例如，在星形工具中，可以改变星的角数，边长，角度，圆角等，但仍然保持星形。形状不像简单路径那样“自由”，但通常更有趣，也更有用。形状总是可以转为路径(Ctrl+Shift+C)，但路径不能转为形状。

形状工具包括：矩形, 椭圆, 星形, 和螺旋线。首先，我们看一下形状工具的基本工作模式，然后再分别详细介绍。

一般性操作

选择相应的形状工具，在画布上拖动drag鼠标可以创建新的形状。新形状创建后(自动被选中)，上面有白色标记状的控制柄，可以通过拖动这些控制柄改变形状的外观。

在任意一种形状工具模式下，或在节点工具模式下(F2)，选择任意类型的形状，都可以使形状的控制柄显现出来。当鼠标悬停在控制柄上时，底部的状态栏将会告诉你拖动或点击将产生的编辑效果。

每一种形状工具在绘图区顶部都会显示一个水平的工具控制栏。通常包括几个数字输入框和恢复默认参数的按钮。如果当前选择的形状工具类型与当前对象的形状类型一致，在控制栏中输入数字也可以改变被选对象的形状。

工具控制栏中参数的修改将应用到新创建的形状上。例如，如果将星形的顶点数改变，新的星形的顶点数将与改变后的数值一致。另外，选择形状后，也会将其具有的特征参数显示在工具控制栏，从而影响以后创建的新形状。

在某一种形状工具模式下，可以通过在形状上单击使其被选中。Ctrl+click (选择群组内对象) 和Alt+Click (选择下层对象)也跟一般的选择模式一样起作用。

矩形

矩形是设计和绘图中最常用、最简单的形状。Inkscape尽力使矩形的创建和编辑简单快捷。

F4或者点击工具列中相应的图标，可以切换到矩形工具模式。请在这个蓝色矩形旁再画一个矩形：

保持在矩形工具模式下，通过单击，交替选择两个矩形。

矩形绘制技巧：

配合Ctrl键，可以绘制正方形或边长为整数比率(2：1，3：1等)的矩形。

刚绘制的矩形(新绘的矩形总是自动被选中)在三个角上显示3个控制柄，但实际上有四个控制柄，如果矩形没有圆角，右上角的两个是重合在一起的。这两个是倒圆控制柄；另外两个(左上和右下)是缩放控制柄。

首先看倒圆控制柄，抓住一个向下拖动，四个角都成圆形，并且第二个倒圆控制柄也显露出来(在角上的原始位置)。现在的圆角是等半径的，如果要使其成为不等半径的圆角，就需要向左移动另一个控制柄。

下图中前两个矩形具有等半径圆角，后两个则是椭圆形圆角：

在矩形工具模式下，点击选择这些矩形，看看他们的倒圆控制柄的区别。

有时，在整个绘图中，我们希望圆角的半径和形状保持不变，而不是随着矩形的大小而变化（例如流程图中，不同大小的矩形具有相同的圆角）。这在 Inkscape中可以轻松实现。切换到选择工具，在工具控制栏的尾部，有四个开关按钮，左数第二个按钮(显示两个同心圆角)，可以控制缩放矩形时圆角是否也同步缩放。

例如，在关闭“缩放圆角”情况下，下图中原始的红色矩形被复制和缩放了数次，各自具有不同的长宽比例：

注意所有矩形的圆角的大小和形状都是相同的，在右上角的位置可以精确地重合在一起。所有的蓝色虚线矩形都是从原始的红色矩形缩放得到，没有对倒圆控制柄做任何修改。

作为对照，下图中是在打开“缩放圆角”时的情形：

矩形的圆角各不相同，在右上角没有一个细微的地方是重合的(可以放大了看)。在将原始矩形转为路径后(Ctrl+Shift+C)，缩放时产生的也是这样的效果。

矩形的倒圆控制柄的操作技巧如下：

按住Ctrl拖动，使两个控制柄的半径相同(等半径圆角)。
Ctrl+click使另一个控制柄的半径与当前相同，不需要拖动。
Shift+click删除圆角。

你可能也已经注意到，矩形工具的控制栏中有水平(Rx)和垂直(Ry)圆角半径，可以精确控制其尺寸。顾名思义，不倒圆 Not rounded按钮可以删除选中矩形中的倒角。

这些操作可以同时对多个矩形起作用。例如，要改变当前图层的多个矩形，只需要Ctrl+A (全部选择)，然后在控制栏中设置需要的参数，而同时被选中的其他类型的形状则不会受到影响。

下面我们来看一下缩放控制柄。你可能会想，既然可以通过选择工具缩放，要这个缩放控制柄不是多余了吗？

选择工具缩放存在的问题是，它总是沿着页面的方向进行水平和竖直缩放，而矩形缩放控制柄是沿边的方向进行缩放，即使矩形经过选择或倾斜。请试一下用选择工具和(矩形工具模式下的)缩放控制柄对下面的矩形进行缩放的区别：

由于有两个缩放控制柄，可以沿任意方向缩放矩形，甚至可以让它沿着一条边的方向移动。这种缩放模式保持圆角半径不变。

缩放控制柄的技巧：

Ctrl拖动沿某一边或对角线方向缩放。换句话所，Ctrl拖动保持宽度、高度或宽高比不变(即使其本身的坐标系经过缩放和倾斜)。

下面是同样的矩形，请用Ctrl拖动，使之沿虚线方向缩放：

通过倾斜和旋转矩形，然后复制和缩放，可以很容易地生成三维效果：

一些其它的矩形构图，其中也用了圆角和渐变填充：

椭圆

椭圆工具(F5) 可以创建椭圆和圆，然后也可以将其分为分割或圆弧段。绘制技巧与矩形工具相同：

Ctrl绘制圆或整数比率(2:1, 3:1,...)椭圆。

让我们详细看一下椭圆的控制柄。首先选择这个椭圆：

与矩形工具相似，可以看到三个控制柄，实际上是四个。最右侧的是两个重合在一起，可以通过它“打开”椭圆。拖动最右侧的控制柄，然后再拖动新出现的一个，可以得到各种扇形的分割和弧线：

要得到分割segment(半径封闭)，拖动时鼠标要保持在椭圆的外面，相反，如果鼠标在椭圆里面，则得到弧线arc。上图中，左侧共有四个分割，右侧有三个弧线。注意，弧线是不封闭的，仅弧线上有轮廓线，两端没有连接在一起。如果去掉填充色，只保留轮廓，看得就更清楚一些：

注意左侧看起来像扇子的很窄的分割，它是通过Ctrl对控制柄进行角度捕捉实现的。下面是弧线和分割控制柄的一些操作技巧：

配合Ctrl键拖动时，控制柄每次改变15度角。

Shift+click使其重新成为完整的椭圆。

捕捉角度的大小可以在Inkscape Preferences(其中的Steps子项)中调节。

另外两个控制柄使椭圆绕中心缩放。操作技巧同矩形的控制柄是相似的：

Ctrl拖动可以使椭圆成为圆(长短轴半径相等)。

Ctrl+单击使椭圆成为圆，不需要鼠标拖动。

与矩形的缩放控制柄类似，椭圆的缩放控制柄也是沿其自身的坐标轴方向调整其长度和宽度。也就是说，旋转或倾斜后仍然可以沿其本身的长短轴方向进行拉伸和挤压。试一下通过缩放控制柄对下面的椭圆进行缩放：

星形

星形是Inkscape中最复杂，最吸引人的形状。如果你想用Inkscape让你的朋友们眼前一亮(Wow！)，只需要给他们展示一下星形工具。它简直是乐趣无穷。

星形工具可以创建两种类似但看起来差别很大的形状：星形和多边形。星形有两个控制柄，决定星形角的长度和形状；多边形只有一个控制柄，拖动时控制旋转和缩放：

在星形工具的控制栏中，第一项是一个选择框，可以将星形转为对应的多边形或者相反。第二项是一个数字框，设置星形或多边形的顶点个数，这个参数只能在这里设置，数值范围为3-1024。如果你的电脑性能不高，最好不要设置过大的数值(超过200)。

当绘制一个新的星形或多边形时，

Ctrl使得其旋转角度为15度的整数倍。

一般来说，星形看起来更吸引人(实际上多边形更实用)。星形的两个控制柄的差异略有不同。第一个初始位置在角的顶端，控制角的长短，但相对于星形的中心旋转时，另一个控制柄也跟着旋转。也就是说，这个控制柄不能让星形的角倾斜。

另外一个控制柄，初始位置在两个角中间的底端，可以沿半径或圆周方向任意移动，而不影响角顶端的控制柄(实际上，它的位置可以超过现有角的顶点而产生新的角)。这个控制柄可以使星形的角倾斜，获得晶体、雪花、曼荼罗以及porcupines等各种形状。

如果你要的只是一个规则的星形，不需要这些装饰性的工艺，你可以抑制倾斜控制柄的倾斜功能：

Ctrl拖动时可以让星的角保持半径方向(不产生倾斜)。

Ctrl+单击取消倾斜，不用拖动鼠标。

作为在画布上拖动控制柄的一种补充方法，控制栏中的轮辐比Spoke ratio可以指定两种控制柄相对于中心距离的比值。

Inkscape兜里还有两个小绝技。从几何上来数，多边形是由直边和锐角组成，但现实中，会出现各种各样的边和圆角，这在Inkscape中也可以实现。倒圆星形和多边形与矩形略有不同，并没有专门的控制柄。但是，

Shift+拖动切向拖动控制柄来倒圆星形或多边形。

Shift+点击控制柄来取消圆角。

“切向”表示垂直于指向中心的半径的方向。如果沿逆时针方向拖动，得到正圆角，顺时针方向拖动，则得到负圆角。（参看下面图中负圆角的例子。）

这里比较一下由矩形工具产生的圆角正方形和由星形工具产生的圆角正四边形：

可以看到，圆角矩形由边和圆角构成，而圆角多边形和星形则没有边，从角的顶点圆滑过渡到底端的点。Inkscape通过在各个顶点产生切线方向一致的Besizer曲线来平滑星形(正多边形)。（将其转为路径后，用节点工具可以看出来。）

控制栏中的圆角Rounded参数是切线相对于多边形/星形的边的长度的比例。为负值时改变切线的方向。其值介于0.2-0.4时是一个“正常”的圆角，其它值将会产生复杂、优美的图案。对于星形，圆角数值过大时曲线沿控制柄延伸很远。下面的例子给出了一些数值对应的图案：

如果希望星形的角上尖锐而底部平滑或相反，可以通过对星形做偏移offset (Ctrl+J)来实现：

Shift+drag拖动控制柄是是调节圆角的最精细的方法之一，可以处理得很好。

为了模仿真实的形状，Inkscape可以对星形和多边形做随机化randomize处理。轻微的随机使形状的规则性下降，更随意，更有趣；而更强的随机则可以产生无法预测的“疯狂”效果。具有圆角的星形随机化后仍然保持圆角。其中的一些技巧：

Alt+drag切向拖动控制柄，产生随机效果。

Alt+click取消随机效果。

当绘制或拖动控制柄随机化星形时，形状会“抖动”，因为控制柄的每个位置都对应一种唯一的随机化效果。不使用Alt，仅移动控制柄，则在当前的随机化水平上重新产生一种随机化效果，使用Alt，拖动控制柄时，保持当前的随机化效果，仅改变随机化的程度。下面的星形具有相同的参数，仅略微移动控制柄，具有不同的随机化效果(随机化程度0.1)。

下面是上图中中间的星形，但随机化程度介于-0.2到0.2：

Alt+drag拖动该行中间星形的控制柄，可以看到它可以变形成为两边的形状。

在你自己的工作中，你将会发现随机化星形具有的各种用途，但我个人尤其喜欢这种像变形虫和凸凹不平的星球表面一样充满幻想的图案：

螺旋线

Inkscape中的螺旋线是一种灵活的形状，尽管不像星形工具那样吸引人，有时也有很大的用途。螺旋线像星形一样，也从中心绘制，编辑也是相对于中心的，

Ctrl+drag保持角度为15的整数倍。

画好后，螺旋线有两个控制柄，分别位于内部的起点和终点。拖动时都可以改变螺旋的角度（增减圈数）。其它的操作技巧：

外部控制柄：

Shift+drag相对于中心缩放/旋转（不改变螺旋的角度）。

Alt+drag改变螺旋角时保持半径不变。

内部控制柄：

Alt+drag竖直拖动改变发散度，向中心汇聚/分散。

Alt+click重设螺旋发散度。

Shift+click将内部控制柄移动到螺旋线的中心。

螺旋线的发散度divergence用来衡量其螺旋的非线性程度。值为1时，为标准的螺旋线；小于1(Alt+drag 向上)，边缘密集，大于1(Alt+drag 向下)，中心密集：

螺旋线的最大圈数为1024。

椭圆工具不仅可以产生椭圆，也可以产生圆弧（具有常曲率），螺旋线则可以产生曲率平滑改变的曲线。与普通的Bezier曲线相比，圆弧和螺旋线更方便，可以随意改变长度（拖动控制柄），而不会影响到其形状。螺旋线默认没有填充，加上填充，取消轮廓线，也可以产生一些有趣的效果。

最有趣的是，轮廓线为点时产生的效果：向内旋转和（点或短线）等间隔分布可以组合成很优美的波纹效果：

小结

Inkscape形状工具的功能非常强大。有时间请注意多加练习，在实际的设计工作中他们不会让你失望的。而且，矢量绘图中形状工具通常比路径工具使用起来更快速，编辑起来更方便。如果你对这些工具有一些改进的建议，请与开发人员联系。

翻译：柳小勤，Liu Xiaoqin, 2008

AnvilStudio:免费的乐谱MIDI编辑打印软件

2008-01-23T19:46:00.000+08:00

AnvilStudio:免费的乐谱MIDI编辑打印软件

在电脑上，可以将乐谱与MIDI文件对应，MIDI文件对乐谱进行标记，调用音频设备产生声音。因此，我们可以简单地认为MIDI就是会发声的乐谱:-)

打谱软件在国内使用最多的似乎是Overture，被国人汉化（破解）后，据说有上千万的下载量。但它实际上是一个商业软件，免费的试用版仅允许15天的使用时间。国内也有人开发了一些类似的软件，但大多针对简谱。
在Linux上，可以选择开源的Rosegarden ，功能相当丰富，可以作为音序器和合成器，也提供了乐谱的编辑、打印功能。不过目前还没有对应的Windows版本。

由Willow Software开发的Anvil Studio 可以作为一个相当不错的MIDI生产器。它的免费版本对MIDI音轨的数目没有限制，可以通过这个软件将乐谱转为MIDI文件。

该软件的运行界面如上图，功能也很全面，乐谱和键盘对照，对初学者应该有很大的帮助。
这个软件的免费版不能打印乐谱，不过我们可以通过其它的免费软件来解决，如Notation Player ,可以将MIDI文件以乐谱的形式显示和打印。Finale Notepad 也可以实现这样的功能，并且也可以编辑和保存乐谱，不过它的编辑功能比较弱，并且软件体积比较大。

Science 2.0: Great New Tool, or Great Risk? 科学2.0：新工具？新风险？

2008-01-18T12:41:00.000+08:00

Science 2.0: Great New Tool, or Great Risk? 科学2.0：新工具？新风险？

科学美国人网站近日刊登了一篇名为《Science 2.0: Great New Tool, or Great Risk?》的文章，对互联网新技术如Wiki和Blog等在科学研究中发挥的作用进行了一些讨论。

科学工作者由于自身工作的需要往往（尤其是信息的获取）与网络都有密切的联系，互联网的发展对科学工作者的工作方式也难免不会产生影响。越来越多的科研人员，不仅是年青人，开始将自己的工作以Blog、Wiki和时下热门的Web2.0社会化网络等方式公开传播。尽管还没有形成普遍的潮流，但星星之火似乎已经表明，这种借助于新兴网络技术的“科学2.0”，比传统的工作方式更有开放性，效率也更高。
“科学并不仅仅是做实验，更重要的是人们对实验的交流。” 在新的工作方式中，研究人员不仅提前展示自己的工作成果，也会记录工作思路的产生、变化，实验中失误、尝试等等。这些不会表现在最终的正式成果中，但对于其他研究者也具有很重要的借鉴意义。
当然，也有很多人对这种完全透明的方式有疑虑，或许你的思路启发了别人，但那人又秘而不宣，抢先发表，断了作者的前途，更糟糕的是，可能会有小人完全剽窃作者的成果。

MIT的两名研究生在2005建立了一个名为OpenWetWare 的维基网站，最初只是用于学生自我展示。但很快，它在交流中的重要意义就展示出来，人们在这里讨论实验室工作的技巧，课程学习中的问题，即将完成学业的学生将他们的经验总结出来，供新来的学生参考。Wiki对实验室管理也有很大的好处，研究成员可以比较系统地整理自己的工作，便于同事之间的相互交流、了解。这种网络工作日志公开后，吸引来的同行可以通过评论等方式方便地与作者开展交流，提出意见，有利于效率的提高。
已经有来自世界各地的很多实验室、数十种课程和兴趣小组，以及数百个讨论组在OpenWetWare上安了家，3000多名注册用户维护着6000多个页面。2007年他们还得到美国国家自然科学基金的资助，以将其建成一个独立于MIT的社区模式的网站。

对于Idea可能被盗的风险，科学2.0的支持者认为，作者在Wiki上的每次更改，发布，都会留下时间标记，这些可以对成果的原创性作出证明，让盗窃者颜面扫地。当然，如果想要申请专利，对成果的公开就要非常小心了。

与Wiki相比，科研工作者介入Blog的进度显得非常缓慢。在2007年四月的北卡罗莱纳科学博客会议上North Carolina Science Blogging Conference，Willard认为，这是由于博客的方式与科研人员受到的训练有天然的矛盾。传统上，他们总是对结果反复验证，发表论文时字斟句酌，因为论文是几乎永久性的，而在博客上发表内容要做到时效性，就难免保证其准确、完整。

有一些人在blog上公布自己的研究工作的细节、进展，以期获得更广泛的交流。也有一些人在上面讨论科研环境中的问题，如科研机构的立足点，对实验室管理的看法，对大boss的期望等等。当然，大多数科研工作者视这些为雷区，讨论时多用假名，以免冒犯了某方面的权威，耽误了自己的前程。也有一些人觉得，在这上面浪费时间还不如酝酿下一篇论文来得实在，毕竟，经过同行评价的正式出版物才是工作和晋升的重要筹码。在博客上，他们无法提升自己的科研资历。

成果的冠名权问题无疑是阻挠科学2.0发展的重要问题。毕竟，论文被看成是科研贡献的一个重要尺度。但大多数人也承认，论文并不是唯一的方面，分享灵感，开展讲座，领导社群，也是很重要的方式。或许，随着信息交流的发展，我们可以找到更全面的度量方式。

当然这也意味着学院内部传统的变革，对科学2.0的倡导者来说，最重要的变化将是从对个人科研名誉的重视，转向真正以科学进展为第一目标。

Inkscape tutorial basic: Inkscape教程：基础

2008-01-13T13:51:00.000+08:00

Inkscape tutorial basic: Inkscape教程：基础

本文翻译自矢量绘图软件Inkscape帮助文件中的基础教程。由于原文档的格式是SVG，可以直接在帮助文件中操作绘图对象，所以请注意文中的叙述方式。
翻译：liuxqsmile (liuxqsmile@gmail.com)

基础

作者：bulia byak, buliabyak@users.sf.net

本教程描述了Inkscape的基本操作方法。本文档以Inkscape的通用文件格式保存，你可以用Inkscape进行查看、复制、编辑、保存等操作。

该教程的主要内容包括：画布浏览、文档管理、形状工具基础，图形选取、变形，群组、填充与轮廓、对齐和叠放。对更复杂的操作，请在帮助菜单中选择其它相关教程。

平移画布

平移画布(卷屏)的方法有很多种。使用Ctrl+arrow键可以用键盘卷屏。(你可以尝试这些按键来卷动本文档。) 也可以通过鼠标中键来拖动画布，或者使用屏幕边缘的滚动条(使用Ctrl+B来显示或隐藏滚动条)。鼠标滚轮wheel可以上下卷动画布，按住Shift键，配合滚轮则可以水平卷动。

放大与缩小

最简单的缩放操作是通过-和+（或=）键。也可以通过Ctrl+middle click 或 Ctrl+right click来放大，Shift+middle click 或 Shift+right click 来缩小画布。也可以用Ctrl键配合鼠标滚轮来缩放。或者在窗口右下角的缩放输入框中输入一个准确的百分比数值。在工具栏中也有缩放按钮，可以缩放到用户选定的区域（对象）。

Inkscape还会记录当前工作会话中使用的缩放历史，按`键回到上一次的缩放比例，Shift+`键来恢复撤销的缩放比例。

Inkscape工具列

Inkscape中的绘图和修改工具集中在左侧的竖直工具列中。在窗口的上方，菜单下面是命令栏(Commands bar)，提供了通用的一些控制命令，下面的工具控制栏(Tool Controls bar)则跟具体的绘图工具有关。窗口底部的状态栏(status bar)则实时显示一些操作提示和信息。

很多操作都有对应的快捷键,在帮助菜单中选择鼠标与快捷键(Help > Keys and Mouse)获取详细的说明。

创建和管理文档

选择菜单文件File > 新建New，或使用快捷键Ctrl+N新建文档。选择菜单文件File > 打开Open (Ctrl+O)打开已有文档。选择菜单文件File > 保存Save (Ctrl+S)来保存文件。或者选择菜单文件File > 另存为Save As (Shift+Ctrl+S)将当前文件以不同的文件名保存。（Inkscape可能有时不够稳定，切记经常保存！）

Inkscape使用SVG(Scalable Vector Graphics可缩放矢量图形)文件格式。 SVG是一种被各种绘图软件广泛支持的开放文件标准。SVG文件是基于XML的，可以用任何文本和XML编辑器来编辑(Inkscape不属于这种文本编辑器)。除SVG外，Inkscape也可以导入和导出其它一些文件格式(EPS,PNG等)。

Inkscape为每一个文档打开一个独立的窗口。你可以用操作系统中的窗口管理器来在各个窗口间切换(例如Alt+Tab键)，也可以使用Inkscape中内置的快捷键Ctrl+Tab在文档间循环切换。(现在可以新建一个文档，尝试在本文档和新文档间切换。)

创建形状

下面我们开始创建一些很漂亮的图形！在工具列中选择矩形工具(Rectangle)(快捷键F4),在(本文档或新文档的)绘图区中点击、拖动：

如上图，默认的矩形内部是蓝色、半透明，轮廓stroke是黑色的。下面将要介绍怎样改变这些特征。在其它工具中，你可以选择创建椭圆、星形、和螺旋线：

这些工具统称为形状工具shape tools。新创建的每一个形状上都有一个或更多四边形的控制器(handles)；试一下拖动这些控制器会产生什么样的效果。在工具控制栏中也可以对形状进行修改。工具控制栏只对当前选中的形状有效(显示出四边形控制器的)，但同时也会成为当前形状工具的缺省参数，影响下次创建的图形。

按键Ctrl+Z可以撤销(undo)上一次操作。(如果你又改变注意了，可以用Shift+Ctrl+Z来恢复(redo)撤销的操作。)

移动、缩放和旋转

Inkscape中最常用的工具是拾取器(Selector)，位于工具列的顶端(箭头形状)，对应快捷键F1 或者空格(Space)。现在你可以选择当前画布上的任何对象。请点击下面的矩形。

可以看到，选择对象的周围出现八个带箭头的控制器。下面你可以：

通过拖动来移动对象。(按下Ctrl来进行水平或竖直移动。)

通过拖动任意的控制器来缩放。(按下Ctrl以保持原始的宽度-高度比例。)

再次在矩形上单击，控制器会发生变化，现在你可以：

拖动对象角落上的控制器来旋转。(按下Ctrl以保持旋转的角度为15度的整数倍。)

拖动中间的控制器来扭曲(倾斜)对象。按下Ctrl以保持扭曲的角度为15度的整数倍。)

在选择状态，也可以在工具控制栏(画布的上方)的输入框中输入数字，精确地控制对象的位置坐标(x,y)和尺寸(宽度W,高度H)。

通过键盘变换

Inkscape区别于大多数其它矢量绘图软件的一个特征是键盘操作的便捷性。几乎所有的命令都可以通过键盘实现，变换操作也不例外。

你可以用键盘来对编辑对象进行移动(arrow 光标键)，缩放(< 和 > 键)，以及旋转([ 和 ] keys)。缺省情况下，每次移动和缩放2px，按下Shift键时每次移动则为20px。而Ctrl+> 和 Ctrl+< 对应的缩放比例则分别为200%和50%。缺省每次旋转 15 度，通过Ctrl键，每次可以旋转90度。

可能更有用的是像素级别的变换(pixel-size transformations), 实现的方法是，在上面的快捷键基础上配合Alt键。例如Alt+arrows可以在当前的页面视图层次上每次移动一个像素(这里是指一个屏幕像素的距离，而不是SVG中的与视图缩放级别无关的长度单位px)。这意味着，如果放大视图，Alt+arrow移动一个像素的绝对距离将缩短。这样，通过缩放视图，就可以任意控制对象的定位精度。

与此类似，Alt+> 和 Alt+< 将选择对象每次缩放一个像素， Alt+[ 和 Alt+] 旋转对象时，距离旋转中心最远的位置每次移动一个像素。

注意：在Linux操作系统中，这些组合键可能在窗口管理器中被指定了其它的用途，执行上述操作时可能不能获得预期的结果。解决的方法当然是相应地修改窗口管理器的配置。

多选

通过Shift+click，可以连续选择多个绘图对象，或者，用鼠标左键拖出一个框来选中框内所有对象，这个也称为弹性区选(rubberband selection)。(从空白处开始拖动时将创建弹性选区，如果在拖动之前先按下Shift，则总是创建弹性选区。) 请尝试选择下面的三个形状：

你可以使用弹性选区选择下面两个椭圆，但不包括矩形：

被选择的对象上会出现一个选择标识(selection cue)，默认情况下是一个虚线矩形框，它可以标识出哪些对象被选中，哪些没有选中。例如，同时选中两个椭圆和矩形时，如果没有矩形标识框，椭圆的选中与否就难以判断。

在已经选择的对象上Shift+click可以取消选择。选中上面的三个对象，然后用Shift+click取消对两个椭圆的选择，只选中矩形。

按Esc取消所有选择，Ctrl+A选择当前图层上的所有对象(如果没有定义图层，则等价于选中文档中的所有绘图对象)。

群组

若干个绘图对象可以组合为一个群组group

。群组可以像普通绘图对象一样进行移动或变换。下图中，左边的三个图形是互相独立的，而右边的三个图形是组合在一起的。试着拖动这个群组看看。

选择一个或多个对象后，按Ctrl+G可以将它们组合在一起。选中一个或多个群组后，按Ctrl+U可以解散组合。群组也可以再次组合，并且群组的嵌套层数没有限制。不过Ctrl+U只能打开最顶层的群组，对于嵌套群组需要多次Ctrl+U才能完全打散组合。

实际上你可以直接修改群组内的对象而不用取消组合。使用Ctrl+click就可以单独选中群组内的一个对象，进行编辑；使用Shift+Ctrl+click则可以选中群组内或群组外的多个对象。不需要解散群组，请试着对上图右面群组中的形状进行单独的移动、变换。然后再选中群组，可以看到这种组合关系仍然存在。

填充与轮廓

Inkscape中的许多功能都借助于对话框的形式。为绘图添加一些色彩的最简单的方法是打开视图View菜单中的调色板Swatches对话框(快捷键Shift+Ctrl+W)，然后为对象选择一种(填充)颜色。

更强大的工具是对象Object

菜单中的填充与轮廓Fill and Stroke对话框。选中下面的形状，然后打开填充与轮廓。

这个对话框中有三个标签面板：填充Fill、轮廓色彩(Stroke paint)和轮廓样式(Stroke style)。填充属性可以修改对象的内部fill 。下面的按钮可以设置填充的类型，包括不填充(图标X)，单色flat color填充，以及渐变(gradients,线性或圆周)填充。对于上面的椭圆，单色填充的按钮是激活的。

这些按钮的下面，是色彩拾取器color pickers，有四种不同的方式：RGB, CMYK, HSL,色盘Wheel。可能最方便是通过色盘来选择，旋转其中的三角形来选择色调，在三角形内可以拾取不同的明暗度。四中拾取方式中都包含一个滑动条来设置对象的透明度(opacity)，即alpha值。

选择不同的对象时，色彩拾取器总是自动更新，对应到当前对象。(选择多个对象时，将显示色彩的平均值。) 在下面的例子上做一下练习，你也可以自己创建这些图形：

在轮廓色彩Stroke paint标签中，可以删除轮廓线stroke，也可以任意为其指定颜色和透明度：

最后一个标签面板，轮廓样式(Stroke style)中，可以设置轮廓的宽度以及其它参数：

最后，除了单色填充之外，可以选择梯度(gradients)模式来填充图形内部和轮廓：

当从单色填充切换到梯度填充时，颜色仍然是前面单色填充时的颜色，不同的是透明度从不透明渐变到完全透明。选择工具列中的渐变工具(Gradient tool,Ctrl+F1)，对象上将会显示出(用线连接在一起的)渐变控制器，拖动渐变控制器(gradient handles)，可以改变色彩梯度的方向和范围。选中某个控制器时(该控制器呈现蓝色)，可以在填充和轮廓中为该控制器单独设置色彩，实现从一种颜色到另一种颜色的渐变。

还有一种改变对象色彩的简便方法是使用滴管工具(Dropper tool,F7)。选择对象后，再选择该工具，然后可以在绘图中单击click任意拾取色彩，这种色彩将自动指定给被选择对象的填充属性(使用Shift+click

再制、对齐和分布

一个常用的操作是生成对象的一个副本，即再制duplicating(Ctrl+D)。新生成的副本与原对象重合(垂直于纸面方向)，并且已经被选中。可以用鼠标或光标键把它移走。想练习一下？请将下面一行用这个黑方块填满：

移动后，新的方块的位置难免些不够整齐，这时对齐对话框(Align dialog,Ctrl+Shift+A)就派上用场了。选择所有的方块(Shift+click或者拖出一个弹性选区)，打开这个对话框，选择“中心水平对齐(Center on horizontal axis)”，再选择“水平等间距分布(Make horizontal gaps between objects equal)”，这些方块的位置就很整齐匀称了。下面是一些利用对齐和分布工具生成的图案：

叠放次序Z-order

z-order指的是绘图中对象的叠放次序，例如，某个对象在最上层，盖住了其它的对象。对象(Object)菜单中的两个命令，置顶(Raise to Top,对应Home键)与置底(Lower to Botton,End键)，将使所选对象置于当前图层叠放次序(Z方向)的顶部或底部。另外两个命令上升(Raise，PgUp键)与下降(Lower，PgDn键)，将使被选择对象上升或下将一个位次，例如，可以将当前对象移动到它上面一个图形的上面。(如果所选对象与其它对象都不重叠，上移和下移分别等同于置顶和置底。)

可以在下面的图形上练习改变叠放次序，让最左边的椭圆位于最上层，而最右边的椭圆位于最下层：

选择叠放的对象时，一个很方便的快捷键是Tab。如果没有选择任何对象，按Tab将会选择最底层的对象；有对象被选中时，将选择其上的对象。Shift+Tab的选择方向则相反，从最顶层开始，往底层逐次选择。默认的叠放次序与图形创建的次序是一样的，所以没有选择对象时，Shift+Tab总是选择刚创建的图形。在上面的叠放椭圆中可以练习一下Tab和Shift+Tab的选择。

选择下面的对象并移动

如果一个对象完全被另一个对象盖住了，该怎么选择呢？如果上面的图形是透明的，你虽然可以看到下面的对象，但点击时选中的却是上面的图形。

这就是Alt+click要干的活。首先在上面的图形上Alt+click，这将选中它，然后在相同的位置上再次Alt+click，这次将选择该位置处，顶层图形下面的对象。对于多层叠放，多次Alt+click实现从顶层到底层的循环选择。

[如果你在Linux系统中工作，Alt+click可能不会像前面描述的那样工作，反而可能会移动整个Inkscape窗口。这是因为窗口管理器为Alt+click指定了其它用途。你需要找到窗口管理器的相应配置，把其中的这个快捷键关掉，或选择其它的组合键。]

选中了被盖住的图形，你又可以做什么呢？可以用光标键移动，可以用鼠标拖动控制器。但是，如果拖动整个对象，则会重新选择顶部的图形(这是点击-拖动的工作模式，总是选中顶部的对象然后拖动)。要让Inkscape拖动当前选择的对象，而不是顶部的对象，需要借助于Alt+drag，这将拖动当前选择的对象，而不论你的鼠标在哪里。

请用Alt+click和Alt+drag选择并拖动绿色透明矩形下的棕色形状：

总结

好了，最后做一下小结。对于Inkscape，这仅仅是开始，但靠这几招，你已经可以做一些简单但不失实用的图形了。更高级的复杂操作，请参照Help > Tutorials中的其它教程。

LiuXiaoqin(liuxqsmile@gmail.com),2008

我的电脑中不显示U盘的问题

2008-01-11T18:49:00.000+08:00

Windows XP USB disk我的电脑中不显示U盘的问题

U盘或移动硬盘插到电脑上以后，并不总是跳出烦人的对话框，有时候也很安静，只在右下角的系统托盘上显示一个“安全移除硬件”的图标就完了。你打开我的电脑，可移动磁盘呢？怎么没有呢？刷新！没有，关掉再重新打开我的电脑，还没有！难道要重启吗？
你点击“安全移除硬件”的图标，打出来的白条上明明写着：安全删除USB存储设备(H:,I:)云云...
恭喜你，有了这个白条意味着U盘已经被操作系统识别并加载了，U盘应该还没有挂掉。怎么方便地操作其中的文件呢，重启不仅麻烦而且也不一定顶事。
你可以：
打开一个应用程序，如记事本Notepad、Word之类。在“文件”菜单中选择“打开”，这将弹出一个打开文件对话框，点击左侧的“我的电脑”，将列出所有磁盘，其中有U盘吗？如果有，你可以接着看下去，如果没有，不好意思浪费你时间了，你可以狠狠地关闭这个页面以泄愤。
在“文件类型”中选择“所有”，这样U盘的文件都显示出来了，现在可以在右键菜单中操作你需要的文件了，复制、粘贴、以及打开等。如果你前面用的是记事本，想要打开U盘中的图片，一定要在图片的右键菜单中选择打开，千万不能直接双击。

也可以试一下在任务管理器中结束explorer.exe进程，再重新打开。
还遇到过一次怪事情，接到电脑上的U盘，第一个不能显示，第二个却可以显示，不妨再找一个U盘来垫底。

这个问题的根源我不知道，我的系统曾经出现过几次这样的情况，或许只是Windows调皮了一下而已。

Ruler with slot带凹槽的塑料直尺

2008-01-07T20:15:00.000+08:00

Ruler with slot带凹槽的塑料直尺

中午去超市买东西，顺便买了一把塑料直尺。回来后想起来以前经常遇到的一个问题，将尺子放到桌面上就很难拿起来。拿新买的尺子试了一下，果然是着这样。尺子本身倒是明晃晃，很漂亮的。就是边上都处理得太光滑了，重新拿起来得费半天的劲，只用一只手的话根本就拿不起来。
如果在没有刻度的一边开一条小凹槽，是不是就能解决这个问题呢？

Baldur's Gate重玩博德之门

2008-01-06T19:50:00.001+08:00

Baldur's Gate重玩博德之门

博德之门和博德之门2算是我唯一付过钱的商业软件了，只玩过一遍就扔在箱底想来实在有些可惜。前不久偶然看到一个关于博德三部曲的Google Group,可以用博德2的引擎将博德1、剑湾传奇、安姆的阴影、巴尔的王座放在一起玩，但是占用的硬盘空间实在是大了一些。于是就只安装了博德1和剑湾传奇。
博德1的设定跟2里有一些差异，比如升级。升级时获得的生命点数也是扔骰子，升级前可以先保存一下，一次不满意可以重新读进度再扔，而在博德2里似乎就不行了，经验值满时，扔出来的骰子就确定了。这个对战士类的角色影响应该是很大的，在剑湾传奇中战士可以升到8级，每次d10，不算其它修正，可以获得将近80的生命值。
如果要相信“命运天注定”的话，8d10，八次平均分布的和，根据中心极限定理，概率分布为均值为40的高斯分布，方差是...，（等我有时间算一下;-)），总之不会离40太远，这个差距也太大了。
以前用圣骑士，现在改用法师了，只能缩在队伍的后面，扔扔石子。与戴娜黑的塑能系专精法师相比还弱了一些，法术专精每等级可以额外记忆一个法术，看来是个不错的选择。我的法师是弄杂了，居然将智慧也弄到17，敏捷狂差，我可没想过转职到牧师。

这个游戏还是太粘人，我已经玩了不下四十个小时，才刚刚进到博德之门。战斗已经不像印象中的那样让人抓狂了，只有石化蜥蜴让我加深了一下印象，嗯，还有贾希拉在卡立德中招时的叫声。
如果将天之痕的剧情画面用博德的引擎、战斗方式来表现.....，只能随便想想了。

update：
春节初一到初四又集中了四天时间终于再次通关，整个游戏大概玩了90个小时。杜拉格之塔花了太多的时间，以至于对最后的决战有些不耐烦，砍倒了沙洛佛克后就草草收场了。由于不用太顾忌队员的阵亡，最后一战的难度实际上并不是太高。
比较郁闷的是，戴娜黑升到顶级后，居然找不到法术卷轴可学，最高级别的两个法术栏位一直空着。

Tips：2007.12工作点滴

2007-12-27T20:15:00.000+08:00

1.MATLAB官方博客：http://blogs.mathworks.com/
其中的Inside the MATLAB Desktop 主要介绍工作环境和绘图中的新功能和技巧。推荐订阅。
最近的一篇"Focusing on Zooming" 一文中介绍了使用命令来控制figure的缩放。通过linkaxes命令，可以将一个figure中的各个axes指定相同的坐标范围，实现各个绘图同步放大、缩小。在缩放一个绘图时，其它也会响应地同步缩放。可以同时关联x、y轴，也可以只关联一个坐标轴。
Giving your code some privacy一节介绍了MATLAB2007b中引入的代码折叠功能，对编写大型的程序应该很有帮助。

2.MATLAB Help中的Interactive Plotting一文提供了一种在figure中用鼠标绘图的方法，主要使用ginput函数。

3.判断点在离散的封闭曲线（多边形）内的方法，基本原理是，多边形内一点发出的任意射线与多边形仅有奇数个交点；

1). 求出封闭曲线的外接矩形的范围。即曲线中所有顶点的最大、最小坐标范围。
2). 判断点是否在该矩形内，否，则结束；是，则执行3)
3). 判断由该点发出的水平（或垂直）射线与多边形的交点个数。该射线与多边形的边有交点的判断（以水平正向为例）：
边的两个端点的纵坐标不能同时大于或小于该点纵坐标，否则，该射线与该条边没有交点。
如果两个端点的横坐标同时小于该点的横坐标，则没有交点，如果同时大于该点的横坐标，则有交点。
如果一大一小，则需要求出该射线与边的交点的横坐标（纵坐标已知），如果交点横坐标大于点的横坐标，则有交点。（如果曲线上的点很密集，判断精度要求也不高，这一步也可以忽略）。

用上面的方法，在MATLAB中手动进行锋电位特征的聚类：

4. Word中图片使用嵌入型时，位置自动移动到页面的边缘：可能是样式引起，将图片的样式改为正文或可解决。实际上
可能是图片所在段落的行距太小造成的，可以设为单倍行距试试。

sudo的简单配置

2007-12-21T18:53:00.000+08:00

sudo的简单配置

sudo的配置文件是/etc/sudoers，它有专门的编辑工具visudo，用su切换到root用户，然后执行这个命令。
对于VI，对我这样的初学者，只要记住几个快捷键就可以了：光标移动，j上一行，k下一行，h向左，l向右。i切换到编辑状态，然后再配合End、Delete、和BackSpace键，完成这个位置的编辑，Esc退出编辑状态；移动光标到下一个位置，再切换到编辑状态，再输入，直到修改完成。然后输入冒号“：”，最下面一行会改为接受命令状态，输入w，保存，q退出。
执行visudo以后，将光标下移，可以看到一行：

root ALL=(ALL) ALL

这个命令定义了root用户的sudo权限，我们要添加的配置也具有相应的命令格式即

user host=[(Runas)] [NOPASSWD] commands
被授权用户主机=[目标用户][不需密码] 命令

即“被授权用户”可以在“主机”上以“目标用户”的权限执行“命令”。方括号内是两个可选参数，目标用户缺省为root。命令必须用绝对路径，之间用","分开。例如

liuxq localhost=/sbin/poweroff

表示用户liuxq可以在本机上以root的权限执行"sudo /sbin/poweroff"命令，而不需要root密码（需要liuxq的用户密码）。
如果加上NOPASSWD，则表示不需要输入密码：

liuxq localhost=NOPASSWD: /sbin/poweroff

基本的定义方法就是这样，如果需要为不同的用户、不同的主机配置不同的命令权限，这样逐条配置显然太繁琐，所以这个配置文件中引入了别名Alias的功能，在配置文件的开头部分，有四个别名类型：

User_Alias; Host_Alias; Runas_Alias; Cmnd_Alias;

分别对应该命令形式中的四个部分。别名实际上定义了相同类型的对象的一个集合。例如

User_Alias LIU=liuxq,lxq

定义了一个LIU集合，其中有成员liuxq和lxq，用LIU代替前面命令中的liuxq，则用户liuxq和lxq都可以执行sudo /sbin/poweroff命令。
其它别名的定义也是类似的。这里就不多啰嗦了，推荐看看Linux 系统中的超级权限的控制，讲得比较详细。

Expectation-maximization：高斯混合模型的期望最大化算法

2007-12-13T21:45:00.000+08:00

高斯混合模型的期望最大化算法

本文大部分内容翻译自Answers.com的Expectation-maximization algorithm 词条。

对于一个不完整的样本集，假设y为观测（已知）数据，z为丢失数据，z和y构成了完整的数据。z可以为丢失的观测量，也可以是隐含的变量，如果这些变量已知的话将会简化问题的求解。例如，在混合模型中，如果产生样本点的各个成分的分布为已知，则可以大大简化似然函数的计算。

假设p为完整数据的联合概率密度函数，该函数具有参数θ：，它也可以看成是完整数据的似然函数（即样本集(y,z)具有θ分布函数的相似性），一个关于θ的函数。进一步，根据贝叶斯公式和全概率公式，丢失数据对于观测数据的条件分布为：

这个式子仅要求计算观测数据对于丢失数据的条件概率（似然函数）和丢失数据的分布概率。（分母仅为归一化需要）

EM算法通过迭代方法，构造更好的 θ₁,θ₂,... ，逐步优化初始的 θ₀， θ的递推优化公式通过下式得到：

Q(θ)是对数似然函数ln的期望。换句话说，我们不知道所有数据的值，所以我们不知道似然函数的准确值，但是对于给定的观测数据y，我们可以得到未知的z的概率分布，对于每一个丢失的观测（z中的一个元素），都有一个关于θ的似然值，这样我们就可以计算似然函数的数学期望。这个期望显然受到θ的值的影响。

Q(θ)通过下式得到：

或者，更普遍地，可以写成条件期望的形式

应该清楚，这里关于的条件期望中用到的z的条件分布（前一个式子中的第一项）是在固定θ_n的情况下取得的。用对数似然代替直接的似然函数，可以在保证极值位置不变的情况下简化计算。也就是说，θ_{n + 1}是通过在θ_n 的基础上，最大化所有数据关于已知数据的条件期望E而取得的。

当样本集合的最大似然估计比较容易获得时，EM算法更为有用。如果似然估计具有给定的形式，则M步就非常简单。一个经典的例子是用于有限的高斯混合模型的最大似然估计，如果混合分布已知，则可以很容易得到各个成分的分布。

高斯混合模型

假设n个高斯分布，从中产生m个样本,每一个 y_j 所从属的分布用 z_j 表示，z_j 的取值范围为1到n。对于任意的y，其来自第i个高斯分布的概率为

我们的任务是估计未知的参数，即每一个高斯分布的均值、方差和先验概率。

首先我们必须设定初始值。

对于无监督的聚类，我们首先必须为每一个样本指定一个类别，这一步可以通过其他的聚类方法实现，如K-means方法，求出各个类别的中心和每一个样本的类别。然后求出各个类别中样本的协方差阵，可以用每个类别中样本的个数来表示该类别的权重（先验概率）。

E-step

使用上一次M步中得到的参数θ，计算未知的z（即每一个分类）对于观测数据的条件分布。即，在当前参数θ_n 下，对每一个样本，计算各个类别的后验概率。下式中第三个等号后面的分子的第一项，为每一个样本在类别i中的分布概率，分子中的第二项为每一个类别的先验概率。
分母为归一化，即对每一个样本，其在各个类别中的分布概率的和为1.

for ti = 1:Ncluster
PyatCi = mvnpdf(featurespace,Mu(ti,:),Sigma(:,:,ti));
E(:,ti) = PyatCi*Weight(ti);
end
Esum = sum(E,2);
for ti = 1:Ncluster
E(:,ti) = E(:,ti)./Esum;
end

M-step

为了最大化联合分布的对数似然函数的期望

上式中最后一项的联合分布可以写成条件分布和边缘分布的乘积的形式，得到

上式需要满足先验概率归一化的要求：

为了求得Q(θ)的极大值的位置，引入拉格朗日算子，然后带入概率密度函数，

为了求得新的 θ_{t + 1},我们需要找出上式取得极值时对应的θ，即其一阶导数。

下面推导满足上式的均值，

for ti = 1:Ncluster
muti = E(:,ti)'*featurespace;
muti = muti/sum(E(:,ti));
Mu(ti,:) = muti;
end

产生新的协方差矩阵

新的各个分类的权重

进行归一化处理

Inserting λ into our estimate:

将新的参数，均值、方差、分类的先验概率，带入E步，重新计算，直到样本集合对各个分类的似然函数不再有明显的变化为止。

下面的程序提供了一种快速的似然函数计算方法，摘自Patrick Tsui的EM_GM样例程序：
http://www.mathworks.com/matlabcentral/fileexchange/loadFile.do?objectId=8636

function L = Likelihood(X,k,W,M,V)
% Compute L based on K. V. Mardia, "Multivariate Analysis", Academic Press, 1979, PP. 96-97
% to enchance computational speed
% a subfunction from Patrick P. C. Tsui's EM_GM code.
[n,d] = size(X);
U = mean(X)';
S = cov(X);
L = 0;
for i=1:k,
iV = inv(V(:,:,i));
L = L + W(i)*(-0.5*n*log(det(2*pi*V(:,:,i))) ...
-0.5*(n-1)*(trace(iV*S)+(U-M(:,i))'*iV*(U-M(:,i))));
end

MAXScript:利用3DSMAX脚本产生原子运动动画

2007-12-10T19:13:00.000+08:00

利用3DSMAX脚本产生原子运动动画

很久，很久以前，曾经做过一些关于分子动力学的工作，最后无疾而终，几乎什么都没有留下。分子动力学（Molecular Dynamic Simulation）计算微观原子间的相互作用，产生其运动轨迹，所有原子的运动形态最终表现为物体宏观的运动和变形等。

如果将这些原子的运动以动画的形式渲染出来，无疑是非常形象的。我曾经用3DSMAX的脚本MAXScript编写了一段程序来完成这个功能，在这里贴出来，供有需要的参考。
以下图中的模型为例，上面是Au探针，压入下面的Ni平板。原子的初始坐标值（x,y,z）分别放在文本文件posau0.txt和posni0.txt;
计算过程中产生一系列的原子坐标文件，这些文件以编号区别，如posau1.txt～posau100.txt;

在脚本文件中，首先读入原子的坐标信息，建立模型，每个原子对应一个球模型；
然后，在时间轴的不同时间坐标上，读入并指定每个原子的新坐标。
脚本的执行我也记不清了。可以在网上搜一下。
完整的程序如下（“--”的注释是把程序粘贴过来后补上去的，不知道有没有什么问题。我是在3DSMAX 5上运行的。）：

--指定数组变量；
posarray = #()
postemp = #()
spAu = #()
spNi = #()
lay = #()
--字符串变量
strti = " "
--Au模型中原子个数，坐标数据量
atomau = 768
dataau = atomau*3
--Ni模型中原子个数，坐标数据量
atomni = 616
datani = atomni*3
--位置坐标文件的个数，不含初始坐标文件
framenum = 14
--原子半径
Rau = 1.44
Rni = 1.3
--指定材质，Ni平板每一层具有不同的颜色
--textMatAu = standardMaterial diffuse:[125, 125, 255]
textMatNi = standardMaterial diffuse:[92, 216, 92]
textMat1 = standardMaterial diffuse:[126, 173, 126]
textMat2 = standardMaterial diffuse:[233, 114, 67]
textMat3 = standardMaterial diffuse:[236, 224, 92]
textMat4 = standardMaterial diffuse:[92, 92, 216]
textMat5 = standardMaterial diffuse:[60, 151, 159]
textMat6 = standardMaterial diffuse:[179, 121, 216]
--读入Au的初始坐标
filename = "F:\data\posau0.txt"
in_name = (filename)
in_file = openFile in_name

if in_file != undefined then
(
for v=1 to dataau do
(
vert = readValue in_file
posarray[v]=vert
)
close in_file
--根据原子坐标建立每个原子的球模型
for v=1 to atomau do
(
ti=(v-1)*3+1
px = posarray[ti]
ti=ti+1
py = posarray[ti]
pz = posarray[ti+1]

spAu[v] = sphere radius:Rau position:[px, py, pz] segs:12 smooth:true
--判断原子所在的层，指定不同的材质
lay = (v-1)/288 as integer
lay = lay+1
if lay == 1 then
(spAu[v].mat = textMat1 )
if lay ==2 then
(spAu[v].mat = textMat2 )
if lay ==3 then
(spAu[v].mat = textMat3 )
if lay ==4 then
(spAu[v].mat = textMat4 )
if lay ==5 then
(spAu[v].mat = textMat5 )
if lay ==6 then
(spAu[v].mat = textMat6 )

)

)

--读入Ni的原子位置，并建立模型
filename = "F:data\pos\ni0.txt"
in_name = (filename)
in_file = openFile in_name
if in_file != undefined then
(
for v=1 to datani do
(
vert = readValue in_file
posarray[v]=vert
)
close in_file

for v=1 to atomni do
(
ti=(v-1)*3+1
px = posarray[ti]
ti=ti+1
py = posarray[ti]
pz = posarray[ti+1]

spNi[v] = sphere radius:Rni position:[px, py, pz] segs:12 smooth:true
spNi[v].mat = textmatNi
)

)

--开始设置坐标-时间动画
animate on (
at time 0
for tj in 1 to framenum do
( ti=tj*10
--指定新坐标对应的时间
at time ti (
strti = tj as string
filename = "F:\data\posau"+strti+".txt"
in_name = (filename)
in_file = openFile in_name
if in_file != undefined then
(
for v=1 to dataau do
(
vert = readValue in_file
posarray[v]=vert
)
close in_file

for v=1 to atomau do
(
ti=(v-1)*3+1
px = posarray[ti]
ti=ti+1
py = posarray[ti]
pz = posarray[ti+1]
--设置Au原子模型的新位置
spAu[v].position = [px, py, pz]
)
)

filename = "F:\data\posni"+strti+".txt"
in_name = (filename)
in_file = openFile in_name
if in_file != undefined then
(
for v=1 to datani do
(
vert = readValue in_file
posarray[v]=vert
)
close in_file

for v=1 to atomni do
(
ti=(v-1)*3+1
px = posarray[ti]
ti=ti+1
py = posarray[ti]
pz = posarray[ti+1]

spNi[v].position = [px, py, pz]

)

)
)

)
)

--max time play

再贴一个效果图(模拟计算结果本身可能不正确):

Task Manager任务管理器中一些功能的简单介绍

2007-11-30T20:21:00.000+08:00

任务管理器中一些功能的简单介绍

表格中的内容来自：http://www.microsoft.com/technet/prodtechnol/windows2000serv/reskit/core/fneb_mon_oyjs.mspx?mfr=true

在任务管理器的查看->选择列中，可以更改进程面板中显示的性能检测项目。

进程统计项目	Description描述
Base Priority 基本优先级	进程的基本优先级决定了这些进程在CPU的处理队列中的排列顺序。基本优先级由进程本身控制，而不是操作系统指定。操作系统设置、改变进程中线程的动态优先级。可以在任务管理器中设置进程的优先级，更详细的内容参考“处理器活动分析 ”。
CPU Time CPU时间	从进程开始运行起，所占用的CPU总时间，以秒计。
CPU Usage CPU占用	上次跟新到本次更新的时间段内，该进程占用的CPU时间的百分比。
GDI Objects 图形设备接口对象数目	进程使用的图形设备接口（GDI）控件的数目。GDI控件是指，从针对图形输出设备的应用程序接口GDI库中调用的控件。
Handle Count 句柄数	进程的对象列表中对象句柄的数目。
I/O Other I/O其它	进程产生的非读非写的输入输出操作的数目，包括，文件、网络和设备I/O。这种操作的一个例子是行使控制功能。不统计指向控制台的操作（控制台输入）。
I/O Other Bytes I/O其它字节数	其它I/O操作传输的字节数。
I/O Read Bytes I/O读取字节数	进程在输入输出操作中读入的字节数。包括文件，网络和设备I/O。指向终端的操作不在统计之列。
I/O Reads I/O读取	进程在输入输出操作中读操作的数目。包括文件，网络和设备I/O。指向终端的操作不在统计之列。
I/O Write Bytes I/O写入字节数	进程在输入输出操作中写出的字节数。包括文件，网络和设备I/O。指向终端的操作不在统计之列。
I/O Writes I/O写入	进程在输入输出操作中的写出操作数。包括文件，网络和设备I/O。指向终端的操作不在统计之列。
Image Name 映像名称	进程名。
Memory Usage 内存使用	进程占用的主内存数，以千字节KB计。
Memory Usage Delta 内存使用增量	两次统计更新之间，内存占用量的变化，以KB计。
Nonpaged Pool 非页面缓冲池	该进程占用的操作系统内存数，以KB计。操作系统内存中从未分页到硬盘的部分。分页是将应用程序内存中的不常用的部分转移到其它存储媒介上，通常是转移到硬盘上。
Page Faults 页面错误	应用程序请求的数据由于不在内存而需要从硬盘上重新读入的次数。从应用程序启动开始计数。
Page Faults Delta 页面错误增量	两次更新间的页面错误增加量。
Paged Pool 页面池	该进程占用的操作系统内存数，以KB计。页面池是可以分页到硬盘的虚拟内存。页面池包括所有的用户内存和系统内存的一部分。
Peak Memory Usage 内存使用峰值	从启动开始，进程占用的物理内存峰值。
PID (Process Identifier) 进程标识	应用程序启动时分配的数字标识。
Thread Count 线程数	应用程序中运行的线程数。
USER Objects 用户对象	进程占用的用户对象数。用户对象是指窗口管理器中的对象，包括窗口、菜单、光标、图标、钩子hooks、加速器、显示器、键盘布局，以及其它内部对象。
Virtual Memory Size 虚拟内存量	为该进程分配的虚拟内存（或地址空间）量。

在任务管理器的性能面板中，通过“查看”菜单中选择“显示核心时间”，可以查看应用程序调用操作系统服务所占用的CPU时间，剩余的时间，即用户模式，是运行应用程序本身的代码所占用的时间。

下表简单地描述了性能面板中一些统计项目的内容。

Task Manager counts	Description
CPU Usage CPU使用
MEM Usage 内存使用	虚拟内存的使用量（虚拟内存（Virtual Memory）是Windows管理所有可用内存的方式）
Total Handles 总句柄数
Total Threads 总线程数
Total Processes 总进程数
Physical Memory: Total 物理内存：总量
Physical Memory: Available 物理内存：可用量	包括清零的，释放的和待分配内存。
Physical Memory: File Cache 物理内存：文件缓存	Amount of physical memory, in kilobytes, released to the file cache on demand. 文件缓存占用的物理内存。
Commit Charge: Total 认可用量：总数	所有进程使用的虚拟内存量
Commit Charge: Limit 认可用量：限制	不扩大页面文件的情况下，可供分配的虚拟内存总量。
Commit Charge: Peak 认可用量：峰值	虚拟内存的峰值用量，如果虚拟内存增大，可能会超过认可用量的限制值。
Kernel Memory: Total 核心内存：总量	页面和非页面内存的总量
Kernel Memory: Paged 核心内存：页面	操作系统占用的页面池容量。
Kernel Memory: Nonpaged 核心内存：非页面	操作系统占用的非页面池容量。

Task Manager counts

Description

CPU Usage

CPU使用

MEM Usage

内存使用

虚拟内存的使用量（虚拟内存（Virtual Memory）是Windows管理所有可用内存的方式）

Total Handles

总句柄数

Total Threads

总线程数

Total Processes

总进程数

Physical Memory: Total

物理内存：总量

Physical Memory: Available

物理内存：可用量

包括清零的，释放的和待分配内存。

Physical Memory: File Cache

物理内存：文件缓存

Amount of physical memory, in kilobytes, released to the file cache on demand.

文件缓存占用的物理内存。

Commit Charge: Total

认可用量：总数

所有进程使用的虚拟内存量

Commit Charge: Limit

认可用量：限制

不扩大页面文件的情况下，可供分配的虚拟内存总量。

Commit Charge: Peak

认可用量：峰值

虚拟内存的峰值用量，如果虚拟内存增大，可能会超过认可用量的限制值。

Kernel Memory: Total

核心内存：总量

页面和非页面内存的总量

Kernel Memory: Paged

核心内存：页面

操作系统占用的页面池容量。

Kernel Memory: Nonpaged

核心内存：非页面

操作系统占用的非页面池容量。

关于内存这部分把我弄晕了，对我来说又没有什么太重要的意义，算了，有心人可以看下面牛人写的：

虚拟内存，物理内存，页面文件，还有任务管理器

[quote]发信人: Hakkk (蓝猫淘气3k问·荷兰篇), 信区: NewSoftware
标题: 虚拟内存，物理内存，页面文件，还有任务管理器
发信站: 水木社区 (Tue Nov 20 21:56:20 2007), 站内

虚拟内存（Virtual Memory）是Windows管理所有可用内存的方式。
对于32位Windows系统，每个进程所用到的虚拟内存地址从0到2^32-1，总容量4GB，
其中2GB是与操作系统以及其他所有进程所共享，
另外2GB分派给进程独占（这就是常说的32位Windows中一个进程最多能用2G内存的由来）。

4GB虚拟地址空间（Virtual Address Space，VAS）中，2GB的进程独占VAS是进程隔离的，
换句话说，每个进程都可以从RAM或者硬盘上映射到属于自己的2GB VAS。
所以虽然32位Windows每个进程最多只能获得4GB可用虚拟内存，
但是所有进程总和可以使用总量超过4GB的虚拟内存。
可惜XP只支持4GB RAM，也就是说超出部分一定得靠页面文件补足；
而比如2003企业版通过PAE支持32GB RAM，可以减少页面文件的用量。
（但单个进程最多仍然只有4GB虚拟内存地址空间，而不是32GB）

所以，一个简单的计算示例：
三个进程最多能用掉多少VAS？
2GB（共享）＋2GB（独占）×3＝8GB

这4GB的虚拟地址空间，按照4KB的大小进行分页（page），
然后以页为单位映射到实际存储单元中，包括：
·物理内存（RAM）
·页面文件（Page File，在Win9x中称为交换文件Swap File，即win386.swp）
·其他文件自身（比如一些长时间未活动的进程的exe文件自身）
可见，通常人们所说的虚拟内存实际上只是指其中的页面/交换文件而已，
这是对虚拟内存的一个错误的理解。
虚拟内存≠页面文件

RAM中除了保存最近读写的文件缓存（File Cache，相当于Win9x中的Vcache）
主要用来存储正在使用的程序代码和数据，
当RAM资源紧张，或者有程序码或数据长时间未使用时，
XP通常会将非活跃程序码所在的地址页映射回程序文件（exe、dll等），
将数据所在的地址页映射到页面文件（pagefile.sys）中并拷贝数据，
然后将它们本来占用的RAM空间释放。
这个过程称为页出（Page Out）。

当系统读取某个虚拟内存地址，而该地址所在的页不在RAM中时，
将产生一个页面错误（Page Fault）中断，
告诉系统从页面文件或者程序文件中取回包含该地址的虚拟内存页，
即将内容拷回到RAM并建立新的虚拟地址映射，并将页面文件中对应部分标记为未使用，
这个过程就是页入（Page In）。
页入成功的话就是一个Valid Page Fault，否则就是Invalid Page Fault。
前者非常普遍，（可以在任务管理器的进程页监视到）
而后者是由程序或硬件错误引起，
如果发生在进程上会导致非法操作，如果是系统自身则很可能蓝屏。

内存用量可以在任务管理器的性能页中看到。
其中物理内存的总量、可用数等是指安装的RAM容量和剩余RAM容量，
而内存使用或者提交更改的总量和限制是指虚拟内存的Commit Limit和Committed Bytes，
可以理解成系统可以使用的虚拟内存总量和当前使用量，
其中总量是由RAM大小＋页面文件大小决定的。

而在任务管理器的进程页中，“内存使用”和“虚拟内存大小”造成的误解很多，
而实际上，这两个值是对应性能监视器中该进程的Working Set和Private Bytes，

Working Set是指一个进程的4GB虚拟地址空间中被映射到RAM中的部分的大小，
通常是该进程的虚拟内存中的活跃部分。
表面看来这个表述和“进程占用的RAM大小”没有太大区别，
但至少有两种情况导致了例外的发生：
第一种是这部分虚拟内存如果是属于2G的系统共享虚拟内存，
那么它映射到的RAM地址可能和其他进程有重复，计算多个进程占用RAM总和就不正确，
使得所有进程的Working Set之和比实际占用RAM要大。
第二种是一个进程中可能有同一段RAM内容的多个引用，比如一个文件被loop多次，
此时虚拟空间中会有多段地址被映射到同一段RAM上，
从而造成该进程的Working Set比实际占用RAM要大。

Private Bytes是一个进程的2GB独占虚拟地址空间中用到的部分的大小，
无论这部分是在RAM中还是在页面文件中，甚至是在exe、dll等文件中。
所以，任务管理器进程页中的虚拟内存绝对不反映pagefile.sys的用量。

所以，任务管理器无论是进程页还是性能页都不能直接反映pagefile的使用情况，
诸如“关掉pagefile还能看到虚拟内存”之类的疑问，应该很清楚了。

下面是一个简单图示，
两个进程，各自拥有2G独占VAS，共同拥有2G共享VAS，
W代表映射到RAM中的部分（Working Set），P代表页出的部分。

0 …………………2GB独占内存……………… 2G ………2GB共享内存……… 4G
进程1 |--PPPP---------------WW-----PPPP--------|
\\ || \\
\\ || \\
\\ || |||| |--PPPPPPPPPPPP--WWWWWWWW--|
exe2 exe1 dll RAMemory pagefile.sys |||||||||||| ||||||||
|||| /// // / || | 共享内核/dll等共享RAM空间
进程2 |-PPPP------PPP---WW--WW------PP----P----|

然后，任务管理器的进程页中的“内存”就是上图中的W，
所以所有进程的内存项加起来肯定比实际用的RAM大，因为有共享部分
此外还有重复引用，比如上图中进程2独占内存中第二和第三个W指向同一段RAM空间。

而任务管理器进程页中的“虚拟内存”就是上图中右半部分，“独占VAS”。
可见pagefile只是其中的一部分，而RAM，只要用于独占VAS，一样计入了该处“虚拟内存”
。

最后，实际的页面文件用量（pagefile.sys中实际使用部分），
可以在性能监视中看到，即Paging File下的% Usage和% Usage Peak
配合当前pagefile.sys的大小就可以计算出字节数，
这个值可以做为设定页面文件最小值的一个参考，
而不是用任务管理器进程页的所有进程的虚拟内存大小相加做为页面文件最小值的参考。

Sitting Too Close to the Computer Screen Can Make You Go Blind

2007-11-28T19:23:00.000+08:00

长时间注视屏幕容易引起视觉疲劳

Sitting Too Close to the Computer Screen Can Make You Go Blind

Eyestrain is a common—and occasionally debilitating—effect of staring at screens

By Molly Webster

原文：http://www.sciam.com/article.cfm?id=is-sitting-too-close-to-screen-making-you-blind&page=1

翻译：liuxqsmile

图片来源http://www.adcock.co.za/LifeStyle.aspx?Page=EyeComponents ；我在上面添加了对应的中文名称。

你晃晃脑袋，松弛一下紧张的脖子，闭上眼睛，揉一揉，然后继续在几英寸之外盯着电脑屏幕。完全不管你的的视线仍然一片朦胧，几难分辨出屏幕上的字符来。这时，你应该会有疑问：再继续下去，键盘上所敲击的每一个按键，YouTube上所观看的每一个视频，以及发送的每一个电子邮件......这些会损坏你的视力吗？

眼科医师，验光师以及其他的眼科专家注意到近视与近距离用眼有一些关系。这些近距离用眼包括，在距眼睛40cm左右读书、注视电脑屏幕等，尽管电脑屏幕一般距离在50cm左右。但是，只有一小部分人会因为近距离用眼而导致近视，不论他们是读书还是读屏幕。“我们还不清楚哪些人会受影响。”俄勒冈州太平洋大学视觉学院教授James Sheedy说。

事实是，近距离用眼不会在所有人中都导致近视，但是，并不是说紧坐在电脑前面就一点问题都没有。尽管对大多数人来说，不会造成永久性的后果，但近距离用电脑却会导致眼睛疲劳综合症。俄亥俄州立大学视觉学院教授Mark Bullimore告诉我们，眼睛疲劳综合症是指，眼睛疲惫不堪：烫痛、干燥和眼肌疼痛，以及所有继续用眼时伴随的不快感受。

理解眼睛为什么会疲劳，并且知道该如何去防止的最简单的方法是了解我们的眼睛能看清细小文字的工作原理。当我们“看”到某个对象的时候，实际上是从物体上反射的光穿过角膜（覆盖在眼球表面的一层透明的穹状结构）。角膜和晶状体（虹膜后面的圆形、透明、弹性结构）对光波进行聚焦，使其准确地到达视杆细胞和视锥细胞，这两种位于视网膜上的光感受器用来收集入射的光线信息。位于最内层的视网膜汇总光线信息，并且通过视神经传递到大脑，形成一幅图像。

Sheedy告诉我们，近距离凝视屏幕会使睫状肌紧张，睫状肌控制晶状体的形状，以准确聚焦，睫状肌必须始终保持张紧，这对这块小肌肉的要求太苛刻了。过紧的盯视也会阻止眨眼。眨眼的重要作用在于将眼泪覆到眼球的表面。如果不眨眼，角膜表面就会干燥，看东西模糊不清。正常情况下，眼睛每分钟眨20次，用电脑时可能会下降到7次，但是，专家们认为不会造成永久性的后果。

眼科教授Dennis Robertson说，如果屏幕周围有耀眼的强光，还可能会导致斜视。斜视可以减弱进入眼球的光线强度（一定程度上缓解眨眼不足），但是斜视让眼球周围的肌肉处于紧张状态，连续九个小时，让它筋疲力尽。

在眼睛放松后，眼睛疲劳通常会在数小时之后消。但只要这些症状一出现，就会严重影响到工作效率，更严重时让人感觉烦躁。这些症状都不是永久性的，不会造成视觉功能的损害。我们也可以通过改善工作环境来避免这些症状的出现。

使用非闪烁的电脑屏幕，改变屏幕的亮度、对比度、和文字的大小，屏幕的位置要略低于眼睛，显示器的上沿与眉毛平齐。如果眼睛的生理有问题，应该赶快配一副矫正眼镜。

关掉周围太亮的光源。用手在眼前遮住灯光，看你的脸部和肩膀是否有一定程度的放松，如果有，说明这些灯光应该做出调整。专家建议，看电视的时候，在比较舒服的位置上跟里面的主人公“休戚与共”。

但是，到目前为止，避免眼部疲劳的最好、最简单的方法仍然是经常性地休息、放松。即使你不渴，也可以到饮水机前走一趟，一定要把你的安乐椅放到电视两英尺（约0.6m）以外，最后，别忘了眨眨眼睛。

译者注：调整Windows桌面字体一般在显示属性->外观中选择大字体，但这样字体反而太大，并且只有桌面中的部分文字受其影响，对应用程序一般没有什么效果。实际上，在设置->高级面板中修改DPI，可以更灵活地控制字体的大小，大部分应用程序中的文字都会相应改变。

Quack medicine: Imagination Goes Wild(3) 伪医学：想象力的疯狂(3)

2007-11-22T19:36:00.000+08:00

Quack medicine: Imagination Goes Wild(3) 伪医学：想象力的疯狂(3)

英文原文出自www.collectmedicalantiques.com ，链接http://www.collectmedicalantiques.com/quack3.html
翻译：liuxqsmile
本文地址：http://liuxqsmile.blogspot.com/2007/11/quack-medicine-imagination-goes-wild3-3.html
（译注：图片有删减）
有很多种产生电的方式，但在这个发展迅速而缺乏规则的行业里，每一种都被逐利的奸商所利用（56）。在交流电应用之前，早期的干电池后来的蓄电池(57-61)，都被用来产生直流电，在诊所(62-64)和病人家里治病。电子腰带可以产生一个“治疗”电磁场，带有一个会发光的灯泡，让人以为它很有效果(65-37)。大型的静电发生器和摩擦生电仪通过摩擦玻璃来产生电荷，存储在莱顿瓶里，放电时就可以用来治疗一切叫得上名字的疾病。不过，这些仪器在早期也为X射线的研究提供电力（68-74）。一些临床医学家甚至在透热疗法的仪器上找到了微波的妙用(75-77)。
在二十世纪，电的应用更为广泛。制造商们也挖空心思，很多的大型设备出现在诊所和医院里，电疗椅(78-82)、电疗浴(83)、电子壁柜(84-86)及橱柜(87)等。

放射性

居里夫人发现的放射性为医学界创造了不可估量的商机。Revigator(88,89)是一个用有弱放射性的粘土制成的罐子，制造商宣称，患者喝了在罐子中过夜的水后可以“水到病除”。不幸的是，直到在那些坚持使用的人中间产生很多的口腔癌患者时，人们才对它产生了怀疑。尽管这个罐子的放射性很低，但氡（译注：参考世界卫生组织|氡与癌症）的含量很高，收藏者要注意把他们放在密闭的容器里。

56)宣传招贴
新奇药品和虚假医疗器械的生产商们通过广告招贴来宣传他们的产品。这是Scott博士的一幅宣传画，用来推销他诸多骗人的电子装置中的一部分。

57-61)家用直流医疗装置
使用电池的“电子组合”是在二十世纪早期很常见的虚假医疗器械，并且在江湖医术的历史中占有重要的位置。正是由于它的销售最终导致了对相似的医疗器械进行立法，对安全和功效进行了标准化。图中展示了诸多附件中的一部分。Williams电磁仪也使用直流电驱动。

62-64)交流电装置
交流电发展起来以后，出现了很多大型装置。图中精美小巧的麦金托什电源盒可以为前面的“电子组合”提供充足的电力，较大的汤普森电子柜则用于耳鼻喉科的诊治。

65-67)治疗腰带
插上电源，打开开关，这个带子会就产生一个磁场，被宣称具有医疗作用，并且很安全（现在看来是很有问题的）。随带子附送一个接有小灯泡的线圈，当线圈接近磁场时，灯泡会发出亮光，“看，不是很有效吗？"。

68-74)大型静电发生器
在二十世纪初，标准化的交流电的应用还比较少，于是“托普勒-霍尔兹Toepler-Holtz”静电发生器就被用来为一些医疗设备提供电源。它同时也为第一部X射线发生器提供能源。它有很多的附件，应用非常普遍，并且在西尔斯商品目录手册上大作广告。图中是非常大的一种静电发生器，由24块玻璃圆盘组成，4个一组，旋转时靠摩擦产生电荷。正在实施治疗的这幅图片选自1908年Franklin Gottschalk博士所著的《实用电疗》。

75-77)透热疗法图中是由匹兹堡透热医疗联合公司生产的仪器，可以产生微波，像这本小册中所述，可以应用于几乎任何疾病的治疗。

78-82)电疗椅: 电疗椅的应用是医学史上最大的讽刺性事件之一。早在1890年，麦金托什电椅就已经开始使用，并且出现在RW St. Clari博士所著的《电学在内科和外科中的实际应用》一书中。另外有很多类似的椅子，其中一种甚至在椅子的两侧装了两个紫色发光管来加强医疗（和眩目）效果。Wappler电子公司在这个领域非常活跃，这里有他们的一种产品，由一个电源柜来供电。还有Davis电子公司1904年生产的精美的瓷质健足盘。

83)电子浴: 早期的医师们彻底为电学的魅力所倾倒，他们甚至让病人进行“电子浴”，完全沉浸在电子里。Max Wocher and Son公司1904年的医疗产品目录中展示了一个用静电发生器供电的电浴装置。
84-86)Wappler电源壁柜: 电源壁柜为电疗椅和其他电子器械供电。图中是1904/1905年，西尔斯医疗商品目录上的产品，并且展示了一些医疗附件。Wappler电源橱柜的具有相同的功能。; .


88-89b) 瓶中的放射性: 放射性的潜力是巨大的，尽管还没有多少严格的验证，它很快就被见缝插针的商人盯上了。Revigator（放射领航者）是一个有低放射性的陶罐，水在其中存放后也就具有放射性。如图中所示，也有几种不同的外观。

<待续>

Normal distribution and Ellipse正态分布和椭圆、椭球

2007-11-19T20:15:00.000+08:00

正态分布和椭圆、椭球

二元正态分布（高斯分布）的等概率曲线是一个椭圆，而三元正态分布的等概率曲面是一个椭球。本文的目的是在已知概率分布参数（均值向量μ和协方差阵Σ）的情况下，在MATLAB中绘出等概率椭圆或椭球。
多元正态分布（multivariate normal distribution ）的概率密度：

很显然，其等概率密度曲面方程为

协方差矩阵作特征值分解：[V,D] = eig(Σ);使得

，即

因为,所以。
代入等概率方程，可以得到

作变量代换，，得到;
二元分布的情况下，y为包含两个元素的列向量，。D为对角阵，两个对角线元素分布为d1，d2；则上式可写为标准的椭圆方程

。
实际上，我们对原始椭圆作一个平移变换（x-u）和一个旋转变换（V*），可以使其中心平移到原点，长短轴的方向与坐标轴重合。而绘等概率椭圆的过程则与此相反，先用标准的椭圆方程产生组成曲线的离散点，然后经过相反的旋转变换和平移，得到原始的椭圆。

举例来说，首先，设置二维正态分布的参数，均值和协方差，并用mvnrnd产生一组符合此分布的随机数。

Mu = [2 3]';
Sigma = [0.9 0.4;0.4 0.2];
p = mvnrnd(Mu,Sigma,100);
plot(p(:,1),p(:,2),'.','MarkerSize',6)

设置半径，进行特征值分解

r =1;
[V,D] = eig(Sigma);

用linspace产生一个坐标轴（y）上的一组等间隔离散坐标值，再根据标准椭圆方程产生对应的x的坐标。

y = linspace(-sqrt(r^2*D(2,2)),sqrt(r^2*D(2,2)),60);
% compute x
x(1,:) = sqrt((r^2-y(:).^2/D(2,2))*D(1,1));
x(1,:) = real(x(1,:));

这只产生了半个椭圆，还要产生另一半（注意两条曲线的坐标旋转方向要一致），然后旋转，平移，画图：

Ellip = [x,-x(1,:)]; % x
Ellip(2,:) = [y,fliplr(y)]; %y
Ellip = Ellip'*inv(V); % rotate
Ellip(:,1) = Ellip(:,1)+Mu(1); %shift
Ellip(:,2) = Ellip(:,2)+Mu(2);
hold on;
plot(Ellip(:,1),Ellip(:,2));

plot(Mu(1),Mu(2),'+'); %Plot center

最终的效果：

另外，在对原始椭圆做旋转变换时，如果在V的前面再乘以一项，改为，则椭圆会变为圆。对多元正态分布的随机变量应用此变换，则其分布在各个方向上也变为均匀的。这就是信号处理中的白化变换。

三维正态分布的等概率曲面为椭球，其绘制过程也是类似的。
产生1/8曲面

xhalf = linspace(sqrt(r^2*D(1,1)),0,Nint);
Ninthalf = round(Nint/2);
zsect = zeros(Nint,Ninthalf);
ysect = zeros(Nint,Ninthalf);
for ti = 1:Nint
r2d = r^2 - xhalf(ti).^2/D(1,1);
ysect(ti,:) = linspace(0,sqrt(r2d*D(2,2)),Ninthalf);

zsect(ti,:) = sqrt((r2d - ysect(ti,:).^2/D(2,2) )*D(3,3));
xsect(ti,1:Ninthalf) = xhalf(ti);
end
zsect = real(zsect);

通过镜像产生1/4

%x>0,Z>0
xsect = [xsect,xsect];
ysect = [ysect,fliplr(ysect)];
zsect = [zsect,-fliplr(zsect)];

1/2

%x>0
xsect = [xsect,xsect];
ysect = [ysect,-fliplr(ysect)];
zsect = [zsect,fliplr(zsect)];

1/1

% make it a whole
xsect = [xsect;-flipdim(xsect,1)];
ysect = [ysect;flipdim(ysect,1)];
zsect = [zsect;flipdim(zsect,1)];

旋转、平移，完成。

% rotate
[lr,lc] = size(xsect);
for ti = 1:lr
for tj = 1:lc
newcodi = [xsect(ti,tj),ysect(ti,tj),zsect(ti,tj)]*inv(V);
xsect(ti,tj) = newcodi(1);
ysect(ti,tj) = newcodi(2);
zsect(ti,tj) = newcodi(3);
end
end
% shift
xsect = xsect+Mu(1);
ysect = ysect+Mu(2);
zsect = zsect+Mu(3);
surf(xsect,ysect,zsect);

Quack medicine: Imagination Goes Wild (2)伪医学：想象力的疯狂(2)

2007-11-13T18:49:00.000+08:00

Quack medicine: Imagination Goes Wild 伪医学：想象力的疯狂(2)

英文原文出自www.collectmedicalantiques.com ，链接http://www.collectmedicalantiques.com/quack2.html
翻译：liuxqsmile
本文地址：http://liuxqsmile.blogspot.com/2007/11/quack-medicine-imagination-goes-wild-2.html

磁

人们对磁现象的好奇心由来已久，而医疗者最早通过琥珀石相互摩擦时产生的磁性注意到这种现象。（23）。 “电子”（有些实际上是磁性）腰带（24-26），刷子（27-30）和梳子（31）在整个十九世纪都非常盛行，一直延续到现在。虽然这些东西大都附有“用了都说好”之类的推荐书，然而严谨的科学证明却从来没有出现过。

23)自然病理疗法
Darrin医生隐藏了他用磁石来治疗的事实（1900年）。

24-26）电子腰带
当金属扣浸入含有硫酸、醋的溶液时，会产生一个微弱的电流。Pulvermacher是这种带子的主要销售商，他们还印刷了大量宣传册子，里面满是使用者的溢美之词。图片中的两个腰带是由Mioxrl制作的，左边是男用的，有一个用于阴部的袢，右边是女用的。

27-30）“电子刷”
这个由Scott医生制作的刷子里面有磁石，除了脱发之外，还用来治疗其它很多疾病。里面还附有一个指南针来表明这个东西是有效力的。“射线”刷需要电池，而“奇迹(Wonder)”刷在按几下上面的平板后则可以感受到一次电击。

31）电子梳子
White牌电子梳由一个小电池充电，除了梳头时可以尝一下电击的滋味之外，毫无用处。

震动

震动施加于身体来治疗也有久远的历史了。它的理论是，在适当的频率下，你可以逐步将疾病晃悠出来，或者加快循环，使身体的虚弱部位得到改善。小件的家用器具，手动（32-37）或电动（38-47）的都曾经很流行。而在医院和疗养院里，则可以看到大的震动治疗器具，震动床、椅等（48-53）。

32－37）手动家用震动按摩器
VeeDee和Macura的装置使患者可以在家里随时对任何疾病进行局部震动治疗。患者的感受与低电压的电击很相似。

38-47）电动震动按摩器
在二十世纪初，更多的企业开始生产方便的电动家用医疗器械。White Cross按摩器附带很多配件，各自针对不同的疾病。从图上也可以看出，这些东西在专业的医疗机构里也能见到，并且老幼咸宜。震动带在二十世纪中期被用来减肥，而在此之前，它已经被用在无数的其它疾病上。

48-53）Kellogg疗养院的震动疗法
在密歇根州的巴特格里克的Kellogg疗养院，病人可能会有机会接受这两种装置的治疗。这种疗法的理论是，震动可以向患处提供更多的“治疗氧”。震动椅非常罕见，并且保存状况良好。

电
十八世纪的后叶，本.富兰克林(Ben Franklin)的风筝实验验证了闪电和电的关系。第一个将电应用在医疗上的是Christian A. Kratzenstein (1723-1795)，他用“好电”来置换风湿、瘟疫等疾病中的“坏电”。1854年，Davis和Kidder发明了利用磁场的发电机，旋转磁场可以产生电流，电流对人的刺激作用则被用来治疗一些神经疾病（54，55）。

54-55）电磁疗法
Davis&Kidder磁电仪是古董展览和拍卖场所的常客。两个相邻的磁石旋转产生电流，通过电缆和金属柄传递给病人。这个装置宣传可以治愈很多的精神紊乱疾病，产生的电击感觉让人觉得似乎有效。盒子很普通，但有海豚装饰的轮子则很少见。

<待续>
Quack medicine: Imagination Goes Wild 伪医学：想象力的疯狂(1)

Quack medicine: Imagination Goes Wild 伪医学：想象力的疯狂(1)

2007-11-07T20:51:00.000+08:00

Quack medicine: Imagination Goes Wild 伪医学：想象力的疯狂(1)

英文原文出自www.collectmedicalantiques.com ，链接http://www.collectmedicalantiques.com/quack.html
翻译：liuxqsmile
本文地址：http://liuxqsmile.blogspot.com/2007/11/quack-medicine-imagination-goes-wild-1.html

只要有疾病的存在，患者总是会寻找医治的方法，草药、咒语或者最新的遗传学方法和经科学证明的手段。但是，由于医疗手段的效果和后果并不总是立竿见影，骗术并不总是会被识破，所以就总是会有庸医和江湖骗子的市场(1)。许多的所谓治疗手段出现，又消失，都需要科学研究进行最终的判定。这种观点最早在1628年才由William Harvey（威廉.哈维）在其伟大著作《动物心血运动解剖论》中提出。骗术和医术的界限总是很模糊的，甚至到今天仍然如此。

1)Kidd医生的惊人宣言
骗子宣扬他们的玩意儿和医术时总是口无遮拦，夸大其词。图中是1900年，《二十世纪评论The Twentieth Century Review》上的一个整版广告：包治百病

骗术登堂入室

第一个获得专利的伪医学专业器械由一些金属探针组成，用来“拽出”疾病，美其名曰“伯金斯牵引器”，于1801年获得专利(2,2a)。从那时起，医疗手段的花样变幻极尽想象之能事，包括水疗，光和色彩疗法，热、磁、振、电和放射性，数不胜数。

2,2a)Gillray的伯金斯牵引器
伯金斯牵引器是医疗骗子所用器材的鼻祖，收藏者的首选藏品。这套特殊的装置由Sandra Louise Perkins销售，他是发明者Elisha Perkins医生的直系后裔。其中的一种名为ìPerkins的探针获得了tractorí专利，用镶金的大红皮革包装。在十九世纪初，讽刺画在英国非常流行的时候，当时著名的讽刺画家James Gillray制作了一幅彩色版画，名为“ìmetallic tractorsî”，来讽刺当时的医疗界。

颅相学

通过身体相貌来评价人的品格和智力是荒谬的，同时也是种族主义和男性至上（大男子主义）的通称表现。但是这门“科学”在十九世纪以颅相学和人相学（相面术，见颅相）的面目获得极大的重视。颅相学由维也纳医师Franz Joseph Gall(弗郎兹·约瑟夫·加尔，1756-1828)建立，由Joseph Spurzheim在1832年介绍到美国。这种理论认为，大脑的形状和体积决定了颅骨的外观，后者进而反应了人的个性(3-5)。依据人相学的观点，体格特征影响个性特征。例如，一个胸膛宽阔的人倾向于追求感官享受，而有强健肌肉的个体则性格坚强。体形协调匀称的人是有才能的，而瘦小的则是冲动的（6）。比较颅相学，通过比较人与动物体形特征的相似性来试图验证颅相学的观点，在这个世纪之交(十九到二十世纪)也非常流行(7-10)。

3-5)人相与人性的研究
OS Fowler1869年在《实践颅相学家》一书中标出了特定的性格特征所对应的位置。这个由LN Fowler所制作的瓷制半身像，在十九世纪很常见，上面绘有相似的标记，。收藏者应该注意这个东西有很多现代的复制品。

6)亨利八世和他的个性
Mary Stanton1913年的著作《面相术百科》（"The Encyclopedia of Face and Form Reading"），对许多名人“以貌取人”，特别请注意那些看起来像某位国王的面孔。

“脸部的综合特征揭示了其性格的自私、专制和荒淫，再加上惊人的野心，其品性令人憎恶。”

7-10)动物形貌与人的个性
许多的出版物认为相貌可以预测人的个性，连动物的性情在这种联系中也没有被忽略。最迟在十九世纪，人们已经注意到了宠物和它主人的相似性。这些图片摘自Samuel Revells 1866年出版的《新人相学》（"New Physiognomy or Signs of Character" ），和医学博士James Redfield 1887所著的《比较人相学》（"Redfield's Comparative Physiognomy"）。

治疗及其手段

在技术手段能够为这些术士提供迷惑人心的嗡嗡声(buzzers)和铃铛之前，他们大多都满足于借用随手可得的东西，包括水 (11,12)，光(13-16)和温度(17,18)。一般的目的是通过加快或减缓局部血液循环，根据疾病的情形，将氧供应到相应的身体部位。古埃及人崇拜太阳和阳光的治疗能力，而罗马人则通过洗浴，不同的病痛使用不同的水温。按摩疗法古代就开始使用，到今天仍然流行，因为从理论上，可以改善患病部位的血液循环和供氧。二十世纪早期，出现了很多宣称可以向患处提供更多“治疗氧”的装置，但实际上一点帮助都没有（21，22）。

11-12)水疗
在很多的疗养院，通过水疗来试图减缓一些慢性疾病的症状。这两幅图片选自医学博士JH Kellogg 1902年所著的《合理水疗》（"Rational Hydrotherapy"），一幅是一个肝病患者正在进行水疗，另一幅是一所位于巴特尔克里克（Battle Creek）的水疗研究所。

13-16) 光疗
戴上眼镜保护罩后，利用弧形灯产生的强光进行光疗。这幅图片选自JH Kellogg 1927年所著的《光疗》，表现了这种备受重视的治疗手段，光疗柜的图片也是出自这本书。铬光谱仪（Spectro-Chrome devices）在二十世纪初也非常流行，拿来治疗一切疾病。它由Dinshah Ghadali设计，将一束强光通过一些不同颜色的滤镜，针对不同的疾病，可以产生各种各样的“Dinshah调制光波”。

17-18) 暖帽
联合默克研究所设计了这种暖帽来促进头发生长。然而，没有任何证据能体现这个加热电灯泡的疗效。

19-20) 眼睛按摩器
锻炼可以使某些肌肉更强韧，为什么视力不可以呢？这就是眼睛按摩器的理论依据。第一种按摩器在伦敦出现，称为“眼科神经活力专家”，第二种出现在纽约，名曰“完美视力恢复器”。

21-22) 电子垫（Electropoise）
类似的装置也曾经非常流行，用来向患处提供更多的治疗氧。当然，拆开来看，是找不到任何可能起作用的部件的。

<未完，待续>

Insert link into 163blog:在网易博客中添加链接

2007-11-01T20:04:00.000+08:00

163的blog中没有内置这项功能，虽然实现的方法很简单，添加几行HTML代码就可以，但对我这样的菜鸟来说，还是很难想到。
在设置中选择装扮博客，然后选择设置首页内容。单击“添加自定义模块”，在模块名称中输入名称，如“友情链接”。然后在下面的输入框中输入HTML代码：

...... ......

<u>表示下划线标记；<br>表示换行，以留出间距。H3表示对中间的文字（即“网站名称”）应用三级标题的格式。可以选择H1，H2，...H6，以选择不同的字体大小。H1最大，H6最小。

选择确定后，这个外部链接模块就出现在首页上了，拖放到合适的地方就可以了。

本来,这个简单的事情到这里就应该结束了。但是,一预览，发现中间的两行代码被解释了，原始的html代码没有了，只有几个带链接的文字。怎么在网页中显示html代码呢？

一种方法是使用<textarea>和</textarea>，强制中间的部分以文字形式显示。另外对于较短的标记，如<br>，则可以用'&'+'lt;'来代替前面的<(不要中间的+)。

参考：http://www.01on.com/a/404.html
http://lxq0xxx0.blog.163.com/blog/static/2587419420071018618331

可怜的Blog

2007-10-31T16:24:00.000+08:00

今天又把MSN Space重新翻出来。现在已经有了三个Blog。Blogspot总是被和谐，而网易的又不知道为什么Google和Baidu总是索引不到其中的内容，只能链接到blog的标题。看着不知道现在的Live Space又怎么样。

我可能会把前面的一些内容断断续续地粘贴到这边。

Blogspot(blogger)有一个优点是可以自己指定日志的发表日期，而Live Space和163的blog都自动指定为当前时间，不可修改，未免显得有些死板了。

Graeme Obree: Homegrown Hero 格莱姆.奥波利：平民英雄

2007-10-17T20:17:00.001+08:00

刚看了一部电影《飞翔的苏格兰人》("The Flying Scotsman") ,对里面的主人公Graeme Obree有一些好奇。正好的BBC的网站上看到一篇介绍他的文章，就打算翻译过来。水平所限，翻译时并不强调词句的精确，只求达意。

作者：Gordon Cairns
原文地址：http://www.bbc.co.uk/scotland/sportscotland/asportingnation/article/0037/index.shtml；
翻译：liuxqsmile

曾经两次创造自行车一小时记时赛世界纪录的苏格兰人Graeme Obree的运动生涯与一件普通的家用电器有着不解之缘。而正是这个家用电器吸引了人们太多的注意，忽视了他所取得的巨大成就。另一方面，尽管只有极少的人曾经打破一小时记时赛的世界纪录，但Obree的公众形象却是一个体育怪人，而不是一个冲刺极限的英雄。

作为一个来自自行车运动的穷乡僻壤的业余运动员，能够在1994年打破世界纪录确实令人难以置信。Eddy Merckx，或许是有史以来最优秀的自行车运动员，在1972年创造了一小时记时赛的纪录，赛后，他说：“绝不会有下次！” 这就是他在这项比赛中所遭受的痛苦。对Obree来说，在自己制造的自行车上打破纪录更是需要勇气。

Obree说，“我最大的遗憾就是曾经跟一个记者说起，我自行车上的一个小零件是从洗衣机上拆下来的。现在，所有人都会把我跟洗衣机扯上关系，再也摆脱不掉。”

Obree承认，和洗衣机的联系减少了他赢得大赞助商的机会。“大赞助商是不会跟使用拼装车的运动员签约的。”

在十年之后，Obree的成就重新获得了认可。在他非常坦诚的畅销自传中，清楚地刻画了他冲破权威的阻挠，战胜自我，最终取得成功的过程。现在Obree的经历已经成为苏格兰场地自行车新星Craig McLean 和Chris Hoy的一个奋斗动力。

Chris Hoy曾经说，“当我刚出道时，我就把Graeme Obree作为我的精神偶像。他的经历就像是一部好莱坞电影，但是，如果你是电影院里看到一个像他那样的故事，一定会觉得那是不真实的。”

Obree很感激这位奥运会冠军的高度评价。他说：“我能够成为像Chris这样有如此成就的人们的一个鞭策，实在是我的荣幸。”

实际上，现在正有人准备根据Obree的经历制作一部电影，这样，继苏格兰传奇体育冠军Eric Liddell（译注：埃里克.利德尔，奥斯卡获奖影片《火战车Chariots of Fire》的原型人物，二战中死在中国）之后，Obree的故事也将走上大屏幕。

Obree早期的焦点在于1990年代中期同英格兰人Chris Boardman的竞争。但是对于Obree来说，能够挑战奥林匹克金牌本身已经是勇气战胜逆境的杰出表现。

最终，Boardman获得了1992年巴塞罗那奥运会的金牌，Boardman的胜利让Obree认识到，在世界自行车运动的舞台上，这是一位有力的竞争者。正是和这位英国“骑士”在80年代末90年代初的对决让他倍感欣慰。

“Chirs赢得金牌让我认识到，我们都有很高的水平，我们针锋相对过数次，互有胜负。”

当Obree正在从车库的杂货堆里寻找零件来制作自己的自行车时，当他为微薄的薪水挣扎，甚至付不起电话费时，Boardman的身后却有一整队的运动专家，心理学专家，以及自行车厂商。Obree并没有因为这些悬殊的差距而抱怨。

"Chris的优势在于拥有专业的团队，他的成功路线组织得非常好，以自己为中心建立团队，逐步取得更好的成绩，扩大影响从而赢得赞助，良性循环，这样使他能够完成环法大赛。

“但对于我来说，只想着打破一小时记时赛的世界纪录，当我达到这个目标后，我不知道下一步该做什么。我根本没有什么职业规划。

“是Chris Boardman掌握着他自己的命运，对我来说，没有什么好抱怨的。”

当然，Obree同样掌握着自己的命运，去冲破周围的阻力，争取跻身于世界纪录创造者的行列。

他是，他似乎并没有从他的成就中获得多大的欢乐。

Obree承认，对失败的恐惧，而不是对成功的渴望，才是他成功后继续前进的最大动力。退役以后，他被诊断出患有燥郁症（bi-polar disorder),这种病的患者经常陷入烦躁和忧郁情绪，精神状态受到严重影响。

“作为一个人，我要尽力达到自己所期望的价值，正是这点激励我去创造纪录。当我听到意味着我超越Moser的枪声时，我倍感轻松。

“即使是天花板要塌了，无论发生什么，在这之后，我已经成为新的纪录，没有人可以把它抢走。”

Obree认为，1995年在哥伦比亚获得的世锦赛4000米追逐赛冠军是最有趣的经历。并不是简单地因为获得了梦寐以求的彩虹衫，而是因为和UCI（国际自行车联盟）主席Hein Verbruggen的积怨。后者在前一年的比赛的前一晚，通过了一项对Obree的骑行姿势极为不利的规则。使Obree未能完成比赛，而这位始作俑者当时正是塞道旁边的旁观者。

“我最兴奋的时候可能是在哥伦比亚获得世界冠军，再加上Verbruggen正好给我颁奖，更让人觉得惬意了。”

令人丧气的是，这种骑行姿势被禁的事情共发生了两次。令Obree遗憾的是，没有像Boardman那样有职业团队支持，错过了参加环法大赛的机会。1994年底，他曾经加入法国Le Groupement车队，但两天后车队就以奇怪的借口将他开除：Obree从没有为团队而骑车。

“我职业生涯中最大的遗憾就是没能参加环法比赛。我喜欢这项赛事，我觉得我有很好的机会来获得冠军。”

Obree现在正在忙着写一本训练手册，并且让自己的腿保持良好的状态。在描述他非凡人生的电影中，这双腿将扮演它自己的角色。

本文在163上的镜像：http://lxq0xxx0.blog.163.com/blog/static/25874194200791983341792

删除系统托盘区的无效图标

2007-10-12T21:06:00.000+08:00

Windows XP用的时间长了，系统托盘区可能会有许多无效的旧图标，长长的一串，看起来挺烦人的。怎么样删除呢？
删除注册表项：
HKEY_CURRENT_USER\Software\Microsoft\Windows\CurrentVersion\Explorer\TrayNotify
中的子键IconStreams 和 PastIconsStream是网上很多人都提到的方法，但对XP SP2似乎不起作用，实际上再加上一步就可以了：
打开任务管理器，强制结束进程：explorer.exe，然后新建任务，选择windows目录下的explorer.exe，重新打开桌面，那些旧图标就全都不见了。
参考：http://www.theeldergeek.com/notification_area.htm

连续小波变换应用于密集模态参数识别

2007-09-27T17:31:00.000+08:00

我的第一篇文章：
柳小勤，岳林，朱如鹏. 连续小波变换应用于密集模态参数识别. 南京航空航天大学学报，2007，39（4），496-500.
并没有很好地解决问题，惭愧。参数拟合区间始终是个大问题。文中没有提到，因为我不知道该怎么比较好地解决。

又见：
http://lxq0xxx0.blog.163.com/blog/static/2587419420078275315470

MATLAB中的振型动画工具：ModeVibra

2007-09-12T15:32:00.000+08:00

ModeVibra是一个简单的模态分析工具，用来查看结构的振型动画。还可以对节点重新排序和生成测量自由度表。
本工具使用MATLAB的GUI平台开发，提供源代码，用户可以自由使用并可以根据自己的需要修改代码。
详细说明说明及下载：http://groups.google.com/group/liuxq-works/web/modevibra?hl=zh-CN

DVR中MP4到AVI的转换

2007-09-07T14:36:00.000+08:00

实验室最近买了一个摄像头、一块视频采集卡，用一个DVR软件来录像，将视频以mp4文件格式保存（昕鼎公司的产品）。由视频压缩卡“硬压缩”出来的文件格式跟一般的mp4文件不同，用一般的播放器都无法播放，更不用说来作一些其它的分析、处理了。所以想要转换成更一般的格式。
通用的播放器中，暴风影音和KMPlayer都不能播放，只有MPlayer可以播放，只是不能正常显示时间长度。试验了一些MP4文件的转换工具，都不能读入文件。
后来想到，既然MPlayer可以播放，那MEncoder应该也能处理这种的视频。在MPlayer的网站上下载了Windows版本（文件夹下包含MEncoder.exe，不需要安装），然后下载一个GUI工具，MeWiG，放到MPlayer的目录下，运行，第一次运行时会先弹出配置窗口，一定要启用测试中的功能。然后就可以转换了。在我的系统上，只能用xvid压缩成divx，才能被其它的播放器识别。

这种方法也只能是对某些DVR系统的MP4文件起作用，我们另一套DVR系统的文件就不能用MEncoder转换。

申请Google Group

2007-09-03T13:42:00.000+08:00

今天申请了一个Google Group，地址是：
http://groups.google.com/group/liuxq-works
提供100MB空间，可以把一些文档放在上面，在国内也可以访问到。就是不知道将来会不会又被封，先用着吧，如果到别的网站上申请空间，又要重新注册，用户名、密码一大堆，太繁琐了。

Linux下内网BT工具：benliud奔流

2007-08-22T14:23:00.000+08:00

Linux下的BT软件有很多，但能在内网使用的却几乎没有。

一种方法是用Wine模拟运行Windows下的BT软件。可以下载一个免安装的BitSpirit压缩包，解压到Wine的Program Files目录，然后用Wine启动就可以了。Wine编译很耗时间，而且对中文的支持有问题，在我的FC6上，在启动命令前加上：lang=zh_CN.GB2312，可以使BitSpirit任务列表中的中文文件名正常显示，但下载到硬盘的实际文件名仍然是乱码。

不过，最近国内有高手做了新的Linux BT客户端：Benliud奔流。Wine和BitSpirit都可以省了，直接下载软件包运行就可以了。作者主要针对Ubuntu7.04，但在我的FC6上也可以正常运行。速度还是很快的。
讨论区见这里：http://groups.google.com/group/benliud

爱国者MP3的录音文件导入MATLAB时一点问题

2007-08-13T13:11:00.000+08:00

想对动物的叫声做一下频谱分析，正好手边有爱国者的一款MP3可以录音，就录了下来。但是分析时出了问题，这个MP3的wav文件不能直接导入到MATLAB中，而Windows中的Sound Recorder录的却可以。后来发现这两个wav的文件头不一样。我用以下方法解决这个问题：

(1)用Windows附件中的Sound Recorder打开爱国者的录音WAV文件，选择另存为，在Format中选择参数PCM，8KHz，8bit，单声道Mono，假定文件名为aigo.wav;

(2) 用Sound Recorder录一个文件，以相同的格式保存，如pc.wav；

(3)用十六进制编辑软件（如PSPad，A.X.E.等）打开这两个文件，在pc.wav中找到fact这四个字符，注意在aigo.wav中也有这几个字符，只是位置不一样。用pc.wav中fact前面的内容覆盖aigo.wav中fact前面的内容。保存之后就可以了。

我用的爱国者mp3录音的精度是4bit，而第一步中转换到8bit，对数据的影响应该不大。

Ubuntu on USB-HDD: 移动硬盘上的ubuntu操作系统

2007-08-06T10:06:00.000+08:00

实际上在USB移动硬盘上安装Ubuntu还是挺简单的，只是我起点太低，折腾了好几天，才算终于启动起来。

我是从硬盘来安装的，从网上下载了安装镜像和引导文件，放在一个FAT32分区上，然后用GRUB4DOS来引导安装。这一步是我以前装Fedora Core时就弄过了的，所以比较顺利就进入了安装界面。一开始我用的是Ubuntu7.04的desktop镜像，安装时总是报错，大概是无法安装光盘文件系统，找不到安装文件。到网上查了一下，很多人都有这个问题，需要改用alternate来安装。

设置磁盘分区时，把我的IDE硬盘认成了sda，而移动硬盘则成了sdb。如果你需要使用USB-HDD模式来启动，最好将移动硬盘的第一个分区mount成boot，并且容量不要设置得太大（网上有人设置为512MB，我设置为256MB，还没来得及验证到底可不可以，因为我的机子不支持USB-HDD。这个可能是因为有些BIOS的磁盘寻址范围有限。）

安装过程没有什么可说的，软件包没有选择，自动安装。然后就是引导的问题。我的机子支持USB-ZIP的模式启动，将一个U盘用GRUB4DOS作成启动盘，然后将boot目录下的initrd和vmlinuz拷到U盘根目录下，在menu.lst中设置

title Ubuntu
root (fd0)
kernel (fd0)/vmlinuz-2.6.20-15-generic rw root=/dev/sdb2 rootdelay=10
initrd (fd0)/initrd.img-2.6.20-15-generic

sdb2是根文件目录。实际上根本不需要对initrd作出任何修改，它已经默认支持USB磁盘了。如果是从USB-HDD模式引导，应该也没有什么大的差别。

我花了不少时间在怎么制作支持USB的启动盘initrd这一步上，虽然最后没有成功，但对linux的文件系统有一些了解，还是有些收获的。

在别的机子上也可以运行，不过图形界面要重新配置：

sudo dpkg-reconfigure -phigh xserver-xorg

还是会把IDE硬盘认成sd*，别的机子上也一样，又在网上查了一下，说是新内核把PATA和SATA统一识别为SATA。

参考文献：