Skip to main content

Posts

Showing posts from 2016

特殊曲线拟合

特殊曲线拟合 S型 \( y = 15 + \frac{10}{(1 + 2e^{-x})^{30}} \) 反S型 \( y = \frac{20}{10 + 2x^8}\) 反S型 \( y = \frac{0.5}{0.4 + e^{5.9x - 8.4}} \) 平滑上升型 \( y = \frac{2.559x}{0.582 + 6.513e^{-2.632x}} \)

正交试验设计

正交试验设计 概念 正交试验设计法(简称正交法)是统计数学的重要分支。它是以概率论数理统计、专业技术知识和实践经验为基础!C 充分利用标准化的正交表来安排试验方案,并对试验结果进行计算分析,最终达到减少试验次数,缩短试验周期,迅速找到优化方案的一种科学计算方法。它是产品设计过程和质量管理的重要工具和方法。 优点 采用正交法安排试验方案可以解决以下问题: 可以节省大量人力、物力、财力和时间。 能够明确影响试验指标各因素的主次顺序,即了解哪些因素重要,哪些因素次要。 可以迅速找到优化方案,在产品开发设计中,迅速找到优化方案,可以大大缩短产品开发设计周期;在生产过程中很快找到优化方案,可以尽快使生产工艺按最佳工艺条件运行,早日实现高效益。 通过试验结果分析,可以进一步指明试验的方向,克服盲目性。 常用名词 指标 在试验中需要考查的效果的特性值,简称为指标。指标与试验目的是相对应的。指标一般分为定量指标和定性指标,正交试验需要通过量化指标以提高可比性,所以,通常把定性指标通过评分定级等方法转化为定量指标。 因素 因素也称因子,是试验中考查对试验指标可能有影响的原因或要素,它是试验当中重点要考查的内容。通常用大写英文字母A、B、C等来表示。一个字母表示一个因素,因素又分为可控因素和不可控因素。可控因素指在现有科学技术条件下,能人为控制调节的因素;不可控因素指在现有科学技术条件下,暂时还无法控制和调节的因素。 正交试验中,首先要选择可控因素列入到试验当中,而对不可控因素,要尽量保持一致,即在每个方案中,要对试验指标可能有影响的不可控因素,尽量要保持相同状态。这样,在进行试验结果数据的处理过程中,就可以忽略不可控因素对试验造成的影响。 水平 试验中选定的因素所处的状态和条件称为水平或位级。例如:加热温度为\(70^\circ C\)、\(80^\circ C\)、\(90^\circ C\),这3个状态,可分别用“1”、“2”、“3”来表示。 正交表 正交法的基本工具是正交表。它是一种依据数理统计原理而制定的具有某种数字性质的标准化表格,以下面的符号表示 $$L_N(M^K)$$ 其中K表示因素数,是正交表中的纵列数;M表示水平数;N表示方案数,是正交表中的横行数;L是正交表代号,是Latin Squ

二维微可压缩流体流动问题

MathJax TeX Test Page 油藏描述 本次算例计算一个二维油藏模型,网格数20*20,油藏大小100ft*100ft*5ft,四口井分布于(5,5)、(5,15)、(15,5)(15,15),四口井分别以100\(ft^3/d\)的产量生产。油藏边界为无流动边界,初始压力5500psi,油藏压缩系数为\(1*10^{-7} 1/Pa\),初始孔隙度在0.3-0.4之间随机生成。x和y方向渗透率在0.6-0.4D之间随机生成。 流体描述 压力和粘度关系如下表所示 4000 5000 5500 6000 6500 7000 7500 8000 0.9100 0.9200 0.9243 0.9372 0.9650 0.9494 0.9812 1.0019 当计算压力值超出插值范围时使用表格中的最小值或者最大值进行替换。 初始油藏状态下的油体积系数为0.9,油压缩系数为\(1*10^{-9} 1/Pa\)。 方程求解 根据物质守恒方程、状态方程和达西定律,单相流动基本方程为 $$\frac{\partial }{\partial x}\left ( \beta_c \frac{A_x k_x}{\mu B} \frac{\partial p}{\partial x} \right )\Delta x + \frac{\partial }{\partial y}\left ( \beta_c \frac{A_y k_y}{\mu B} \frac{\partial p}{\partial y} \right )\Delta y - q_{sc }=\frac{V_b}{\alpha_c}\frac{\partial }{\partial t}\left ( \frac{\varphi }{B} \right )$$ \(\beta_c\)为传到率转换因子,在油田单位制下为1.127,在国际单位制下为1. \(\alpha_c\)为体积转换系数,在油田单位制下为5.614583,在国际单位制下为1. \(V_b\)为网格体积,\(A\)为截面积。 $$A_x = \Delta y * \Delta z$$ $$A_y = \Delta x * \Delta z$$ \(\mu\)是流体粘度,使用插值进行

课题追踪的两种方法

RSS订阅 1 什么是RSS RSS(简易信息聚合)是一种消息来源格式规范,用以聚合经常发布更新数据的网站,例如博客文章、新闻、音频或视频的网摘。RSS文件(或称做摘要、网络摘要、或频更新,提供到频道)包含全文或是节录的文字,再加上发布者所订阅之网摘数据和授权的元数据。 Really Simple Syndication“聚合真的很简单”就是RSS的英文原意。把新闻标题、摘要(Feed)、内容按照用户的要求,“送”到用户的桌面就是RSS的目的。RSS一词有时候大体上意为社会性书签,包括各种RSS的不同格式。例如,Blogspace对使用网摘于一集成器内之动作标为RSS info和RSS reader。虽然它的第一个句子就包含明确的 Atom 格式:“RSS和Atom文件能够用简单的格式从网站更新消息至你的电脑!”( 摘自维基百科 ) 另外XML也是一种常见的格式。 2 用什么工具 获取或订阅RSS信息的工具可以称作RSS阅读器。我现在使用的RSS阅读是Feedly和inoreader,其他别人推荐的RSS阅读器还有很多,请参考知乎问答: 2015 国内最好的 RSS 阅读器是什么? 个人认为一个好的RSS阅读器要能够做到: (1)即时抓取订阅的信息; (2)友好的阅读界面; (3)方面的保存和稍后阅读功能; (4)友好的分享功能; (5)多平台应用支持。 基于以上原因,在此推荐inoreader。 3 如何使用 所有RSS阅读器的订阅步骤都大同小异,下面以inreader为例,详细介绍如果订阅Nature期刊中的感兴趣主题。 (1)注册账户,登录 inoreader网站 。 (2)寻找RSS订阅源 打开 Nature期刊主页 ,在搜索框中搜索感兴趣的话题,比如“Fluid”,结果页面会出现相关文章,并且右边会出现相关的主题,如图中的Fluid dynamics。 点击进入主题,就会看到有RSS图标样式(几乎所有网站都使用此图标)的链接,右键此链接-复制链接地址,这个RSS Feed的地址就是要寻找的源。 (3)将订阅源添加到阅读器 将上一步得到的源地址粘贴到inoreader中的搜索/订阅框中,阅读器会自动识别此订阅源,识别出后点击订阅即可订阅成功。以后当此网站的主题更新后会立刻被

当我跑步时我在想些什么

从《挪威的森林》开始接触村上春树的作品,虽然去年也在“北京不止读书会”听别人聊过村上的《当我谈跑步时我谈些什么》,但其实那时并没有读,直到最近才真正读完这本写他自己的跑步随笔。给我的第一个感觉是:他跑的好多。每天十公里的跑量我是望尘莫及,几乎每次都是四小时以内的马拉松成绩更是让我这个二十几的小伙自愧不如。 曾经有很多人问我,当我跑步时我在想些什么。其实大部分时候真的什么都不想,只是戴着耳机听歌,途中看到和我速度差不多的人就一起迈步向前,跟不上了就放弃追赶,嫌他速度慢就大步超越;当确实有事情想的时候,跑着跑着也不再去想了,或者还是事情的原本状态而没有进一步思考。很多人觉得跑步很枯燥,除了跑以外就没有了其他动作,但其实跑步时所看到的、听到的、以及引发的联想,这些东西都会使跑步更加有趣。每次身边一起跑的人会不一样,看到的事情会不一样,或者偶尔换一下跑道欣赏路旁不同的景色,将这些当作跑步时的小确幸也着实不为过。 跑步也是一件很有哲学韵味的运动。有句话叫“人生就像一场马拉松”,体会这句话最深的时候是在2014年的北马。我和同校的小伙伴们一起出发,然而在十公里饮水处走散。从十公里到半程终点,一直跟着一个在外放音乐的大叔,边听歌边跑感觉身上充满了力量。可是到25公里后渐渐赶不上大叔的脚步,随后一直到全程的终点都是自己一个人尽全力地跑,被很多人超越,也超越了很多人。联想到自己的经历,从小到现在有的人只是在一段时间内相识就再也没见过,有的人却是一路陪伴,更多的只是身边的人来来去去,最终也将只剩下自己迎接终点。 村上说跑步和写小说一样需要耐心,对我来将跑步和写文章一样都要扎实稳进。一点一点地想一段一段地写才能完成预想的目标,一步一步地跑才能达到的运动量。跑步带给我的不仅仅是一个还算强壮的身体,还有一份承受寂寞脚踏实地向前走的心态。 以下是从书中摘抄的句子。 说起来,我这个人是那种喜爱独处的性情,表达得准确一点,是那种不太以独处为苦的性情。每天有一两个小时跟谁都不交谈,独自一人默默地跑步也罢,四五个小时伏案独坐,默默地写文章也罢,我都不觉得难熬,也不感到无聊。这种倾向从年轻时起便一以贯之,始终存在于我的身上。和同什么人一起做什么事相比,我更喜欢一人默不作声地读书,或是全神贯注地听音乐。只需一个人做的事情,我可以想出许多许多来。 当受到某人无缘无故(至少我看来是如此)

telegram中的Sci-Hub机器人,又一文献下载利器

或许你看到标题会问什么是telegram,什么是Sci-Hub?请听我一一道来。 什么是Sci-Hub Sci-Hub是一个线上 数据库 ,其上提供48,000,000篇科学学术论文和文章。网站透过“.edu”代理服务器访问相关页面,每天会上传新的论文文章。2011年,哈萨克研究生亚历珊卓·艾尔巴金(Alexandra  Elbakyan)因为研究论文成本过高,每篇论文在付费墙机制下通常需要花费30美元,而决定成立Sci-Hub。2014年,学术界开始预测网站将会发展为类似Napster的服务。不过到了2015年,学术出版社爱思唯尔向纽约地方法院提交诉讼,指控Sci-Hub已经侵犯版权。纽约地方法院在2015年10月28日仍下令Sci-Hub原本使用的网域名称“Sci-Hub.org”必须终止。爱思唯尔在法院上获得胜诉后,一群研究人员、作家和艺术家则连署一封表态支持Sci-Hub和创世纪图书馆的公开信,声称这次诉讼对于世界各地的研究人员是“重大打击”,并指出:“它同样贬低我们、作者、编辑和读者。它寄生于我们的劳动,它阻挠我们为大众服务,它阻拦我们进入。”而该计划于11月因法院命令中止后,在同一个月内便改用网域名称“.io”重新上线,并开放使用Tor浏览。2016年1月时,Sci-Hub平均每天约有200,000人访问,Sci-Hub则声称网站服务每天平均有数十万次档案请求。  Sci-Hub是目前已知第一个提供大量自动且免费的付费学术论文的网站,使用者不需要事前订阅或付款,就能够使用原本存放在付费数据库的论文文章,并提供搜寻原先出版社网站内的文件档案服务。 以上介绍来源于维基百科词条 Sci-Hub Sci-Hub网站被屡次下线,但是又通过更换域名重新上线。以下三个网址经测试可以使用:  http://www.sci-hub.bz/   http://www.sci-hub.ac/   http://www.sci-hub.cc/   广大学者将自己的文章发表至学术期刊(免费或者支付版面费),然而当需要查看其他学者的文章时还需要向出版商付费,你是不是也觉得这完全阻碍了科学文化的传播。艾尔巴金在为自己辩护时援引联合国《世界人权宣言》第二十七条所提的:“人人有权自由参加社会之文化生活,欣赏艺术,并共同襄享科学进步及其利益。”关于这个问题,另一位被称

梦想、现实与生命的尊严哪一个更重要

以前读完一本书习惯在豆瓣标记一下,不会写下些什么,最多只是回想一下整本书讲了什么、还记得什么,但是时间久了最后留在记忆里的也只剩自己曾经看过这本书这件事了。所以决定以后看完书尽量记录下自己的想法,抄下一些印象深的句子,哪怕仅作为回首时还能够追忆的一丝痕迹。 上午看完了《月亮与六便士》,毛姆通过叙述主人公追寻绘画艺术的经历以及和主人公有接触的人的生活,来阐述人在梦想和现实面前的抉择。 当我在写题目时,想到的就是《我是一只小小鸟》中的歌词:“生活的压力和生命的尊严哪一个重要”。每个人在做一个对自己很重要的选择时,究竟以怎样的缘由作为目标。经历和环境的不同导致选择自然不同,但是无论选择什么都没有对错之分,旁人也无需过多评论。 一些摘抄到的句子 : "If you look on the ground in search of a sixpence, you don't look up, and so miss the moon."   我想他这一段生活是很富于浪漫情调的,但是他自己却绝对没有看到任何浪漫的色彩。或许一个人如果想体会到生活中的浪漫情调就必须在某种程度上是一个演员;而要想跳出自身之外,则必须能够对自己的行动抱着一种既超然物外又沉浸于其中的兴趣。   我很怀疑,阿伯拉罕是否真的糟蹋了自己。做自己最想做的事,生活在自己喜爱的环境里,淡泊宁静、与世无争,这难道是糟蹋自己吗?与此相反,做一个著名的外科医生,年薪一万镑,娶一位美丽的妻子,就是成功吗?我想,这一切都取决于一个人如何看待生活的意义,取决于他认为对社会应尽什么义务,对自己有什么要求。但是我还是没有说什么,我有什么资格同一位爵士争辩呢?   轮船缓缓驶出咸水湖,从珊瑚礁的一个通道小心谨慎地开到广阔的海面上,这时,一阵忧伤突然袭上我的心头。空气里仍然弥漫着从陆地飘来的令人心醉的香气,塔希提离我却已经非常遥远了。我知道我再也不会看到它了。我的生命史又翻过了一页;我觉得自己距离那谁也逃脱不掉的死亡又迈近了一步。 上帝的磨盘转动很慢,但是却磨得很细。   “我不需要爱情。我没有时间搞恋爱。这是人性的弱点。我是个男人,有时候我需要一个女性。但是一旦我的情欲得到了满足,我就准备做别的事了。我无法克服自己的欲望,我恨它,它囚禁着我的精神。我希望将来能有

MatLab中patch函数的基本用法

patch是用来构建多边形的一个基本函数。 用法一 patch(X,Y,C) patch(X,Y,Z,C) patch( 'XData' ,X, 'YData' ,Y) patch( 'XData' ,X, 'YData' ,Y, 'ZData' ,Z) 1.1 说明 patch(X,Y,C)用来构建一个或者多个可填充的多边形,其使用X和Y作为每个点的坐标值,patch将会按顺序连接每个点。如果要得到一个多边形,将X和Y设置为向量;如果要得到多个多边形,将X和Y设置为矩阵,没一列对应一个多边形。C决定多边形的颜色,可以是系统认定的字符,也可以是一个数值,也可以是RGB向量。 patch(X,Y,Z,C)用来构建三维坐标下的多边形。 patch(‘XData’,X,’YData’,Y)和patch(‘XData’,X,’YData’,Y,’ZData’,Z)的用法与patch(X,Y,C)和patch(X,Y,Z,C)的用法类似,只是不设定颜色。 1.2 例子 1.2.1 x = [ 0 1 1 0 ] ; y = [ 0 0 1 1 ] ; patch(x,y, 'red' ) x和y都是1*4的向量,表示将四个点(0,0)、(1,0)、(1,1)和(0,1)依次连接,最后闭合形成一个四边形,设定颜色为红色。 1.2.2 x2 = [ 2 5 ; 2 5 ; 8 8 ] ; y2 = [ 4 0 ; 8 2 ; 4 0 ] ; patch(x2,y2, 'green' ) x2和y2都是3*2的向量,两列表示画两个多边形。第一个多边形连接的点依次是(2,4)、(2,8)和(8,4),第二个多边形连接的点依次是(5,0)、(5,2)和(8,0),颜色设定为绿色。 1.2.3 如果上例的三角形第一个是红色,第二个是绿色,那么patch代码修改为 x2 = [ 2 5 ; 2 5 ; 8 8 ] ; y2 = [ 4 0 ; 8 2 ; 4 0 ] ; patch(x2(:, 1 ),y2(:, 1 ), 'red' ) pat