代码改变世界

[转]科研魔鬼训练术

2012-12-13 20:58 by H_OU, 241 阅读, 0 推荐, 收藏, 编辑
摘要:如果你愿意在科研上有所成就,或者说靠科研这个来生存,希望能有超过常人的能力和见识,除了天资以外,可以采用以下这三种自我训练方法,暂时说是魔鬼训练。当然作为导师也可以利用这两种方法来训练学生,据说效果很好。这些方法不是从一个途径听到的,但都有人在使用,在使用过程中,可以选择其中一种或两种,不一定全部采用。一日一文:让学生每天阅读一篇文献,无论是否与本领域相关,要求学生把文章的假说、研究方法、研究设计、研究结果和研究说明的问题等总结出来,找出文章的亮点和不足。一天一篇文献,任务看上去不重,但要把所有的内容都吃透,而且用自己的语言把文献总结出来。没有相关文献的汇总和大量阅读,几乎是不可能完成的任务。 阅读全文

[转]grads中站点处理和oacres&maskout函数应用心得

2012-12-04 11:43 by H_OU, 5668 阅读, 0 推荐, 收藏, 编辑
摘要:SJ 发表于: 2011-1-05 13:22 来源: MNMUC空间昨儿刚刚帮同事画了两张降水图(不同时间段的),很是费劲,中间碰到若干问题,但好在最后基本达到了满意的效果,下面把这次画图中遇到的问题及解决办法贴出来,供大家参考或指正:1,数据处理的问题 由于micaps的24小时降水资料中,并不是每个时次的文件里的站点个数都是相同的,就更不用提两个文件里有完全相同的站点了,这样就给各个站点降水的累加处理造成了极大的困难(如果一开始就有人或专门的机构来做这件事,将大大提高科研人员的工作效率),但本帖不打算过多描述数据处理的情形,因此这个问题就到此为止。2,站点数据中打网格和maskout的联 阅读全文

【转】Courant–Friedrichs–Lewy condition(CFL条件,柯朗-弗里德里希斯-列维条件)

2012-11-26 17:30 by H_OU, 6703 阅读, 0 推荐, 收藏, 编辑
摘要:In mathematics, the Courant–Friedrichs–Lewy condition (CFL condition) isa necessary condition for convergence while solving certain partial differential equations (usually hyperbolic PDEs) numerically by the method of finite differences.[1]虽然在Citcom解地幔对流中,每一个时间步是有限元法,但是在时间演化上,仍然是有限差分。我现在,还没见到过有限元用在时 阅读全文

(转)u检验、t检验、F检验、X2检验

2012-11-22 00:59 by H_OU, 4845 阅读, 0 推荐, 收藏, 编辑
摘要:常用显著性检验 1.t检验 适用于计量资料、正态分布、方差具有齐性的两组间小样本比较。包括配对资料间、样本与均数间、两样本均数间比较三种,三者的计算公式不能混淆。 2.t'检验 应用条件与t检验大致相同,但t′检验用于两组间方差不齐时,t′检验的计算公式实际上是方差不齐时t检验的校正公式。 3.U检验 应用条件与t检验基本一致,只是当大样本时用U检验,而小样本时则用t检验,t检验可以代替U检验。 4.方差分析 用于正态分布、方差齐性的多组间计量比较。常见的有单因素分组的多样本均数比较及双因素分组的多个样本均数的比较,方差分析首先是比较各组间总的差异,如总差异有显著性,再进行组间... 阅读全文

contains

2012-11-20 11:16 by H_OU, 296 阅读, 0 推荐, 收藏, 编辑
摘要:prigram main...... ......call p1()call p2().......... containssubroutine p1() ............end subroutine p1subroutine p2() ............end subroutine p2end! 子程序中的变量无需重新说明,变量值无需通过哑实结合传递。elemental代表的是你可以传数组给这个函数,而这个函数会一个一个地处理这个数组的,并且是相互独立的(可并行)external代表这后面声名的子程序是在别处定义的,不属于当前子程序的内部子程序contains代表下... 阅读全文

[转]在linux下如何使用Makefile对fortran程序进行编译

2012-11-14 20:38 by H_OU, 6473 阅读, 0 推荐, 收藏, 编辑
摘要:-在linux下写过不少的fortran程序,但很少写makefile文件,通常情况下我是手动用ifort或pgf90进行编译,虽然也在大型程序中写过一些代码,但因为都有现成的makefile文件,我也只是修改一下配置参数,先make clean,再make,只要生成可执行程序就算搞定了。今天心血来潮,自己写了几个Makefile,供大家参考。一、不用Makefile的直接编译比如有2个简单程序,主程序为main.F90,子程序为sub.F90. 如果直接编译的话,键入ifort main.F90 sub.F90就可以生成一个a.out文件,即为可执行程序了。你也可以指定一个可执行程序的名字, 阅读全文

【转】一些解决变态数学公式的算法地址

2012-11-14 20:26 by H_OU, 993 阅读, 0 推荐, 收藏, 编辑
摘要:一些解决变态数学公式的算法网址2011-09-26 09:48:16|分类: 默认分类 |标签:矩阵分析 |字号大中小订阅http://www.acm.org/calgo/contents/各种算法的程序。http://www.kinjo.ac.jp/~kawabe/qc/mig.htmlMIG,分子积分生成程序。直接下载,可运行于Windows/UNIX。 http://forcal.go.nease.net/Forcal数学软件,用于数学计算的库函数。http://www.cnvf.com/Fortran 博物馆。http://www.matlab-world.com/MATLAB大观园。 阅读全文

【转】fortran 广义逆 子程序

2012-11-14 16:56 by H_OU, 332 阅读, 0 推荐, 收藏, 编辑
摘要:SUBROUTINE BGINV(M,N,A,AA,L,EPS,U,V,KA,S,E,WORK)DIMENSION A(M,N),U(M,M),V(N,N),AA(N,M)DIMENSION S(KA),E(KA),WORK(KA)DOUBLE PRECISION A,U,V,AA,S,E,WOEKCALL BMUAV(A,M,N,U,V,L,EPS,KA,S,E,WORK) !一般实矩阵奇异值分解IF (L.EQ.0) THENK=110 IF(A(K,K).NE.0.0) THENK=K+1IF (K.LE.MIN(M,N)) GOTO 10ENDIFK=K-1IF(K.NE.0) THE 阅读全文

【转】SVD分解

2012-11-14 11:13 by H_OU, 495 阅读, 0 推荐, 收藏, 编辑
摘要:分类的关键是计算相关性。我们首先对两个文本计算出它们的内容词,或者说实词的向量,然后求这两个向量的夹角。当这两个向量夹角为零时,新闻就相关;当它们垂直或者说正交时,新闻则无关。当然,夹角的余弦等同于向量的内积。从理论上讲,这种算法非常好。但是计算时间特别长。通常,我们要处理的文章的数量都很大,至少在百万篇以上,二次回标有非常长,比如说有五十万个词(包括人名地名产品名称等等)。如果想通过对一百万篇文章两篇两篇地成对比较,来找出所有共同主题的文章,就要比较五千亿对文章。现在的计算机一秒钟最多可以比较一千对文章,完成这一百万篇文章相关性比较就需要十五年时间。注意,要真正完成文章的分类还要反复重复上述 阅读全文

【转】grads遇到的一些问题

2012-11-07 15:43 by H_OU, 2415 阅读, 0 推荐, 收藏, 编辑
摘要:1.【转】 grads中同时使用多个地图叠加的方法使用GrADS作图时,有时候会用到不同分辨率的地图叠加在一起,比如说,世界地图+省地图+区县地图,以前笨笨总是多次设置,然后画地图,今天偶然发现可以一次设置多个地图,然后一起画的。'set mpdset cnworld jiangsu jiangsu_q''draw map'可以将地图放在grads-1.9b4/data目录下2.Error opening stroke character data set 错误需要加两个系统变量 GADDIR=d:\pcgrads\dat GASCRP=d:\pcgrads\l 阅读全文