MATLAB 函数编写方法与应用

所谓 MATLAB 程序，大致分为两类：M 脚本文件 (M-Script) 和 M 函数 (M-function), 它们均是普通的 ASCII 码构成的文件。M 脚本文件中包含一族由 MATLAB 语言所支持的语句，它类似于 DOS 下的批处理文件，它的执行方式很简单，用户只需在 MATLAB 的提示符 >> 下键入该 M 文件的文件名，这样 MATLAB 就会自动执行该 M 文件中的各条语句，并将结果直接返回到 MATLAB 的工作空间。M 函数格式是 MATLAB 程序设计的主流，一般情况下， 不建议您使用 M 脚本文件格式编程。
MATLAB 的 M 函数是由 function 语句引导的，其基本格式如下：

function [返回变量列表] = 函数名 (输入变量列表)
注释说明语句段, 由 % 引导
输入、返回变量格式的检测
函数体语句

这里输入和返回变量的实际个数分别由 nargin 和 nargout 两个 MATLAB 保留变量来给出，只要进入该函数，MATLAB 就将自动生成这两个变量，不论您是否直接使用这两个变量。返回变量如果多于 1 个，则应该用方括号将它们括起来，否则可以省去方括号。输入变量和返回变量之间用逗号来分割。注释语句段的每行语句都应该由百分号 % 引导，百分号后面的内容不执行，只起注释作用。用户采用 help 命令则可以显示出来注释语句段的内容。此外，正规的变量个数检测也是必要的。如果输入或返回变量格式不正确，则应该给出相应的提示。我们将通过下面的例子来演示函数编程的格式与方法。


[例子]：假设我们想生成一个 n*m 阶的 Hilbert 矩阵, 它的第 i 行第 j 列的元素值为 1/(i+j-1)。我们想在编写的函数中实现下面几点：

1）、如果只给出一个输入参数，则会自动生成一个方阵，即令 m=n
2）、在函数中给出合适的帮助信息，包括基本功能、调用方式和参数说明
3）、检测输入和返回变量的个数，如果有错误则给出错误信息
如果调用时不要求返回变量，则将显示结果矩阵。其实在编写程序时养成一个好的习惯，无论对程序设计者还是对程序的维护者、使用者都是大有裨益的。

采用 MATLAB 函数编写格式和上述要求，我们可以编写出一个函数

function A=myhilb(n, m)
%MYHILB a demonstrative M-function.
% A=MYHILB(N, M) generates an N by M Hilbert matrix A.
% A=MYHILB(N) generates an N by N square Hilbert matrix.
% MYHILB(N,M) displays ONLY the Hilbert matrix, but do not return any
% matrix back to the calling function.
%
%See als HILB.

% Designed by Professor Dingyu XUE, Northeastern University, PRC
% 5 April, 1995, Last modified by DYX at 21 March, 2000
if nargout>1, error('Too many output arguments.'); end
if nargin==1, m=n;
elseif nargin==0 | nargin>2
error('Wrong number of iutput arguments.');
end
A1=zeros(n,m);
for i=1: n
for j=1:m
A1(i,j)=1/(i+j-1);
end, end
if nargout==1, A=A1; elseif nargout==0, disp(A1); end


这样规范编写的函数用 help 命令可以显示出其帮助信息：

>> help myhilb
MYHILB a demonstrative M-function.
A=MYHILB(N, M) generates an N by M Hilbert matrix A.
A=MYHILB(N) generates an N by N square Hilbert matrix.
MYHILB(N,M) displays ONLY the Hilbert matrix, but do not return any
matrix back to the calling function.

See als HILB.


有了函数之后，可以采用下面的各种方法来调用它，并产生出所需的结果。

>> A=myhilb(3,4)
A =
1.0000 0.5000 0.3333 0.2500
0.5000 0.3333 0.2500 0.2000
0.3333 0.2500 0.2000 0.1667
>> A=myhilb(4)
A =
1.0000 0.5000 0.3333 0.2500
0.5000 0.3333 0.2500 0.2000
0.3333 0.2500 0.2000 0.1667
0.2500 0.2000 0.1667 0.1429
>> myhilb(4)
1.0000 0.5000 0.3333 0.2500
0.5000 0.3333 0.2500 0.2000
0.3333 0.2500 0.2000 0.1667
0.2500 0.2000 0.1667 0.1429

另外，因为 MATLAB是一种解释性语言，所以即使在某个或某些函数中存在语法错误，但如果没执行到该语句时可能就不会发现该错误，这在一个成功的程序设计中是不能容许的。要查出某目录中所有的M函数语法错误，首先应该用 cd 命令 进入该目录，然后运行 pcode * 命令进行伪代码转换。因为该命令会将 MATLAB 函数转换成伪代码，而在转换过程中该程序将自动翻译每一条语句，所以一旦发现有语法错误，将会停止翻译，给出错误信息。改正了该语法错误后，再重新执行 pcode 命令，直到没有错误为止。至少这样会保证目录下所有的程序不含有语法错误。

colloquium:2011-2-16

TVA-NGDC Activities

Wally Justice

site manage engineering

美国之音 中国政府的“自律反腐说”令人失望

中共中央总书记胡锦涛日前发表讲话，表示要加强打击腐败的斗争。与此同时，评论人士说，在中国很多共产党和政府官员凭借不受制约的权力肆意贪污腐败导致民怨的今天，胡锦涛再次强调通过领导干部廉洁自律、加强思想道德教育来制止腐败，令人感到失望。

有关腐败的最著名的言论来自19世纪的英国阿克顿勋爵。他说："权力导致腐败，绝对的权力导致绝对的腐败。"领导共产党夺取中国大陆政权的中共已故领袖毛泽东虽然没有明确表示接受西方民主制度的权力制衡的思想，但是他在60多年前也明确表示怀疑绝对的权力会进行自我监督，进行有效的自律，保证自身不腐败。

毛泽东1945年夏天表示，中共已经找到了一条新路，可以摆脱过去历代王朝在腐败中倒台的周期律。他说："我们已经找到了新路，我们能跳出这周期律。这条新路就是民主。只有让人民来监督政府，政府才能不敢松懈。"

观察人士普遍认为，毛泽东所说的民主只是口惠而实不至。共产党实行了将近60年的一党独裁，结果是发展到今天，中国进入了不少中国民众所说的几乎"无官不贪"时代。

在这种情况下，中共新领袖胡锦涛发表了官方媒体所说的有关反贪污腐败的重要讲话，突出强调"要从思想道德教育这个基础抓起，不断夯实廉洁从政的思想道德基础、筑牢拒腐防变的道德防线"。

有评论人士认为，这种讲话是空话、废话，了无新意，因为毛泽东还敢从人民监督政府的角度谈民主，而如今提到民主的时候则是讲"严格执行民主集中制"，也就是强调专权。

在中国公众对党政官员凭借不受制约的权力肆意贪污腐败的今天，胡锦涛最新的有关反贪污腐败的讲话的确令人失望，因为他的观念和说法至少是落后的。

中共政府官员的肆无忌惮的贪污腐败，根基是中共权力没有公众的制约，只要没有制度性的改变，党政官员的腐败只能是一如既往，愈演愈烈。

Technorati : 反腐 民主

中国“大屋顶”政治如何完善？

去年11月到新加坡一游，发现在政府自上而下的治理结构上，中国和新加坡有很多相似之处。新加坡也是只有一个政党当家，但有直接选举。执政党能保证每次大选中百无一失。人们不禁要问，为什么中国不学一学呢？新加坡地方小，要藏拙是有困难的。我们在新加坡市容上看不到贫民窟和穷人，说明政府管理是到家了。这又是值得中国学习的地方。

持续两千年的大国架构

　　话又得说回来，新加坡那么小，中国那么大，不可同日而语。中国是世界上最最特殊的大国，有两大特点值得注意。第一，世界上有十条大河，却只有长江与黄河不是国际公河。中国是几千年来举世无双的大型灌溉农业高度发达的国家，被二十世纪美国著名汉学家、前德国共产党领袖魏复光（Karl August Wittfogel）誉为超大型水利社会。

　　第二，正因为中国具备这种水利独特享受，欧亚大陆东端不同人种、民族、部落、群体纷纷加入中国，使中国变成容纳人类一大部分人口的政治"大屋顶"。根据英国经济历史学家麦迪逊（Angus Maddison）统计：从公元50至960年，中国人口是全球的六分之一，1280至1700年上升到四分之一，1820年更是全球三分之一，1995年是全球五分之一。哈佛大学中国研究泰斗费正清在1980年代形象地比较说：世界上有两个10亿人口，一个分散在欧洲和美洲的五十多个国家内，另一个却收容在中国一个政治实体中。一个"大屋顶"下容纳了欧美两洲人口的总数。

　　从秦始皇公元前221年吞灭六国开始，中国和"大屋顶"政治结缘已有22个世纪。在这漫长时期中，地球上许多"大屋顶"政治都垮掉，中国也饱受外力侵略。今天形势不同，在紧邻出现了另外一个"大屋顶"政治--印度。还有，大洋彼岸的美国也算是"大屋顶"政治，虽然人口比中国、印度少得多。中、印、美这三个"大屋顶"政治还有一个共同特点，就是在现代化的进程中调整了人民的空间，不但重新分配了财富，而且重新分配了尊严（例如中国农民翻身、印度"贱民"当上总统、美国黑奴后裔竞选总统）；后一点比前一点做得更好。中国应该从美国与印度的模式中取长补短，然后走出自己的道路。

　　中国不断把行政官员送到美国各大学学习，有的学完回国，听说已经进入中央和地方的决策层了。许多人批评中国也多半援引美国例子作正面教材。我认为中国可以学习美国的现代管理、教育与科研的卓有成效的环境（美国各大学好像在包办科技与社会科学的诺贝尔奖金），学习它的政府体制把"官"隐藏或埋没，学习它公私机关办公人员接待群众的礼貌与谦虚，学习它全社会对生命的重视以及救死扶伤的认真，学习它知识精英、特别大学校园内"天下一家"而较少肤色、人种与宗教信仰的成见、言论与活动的自由等等。

美国和印度模式

　　但是美国的"大屋顶"政治也有很多弊病。可以从四方面来看。第一，美国最吃香的人是表演明星、包括会做广告的企业家和善于蛊惑人心的政客。他们的共同特点是戏剧表演，观众欣赏时喜欢大喊大叫。美国许多领导人如果以中印文明的"真善美（satyam shivam sundaram）的标准来打分是不及格的。如果拿孙中山的"圣、贤、才、智、平、庸、愚、劣"八个等级标准来衡量，美国最近两次赢得总统大选的人肯定不属于前四类。美国民意测验中对总统的支持率一直偏低，最近几个月落到30%以下。在任何别的制度下，这样的政府早就垮台了。总的来说，美国是个奇怪的提倡"平庸"、甚至"愚劣"政府的国家。中国绝不能学这种模式。

　　第二点，美国的明星式的政治人物，都是伪善的个人野心家。美国舆论评论起来对这点也不讳言。美国政治的最大毛病叫做partisanship（党派行为），小布什执政期间达到新的高潮。美国内政外交的很多问题都在共和、民主两党"党派行为"的拉锯争夺战中变得无法解决。

　　第三点，美国政治的暗点是公开性的lobbyism（游说疏通），有登记的、合法的lobbyist（说客）机构和人员可以自由在国会山活动。说穿了，这就是"钱权交易"。美国国会基本上是受大军火商、石油公司、制药公司、农场主组织等控制的。

　　第四点，美国是国际上最逞强、最富侵略性、最容不得匹敌、名气最臭的大国，一个"民主"国家怎么会变成这样，值得人们深思。中国千万不能朝这个方向发展。

　　再看印度的"大屋顶"模式，和中国相比有"民主"的优势，在两方面值得称道：一是政府如果不为老百姓办事，选民就一定会把它撤换（不要说中国，连美国也做不到）；二是民间的组织强大，政府必须倾听他们的诉求，甚至在他们强大的抗议面前让步。

　　但是，印度模式的明显缺点有三：一、政局不稳定，经常举行中期选举，难有政策连贯性，昙花一现的掌权政府没有远景发展计划；二、离心力大，特别是最近十年来进入联合政府时期以后，小党乘机在两大阵线之间讨价还价，朝野之间与政权中的党派之间相互牵制，很难贯彻重大改革；三、统治、管理质量与效率不高，最近首都新德里发现的一富家的主仆两人拐骗、强奸了十余女孩，然后把她们杀死，割掉器官出售，由于受害者都是附近贫民窟住的、从孟加拉国逃来的难民，警察听到失踪报警后不但不调查，反而讥讽说："你们生这么多孩子干嘛？"这种丑闻使得国际舆论震动。

中国的强势与弱势

　　在三个"大屋顶"政治中，中国显出三大优点：一、政府强大、自上而下有力贯彻各种改革政策，人民向心力强，全国一盘棋，因此出现发展奇迹；二、强调崇高理想，与时俱进，不时推出新的精神，虽然有点口号式、教条主义、言过其实，但具有精神文明特色，堪称当今世界一面文明旗帜；三、众志成城，几亿人民发奋图强，全国各地你追我赶，经济发展全面开花。

　　中国也有不如别人的地方，说来话长。主要一点是自上而下的压力过大，自下而上的动力不足。所谓"世上好言佛说尽"，只有"最高"脑子里出的智慧才有价值，这样的体制就违反人类进化的规律。君不见历史上所有重大发明、改革都出自民间！

　　还有一点，中国政治"官"气发酵，民间文化活动似乎没有"官"来题字、剪彩、致辞就不光彩。中国锦绣山河、旅游胜地总有些比小学生还差的书法高高挂，这又和美国"治人者"的"平庸愚劣"相辉映了。可是在美国，人们很少拜"官"。中国则不但拜官，而且买官、卖官、嫁官、贿官、怕官、骂官、告官，官势、官架、官文、官腔、官伞、官私、官帮等不一而足，使社会多一层污染。中国不改掉"官本主义"就不会有完美的现代化。

　　共产党人应该和中国两千多年的"官本位"水火不相容的，为什么却让这样的社会污染在中国散布呢？美国是个人英雄主义泛滥的国家，但总统在记者招待会上用尊敬的口气叫别人sir，一般人对领导都直呼其名（不加称号）。中国是提倡集体主义的"社会主义"国家，芝麻大的领导也必须在称呼中表现出来：什么"黄总"、"白队（长）"的，比资本主义还肉麻！听说"同志"已经变成贬词、忌词--代表"同性恋"了，这是什么社会进步呀？中共要移风易俗，先要把社会语言的质量提高一番才行。

　　怎么完善中国"大屋顶"政治是知也难，行也不易。但是应该及早着手，一方面为了"大屋顶"下的13亿人改善生活着想，另一方面也为了改善中国的形象。当中国人民自己为"崛起"沾沾自喜时，却看不到别国对中国的批评与不满。

　　美国越来越多的工厂关闭、中产阶级破产、工人失业，都怪罪于中国。在这次11月国会改选的大辩论中，骂娘的声音很响；主要咒骂对象是布什政权，其中也夹杂着中国的名字。有人说白宫被北京的"独裁者"摆布，还有人甚至说中国产品价廉是因为中国工人的非人道"奴隶"待遇。许多说过类似的话的民主党人都进了新国会，即将向布什政权发起"秋后算账"，要把free trade（自由贸易）改为fair trade（公正贸易）。中国"大屋顶"政治的缺点将会成为民主党人攻击"不公正贸易"的主要目标之一，中国不能掉以轻心。

　　·作者是从印度退休的学者，现旅居芝加哥

Technorati : 中国 政治

中国军机再度失事：北海舰队军机青岛市区内爆炸

解放军北海舰队一架小型军用飞机11日晚在山东省青岛市区株洲路附近迫降后爆炸，飞行员提前跳伞，没有受到爆炸的伤害。另外，由于机师选择迫降的路面较宽，没有对群众和财产造成损失。这已经是解放军进入2007年以来的第三次军机坠毁事件。

　　香港文汇报引述现场群众介绍，军机大约在11日17时40分许掉下来，军机坠地后发出巨大的爆炸声。然后出现烟雾和火光，附近工厂忙打110和119，火很快被扑灭。据了解，由于机师提前跳伞，并没有受到爆炸的伤害。医护人员赶到现场后，对受轻伤的飞行员进行包扎。

　　事故发生后，方圆两公里都被穿蓝色海军服装的官兵拉上警戒线，一看到有车辆过来，连忙挥手倒车。军机残骸随后被军方用集装箱装走。由于事故发生的这条马路有二十多米宽，并没有对群众和财产造成损失。据悉，这是北海舰队军机首次在市区失事。目前，事故原因正调查中。

　　这已经是解放军进入2007年以来第二次军机坠毁事件，此前广州军区一架军机在广东揭阳市揭西县失事，飞机坠毁并发生爆炸，至少一名机师跳伞逃生时受伤，另一名飞行员牺牲。

　　据报道，由于涉及军方事故，现场被大批公安及武警严密封锁。据当地居民称，出事的是一架教练机，型号和失事原因不详，没有当地居民伤亡的消息。

　　此外，兰州空军郑州籍飞行员李剑英驾驶的一架歼击机于2006年11月14日发生意外：在训练中撞上飞鸟群，由于机上弹药与燃油充足，且靠近村庄，飞行员担心跳伞后飞机失控坠入村庄，最终迫降失败殉难。

　　空军某部领导李少华表示："李剑英同志是为了保护国家财产和人民生命安全，献出了年轻的生命。"2006年12月21日，空军党委给某部飞行员李剑英追记一等功，并追授他功勋飞行员金质奖章。

　　近来解放军的各种演习不断，并且多为实战状态下军事演习，尤其是解放军重视空军的发展，因此解放军空军的实战训练也空前增多，这从一定程度上增加了发生军事事故的可能。

　　军方多次的坠机事件引起高层的重视，解放军高层在2006年结束之前曾要求各部队在训练时要注重军事安全。中央军委副主席郭伯雄在沈阳战区部队调研指出，解放军职能任务的拓展、武器装备的更新、人员成分的变化等，给部队的安全工作带来许多问题，必须更加重视安全发展，确保部队安全稳定。

Technorati : 军事

十七大临近政改思潮涌动 渐进式改革方略渐清晰

随着中共"十七大"临近，中央以近期地方党委换届和十七大代表选举为契机，正推行一系列党内民主措施，为十七大政治改革的新局面奠定基础。政情观察家指出，中央积极推进民主政治，秉承渐进式改革的方略已越来越清晰。

　　香港《文汇报》报道，中央党校教授谢春涛指出，在地方党委换届中，对新领导班子人选进行全额定向民主推荐，对现任领导班子成员进行民主测评；适当扩大会议推荐范围，在会议推荐的基础上，根据实际情况，可按一定差额比例进行二次会议推荐；推荐党委委员、纪委委员人选时，要适当扩大差额比例，等等。所有这些，都显示出这次换届工作的全过程，都贯穿着中央扩大党内民主的明确要求。

　　在十七大党代表选举工作中，扩大民主的特点更为突出。例如，扩大差额选举的比例，选举单位的代表候选人一般应多于应选名额的15%，差额比例比十六大时增加了5%，扩大了选举人的选择范围。再如，这次选举工作首次要求省区市在确定代表候选人预备人选之前，要向民主党派、工商联和无党派人士通报并听取意见。

　　中央党校副校长李君如指出，中国共产党的党内民主建设正在改革和完善党内选举制度、决策机制、监督制度，以及在市、县扩大党代表大会常任制度试点这四个方面积极推进。李君如表示，中共党内民主化的程度直接影响着国家政治生活民主化的程度。

　　然而，中共中央党校最近对省部级高官进行的一项调查显示，受访的很多高官把政治体制改革的必要性排在非常靠后的位置。中共中央编译局当代所所长何增科1月10日在《2006中国法治蓝皮书》发布会上对这一调查结果表示忧虑，他认为"我们患上民主的恐惧症，只有坚定不移地进行政治体制改革，才能在反腐败的过程中掌握主动权，我们也期待着2007年能够收到更多的好消息。"

　　中共中央编译局这位专家的大胆陈词，不禁令人联想起中央编译局副局长俞可平目前广为传播的文章《民主是个好东西》。2006年10月23日，中共北京市委机关报《北京日报》刊登了俞可平的文章《关于"民主是个好东西"的辨正》，中共中央党校报纸《学习时报》随后转载该文。俞可平在此篇文章中指出，"即使有再好的衣食住行，如果没有民主的权利，人类的人格就不完整"；但同时他也提醒民主的局限，"可能要付出痛苦代价"，所以，他倡导渐进式、增量式的民主化。

　　俞可平的文章被大陆官方报章和海外媒体广为转载，引起巨大反响。在十七大日益临近的背景下，俞可平的这篇文章被海外反复解读，甚至被视为中共未来政治改革的一个风向标。

　　就这篇文章及其引发的话题，俞可平目前未做出回应，故有媒体指出，《民主是个好东西》到底是俞可平的专家之言，还是有决策层授意背景，目前无法确认。

　　不过，结合中央编译局当代所所长何增科公开引述高官调查内容，忧虑"民主恐惧症"的言论，不难看出官方体制内部学者的确是有感而发的。而专家们如此公开倡言民主，呼吁加快政治改革，自然是理论界思想论争渐入活跃期的一个表现。

Technorati : 政改