发表于
2012-10-31 00:00:00
调试代码(Python)
调试代码翻译自：http://scipy-lectures.github.com/advanced/debugging/index.html作者：Gaël Varoquaux本教程探索工具，以更好地理解你的编程基础：调试，找到并修复错误。这并不特定适合Python科学计算社区，但是我们采取的策略量身定制。先决条件NumpyIPythonnosetests (http://readthedocs.org/docs/nose/en/latest/)line_profiler (http://packages.python.org/line_profiler/)pyflakes (http://pypi.python.org/pypi/pyflakes)gdb 用于C调试部分目录避免错误避免麻烦最佳编码法pyflakes:快速静态分析在当前编辑文件运行pyflakes一个在线拼写检查集成调试工作流使用Python调试器调用调试器尸检单步执行其它开始调试的方法调试器命令和交互使用gdb调试段错误总结练习Footnotes避免错误避免麻烦最佳编码法我们都写有bug的代码。接受它，解决它。写程序时，在心里测试和调试。Keep it Simple, Stupid(KISS).可能有效的最简单方法Don‘t Repeat Yourself (DRY).每个知识都必须有个单独、明确、在一个系统中的正
...继续阅读
(15)
发表于
2012-10-30 00:00:00
斯普罗站最大风力预测
发表于
2012-10-22 00:00:00
Scipy:高端科学计算
发表于
2012-10-19 00:00:00
优化代码(Python)
优化代码翻译自：http://scipy-lectures.github.com/advanced/optimizing/index.html作者:Gaël VaroquauxLicense:Creative Commons Attribution 3.0 United States License (CC-by) http://creativecommons.org/licenses/by/3.0/us过早的优化是罪恶的根源。——Donald. Knuth这个章节涉及使Python代码运行更快的策略。先决条件line_profiler(http://packages.python.org/line_profiler/)目录优化工作流剖析(profile)Python代码Timeit分析器(Profiler)Line-profiler让代码更快算法优化SVD的示例写更快的数值代码优化工作流让它工作：以简单_清晰_的方式书写代码。让它可靠的动作：书写自动化的测试实例，确认你的算法是正确的。如果你中止它，测试将捕捉到中断。优化代码：通过剖析(profile)简单的用例来发现瓶颈，并且加速这些瓶颈，找到更好的算法或实现。记住在剖析一个现实的实例和代码的简洁与执行速度之间权衡。对于有效率的工作，最好让剖析运行约十秒。剖析(profile)Python代码非测量，不优化测量：剖析(profi
...继续阅读
(27)
发表于
2012-10-14 00:00:00
Scikit Learn: 在python中机器学习
Scikit Learn: 在python中机器学习Warning警告：渣翻译，很多我自己都没看懂.但我会慢慢看慢慢修改，此文档维护中……翻译自：Scikit Learn:Machine Learning in Python作者: Fabian Pedregosa, Gael Varoquaux先决条件Numpy, ScipyIPythonmatplotlibscikit-learn目录载入示例数据一个改变数据集大小的示例：数码数据集(digits datasets)学习和预测分类K最近邻(KNN)分类器训练集和测试集分类支持向量机(SVMs)线性支持向量机使用核聚类：将观测值聚合k均值聚类应用到图像压缩用主成分分析降维将一切放在一起：人脸识别线性模型：从回归到稀疏稀疏模型同一问题的不同算法模型选择：选择估计器和它们的参数格点搜索和交叉验证估计器格点搜索交叉验证估计器Footnotes警告：在0.9版中(2011年9月发行)，scikit-learn的导入路径从scikits.learn更改为sklearn载入示例数据首先我们载入一些用来玩耍的数据。我们将使用的数据是非常简单的著名的花朵数据——安德森鸢尾花卉数据集。我们有一百五十个鸢尾花的一些尺寸的观测值：萼片长度、宽度，花瓣长度和宽度。还有它们的亚属：山鸢尾（Iris setosa）、变色鸢尾（Iris versicolor）和
...继续阅读
(39)
发表于
2012-10-14 00:00:00
Always Look on the Bright Side of Life
blahblahblah,这几天看了三部电影，随便感慨下。Always Look on the Bright Side of Life也许执念太深，从一年前开始看到无所事事的人就极度不爽，也无法忍受自己浪费一分一秒，然而事情总是这样：你越想抓住什么，什么就越从你的指缝间溜走得越快。一年来总是给自己一些莫名奇妙的压力，天天为各种各样的琐事烦得一塌糊涂。结果身体一直也不好，要不要感冒每个月就像大姨妈一样来一次？1有很多人和事，只觉得他们二逼，然而搅扰的自己也无法静下心。也许只是我自己太在意。我想黑我们学校，黑自己专业，黑这学校很多人很久了。然而最终还是忍住，也许世界就该这样，或者你爱怎样就怎样。没必要因此搅扰得自己不安。也许我就该笑笑，依旧回头走我的路，如果路被人堵上，换另一条继续走，条条大路通罗马。也许我就该放下许多执念，在心灵上解放自我，在行动上自我解放。也许我就该随波逐流，无所谓得失，让该来的来，让该去的去，万事随缘。每天入夜，当我怨恨时，我不断暗示自己，神赐你一颗平和的心，不要怨恨。当我为周围的风吹草动胆颤心惊时，不断鼓舞自己，神赐给你的是一个刚强的心，无所畏惧。让心平和下来吧，一直想考研，跨掉自己的专业，离开熟悉的学校，寻找更好的环境和人2。但没必要逼自己，让心平和下来吧好好锻炼身体，考不上明年再来。再考不上再做打算。留得青山在，不怕没柴烧。总是看见生活光明的一面。some
...继续阅读
(16)
发表于
2012-09-20 00:00:00
创新项目要结束了
发表于
2012-09-12 00:00:00
科学计算工具和工作流
Scientific computing with tools and workflow好长的标题……从scipy lecture notes引言翻译。这是篇python 做科学计算的水文……作者:Fernando Perez, Emmanuelle Gouillart, Gaël Varoquaux。感谢他们贡献了如此棒的教材。目录为什么Python？科学家的需求要求已经存在的解决方案python科学工具交互工作流:IPython和文本编辑器命令行交互在编辑器中拟定算法为什么Python？科学家的需求获取数据(模拟，试验操控)操作和处理数据可视化结果……让我们理解我们在做什么！交流结果：产生用来做报告或出版的图像，写演示文稿。要求对相应经典数值方法和基本作用已经存在丰富的程序块：我们不想再次为绘制曲线、傅利叶变换或者拟合算法重新编程。不要重复发明轮子。简单易学：计算机科学不是我们的课程也不是工作。我们想要在几分钟内绘制曲线、平滑信号和做傅利叶变换。易于和合作者、学生、顾客交流，让代码在试验室或公司内部都生动起来：代码应该像本书一样易读。因此，语言应该包含尽量少的语法符号或不必要的规范，使数学家和科学家读者能够理解代码。有效率地编码和快速执行……但是不必说如果我们花太多时间来编写它，非常快的代码也变得无用。所以，我们需要同样短的开发时间和短的执行时间。一个单独的环境/语言解决所有事
...继续阅读
(50)
发表于
2012-09-07 00:00:00
Matplotlib Tutorial(译)
发表于
2012-08-30 00:00:00
SymPy Tutorial(译)
发表于
2012-08-24 00:00:00
SciPy入门指南(译)
发表于
2012-08-22 00:00:00
试验性NumPy教程(译)
试验性的Numpy教程原文来自Tentative NumPy Tutorial目录先决条件基础篇一个例子创建数组打印数组基本运算通用函数(ufunc)索引，切片和迭代形状操作更改数组的形状组合(stack)不同的数组将一个数组分割(split)成几个小数组复制和视图完全不拷贝视图(view)和浅复制深复制函数和方法(method)总览进阶广播法则(rule)花哨的索引和索引技巧通过数组索引通过布尔数组索引ix_()函数用字符串索引线性代数简单数组运算矩阵类索引：比较矩阵和二维数组技巧和提示“自动”改变形状向量组合(stacking)直方图(histogram)参考文献先决条件在阅读这个教程之前，你多少需要知道点python。如果你想重新回忆下，请看看Python Tutorial.如果你想要运行教程中的示例，你至少需要在你的电脑上安装了以下一些软件:PythonNumPy这些是可能对你有帮助的:ipython是一个净强化的交互Python Shell，对探索NumPy的特性非常方便。matplotlib将允许你绘图Scipy在NumPy的基础上提供了很多科学模块基础篇NumPy的主要对象是同种元素的多维数组。这是一个所有的元素都是一种类型、通过一个正整数元组索引的元素表格(通常是元素是数字)。在NumPy中维度(dimensions)叫做轴(axes)，轴的个数叫做秩(rank)。
...继续阅读
(25)
发表于
2012-08-17 00:00:00
To Be An Epicurean
胡言乱语，请随意喷。起晚上我正对着一堆python的东西不知所措时，室友自习回来，怒斥着考研太难过，表示不想考了。瞬间我们寝室便从四人考研将要沦落到我一个人了，这不是什么问题，也就是将来一个人困难更多些而已。然后他补充了一句：人活了为了什么？不就是为了爽吗？我忽然感觉语塞，不知如何应答。说不出有什么问题，但心中拼命有个声音在问着：是吗是吗是这样吗？感觉不对劲，很不对劲。我不是一个反享乐主义者，事实上，我觉得每一个人都应该尽力享受他们的生命。若干月前我在图书馆看到一本享乐主义手册，书中虽然罗列了史上臭名昭著的享乐主义者，但极力鼓励人们享受自己的生活，过一个轻松愉快的人生。但总感觉有哪里不对劲，但又理不清究竟是哪里不对。承于是心神不宁地去操场跑步，仰头看着久违的美丽星空。矛盾在心中一直萦绕，这无关紧要的问题为什么让我如此心神不宁？我也许不是个纯粹的享乐主义者，甚至在很多方面我都倾向于苦修主义者，希望借以磨砺自己的意志。但我只是跟他人方向不同罢了，互联网和文字使那些与自我时空千里之遥的思想在我心之土壤生根发芽，比周边的人更多塑造了我。我只是把自己的执迷不悟和强烈的欲望发泄在了我因此所感兴趣的领域上——数学、艺术、计算机。但我依然在他人不关心的方面贪嗔痴怪，不觉不慧。我寻找、追寻、计划、妄想，依然也是使我感觉更“爽”，像在技术上的成就感、hack的乐趣，艺术的美和数学的美和精致。既然我亦如
...继续阅读
(26)
发表于
2012-08-15 00:00:00
GALGAME
发表于
2012-08-11 00:00:00
随笔——关于学习
发表于
2012-08-04 00:00:00
这只是篇随笔
2012年08月04日随笔其实是好久没在这写东西了。^o^无锡的夜晚上去跑步时，JNU的夜shocked me.当站在宽阔的操场上，凝视着安详高远的天空，宁静平和的初升明月,南边一溜美丽祥和的云朵,东边若隐若现的城市灯火，黑暗中飒飒作响的树影，西边绵延而沉默的山峦，头顶点点的星辰。夏夜的晚风穿过身旁，撩拨着我的衣裳，也撩拨着我的神经。好像那片云云千奇百怪，诡奇壮丽，它们是天空最亮丽的风景线。它洁白如雪过，它压抑阴沉的想要坠落过,它曾排山倒海压来，也曾宁静的好像不变的山峦。然而它最终要散，无影无踪。就像我所追求的，我所喜爱的，我所经历的，我所拥有的。在生命经历一个个奇幻的形态后烟消云散。曾经的美丽只留存在他人的记忆之中。跑完步时，南边的云散得无影无踪，好像从未存在。linuxer周年记入门是件很奇特的事，我视之为缘分,思之为注定。一切在阴差阳错间成就，也许像从前曾有缘过的一样缘尽于某日然而就这么改变了我一年的人生之前生命中点点阴差阳错点点成线，汇集成今那个轨迹，就是人们所说的命运？是否在他人眼里是彻头彻尾的宅loser？是否被看成不通世情的sb？It dosen’t matter.我喜欢开源社区的氛围，开放、共享、自由与进步。我认同黑客的精神，自强、求索、执着、包容。我相信人生有很多值得去体验和经历的事，有很多有趣的东西：哲学、文学、艺术、数学、计算机、还有风景，还有girl……我相
...继续阅读
(14)
发表于
2012-07-13 00:00:00
回忆
回忆我向来厌恶回忆，当被它的触角在夜幕中拉住，你永远别指望从中抽身，回忆越陷越深，一层一层向过去延伸，好像永远没有终点。我向来厌恶面对过去，过去的哀伤与委屈，过去的快乐与迷恋，都在回忆中幻化成不可抑制的情愫与遗憾。回忆越陷越深，它在你的心头缠绕了一圈又一圈，当回忆越行越远，往事便像猴子头上的紧箍咒一样越来越紧，让人压抑的喘不过气。不知为什么，感觉往事不堪回首。虽然它普通的不能再普通，虽然它平淡的不能再平淡，但我不愿回头.让记忆随着光阴越流越远吧，我希望永远不要再见到它。可是它总是找上我。终有一天，还是要面对自己的懦弱，愚蠢与自私。有一种被人抛诸脑后的过去，包含一些不该遗忘的点点滴滴。欢笑之余，我将故旧拾起，才刹见自己的薄情。人世间似乎存在着某种公平，每个人都曾被抛掷在他人的遗忘中；而今你所有哽咽难宣的话语，也正是你曾为人所期盼，却永远听不见的字句。赫曼· 赫塞—-大学之前不想多说，封尘的不想重启。简单的写下，简单地卸下，轻装面对未来。大学所谓大学者,非谓有大楼之谓也,有大师之谓也。梅贻琦1931年在清华的就职演说—-说什么好呢，本来没想到来这里。高考之前想去南京大学，考完试知道毁了，打点好好准备了一堆要报的二本大学，一本随手填了个江南大学，只因为作为革命老区的朴素人民，觉得无锡这地方比较发达。报专业的时候报的教育技术，因为自己虽然是学理科的，但一直想搞文科1。但老爸是搞教育的，想这
...继续阅读
(21)
发表于
2012-07-11 00:00:00
这两天的二逼日子
发表于
2012-07-09 00:00:00
land of lisp读书笔记
Happy Hacking我第一次读这本书，便被其浓郁的黑客风格与幽默所吸引，它向我展示了common lisp这门语言强大与趣味的一面。之前从来没有试着度过原版图书，但看到这么可爱一本书，竟然吸引我一边用着取词软件一边用着google翻译给看完了。说是看完了，其实有很多地方都是略略而过，与其说我是学习了一遍common lisp，不如说是大致了解了一下它。很多东西像AI，程序设计，算法的东西都是云里雾里，好在作者假设读者什么也不知道，详细而幽默的解说引人入胜，看起来虽然很累但趣味性十足。所以，即使是lisp菜鸟，这本书也是很推荐。1说说我学到的：学到了很多优秀而强大的common lisp技术，在有些地方还进行了自己的hack，甚至学到了服务器、http协议、svg和html、博弈算法的一些东西。也说说从这本书中衍生的新想法：lisp做web应用还是蛮好的，正好弥补lisp缺乏图形界面库的短版，common lisp强大但复杂，也许scheme才是我心仪的语言2语言障碍仍然是个很大的问题，要好好学习英语这门世界性语言。将我4月26日至7月7日的读书笔记集合在这里,虽然水平很烂但斗胆放在这里，希望对他人有帮助land of lisp读书笔记忽然发现这个还在草稿里……2012年09月22日 星期六 13时26分44秒Footnotes好像此书就是为了初学者设计的……↩对一
...继续阅读
(16)
发表于
2012-07-07 00:00:00
Dice of Doom v4
最终版的Dice of Doom这个游戏终于结束了，基本上这本书的主要内容也算看完了。实际上，我根本就写不出这些东西，源码也读的云里雾里。但这是一本很有趣的书，我想看完htdp之后用htdp推荐风格重新处理下这本书的代码，顺便也复习一下。让Dice of Doom 更加的有趣我们上回完成的版本虽然很棒，但总缺乏些什么。因为之前程序的设计，我们可以很轻易的扩展它，让它变得更有趣先加载我们之前的游戏1(load"dice-v3");增加玩家数因为我们之前的架构，增加玩家数很简单。为了增加玩家数后计算机能反应敏捷，让它变“傻”一点。(defparameter*num-players*4)(defparameter*die-colors*'((2556363)(6363255)(6325563)(25563255)))(defparameter*max-dice*5)(defparameter*ai-level*2)建立概率节点我们本质上是向游戏树中加入几率节点，先升级attacking-moves函数，增加个board-attack-fail分支。(defunattacking-moves(boardcur-playerspare-dice)(labels((player(pos)(car(arefboardpos)))(dice(pos)(cadr(arefboardpos))))(l
...继续阅读
(19)
发表于
2012-07-05 00:00:00
概率论与数理统计基础
发表于
2012-07-02 00:00:00
基于web的Dice of Doom
发表于
2012-06-29 00:00:00
Dice of Doom v2
Dice of Doom v2我已经不怎么看得懂在讲什么了，所以胡言乱语请见谅,有几个问题:记不起之前的游戏树什么样了，我又懒得回头去看。是不是某种决策树？不停的函数接着函数扩展看起来很方便，但又要求你时刻记得之前实现的细节，以实现正确的改进作者的编程风格，以后还是按htdp里说的玩去吧，起码不会像现在这样晕。英语还是障碍啊，有时候看完一大段文字后不知道在讲什么，囧了。Dice of Doom v2首先加载之前我们完成的Dice of Doom源码和为惰性求值所实现的一系列库文件。没有的话就从land of lisp网站上下载吧，话说这本书真贵要40刀左右，豆瓣上竟然380元人民币……(load"dice_of_doom_v1.lisp")(load"lazy.lisp")然后先把游戏board扩大到4x4。(defparameter*board-size*4)(defparameter*board-hexnum*(**board-size**board-size*))下面把Dice of Doom游戏全面lazy化吧(defunadd-passing-move(boardplayerspare-dicefirst-movemoves)(iffirst-movemoves(lazy-cons(listnil(game-tree(add-new-diceboardplayer(1-s
...继续阅读
(14)
发表于
2012-06-28 00:00:00
Google glasses
最近耳闻google I/O大会，看到个消息：开发版 Google 眼镜开始接受限量预定，1500 刀超高价格明年年初发货然后随便上网查查，奇怪为啥这么多人吹捧。先看看google的宣传片吧，有种未来触手可及的感觉。下一步就该把眼镜和人的视觉系统整合了吧，攻壳机动队全义体化的假想要实现了么？还有个国内吐槽版：真欢乐……
...继续阅读
(21)
发表于
2012-06-26 00:00:00
Common lisp菜鸟指南(译)
Common Lisp菜鸟指南翻译自:http://ghostopera.org/blog/2012/06/24/the-newbie-guide-to-common-lisp/进入Common Lisp的世界可能非常让人退缩，甚至对有经验的程序猿也是这样。这个语言兼有力与美，但也有许多清晰的边界和没有真正被普遍理解的开端。这个指南可能将持续一段时间，但是希望它将使踏入Common Lisp的任务变得不那么可怕。理解Common LispAriel Networks有个适当形式的 Common Lisp 指南。TheHyperSpec有完整的 Common Lisp 文档。CLQR是一个好的课下载的 Common Lisp 快速参考手册。了解你的工具Emacs1Emacs是一个可实用lisp高度定制的编辑器2。它可以为lisp提供一个非常独特的实时开发流程。查看David O’Toole在通过Emacs和他的 “Blocky.io” 游戏开发系统在一个游戏中的lightning talk，来查看实时开发实例。很酷不是吗？Steel Bank Common Lisp当这世界上有很多可用的 Common Lisp 实现时，"SBCL":sbcl 是非常标准的实现。它很积极的被开发着，有适合的调试器、原生线程、多平台支持。QuicklispQuickl
...继续阅读
(24)
发表于
2012-06-23 00:00:00
宽容
发表于
2012-06-22 00:00:00
A Mathjax Test
TeX Samples以下公式以HTML源代码呈现LaTeX表达式。洛伦兹公式$$\begin{aligned}\dot{x} & = \sigma(y-x) \\\dot{y} & = \rho x – y – xz \\\dot{z} & = -\beta z + xy\end{aligned}$$柯西-施瓦茨不等式$$\left( \sum_{k=1}^n a_k b_k \right)^2 \leq \left( \sum_{k=1}^n a_k^2 \right) \left( \sum_{k=1}^n b_k^2 \right)$$叉积公式$$\mathbf{V}_1 \times \mathbf{V}_2 = \begin{vmatrix}\mathbf{i} & \mathbf{j} & \mathbf{k} \\\frac{\partial X}{\partial u} & \frac{\partial Y}{\partial u} & 0 \\\frac{\partial X}{\partial v} & \frac{\partial Y}{\partial v} & 0\end{vmatrix}$$n次实验k次成功的概率$$P(E) = {n \choose k} p^k (1-p)^{ n-k}$$拉马努金恒等式$$\frac{
...继续阅读
(52)
发表于
2012-06-21 00:00:00
罗列下linux的好处
是什么把你留在了linux下不懂者请自行google，不喜者请自行绕道和windows和mac OS相比，linux桌面可谓默默无闻，小众的一塌糊涂。经常能在各种相关论坛看到各种抱怨。驱动、flash、游戏、不稳定、不用户友好……千宗之罪，不可胜数。不过既然linux桌面这么“不堪”，是什么把你留在了这里？它是黑客友好的，透明性无与伦比。它是靠无数有热情的人们用智慧和激情维系的世界，它的价值观是开放、自由、进步和共享的。它追求平等，赋予每个人以平等的机会和权利1。它不会莫名弹窗，不会各种流氓。事实上，在firefox adblock的护佑下我都感觉不到网络世界广告的存在2。这里是天堂它适于DIY。你要愿意，可以搞得天花乱坠秒杀苹果，可以搞得特立独行到让人认不出这是个人电脑，也可以弄得朴素平实得像win98，也可以搞得非常大众3。有兴趣看截图的请自行google4。各种配置都可以自己DIY或选择发行版轻快的使用。不满意还可以自己手动hack，当然只要你有能力和精力……懒人的包管理器。曾经linux程序的依赖关系完全可以把一个正常人逼疯，时过境迁，现在回到windows装软件人完全会疯，虽然windows有那啥360来着5。百毒不侵。可怜的windows……不知道mac怎么样。稳定。稳定在于你想稳定，天天像作者一样更新可不好。想稳定，用debian、用bsd。服务器系统的稳定不是盖的。
...继续阅读
(10)
发表于
2012-06-19 00:00:00
惰性求值
惰性求值啥叫惰性求值具体的含义可以参看这里惰性求值。我的看法就是一种延迟计算方法，好像一个懒汉，只计算需要计算的东西。为什么要惰性求值因为惰性求值能极大优化程序的性能，想想那个Dice of Doom的游戏，如果不需要每次都生成不需要的游戏树，将会节省多少时间和电力。同时，惰性求值能做到一些几乎不可能的事，构造无穷的东西。在common lisp中实现一整套惰性求值函数库据作者说，像Haskell和Clojure这样的预言直接把惰性求值作为语言默认的一部分，common lisp则没有这样，不过实现一整套惰性求值库也没有什么难度延迟求值我们想要得到这样的东西：lazy和force.lazy只返回函数。(lazy(+12))#(force(lazy(+12)))3lazy可以用宏这么实现(defmacrolazy(&body;body)(let((forced(gensym))(value(gensym)))`(let((,forcednil)(,valuenil))(lambda()(unless,forced;unless，when后是并行结构(setf,value(progn,@body))(setf,forcedt)),value))))force可以直接用简单的函数实现(defunforce(lazy-value)(funcalllazy-value))但我们需要的不只这些
...继续阅读
(19)
发表于
2012-06-16 00:00:00
传说中的DSL
lisp中的dsl什么是dsl，为什么要dsl显然对于特定的领域，特定的程序是更好的解决方案，这个程序就是DSL1.lisp非常方便来实现DSL。在land of lisp一书中进行了两种dsl，一种是生成xml或html，另一种是对我们之前wizard游戏的中的游戏命令进行dsl。生成svg文件关于svg文件不加以详述，读者可以参考这里。这非常重要，只有理解了xml文件的格式才能继续我们的dsl。svg采用xml格式，顺便我们也说说生成网页的dsl。首先写一个辅助函数，;Creating XML and HTML with the tag Macro;Writing a Macro Helper Function(defunprint-tag(namealstclosingp)(princ#\<)(whenclosingp(princ#\/))(princ(string-downcasename))(mapc(lambda(att);使用mapc是为了副作用。它总返回第一个参数(formatt" ~a=\"~a\""(string-downcase(caratt))(cdratt)))alst)(princ#\>))可以试试(print-tag'mytag'((color.blue)(height.9))nil)然后创造生成标签的宏，宏能在以下几个方面提供不可替代的改善：标签总是
...继续阅读
(15)