科学网

 找回密码
  注册

tag 标签: 实例

相关帖子

版块 作者 回复/查看 最后发表

没有相关内容

相关日志

房价报道:线性思维举例
热度 1 Mech 2016-4-11 11:37
媒体上不时看到房价太高的报道。太高的依据就是所谓不吃不喝要多少年才能买多大的房。房价是不是太高我不知道,但这种推理方式是典型的线性思维。 这类房价报道中,最早的可能是网上看到 1989 年 2 月 10 日 《人民日报》第二版一个截图。标题的四六句说得过去,但内容真是不行。图片如下,全文附在后面。 其中所谓“一名大学生从参加工作起就日日缩衣节食,每月存储 50 元,已是最高极限,需 100 年才能买上两居室。” 这就是典型的线性思维,我在通识课上,常以此举例,虽然在校学生可能并不太关心房价。媒体仍不断出现类似言论。有兴趣的学生,不难收集这类言论,来课堂上分析。 这种分析与房价是否合理没有关系。只是说这种论证方式基于线性思维。众所周知,论证的错误,在逻辑上并不意味着结论的错误。线性思维的预期与实际不符。 当年的大学毕业生不知道现在有多少人买了房,才过去 27 年,远没有到 100 年。这到让我想起中学古文课文中的话, 天下事有难易乎?为之,则难者亦易矣;不为,则易者亦难矣。 后面的类比也适合对年轻人讲, 人之为学有难易乎?学之,则难者亦易矣;不学,则易者亦难矣。 至于说要立法解决问题,那就是另一种线性思维。如果学生在我后面的课中听到“日凿一窍,七日浑沌死。”的寓言时,能联想到这个立法的说法。那就有些“通”的意思了。 附:报道全文 房价猛涨 百姓望楼兴叹 势在必控 国家正拟法规 本报讯  据今天的《市场报》报道:记者最近从建设部城市住宅局了解到,针对目前各地住宅市场出现的价高、滞销等问题,建设部、国家物价局、建设银行正在会同制定《商品房价管理暂行办法》、《关于住宅市场的若干规定》等6个法规,将明确规定经营者虚报成本、图谋厚利属非法行为,商品房价混乱猛涨等现象可望依法纠正。 据最新统计,全国储蓄存款加上私人口袋里的现钞,总数达 5500 亿元。如何把这些消费资金变为生产基金,拿出大量的商品住宅组织货币回笼是一个必要办法。 但是,买房对普通干部来说,有人打比喻,犹如从夜空里摘颗星星。今年 1 月,上海住宅市场年初出售 11000 多套住房,成交不到 1/10 ,幢幢楼房无人买,几十万平方米新住宅空空荡荡。原因是价格昂贵,每平方米最高价已达 2300 元,令人咋舌。北京最近提供 2 万多平方米住房,每平方米为 1600 元 —1900 元。若买两居室,少说也要 6 万多元。一名大学生从参加工作起就日日节衣缩食,每月存储 50 元,已是最高极限,需 100 年才能买上两居室。 商品房价格猛涨,追寻原因,一是经营房地产被看作是最肥的行当。几年来,各类房产公司蜂拥而上,全国陡然间出现 2700 多家。当然,不少房产公司为发展商品房建设作出了贡献,功不可没,但也有贪图享利、干坑骗消费者勾当的,遍地拉线,到处拆建,违法违纪严重。二是制订价格构成缺乏相应的政策法规,各自为政,极不统一,有的规定为 10 项,有的则十几项。基本造价之外,还有居委会、派出所、商店、中小学、幼儿园、公园、卫生所、道路等等的配套设施费用竟相当于基本造价的两倍甚至 3 倍多。今后商品房的成本计算、价格、利润将受到法律约束。
个人分类: 教学行思|3471 次阅读|2 个评论
NP是可计算的吗?- “算法”的二个层次
热度 3 liuyu2205 2015-12-30 17:26
于诸博文(注)我们解读流行观念“NP是可计算,但是难计算”,认为存在着认知错误,其根源在于人们未深究“算法”的本质(可计算性)。实际上,“算法”这一概念涉及到二个不同的层次:实例和问题,人们混淆了这二个层次,导致对“算法”概念的模糊。 这里,我们用下面的图帮助说明“算法”涉及到的二个层次: 图灵机 算法 ------- 实例 人 算法 ------- 问题 面对“实例”,对“图灵机”而言,比如说可通过基于枚举的“穷举法”计算实例,即涉及“机器”层次;面对“问题”,因问题是实例的抽象,对“人”而言,能否将计算实例的算法推广到计算问题中的任何一个实例,需考虑该算法的“复杂度”、“可计算性”,即涉及“人”的层次。 对于一个问题,只要人能举出实例,机器总可以计算该实例,比如用“穷举法”;真正的问题在于,NP的搜索空间是指数型的,人根本无法用枚举法举出指数增长规模的空间实例,因此机器也就无法计算任何一个实例了,说“NP问题实例是难计算”,实际上是“NP问题难到不可计算”,所以NP才是“不确定性问题(Nondeterministic Problem) ”,。。。 岁末感言: 当初 开博客的心愿就是来与大家分享NP理论研究的心路历程、同参共学的, 如今博客满周岁, 一篇篇博文见证了与大家的互动和共同的进步,感谢在后面支持的师友!感谢在这里相遇的网友!。。。在新的一年里,我们将把“不确定性问题(NP)”的讨论延伸到英语博客,与国际同行一起讨论,希望大家同行! 辞旧迎新之际,祝大家新年愉快!
个人分类: 不确定性问题和算法讨论|3596 次阅读|6 个评论
3-SAT分段子句消去法求唯一解例题
热度 1 accsys 2015-10-6 08:49
3-SAT 分段子句消去法求唯一解例题 姜咏江 利用子句块唯一解和动态块唯一解来处理 3-SAT 求解,最方便的是 3-CNF=1 有唯一解的情况。这也是用其他方法求解 3-SAT 问题最不易的事情。 在用关联段子句消去的过程中,如果出现动态块消去呈现多解的情况,如果选择的变量值不合适,有可能出现后续无解的情况。正确的做法是暂不设值,而是沿着一元关联的子句块寻找唯一解的剩余子句块。只要找到关联的满足唯一解的子句块,就可以先设定这个子句块的变量值,然后反方向去确定前面多解的选值,这样可以最终正确地求出 3-SAT 的解。如果最终都是多解的子句块,那么这个3-SAT的问题肯定是多解了,此时确定最后一个子句块解,反向可得解。 3-SAT 有唯一解的情况,用分段子句消去法既方便又准确。有兴趣的朋友不妨试试。我还没有编写程序处理。不过手工操作弄个百八十个变量也不是难事。就是上千个变量也未必是难事,操作过程中随着子句的消去,工作量会迅速减少,时间自然会加速减少。 本文给朋友一个唯一解的例子,希望有心人能够自己做做看,体会一下求解3-SAT问题的关键点在何处。表下面一行是实际消去子句后的结果。 唯一解实例 : 3-SAT分段子句消去法唯一解例题.xls 参考: http://blog.sciencenet.cn/blog-340399-925143.html 2015-10-6
个人分类: 机器计算|2648 次阅读|2 个评论
[转载][《Flash8完全自学手册》][图形图像][PDF/107M]
lcj2212916 2014-12-7 13:47
这是一本教用户如何创建Flash动画的资料。本书并不是单纯地介绍Flash技术层面的东西,旨在提供以某种技术为核心的整个开发思路和流程。本书采用鼠标加大脑的学习模式,通过精心设计的实例深入细致地介绍了Flash 8完整的基础动作指令及进阶技巧。通过书中实用动画实例的学习,可以在较短时间内学习和掌握如何制作精彩的动画及制作技巧。 本书使用大量实例对中文版Flash Professional 8的常用功能和使用方法进行了详细讲解,特别是把Flash 8的新功能融入到实际的动画制作过程,另外对ActionScript进行了深入介绍,引导读者学会使用ActionScript编程制作复杂的动画。内容主要包括FlaSh的基本操作、常用工具的使用、各种Flash动画的创作方法、组件的使用和创建、ActionScript语言等。本书适合从事动画制作及网页编程的用户阅读,也作为各类大中专院校及相关职业培训学校的教材。 下载地址: http://www.400gb.com/file/79953447
806 次阅读|0 个评论
GAMIT-GLOBK处理实例——从零开始,到解算结束
热度 2 POPOPOLAR 2014-8-28 09:55
确认 GAMIT/GLOBK 软件安装成功之后,为了能正确处理数据和绘图,还需要做两件事情。第一,更新表文件。在 ftp 服务器上下载最新的 tables 表文件( ftp://garner.ucsd.edu/pub/gamit/tables/ ),然后复制到安装目录下的 tables 中,遇到相同文件选择覆盖即可;第二,安装并配置好绘图工具 GMT ( The Generic Mapping Tools )。 友情链接: Linux/Ubuntu 环境下轻松安装 GAMIT/GLOBK—— http://blog.sciencenet.cn/blog-1217335-814640.html Linux/Ubuntu 环境下轻松安装配置 GMT—— http://blog.sciencenet.cn/blog-1217335-758694.html 下面的实例中,处理的数据为 2010 年第 56 到 60 共 5 天的 3 个 IGS 跟踪站( cas1 、 dav1 和 mac1 )和 2 个待求站( grw1 、 zhn1 )的数据。 一、处理前的准备 1 、在主文件夹内新建 nanj 项目文件夹,项目内新建 brdc 、 igs 、 rinex 三个文件夹,分别存放当天的广播星历、精密星历和观测值文件。 注意:如果文件采用 .z 压缩格式,则应当使用 gunzip 命令进行解压;如果解压出的观测值文件仍采用 .d 的压缩格式,则应使用 crx2rnx 命令将其转化为 .o 的标准 RINEX 格式( sh_crx2rnx -f *.*d 可用于 .d 到 .o 文件的批量转换)。 2 、终端进入 nanj 项目文件夹内链接 tables ,运行: sh_setup -yr 2010 3 、生成 station.info 文件。将 nanj/tables 下的 station.info 文件拷贝到 rinex 文件夹下,打开并编辑,仅保留以 # 或 * 开头的前几行,保存并关闭。 打开终端并进入 rinex 文件夹,运行 sh_upd_stnfo -files *.10o 运行成功以后 station.info 文件便存放了此项目内的开始、结束时间,站名和接收机、天线类型等信息。 4 、建立 lfile. 文件。 lfile. 是测站的先验坐标文件。打开终端并进入 rinex 文件夹,用批处理的方式生成 lfile. 文件可以分为三步: a )提取观测值 .o 文件的先验 XYZ 坐标 grep POSITION *.10o lfile.rnx b )将 .rnx 文件转化为 .apr 文件 rx2apr lfile.rnx 2010 056 c )由 .apr 文件生成 lfile. 文件 gapr_to_l lfile.rnx.apr lfile. 2010 056 将 station.info 和 lfile. 两个文件拷贝到 nanj/tables 文件夹,覆盖原文件。 5 、 sestbl. 的配置。 nanj/tables 下的 sestbl. 文件是测段分析策略文件,该文件内部有详细的说明。一般来说采用默认配置即可,通常需要修改的三个地方是: Choiceof Experiment 选择处理方式 Choice of Observable 选择观测值类型和模糊度解算 Use otl.grid 选择是否使用潮汐文件 对于 Choice of Experiment ,选择 BASELINE 时将固定轨道并在 GAMIT 处理中和输出 h- 文件时忽略轨道参数;选择 RELAX 时将采用松弛解,合并全球 IGS h- 文件时需要。要想点位置精度高用 RELAX ;若目的是求基线后面平差则用 BASELINE 。在此实例中采用默认的 BASELINE 。 对于 Choice of Observable ,选择 LC_AUTCLN 为采用宽巷模糊度值并用伪距在 autcln 中解算;对于小于几公里的基线,用 L1 和 L2 独立载波相位观测值( L1,L2_INDEPENDENT )或者仅用 LI ( L1_ONLY ),相比用无电离层组合( LC_HELP )可以减少噪声水平。 对于 Use otl.grid ,由于这里已在 ftp 上更新 tables ,有了最新的 otl.grid 文件,所以这里选择 Y 。 6 、 sittbl. 的配置。 nanj/tables 下的 sittbl. 文件对各个测站的先验坐标(或钟差、大气模型等)进行约束。对高精度的已知坐标采取强约束,而对待求点采用松弛约束。如 IGS 站的坐标分量约束在较小的 1 ~ 75px ,对未知点的约束可以到 5 ~ 10m 。 7 、 sites.defaults 和 process.defaults 的配置。 sites.defaults 文件用来控制需要参与解算的测站。在文档末尾可以根据提示编辑,来给定那些不参与解算的测站或是测站的某些天。 process.defaults 文件用来控制处理过程中的很多细节,比如 sampling interval, number of epochs, start time for processing , default globk .apr file 等等,根据需求和提示进行编辑。 在这里均使用它们的默认值。 二、利用 GAMIT 解算基线 用终端进入 nanj 的项目文件夹,输入批处理命令进行解算: sh_gamit -expt nanj -s 2010 056 060 -orbitIGSF -yrext -noftp -dopt D ao c x sh_gamit.log 参数说明: -expt :指定四个字符的项目名称 -d :指定需要处理的指定日期,例如 -d 2010 56 60 ,指的是处理 2010 年第 56 和 60 天。 -s :指定需要处理的时间序列,例如 -s 2010 56 60 ,指的是处理 2010 年第 56 到 60 天。 -orbit :卫星轨道类型。 -yrext :给日目录前添加年前缀,例如 2010_006 。 -noftp :处理过程中不连接 ftp 下载数据。 -copt :数据处理完成后待压缩的文件类型,例如 -copt o q m k x 。 -dopt :数据处理完成后待删除的文件类型,例如 -dopt D ao c x 。 结果文件将存放在名称为年积日的文件夹内,此例中,可供参考的结果文件为: sh_gamit_2010_006.summary 解算总结 qnanja.006 解算记录 onanja.006 解算记录的简略版,一般关注此文件 hnanja.10006 协防差矩阵、参数平差值 在上述文件中,可根据描述来判断解算结果是否符合相应的需求。其中,基线解算结果( O 文件,即这里的 onanja.006 )中的 postfit_nrms 项优于 0.3 左右时最佳;如果大于 1.0, 则表示此解存在问题。 三、利用 GLOBK 进行平差处理 用终端进入 nanj 的项目文件夹内,运行: sh_glred -expt nanj -s 2010 050 2010 065 -yrext -opt H G E sh_glred.log 参数说明: H :运行 htoglb ,把文件转换为二进制文件 G :运行 glred ,合并文件 E :绘图 运行成功后,进入 gslon 文件夹,以 psbase 开头的文件即为各个站三维坐标时间序列的图形表示。 用终端进入 gsoln 目录下,依次运行: ls../glbf/h*glx nanj.gdl glred6 globk_comb.prt globk_comb.log nanj.gdl globk_comb.cmd 生成的 globk_comb.org 文件中就包含了解算点的三维坐标和相关参数。 至此,利用 GAMIT-GLOBK 解算 GPS 基线并进行平差的实例叙述完毕。 本文出于帮助初学者利用 GAMIT-GLOBK 软件解算 GPS 基线并平差的目的撰写,但受笔者水平限制,难免存在疏漏和错误,欢迎大家批评指正。
个人分类: GAMIT|17428 次阅读|2 个评论
圣经给我智慧的一个实例
热度 4 jiangjiping 2014-8-24 00:46
圣经给我智慧的一个实例 蒋继平 2014 年 8 月 24 日 圣经中充满着人类文明的智慧, 其中在我记忆中影响最深是一个年轻国王的智慧。 现在写出来与大家一起分享。 故事的情节是这样的:从前某国有一个年轻人被上帝选为该国的国王。有一天深夜, 上帝来到这位年轻人的面前, 让他许愿, 只要他提出要求, 不管多少财富资产, 或者其他要求,上帝一定会满足他的要求。这位年轻人确信上帝有这个能力, 可是,他没有要求上帝给他财富。 他对上帝说:“仆人住在你所拣选的民中,这民众多得不可胜数。所以求你赐我智慧,可以判断你的民,能明辨是非。不然,谁能判断这众多的民呢?”上帝当然很乐意地满足了这个年轻国王的要求。 这个年轻国王上台后, 应用他的智慧才能,尽力治理国家, 体察民意, 解决实际问题,深受人民的信任和爱戴。 因而,人民遇到疑难问题时,总是找他亲自判定。 有一天,他的部下带来了两个妇女和一个小孩。 部下告许国王,这两个妇女在一个无人的地方生下了两个男婴,其中的一个生下来三天后就被他自己的妈妈压死了。 这位妇人就乘另一个妈妈睡熟了的时候将另一个婴孩抢走了。现在她们两人都说这个婴孩是她们自己生的,我们没有任何证人可以证明当时的实际情景, 所以, 我们不能做出任何判定, 请国王亲自裁决。 各位读者, 当您读到这里的时候, 请您理解当时的历史条件下, 科技水平还很低下, 没有验血这样的亲子鉴定测试手段可以利用。 在您继续看下文之前, 要是您没有读过这段圣经,您自己如何来作出判定呢? 您将采取何种方式来使两位妇女露出破碇呢? 话说这位年轻但是充满着智慧的国王, 冷静地审视了面前的两位妇人,要求她们中抢走婴孩的那位将孩子归还给亲生妈妈。 可是,两位妇人坚持说孩子是她们亲生的, 没有一个愿意让步。在这样的僵持局面下, 国王叫他的部下拿一把剑来,他拿着利剑对两位妇人说, 既然如此, 我只能将孩子砍成两半, 你们每人得到一半。 听到国王的话后, 两个妇人的反应截然不同。 其中的一个非常同意国王的裁决, 而另一个非常惊恐, 要求国王千万不要如此做,请国王保全孩子的性命, 即使判给另一个妇人她也认了。 看到两个妇人的不同反应, 国王立即命令部下将孩子交给那位要求保全孩子性命的妇人。他确信这位妇人才是孩子的亲生妈妈。 读到这里, 各位读者, 您能理解这位国王的智慧吗? 当时, 我读到这个故事的时候, 我没有想到这个举措。 所以,我觉得我跟这位国王在智慧上还存在着一大差距。因而,我对这位国王产生了由衷的佩服和敬意。 这位年轻国王叫所罗门 ( Solomon )。
个人分类: 故事感想|3321 次阅读|4 个评论
[转载]【网络编程】【W2C标准asp入门教程】
lcj2212916 2014-6-25 11:11
【名稱】: W2C标准asp入门教程 【作者】: --- 【大小】:1.36MB 【格式】:PDF 【語言】:簡體中文 【內容簡介】: 在我们的 ASP 教程中,您将学到 ASP 的相关知识,以及如何在服务器端执行脚本。 您会发现,ASP 是一种生成动态交互性网页的强有力工具。 开始学习 ASP ! ASP 参考手册 在 W3School,我们为您提供完整的 ASP 参考手册,其中包括内建对象和组件,以及它们的属性和方法。 ASP 参考手册 ASP 实例 通过实例来学习!因为 ASP 脚本只能在服务器端执行,所以你无法在浏览器中查看 ASP 代码,你能看到的仅仅是由 ASP 输出的纯粹的 HTML 代码。在 W3School,每个实例均可显示出以往被隐藏的 ASP 代码。这样,您就可以更容易地理解它们的工作原理。 .............. 【下載載點】: http://www.400gb.com/file/67154533
1038 次阅读|0 个评论
[转载]【网络编程】【PHP4.0 时尚编程百例】
lcj2212916 2014-6-6 20:01
【名稱】:PHP4.0 时尚编程百例 【作者】:-- 【大小】:7.4MB 【格式】:pdg 【語言】:簡體中文 【內容簡介】: PHP 4.0是源代码公开的服务器端脚本语言,用于商业动态网页的制作以及其他网络需求。它具有丰富的函数库和多种数据库支持函数,并且是完全免费系统。本书通过100个实例,全面讲解PHP 4.0的开发方法和技巧。内容包括:如何进行一个普通网页程序的开发、如何进行一个数据库应用程序的开发、对数据库进行查询、对数据进行处理、如何使用多种PHP 4.0常用文件操作、数据库操作函数等,并将各种技巧应用于实例程序。\r\n 本书既可作为不同层次培训班的教材,也可以作为PHP 4.0开发人员的代码和命令参考手册使用 【下載載點】: http://www.400gb.com/file/66052153
1322 次阅读|0 个评论
[转载]matlab中SVM实例-Xu Cui
wl2119 2014-5-26 00:36
转载Xu Cui的SVM练习实例:http://www.alivelearn.net/svm-support-vector-machine-with-libsvm/ 代码在下面,拷贝复制到.m文件中运行即可。 cuixu_test_svm1 Linear Nonlinear example (radial basis) Nonlinear, circle Nonlinear, two circles Nonlinear, quadrant 3-class example % SVM with libsvm % Xu Cui, 2009/10/07 %% two classes, linear % data setup: our data contains two classes, each N samples. Thedata is 2D N = 500; l = ; % label d = ; % data % plot original data for visual inspection figure('color','w'); pos = find(l==1); plot(d(pos,1),d(pos,2),'r.'); pos = find(l==-1); hold on; plot(d(pos, 1),d(pos, 2),'b.'); axis equal % SVM with linear kernel (-t 0). We want to find the bestparameter value C % using 2-fold cross validation (meaning use 1/2 data to train,the other % 1/2 to test). Please note the parameter -g (gamma) is uselessfor linear % kernel bestcv = 0; for log2c = -1.1:3.1, for log2g = -4.1:1.1, cmd = ; cv = svmtrain(l, d,cmd); if (cv = bestcv), bestcv = cv; bestc =2^log2c; bestg = 2^log2g; fprintf('%g %g %g(best c=%g, g=%g, rate=%g)\n', log2c, log2g, cv, bestc, bestg, bestcv); end end end % After finding the best parameter value for C, we train theentire data % again using this parameter value cmd = ; tic;model = svmtrain(l, d, cmd);toc % now plot support vectors hold on; sv = full(model.SVs); plot(sv(:,1),sv(:,2),'ko'); % now plot decision area =meshgrid( , ); dd = ; tic; =svmpredict(zeros(size(dd,1),1), dd, model);toc pos = find(predicted_label==1); hold on; redcolor = ; bluecolor = ; h1 =plot(dd(pos,1),dd(pos,2),'s','color',redcolor,'MarkerSize',5,'MarkerEdgeColor',redcolor,'MarkerFaceColor',redcolor); pos = find(predicted_label==-1); hold on; h2 =plot(dd(pos,1),dd(pos,2),'s','color',bluecolor,'MarkerSize',5,'MarkerEdgeColor',bluecolor,'MarkerFaceColor',bluecolor); uistack(h1, 'bottom'); uistack(h2, 'bottom'); %% two classes, non-linear, radial basis function:exp(-gamma*|u-v|^2) % data setup: our data contains two classes, each N samples. Thedata is 2D N = 500; d = (rand(2*N,2)-0.5)*6; l = -1*ones(size(d,1),1); % here are 3 examples, uncomment one of them pos = find((d(:,1).^2 + d(:,2).^2)1); % one circle in themiddle %pos = find(((d(:,1)+1).^2 + d(:,2).^2)1 | ((d(:,1)-2).^2 +d(:,2).^2)1); %two circles %pos = find(d(:,1)0 d(:,2)0 | d(:,1)0 d(:,2)0); % quadrant l(pos) = 1; % normalization, but we don't do it here % d = (d-repmat(min(d, ,1)-min(d, ; cv = svmtrain(l, d,cmd); if (cv = bestcv), bestcv = cv; bestc =2^log2c; bestg = 2^log2g; fprintf('%g %g %g (bestc=%g, g=%g, rate=%g)\n', log2c, log2g, cv, bestc, bestg, bestcv); end end end % After finding the best parameter value for C, we train theentire data % again using this parameter value cmd = ; tic;model = svmtrain(l, d, cmd);toc % now plot support vectors hold on; sv = full(model.SVs); plot(sv(:,1),sv(:,2),'ko'); % now plot decision area =meshgrid( , ); dd = ; tic; =svmpredict(zeros(size(dd,1),1), dd, model);toc pos = find(predicted_label==1); hold on; redcolor = ; bluecolor = ; h1 =plot(dd(pos,1),dd(pos,2),'s','color',redcolor,'MarkerSize',5,'MarkerEdgeColor',redcolor,'MarkerFaceColor',redcolor); pos = find(predicted_label==-1); hold on; h2 =plot(dd(pos,1),dd(pos,2),'s','color',bluecolor,'MarkerSize',5,'MarkerEdgeColor',bluecolor,'MarkerFaceColor',bluecolor); uistack(h1, 'bottom'); uistack(h2, 'bottom'); %% three classes, linear % data setup: our data contains two classes, each N samples. Thedata is 2D N = 500; l = ; % label d = ; % data d(2*N+1:end,:) = ; % data % plot original data for visual inspection figure('color','w'); pos = find(l==1); plot(d(pos,1),d(pos,2),'r.'); pos = find(l==-1); hold on; plot(d(pos, 1),d(pos, 2),'b.'); pos = find(l==-3); hold on; plot(d(pos, 1),d(pos, 2),'k.'); axis equal % SVM with linear kernel (-t 0). We want to find the bestparameter value C % using 2-fold cross validation (meaning use 1/2 data to train,the other % 1/2 to test). Please note the parameter -g (gamma) is uselessfor linear % kernel bestcv = 0; for log2c = -1.1:3.1, for log2g = -4.1:1.1, cmd = ; cv = svmtrain(l, d,cmd); if (cv = bestcv), bestcv = cv; bestc =2^log2c; bestg = 2^log2g; fprintf('%g %g %g(best c=%g, g=%g, rate=%g)\n', log2c, log2g, cv, bestc, bestg, bestcv); end end end % After finding the best parameter value for C, we train theentire data % again using this parameter value cmd = ; tic;model = svmtrain(l, d, cmd);toc % now plot support vectors hold on; sv = full(model.SVs); plot(sv(:,1),sv(:,2),'ko'); % now plot decision area =meshgrid( , ); dd = ; tic; =svmpredict(zeros(size(dd,1),1), dd, model);toc pos = find(predicted_label==1); hold on; redcolor = ; bluecolor = ; blackcolor = ; h1 =plot(dd(pos,1),dd(pos,2),'s','color',redcolor,'MarkerSize',5,'MarkerEdgeColor',redcolor,'MarkerFaceColor',redcolor); pos = find(predicted_label==-1); hold on; h2 = plot(dd(pos,1),dd(pos,2),'s','color',bluecolor,'MarkerSize',5,'MarkerEdgeColor',bluecolor,'MarkerFaceColor',bluecolor); hold on; pos = find(predicted_label==-3); h3 =plot(dd(pos,1),dd(pos,2),'s','color',blackcolor,'MarkerSize',5,'MarkerEdgeColor',blackcolor,'MarkerFaceColor',blackcolor); uistack(h1, 'bottom'); uistack(h2, 'bottom'); uistack(h3, 'bottom'); %% Test libsvm's performance result = for dim = l = ; % label d = ; % data tic;model =svmtrain(l,d,'-t 0');x = toc; tic; = svmpredict(l, d, model);y = toc; result = ]; end end disp(result) %% % calculate w and b w = model.SVs' * model.sv_coef; b = -model.rho; if model.Label(1) == -1 w = -w; b = -b; end disp(w) disp(b) % plot the boundary line % x = ; % y = (-b - w(1)*x ) / w(2); % hold on; % plot(x,y) pause(1) return
个人分类: Matlab|3408 次阅读|0 个评论
R语言函数定义中参数的设定实例
Bearjazz 2012-11-9 18:13
#R语言函数定义中参数的设定实例# 熊荣川 xiong rongchuan 六盘水师范学院生物信息学实验室 xiongrongchuan@126.com http://blog.sciencenet.cn/u/Bearjazz fasextract - function(singlesave = TRUE){ + if (singlesave == TRUE) + return("welcome") + else return("well no") + } #定义一个函数,参数single初始设定为真,即默认为真 fasextract(F) #输入F,False的简写 "well no" fasextract(T) #输入T,True的简写 "welcome"
个人分类: 我的研究|5113 次阅读|0 个评论
公事出差报告的一个实例
热度 8 jiangjiping 2012-8-14 08:03
Five Major Advantages for Seed Production in China Jiping Jiang August 18, 2006 I made my China trip visit from June 18 to 28, 2006. During this trip, I had a meeting with Mr. Mengyu Zhang, who is the general manager of the vegetable section in China national seed group Corporation (CNSGC). The meeting was hold in Nanjing. We discussed many issues related to seed production in China. I also talked to many other professional seed production managers. Also, as an employee of Seminis, I have involved in seed production in China since 1996. I have met many professional people in seed industry in China, Besides Mr. Zhang, they were the CEO of China National Seed Group Corporation, Wei He, and many provincial and local managers in the seed business. I have visited many production sites in China, Which included Hebei, Tianjin, Liaoning, Gansu, and Shanxi. Therefore, I believe that I know the seed production business in China quite well. I think there are five major advantages for seed production in China. China is a big country. It has all weather types, including tropical and subtropical ones. It is good for seed production in 4 seasons year around. So it can ensure a continual seed production at any time of the year. Very favorable weather type for seed production in Northwest region, especially in Jiuquan area, Gansu Province. The vast Northwest region has a weather type with abundant sunlight, few rains, and low humidity during growing season, which is very favorable for the production of high seed quality. Abundant labor source. China has more than 1.3 billion people and more than half of them are farmers. Therefore, it is quite easy to find the farmers for the high labor-consuming job in China. This is an important factor, especially in the hybrid seed production section. Cheaper labor and land, but hard work and experienced people. The cost for the labor is quite low in China, especially in rural northwest area. It is about 1/10 cost compared to the same job in USA. The land is also cheaper in these areas comparing to USA. But the farmers are the people with the hard-working characteristics. They are experienced farmers with necessary skill to manage a good crop. Government supports the business. From Chinese central government to local governments, all level’s policy is favorable to the seed business, especially for the seed production for foreign seed companies. It provides priority for seed production in critical resource management, which includes water, fertilizers and pesticides. Based on these advantages, I strongly suggest that we should consider to make our hybrid tomato seed production in China. Or at least, we should shift some of our current production in other countries to China. 备注:这是我 2006 年到中国出差,回美国后给公司提交的一份报告。括弧中的中文是我的翻译。 在原始报告中没有中文。 顺便说一下,我本来跟公司请假,回国参加外甥女的婚礼, 结果公司让我以公事的身份回国,来回费用全部由公司负责,其交换条件就是在中国期间了解一下中国种子生产的一些相关情况。 自从我提交这份报告后,公司就将原来在南美洲智利生产的一部分计划转移到中国。 现在每年在中国的生产订量越来越大。
个人分类: 科研工作报告|7211 次阅读|18 个评论
加了几个安全标识牌
热度 3 Greg66 2012-8-14 01:06
这种安全标识牌很好 点击查看
3939 次阅读|4 个评论
只懂行业技术的人员不能管好安全的实例
热度 3 Greg66 2012-8-5 23:28
案例(基坑坍塌事故): 一技术人员在未进行细致的全风险(危险)评估的基础上制定了一条实际上无法执行的安全措施,未得到操作工人的执行,结果导致事故 , 造成 4 人死亡。 分析: 这案例里面包含的行业共性的安全方法是:( 1 )制定措施时要进行全面风险 分析和评估 后制定技术措施。( 2 )制定措施后,要 评估 措施执行过程中的困难及其解决办法,分析所制定措施是否能够被执行。( 3 )组织执行的人接到措施后,首先要 分析 措施执行的可能性、困难及其注意事项(而不是只当传话筒),再安排给工人具体执行。( 4 )发布措施后,要 观察、确保 工人的正确执行了措施,而不是发步完了就算了事。 这个共性的方法里面包含有深入研究的问题(措施的制定、执行程序)。这个措施执行过程中,没有安全专业人员或者没有充分安全知识的人,肯定难以保证,所以安全专业人员的作用是独特的。这个案例的详细描述就可以证明我说的是正确的。 我用这个例子仅仅是为了说明,行业间是有共性的安全问题、方法值得深入研究的。安全管理不等于行业技术,只会行业技术不能管理好安全,所以需要安全专业人员时,招聘行业技术人员来是不够的。
4265 次阅读|5 个评论
震群的定义与实例
qsqhopeiggcas 2011-12-18 22:29
震群在地震教科书上的定义并不严谨。根据我们的理解,震群事件可定义为在一个较小的区域范围内(属于同一个地震区)、 1 年左右的时间范围内,发生的一连串不同大小的地震事件。在震群事件中,往往有一个或两个震级显著较大的地震,称之为震群中的标志性地震事件。 通常认为,强震是断层中高强岩石发生宏观破裂的结果。岩石要发生宏观破裂,须经过岩石变形破坏过程的体积膨胀点,因此震群事件必然发生。震群事件是强震前唯一的、人类目前监测技术可识别的前兆现象。 以下是近些天来发生在 PUERTO RICO 的部分地震目录,这是一个典型的震群事件实例。 引自: http://earthquake.usgs.gov/earthquakes/recenteqsww/Quakes/quakes_all.php Update time = Sun Dec 18 13:05:13 UTC 2011 MAG UTC DATE-TIME y/m/d h:m:s LAT deg LON deg DEPTH km Region MAP 2.6 2011/12/18 07:31:33 18.196 -67.354 15.2 MONA PASSAGE, PUERTO RICO MAG UTC DATE-TIME y/m/d h:m:s LAT deg LON deg DEPTH km Region MAP 2.8 2011/12/17 23:08:43 18.202 -67.364 13.8 MONA PASSAGE, PUERTO RICO MAP 3.2 2011/12/17 23:07:20 18.177 -67.367 17.7 MONA PASSAGE, PUERTO RICO MAP 2.5 2011/12/17 19:47:51 18.180 -67.356 15.9 MONA PASSAGE, PUERTO RICO MAP 3.0 2011/12/17 11:43:20 18.166 -67.366 12.4 MONA PASSAGE, PUERTO RICO MAP 2.6 2011/12/17 10:30:42 18.191 -67.307 16.0 MONA PASSAGE, PUERTO RICO MAP 2.9 2011/12/17 10:05:32 18.212 -67.346 18.1 MONA PASSAGE, PUERTO RICO MAP 3.1 2011/12/17 10:05:32 18.165 -67.375 17.8 MONA PASSAGE, PUERTO RICO MAP 2.8 2011/12/17 09:48:34 18.188 -67.372 11.8 MONA PASSAGE, PUERTO RICO MAP 2.5 2011/12/17 09:13:07 18.167 -67.369 19.2 MONA PASSAGE, PUERTO RICO MAP 2.5 2011/12/17 08:00:57 18.194 -67.340 15.6 MONA PASSAGE, PUERTO RICO MAP 2.6 2011/12/17 07:34:48 18.191 -67.342 17.2 MONA PASSAGE, PUERTO RICO MAP 3.1 2011/12/17 07:32:20 18.179 -67.371 16.2 MONA PASSAGE, PUERTO RICO MAP 3.1 2011/12/17 07:30:21 18.195 -67.380 16.8 MONA PASSAGE, PUERTO RICO MAP 2.7 2011/12/17 07:21:57 18.238 -67.355 13.2 MONA PASSAGE, PUERTO RICO MAP 2.8 2011/12/17 07:11:13 18.161 -67.368 16.6 MONA PASSAGE, PUERTO RICO MAP 3.1 2011/12/17 07:03:18 18.178 -67.366 21.5 MONA PASSAGE, PUERTO RICO MAP 2.5 2011/12/17 07:01:08 18.174 -67.382 16.3 MONA PASSAGE, PUERTO RICO MAP 2.7 2011/12/17 06:55:14 18.174 -67.377 24.0 MONA PASSAGE, PUERTO RICO MAP 2.8 2011/12/17 06:49:17 18.173 -67.362 21.2 MONA PASSAGE, PUERTO RICO MAP 3.2 2011/12/17 06:45:56 18.184 -67.394 7.6 MONA PASSAGE, PUERTO RICO MAP 2.5 2011/12/17 06:40:09 18.152 -67.397 17.2 MONA PASSAGE, PUERTO RICO MAP 2.7 2011/12/17 06:38:23 18.160 -67.391 15.1 MONA PASSAGE, PUERTO RICO MAP 3.1 2011/12/17 06:32:18 18.158 -67.365 16.0 MONA PASSAGE, PUERTO RICO MAP 2.6 2011/12/17 06:23:08 18.189 -67.357 13.0 MONA PASSAGE, PUERTO RICO MAP 3.2 2011/12/17 06:21:03 18.185 -67.380 17.0 MONA PASSAGE, PUERTO RICO MAP 2.5 2011/12/17 06:16:39 18.186 -67.355 12.0 MONA PASSAGE, PUERTO RICO MAP 2.6 2011/12/17 06:15:40 18.194 -67.367 15.6 MONA PASSAGE, PUERTO RICO MAP 3.0 2011/12/17 06:14:20 18.186 -67.319 7.2 MONA PASSAGE, PUERTO RICO MAP 3.2 2011/12/17 06:11:03 18.170 -67.165 41.0 PUERTO RICO MAP 5.3 2011/12/17 06:09:09 18.172 -67.371 17.0 MONA PASSAGE, PUERTO RICO MAP 5.1 2011/12/17 06:06:12 18.188 -67.370 14.0 MONA PASSAGE, PUERTO RICO MAG UTC DATE-TIME y/m/d h:m:s LAT deg LON deg DEPTH km Region MAP 3.0 2011/12/16 22:27:44 19.235 -67.582 19.2 PUERTO RICO REGION MAP 2.8 2011/12/16 20:31:03 19.146 -66.805 16.0 PUERTO RICO REGION MAP 3.0 2011/12/16 19:59:30 19.141 -66.850 32.5 PUERTO RICO REGION MAG UTC DATE-TIME y/m/d h:m:s LAT deg LON deg DEPTH km Region MAP 2.8 2011/12/15 08:34:28 18.003 -66.939 13.0 PUERTO RICO MAG UTC DATE-TIME y/m/d h:m:s LAT deg LON deg DEPTH km Region MAP 3.2 2011/12/14 00:48:24 19.101 -66.808 38.0 PUERTO RICO REGION MAG UTC DATE-TIME y/m/d h:m:s LAT deg LON deg DEPTH km Region MAP 2.7 2011/12/13 17:06:22 19.084 -67.402 1.7 PUERTO RICO REGION MAP 3.1 2011/12/13 11:41:38 19.237 -67.488 28.0 PUERTO RICO REGION MAP 3.2 2011/12/13 11:06:12 19.215 -67.462 18.4 PUERTO RICO REGION MAG UTC DATE-TIME y/m/d h:m:s LAT deg LON deg DEPTH km Region MAP 2.5 2011/12/12 18:49:39 18.288 -67.174 31.1 PUERTO RICO MAP 3.1 2011/12/12 15:20:55 18.566 -66.644 102.9 PUERTO RICO REGION MAG UTC DATE-TIME y/m/d h:m:s LAT deg LON deg DEPTH km Region MAP 3.0 2011/12/11 20:34:28 18.304 -68.037 4.5 MONA PASSAGE, DOMINICAN REPUBLIC MAP 2.9 2011/12/11 19:11:17 19.129 -66.792 13.8 PUERTO RICO REGION
个人分类: 科普|3534 次阅读|0 个评论
[转载]晶格优化中ISIF取值的探究(六方MgO为例)
swx0789 2011-7-23 15:54
一、目的:搞清楚用vasp进行晶格优化中ISIF值的重要性 二、方法:采用分步优化和一次性优化两种方法探究 三、重点参数说明: ISIF:2-relax ions?/Y; || change cell shape?/N; || change cell volume?/N 3-relax ions?/Y; || change cell shape?/Y; || change cell volume?/Y 4-relax ions?/Y; || change cell shape?/Y; || change cell volume?/N 5-relax ions?/N; || change cell shape?/Y; || change cell volume?/N 6-relax ions?/N; || change cell shape?/Y; || change cell volume?/Y 7-relax ions?/N; || change cell shape?/N; || change cell volume?/Y POSCAR: a 0.0 -1.0 0.0 0.8660254037844 0.5 0.0 0.0 0.0 1.6230529595 2 Direct 0.6666666666666667 0.3333333333333333 0.750 0.3333333333333333 0.6666666666666667 0.250 四、实验步骤: A、采用分步优化法: 算法原理: a.先进行一次驰豫优化计算 b.将CONTCAR拷贝成POSCAR c.再进行一次静态计算,求出总能和体积,拟合状态方程,得出晶格常 数 I、取ISIF=2时:(#上下对应。上行为系列缩放系数;下行为相应的总能;所取CONTAR&volume为缩放系数取3.71时的结构情况) 2.91 3.01 3.11 3.21 3.31 3.41 3.51 3.61 3.71 -0.09768706 -2.60381322-2.62609278 -2.62612787 -2.62612791 -2.60064488-2.59725747 -2.59725252 -2.597253 volume of cell : 71.78 CONTCAR: 3.710000000000000 0.0000000000000000 -1.0000000000000000 0.0000000000000000 0.8660254037844000 0.5000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 1.6230529595000000 II、取ISIF=3时: 2.91 3.01 3.11 3.21 3.31 3.41 3.51 3.61 3.71 -1.47880875 -3.56749270-3.57829654 -3.57830700 -3.57830701 -3.55534231 -3.55344206 -3.55344601 -3.553446 volume of cell : 43.12 CONTCAR: 3.710000000000000 0.0000000000000000 -0.8480309609621116 0.0000000000000000 0.7344163553888973 0.4240154804810558 0.0000000000000000 0.0000000000000000 0.0000000000000000 1.3556928126902820 III、取ISIF=4时: 2.91 3.01 3.11 3.21 3.31 3.41 3.51 3.61 3.71 -0.24918897 -2.60258815 -2.63108804 -2.63118786 -2.63118806 -2.60762910 -2.60458644 -2.60459724 -2.604597 volume of cell : 71.78 CONTCAR: 3.710000000000000 0.0000000000000001 -1.0553127846071310 0.0000000000000000 0.9139276804085778 0.5276563923035652 0.0000000000000000 0.0000000000000001 0.0000000000000000 1.4573715661012910 IV、取ISIF=5时: 2.91 3.01 3.11 3.21 3.31 3.41 3.51 3.61 3.71 -0.24918897 -2.60258815 -2.63108804 -2.63118786 -2.63118806 -2.60762910 -2.60458644 -2.60459724 -2.604597 volume of cell : 71.78 CONTCAR: 3.710000000000000 0.0000000000000001 -1.0553127846071310 0.0000000000000000 0.9139276804085778 0.5276563923035652 0.0000000000000000 0.0000000000000001 0.0000000000000000 1.4573715661012910 V、取ISIF=6时: 2.91 3.01 3.11 3.21 3.31 3.41 3.51 3.61 3.71 -1.47880875 -3.56749270 -3.57829654 -3.57830700 -3.57830701 -3.55534231 -3.55344206 -3.55344601 -3.553446 volume of cell : 43.12 CONTCAR: 3.710000000000000 0.0000000000000000 -0.8480309609621116 0.0000000000000000 0.7344163553888973 0.4240154804810558 0.0000000000000000 0.0000000000000000 0.0000000000000000 1.3556928126902820 VI、取ISIF=7时: 2.91 3.01 3.11 3.21 3.31 3.41 3.51 3.61 3.71 -1.28296288 -3.56234825 -3.57418319 -3.57419641 -3.57419642 -3.55112710 -3.54922372 -3.54922856 -3.549229 volume of cell : 43.07 CONTCAR: 3.710000000000000 0.0000000000000000 -0.8434465875014347 0.0000000000000000 0.7304461715115043 0.4217232937507173 0.0000000000000000 0.0000000000000000 0.0000000000000000 1.3689584800243790 B:采用一次优化法: 与上面方法不同的是INCAR中取ISIF=3,取定的实验缩放系数a=3.21,计算的精度选得更高,且仅做一次优化。 结果为: TOTEN= -3.549568 eV volume of cell : 43.12 CONTCAR: Mg-hex 3.210000000000000 0.0000000000000000 -0.9801232359458133 0.0000000000000000 0.8488116211684474 0.4900616179729067 0.0000000000000000 0.0000000000000000 0.0000000000000000 1.5669819270266370 2 Direct 0.66666666666666430.33333333333333570.7500000000000000 0.33333333333333570.66666666666666430.2500000000000000 结论: 1、特殊位置的粒子随晶格形状的改变而变化,不论允许驰豫与否。 2、不同ISIF的取值对晶格优化影响不大,本例中缩放系数均在3.21~3.31之范围内 3、采用分步优化法得出的优化晶格比较可靠
个人分类: VASP|7936 次阅读|0 个评论
和谐颂-- 写在全国与全省“两会”前夕
hucs 2011-2-11 15:12
博士研究之人文 和谐颂 -- 写在全国与全省“两会”前夕 南昌大学(医学院) 胡春松 在职博士生 以人为本惠民生, 促进和谐推新政。 务实合作共发展, 韬光养晦谋和平。 给力幸福强军民, 预腐污堵净环境。 盛世颂党迎华诞, 锦绣中华家园春。 注: 这里,转载近期党和政府的两项“爱民为民”务实举措是提升百姓幸福指数或幸福感的两个有力实例。 转载来源: 埃及局势骤紧,滞留埃及中国游客的处境受到党中央、国务院的高度重视。中共中央总书记、国家主席胡锦涛作出重要批示,要求外交部协调有关部门,做好接回部署。 24 小时求助热线随即开通,外交部、国家旅游局、公安部、民航局等部门围绕尽最大努力在节前接回所有滞埃同胞的目标,迅速行动起来。 有国外媒体评论,中国是所有自埃及接返国民的国家中相关工作启动最早、速度最快、效率最高、规模最大的国家之一。 …… ( http://www.china.com.cn/news/txt/2011-02/10/content_21890799.htm ) 在 2011 年底之前,国家级、省级美术馆将全部向公众免费开放,同时全国所有的公共图书馆、文化馆(站)也实现无障碍、零门槛进入,公共空间设施场地全部免费开放,所提供的基本服务项目全部免费。日前,文化部、财政部共同出台《关于推进全国美术馆公共图书馆文化馆(站)免费开放工作的意见》,就全国美术馆、公共图书馆、文化馆(站)免费开放进行总体部署,这是继 2008 年博物馆、纪念馆实现全面免费开放后又一推进公共文化服务体系建设的重要举措。 …… ( http://news.xinhuanet.com/video/2011-02/11/c_121064723.htm ) 继续转载: 地方“一把手”论如何当好领导干部 2011 年 02 月 11 日 09:57:57  来源: 新华时政 汪洋:做维护群众利益"守护神" 俞正声:领导干部都该自问值得群众信赖吗 薄熙来:领导干部要算好人生大账 李鸿忠:饱含深情 提振精神 扑下身子
1140 次阅读|0 个评论
姜罗罗的参赛博文:生态复杂系统中 “石头-剪刀-布”博弈的实例
Fangjinqin 2010-6-7 07:37
新参赛博文: 姜罗罗的参赛博文:生态复杂系统中 石头-剪刀-布博弈的实例 姜罗罗 汪秉宏 1996 年,美国加利福尼亚大学的行为遗传学家 B. Sinervo 和 C. M. Lively 首次发现真实生态系统中的石头 - 剪刀 - 布博弈 。生活在美国加利福尼亚的雄性侧边斑点蜥蜴 (side-blotched lizard) 的喉部呈现出黄、橙、蓝三种不同的颜色,如图 1( 左图 ) ,并且喉部不同颜色的雄性侧边斑点蜥蜴表现出不同的行为。黄喉雄性侧边斑点蜥蜴较为柔弱,行踪诡秘,在与其他亚种蜥蜴争夺配偶的竞争中扮演偷情者的角色;橙喉雄性侧边斑点蜥蜴较为生猛好斗,占据较大的领地,拥有为数众多的雌性配偶,奉行一夫多妻制;蓝喉雄性侧边斑点蜥蜴的武力介于其他两个亚种蜥蜴之间,同伴间相互合作捍卫自己的领地和雌性配偶,奉行一夫一妻制。图 1( 右图 ) 展示了三个亚种蜥蜴活动范围统计:橙喉雄性侧边斑点蜥蜴的活动范围最大,蓝喉雄性侧边斑点蜥蜴次之,黄喉雄性侧边斑点蜥蜴的活动范围最小。 图 1 、 (左图)生活在美国加州的雄性斑点蜥蜴,依据其喉咙部位的颜色依次可分为三个亚种(从左到右):橙喉、蓝喉、黄喉雄性斑点蜥蜴。(右图)三个亚种蜥蜴的活动范围统计图。摘自文献 。 三亚种雄性侧边斑点蜥蜴竞争雌性配偶的策略相互循环制约,形成石头 - 剪刀 - 布博弈。在这个循环博弈中,如图 2( 左图 ) ,橙喉雄性侧边斑点蜥蜴武力争夺蓝喉雄性侧边斑点蜥蜴的雌性配偶(石头抑制剪刀);蓝喉雄性侧边斑点蜥蜴武力捍卫自己的雌性配偶不受黄喉雄性侧边斑点蜥蜴袭扰(剪刀抑制布);由于橙喉雄性侧边斑点蜥蜴领地太大,顾此失彼,黄喉雄性侧边斑点蜥蜴得以偷偷摸摸地溜到没有设防的领地,找到雌蜥蜴成功交配 ( 布抑制石头 ) 。这种雄性个体间繁殖策略的石头 - 剪刀 - 布博弈直接导致了各亚种的种群数量呈现周期振荡,如图 2( 右图 ) 。 B. Sinervo 等人发现,各亚种的种群数量在 1990-1995 六年间呈现周期振荡:当蓝喉雄性侧边斑点蜥蜴的数量占据优势时,橙喉雄性侧边斑点蜥蜴的数量增加,从而抑制蓝喉雄性侧边斑点蜥蜴继续占据优势;随后橙喉雄性侧边斑点蜥蜴占据优势,黄喉雄性侧边斑点蜥蜴的数量增加,从而抑制橙喉雄性侧边斑点蜥蜴继续占据优势;紧接着黄喉雄性侧边斑点蜥蜴优势,蓝喉雄性侧边斑点蜥蜴数量增加,抑制黄喉雄性侧边斑点蜥蜴继续占据优势。如此往复,三亚种雄性蜥蜴互相制约,轮流居于优势地位,从而使蜥蜴总体数量上不会繁殖过快。因此,这种石头 - 剪刀 - 布博弈的繁殖策略被认为是生态系统维持稳定的一个重要途径。此外,蓝喉雄性侧边斑点蜥蜴个体间合作以及利他性也备受关注。当两只蓝喉雄性侧边斑点蜥蜴正在保护其领地不受橙喉雄性蜥蜴的侵犯时,其中一只蓝喉雄性蜥蜴将会忘我地上前迎战入侵者。它的这种利他行为可能会使其丧失与雌性蜥蜴成功交配的机会。蓝喉雄性侧边斑点蜥蜴的这种利他行为与动物将基因继续传递给下一代的生存本能大相径庭。 B. Sinervo 等人深入分析了 1990-2003 年观测到的侧边斑点蜥蜴数据后指出,尽管利他主义行为可能会损害某一只参加战斗的蓝喉雄性蜥蜴自身的繁殖机会,但是它却成全了其他蓝喉雄性蜥蜴与雌性交配,从而保护了它们的基因在下一代个体中能够继续存在 。 图 2 、(左图)三个亚种雄性斑点蜥蜴的繁殖策略相互循环抑制。(右图)依据 1990-1995 六年间各亚种的种群数量做出的各亚种比率振荡图。摘自文献 。 无独有偶, B. Sinervo 等人还发现雄性欧洲普通蜥蜴 (Lacerta vivipara) 在求偶时也会采取类似于石头 - 剪刀 - 布博弈的繁殖策略 。雄性欧洲普通蜥蜴依据腹部颜色的不同可分为橙色、黄色、白色三种。橙色腹部的雄性蜥蜴是暴力征服者,它会侵入别的蜥蜴的领地,并强行与那里的雌性蜥蜴交配;黄色腹部的雄性蜥蜴则是欺骗者,它会趁机潜入橙色腹部雄性蜥蜴空出来的领地与没有防备的雌性交配;白色腹部的雄性蜥蜴则是合作者,它会紧紧守护它的配偶,并与其它白色腹部蜥蜴合作共同阻止黄色腹部雄性蜥蜴侵入。这样就形成了 征服 - 合作 - 背叛 相互循环制约的繁殖策略:征服者 ( 橙色腹部雄性蜥蜴 ) 能战胜合作者 ( 白色腹部雄性蜥蜴 ) ,合作者又能战胜欺骗者 ( 黄色腹部雄性蜥蜴 ) ,欺骗者又能战胜征服者。 石头 - 剪刀 - 布博弈能够非常简洁的描述生态体系中个体间的相互作用,可以帮助我们理解生态群落 (biological communities) 是如何建立的。为了理解生态群落中的物种多样性,很多因素被引入石头 - 剪刀 - 布博弈,如长程捕食、病毒传播等等。然而,生态体系的另一个重要特征――群落结构,却知之甚少。这种群落结构的简单形式是一些规则的空间斑图如螺旋波、靶波等 ,我们对这一类自组织结构有了较好的理解。一方面,这些群落结构在宏观上可以用复杂网络加以描述,如生态食物链网络。另一方面,其微观机制则可以通过动力学加以描述。因此,通过演化博弈理论和复杂网络理论研究生态群落结构的形成以及其对物种多样性的影响将是被关注的热点课题。 参考文献: B. Sinervo and C. M. Livel, The rock-paper-scissors game and the evolution of alternative male strategies, Nature 380 , 240 (1996). B. Sinervo, A. Chaine, J. Clobert, R. Calsbeek, L. Hazard, L. Lancaster, A. G. McAdam, S. Alonzo, G. Corrigan and M. E. Hochberg, Self-recognition, color signals, and cycles of greenbeard mutualism and altruism, Proc. Natl Acad. Sci. USA 103 , 7372 (2006). B. Sinervo, B. Heulin, Y. Surget-Groba, J. Clobert, D. B. Miles, k A. Corl, A. Chaine, and A. Davis, Models of density-dependent genic selection and a new rock-paper-scissors social system, Am. Nat. 170 , 663 (2007). B. Kerr, M. A. Riley, M. W. Feldman, and B. J. M. Bohannan, Local dispersal promotes biodiversity in a real-life game of rock-paper-scissors, Nature 418 , 171 (2002). T. Reichenbach, M. Mobilia , E. Frey, Mobility promotes and jeopardizes biodiversity in rock-paper-scissors games, Nature 448 , 1046 (2007). L.-L. Jiang , T. Zhou, M. Perc, X. Huang and B.-H.Wang, Emergence of target waves in paced populations of cyclically competing species, New J. Phys. 11 , 103001 (2009). 姜罗罗的参赛博文:生态复杂系统中 石头-剪刀-布博弈的实例
个人分类: 博文比赛|6398 次阅读|2 个评论
逻辑、范式、理论、 实例
wns 2010-2-3 20:34
逻辑、规范、模式、 理论、实例
个人分类: 生活|3428 次阅读|2 个评论

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-4-28 00:59

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部