Journey of Learning and Sharing

Anaconda是Python的一个免费增值（freemium）发布，用于处理大规模数据处理、预测式分析和科学计算，目标是简化包的管理和部署。它的包管理系统是conda。Anaconda包含了超过330个Python开源包。所以使用Anaconda平台进行python编程只需专注于内容而不用担心python支持包的问题。 前几天六维开放注册，就有新人问怎样才能快速涨分，我的经验是： 在热门资源的早期就开始做种 。按照这个思路就有两个问题出现： 1、怎么识别热门资源 2、如何在其早期就开始做种 怎么识别热门资源呢？我对电影区比较熟悉，就拿电影区举例。电影区的热门资源有两条规律： （1）院线上映大片的前几部或者导演演员的前续作品都很有可能是热门资源。比如王牌特工2上映前，王牌特工1肯定是热门资源。 （2）院线大片有资源后只要在电影区出现肯定是热门资源。 当然电视剧集区里，微博微信的热门话题剧集就是热门资源了。 如何在热门资源早期就开始做种？由于六维下线了RSS功能，所以如果没有其他提醒方式帮助，只能频繁的手动查看资源区，这样很浪费时间的。 这种重复性的工作就交给脚本吧！ 思路：六维新发布的资源都会显示在相应资源区的第一页，并且这一页的网址不会变。因此用脚本监测这一页的内容，如果出现想要的资源马上提醒。 方案：使用python制作脚本，使用urllib包登录六维并抓取监测页面内容，使用BeautifulSoup包解析抓取的内容、查找想要的资源是否出现，使用email、smtplib包发邮件提醒新资源出现了（如果手机有邮箱客户端就可以很快收到提醒了）。 以下贴上脚本源码，在红色位置填上自己的信息即可。 # 登录neu6抓取内容 import urllib.request  import urllib.parse  import http.cookiejar #设置时钟 import time # 解析抓取内容 from bs4 import BeautifulSoup import io  import sys # 发送邮件提醒 from email import encoders from email.hea...

不知什么原因，近几个月服务器上的shadowsocks总是会莫名地掉线，只能再次登录开启一遍。于是写了一个shell脚本去保持运行。 # keep_ss.sh #! /bin/bash count=$(ps -ef | grep shadowsocks | grep -v "grep" | wc -l) if [ $count = 0 ]; then     # shadowsocks shutdown, run again     setsid /etc/init.d/shadowsocks start     echo "shadowsocks restart at $(date)" >> ss.log fi 然后在crontab中添加例行任务，每10分钟运行一次脚本 01,11,21,31,41,51 * * * * root bash keep_ss.sh

Shadowsocks端口流量统计 添加端口到流量统计规则中 添加单个端口 iptables -I INPUT -p tcp --dport 1234 iptables -I INPUT -p udp --dport 1234 iptables -I OUTPUT -p tcp --sport 1234 iptables -I OUTPUT -p udp --sport 1234 添加多个连续端口 shell脚本 #!/bin/sh for ((i = 1234 ; i <= 1243 ; i++)) do iptables -I INPUT -p tcp --dport $i iptables -I INPUT -p udp --dport $i iptables -I OUTPUT -p tcp --sport $i iptables -I OUTPUT -p udp --sport $i done 将流量数据写入文件 shell脚本： check.sh #!/bin/sh # 分割线 echo "*******************************************************************" >> flow.txt # 时间戳 date >> flow.txt # 写入流量数据 iptables -n -v -t filter -L INPUT >> flow.txt iptables -n -v -t filter -L OUTPUT >> flow.txt 将流量数据清零，记录文件清空 shell脚本： set-zero.sh #!/bin/sh iptables -Z var=`date "+%Y-%m-%d %H:%M:%S" ` echo "Flow set to zero at ${var} ." > flow.txt 添加例行任务 目的：每周二20点记录一次数据，每月1号19点清零 在/etc/crontab文件中添加下面两行 00 20 * * 2 roo...

一、找BotFather拿到bot token     在telegram中私聊BotFather建立自己的bot，给bot取名，名字必须要以bot结尾。建好后自己的bot就有一个唯一的token，类似下面的一串字符 164354723:AAEjT6-IyNoXjt7miD0dwa-P5VmDTtHQC8 二、确认bot响应文件的位置     在写好bot响应文件后，要把bot放在网络上的一个位置，并且这个位置必须要加密的，即以https开头的一串网址。比如响应文件的名称为telbot.php，把它放在下面这个网址的位置： https://my.webhost.com/ 164354723:AAEjT6-IyNoXjt7miD0dwa-P5VmDTtHQC8 /telbot.php 上面网址中的红色设置和bot的token一样是为了确定这个唯一的位置，当然也可以任意设置。 三、告诉Telegram响应文件的位置 Telegram用下面网址的形式来设定webhook响应方式 https://api.telegram.org/bot [myauthorization-token] /setwebhook?url= [myboturl] 按照上面的网址形式，把自己创建的bot的token以及响应文件的位置填入，然后在浏览器中运行一下即可设置成功。比如： https://api.telegram.org/bot164354723:AAEjT6-IyNoXjt7miD0dwa-P5VmDTtHQC8/setwebhook?url=https://my.webhost.com/164354723:AAEjT6-IyNoXjt7miD0dwa-P5VmDTtHQC8/telbot.php 设置成功后，页面会显示下面的内容： {"ok":true,"result":true,"description":"Webhook is already set"} 四、在Telegram中给自己的bot发消息进行验证 php响应文件例子 <?php  define('BOT_TOKEN', 'YOURBOT:TOK...

特殊曲线拟合 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\)是流体粘度，使用插值进行...