2017年12月_Mr-Cat伍可猫

转载 python运算符

转自：http://www.runoob.com/python/python-operators.htmlPython 运算符什么是运算符？本章节主要说明Python的运算符。举个简单的例子 4 +5 = 9 。例子中，4 和 5 被称为操作数，"+" 称为运算符。Python语言支持以下类型的运算符:算术运算符比较（关系）运算符赋值运

2017-12-25 17:45:27 265

原创 python中格式符的应用%s,%d,%f以及format（）的实例以及输出格式

python中用%代表格式符，表示格式化操作，常用的操作有%s,%d,%r等.1.%s,%r,%d分别表示字符串以str（），rper(),以及十进制整数表示，%f表示结果为浮点型，更多见下表：例如：输入：结果为：如果上图中括号内的6和school顺序对换，结果如下：提示说%d要求输入的数字类型，而不是字符串类型。紧接着的两条命令分别

2017-12-25 17:44:14 14024

原创 python 死循环捕捉输入异常

死循环即无法靠自身控制终止的循环。死循环并非一无是处，while True or while 1是死循环的普遍用法，死循环一直执行等待程序中断，然后去处理中断程序。举例1while True: try: x=int(input('输入一个值：')) break except ValueError: print('Oops,tha

2017-12-25 17:43:26 4483

原创 python以及MATLAB终止循环的快捷键

MATLAB：Ctrl+c或者Ctrl+Breakpython:Ctrl+c

2017-12-25 17:42:04 5250

原创 python对字幕的改动

#想学英语，下了英语字幕，发现有很多时间戳，打印下来会很占空间，于是用python将时间戳去掉。def model1(): f1=open('coherence.txt','r')#将下载的字幕复制（或者更名）为coherence.txt（只要是txt文件即可）并放在python根目录下 f2=open('test2.txt','w')#打开（没有则新建）test2文件 a=['

2017-12-25 17:38:34 657

原创 mpirun命令

MPI程序的一般启动方式是mpirun -np 如：在目录下运行的程序：~/programs/CosmoMC/camb>mpirun -np 4 ./cosmomc test.ini

2017-12-25 17:37:18 41201

原创 Xshell Linux常用命令及vim与权限修改等常见问题

Xshell :一个强大的安全终端模拟软件，它支持SSH1, SSH2, 以及Microsoft Windows 平台的TELNET 协议.（百度百科）此文讲述的是xshell下的Linux命令，此命令在help下大都可以找到。在xshell终端下常用的命令有：@[TOC] 目录## 第一部分0.帮助命令与清屏help 查看帮助命令ls --help ...

2017-12-25 17:36:15 7292 1

原创 matlab 2016b更改工作路径

1打开matlab2从打开的matlab窗口中“打开”（或者open）matlab下的toolbox\local\matlabrc.m,在最后门新增cd yourworkpath,如： 3关闭matlab,重新启动即可，此时输入cd ,出来的即是修改的工作路径

2017-12-25 17:33:24 2266

原创 matlab fsolve()函数的使用

fsolve()函数通常用于数值求方程或方程组的解，更常用于求解非线性方程组。其基本结构如下：（可以在matlab命令页面输入help fsolve查询）fsolve()解决的方程形式为F(X)=0; 用法为：X = fsolve(FUN,X0,OPTIONS)或者写为[x,fval,exitflag]=fsolve(fun,x0,options)其中f...

2017-12-25 17:31:45 147336 9

原创 matlab 分子式

对于分子式求解中，因为要求个数是整数，故重新采用循环的结构：clccleari=1;file=importdata('FTICR1.txt');C=12.0107;H=1.00794;O=15.9994;N=14.0067;P=30.973762;S=32.065;fid=fopen('FT1.txt', 'w');for i=1:length(file) minfun=ab

2017-12-25 17:30:08 867

原创 mathematica求定积分和不定积分

参考网址：http://xuxzmail.blog.163.com/blog/static/251319162011784175529/A.定积分：格式：integrate[f(x),{x,a,b}]举例1普通求解： 2.对于含有三角函数等需要将其化为数值的情况：0.0代表是浮点值，因此结果也要化为浮点值。3.对于下限是变量的情况

2017-12-25 17:22:25 16044

原创 CMB中相关仪器设备和术语

DASI （degree angular scale interferometer）度角尺度干涉仪，位于南极NSF （美国）国家科学基金会（National Science FoundationNASA 美国国家航空和宇宙宇航局（National Aeronautics and Space Adiministration）WMAP 威尔金森微波

2017-12-25 17:21:07 384

原创置信度和置信区间

置信度：也称为可靠度，或置信水平、置信系数，即在抽样对总体参数作出估计时，由于样本的随机性，其结论总是不确定的。因此，采用一种概率的陈述方法，也就是数理统计中的区间估计法，即估计值与总体参数在一定允许的误差范围以内，其相应的概率有多大，这个相应的概率称作置信度。简言之即是样本估计总体时的可信度。比如说一百组实验来估计总体均值时，有90组实验可以得到总体均值，则置信度为90%。置信区间：在

2017-12-25 17:19:31 6385

原创 excel实现统计词频功能

1）将文件复制到excel中2）统计A列中分子式的频数，将A列复制到B列中并选中B列3）选择B列4）在C1的单元格中输入命令=COUNTIFS(A:A,B1)并下拉到B列结束的位置：统计B列单元格中的每个词语在A列中出现的次数

2017-12-25 17:17:59 20559

原创解决电脑没有任何图标问题（以及没有底部的任务栏）

此问题是关闭了explorer所致，所以解决办法是运行explorer.exe，方法如下：Ctrl+alt+delete----“启动资源管理器”，在应用程序一栏点击右下角的新任务，输入explorer.exe，点击确定，关闭任务管理器即可。

2017-12-25 17:15:29 1871

原创 office 中解决ppt图标显示不正常问题（word,excel正常）

笔者下载了office2010安装完成后ppt的图标（后缀为pptx）显示异常，显示为空白; （原因是之前使用了默认的wps,之后安装office考虑到与wps可能不兼容的问题卸载了wps,但系统依然默认以wps打开的）笔者采用了网上说的从属性中更改图标的方法，结果图标显示正常，但是打开文件异常，不能正常显示文件内容（可以从最左下角开始中打开Microsoft powerpoi

2017-12-25 17:13:38 15794

转载微波背景辐射的发现(获诺贝尔奖)——彭齐亚斯和威尔逊

转载微波背景辐射的发现——彭齐亚斯和威尔逊　　早在20世纪40 年代，伽莫夫等人根据热大爆炸宇宙学说的观点，预言宇宙空间应该充满着残余辐射，它们的温度已经相当低了，大致为几K 或至多几十K。“K”是开氏温标或者叫热力学温标中的温度单位“开耳芬”的符号，正像我们用“℃”来表示摄氏温标中的温度的道理是一样的。0 摄氏度（0℃）相当于273.15K，3K 相当于-

2017-12-25 17:12:05 5900

转载 MCMC蒙特卡罗方法

有一位博客园的博主对此写得很详细了，链接摘抄如下：MCMC(一)蒙特卡罗方法文章转载于：http://www.cnblogs.com/pinard/p/6625739.htmlMCMC(二)马尔科夫链http://www.cnblogs.com/pinard/p/6632399.htmlMCMC(三)MCMC采样和M-H采样

2017-12-25 17:09:43 1608

原创 CosmoMC命令CAMB文件画功率谱

CosmoMC命令CAMB文件画功率谱画功率谱步骤：1)进入software文件夹，解压CAMB-Nov2016.tar.gz文件 2）进入刚建立的文件夹CAMB-Nov2016文件夹，vim Makefile(进入后将-fast去掉)-------make all(注意：文中make clean是清除刚才编译的文件

2017-12-25 17:05:42 1295

原创笔记本/台式机作为扩展屏

A，笔记本作为主机，外接扩展屏（如投影仪，台式电脑等）1,笔记本开机 2,用VGA线连接笔记本和投影仪 3win7下按win+p键有 “复制”—电脑和投影仪同时显示 “扩展”—将笔记本电脑扩展到投影仪上（相当于加宽笔记本屏幕） “仅投影仪”—-电脑没有显示，只有投影仪有显示 (可以参考百度经验)B，台式电脑作为主机，另一台电脑或者投影仪作为扩展（没有找到将笔记本作为扩展的方法，因为笔记本

2017-12-24 14:13:28 14071

原创 matplotlib画图中文显示问题

matplotlib不支持中文显示，以下以画一个cos函数为例代码如下：import matplotlib.pyplot as pltimport numpy as npx = np.linspace(-3,3,1000)y = np.cos(x)plt.plot(x,y)plt.title ('y和x的关系图')plt.show()图如下：可以看到有中文的地方出现乱码以下给几

2017-12-22 20:58:31 415

原创 win7环境变量设置

以下以google浏览器为例要将google浏览器设置为可以通过终端（ternimal）打开，方法如下： 1）将谷歌浏览器图标（或者exe）右键附加到开始附件，然后在开始中输入google，找到图标，右键“属性”，复制“快捷方式”下的“起始位置”。 2）“计算机”（或者“我的电脑”）右键“属性”—–高级属性设置—–“高级”下的—-环境变量—-系统变量中找到path—-点击编辑—-在变量名最后

2017-12-22 15:53:43 1070

原创反弹模型（bounce model）----adot, H, Hdot变化图

import matplotlib.pyplot as pltimport numpy as np#这里主要以y = x^2+1(y即尺度因子a)来模拟a,H,aH,Hdot,adot,等随x的变化的图def slope_y(x,y):#y的斜率，即adot xs = x[1:] - x[:-1] ys = y[1:] - y[:-1] ydot = ys/xs H

2017-12-21 17:20:57 1095

原创 numpy shape函数

shape函数返回list，array,matrix等的一维和二维长度值。 1）list列表列表只有一维，二维为0>>> a = [1,2,3]>>> np.shape(a)(3,)上面是一维列表，长度是3 值得注意的是，列表是没有shape函数的，如下：>>> a.shape[0]Traceback (most recent call last): File "<pyshell#

2017-12-19 13:43:43 4250

原创 python bar

以下为bar的一个实例import matplotlib.pyplot as pltimport numpy as npdata = [[13,3,2,55,2,51,4,23], [12,34,1,13,55,24,13,13,4], [11,23,13,34,52,13,45,12,32,13], [12,43,24,13],

2017-12-19 13:21:53 4709

转载 '','args','','kwargs'

args即postional arguments,kwargs即keyword arguments 这两个是python中的可变参数。*args 表示任何多个无名参数，它是一个tuple；**kwargs 表示关键字参数，它是一个dict。并且同时使用*args和**kwargs时，必须*args参数列要在**kwargs前. 具体可以参考博文，写得很详细

2017-12-18 22:20:36 137

原创熵扰动和绝热扰动

宇宙中的结构起源于早期的密度扰动，初期的扰动分为绝热扰动和熵扰动。绝热扰动：参加扰动的各种组分之间没有热交换，因此每个组分的熵不变。绝热扰动中的每个组分对平均密度的偏离正比于该组分的平均密度，即平均密度越大，扰动就越大。此类扰动也叫曲率扰动，因为在扰动理论中，总密度扰动（从而每组分的扰动）正比于空间曲率。熵扰动：扰动过程中，熵不再是不变量，而是改变的量。对于每种组分而言，粒子数的相对扰动减去

2017-12-18 22:14:04 771

原创 numpy.loadtxt画功率谱图

（非专业相关可以选择性阅读）从CAMB中可以得到.dat文件，由于第一行有L,TT,EE等字符使得numpy无法正常读取文件；故可以手动将第一行去掉，剩下的以python写代码画图，如下：import matplotlib.pyplot as pltimport numpy as npdata = np.loadtxt('test_scalCls.dat')L = data[:,0]TT

2017-12-17 23:08:55 2425

原创 python读取txt文件并画图

1，使用python读取txt文件已知txt文件内容如下： 0 0 1 1 2 4 3 9 4 16 5 25 6 36请以第一列为x轴，第二列为y轴画图步骤如下： 1）使用readlines读取文件 2）建立两个空列表X,Y，将第一列的数字放入X，第二列的数字放入Y中 3）以X，Y为轴画图实现如下：im

2017-12-17 22:31:41 32221 2

原创 python自定义函数和类并调用

python可以自己定义一个函数，同时在其他地方调用，举例如下：定义一个斜率函数import matplotlib.pyplot as pltdef slope(x,y): xs = x[1:] - x[:-1] ys = y[1:] - y[:-1] plt.plot(x[1:],ys/xs)保存为testgeany.py，保存路径为python根目录：C:\User

2017-12-17 20:16:47 7473

原创 numpy.log(math.log)

此处先介绍log常用的两个底数计算（以10为底和以e为底），以后再补充完毕 1以10为底：import numpy as npnp.log10(x)如：>>> np.log10(100)2.0>>> 2，e为底 log下什么都不写默认是自然对数如：>>> np.log(np.e)1.0>>> np.log(10)2.3025850929940459

2017-12-14 20:34:38 80274

原创将带空格的字符串去掉空格.join()函数

举例：f = " pri ng tshg "f是前后中间都有空格，如果要去掉前面的空格，可以如下：>>> print(f.lstrip())pri ng tshglstrip()清左边空格，rstrip()清右边空格，但不能清除中间空格，可以借助”.join()函数f = " pri ng tshg "a = f.split()b = ''.join(a)print(

2017-12-02 22:10:16 3328

原创设置添加环境变量Path

以geany编译器为例：在终端打开geany时要进入到其安装目录下才能打开，这很麻烦，因此可以将其设置添加到环境变量中，过程如下： 1，点击Geany图标，右键属性，找到目标路径（如果没有快捷键可以在开始中输入geany再右击） 2，复制这个路径 3，打开控制面板—系统与安全—-系统—-高级系统设置—-高级—–环境变量—-在系统变量中找到Path—-点击编辑—–在变量值中移到最右边（

2017-12-02 20:51:30 20749

原创解释和编译的区别

这方面在知乎上回答的很多，也很到位. 先给出参考的文献和链接如下：百度百科，知乎，《c语言程序设计》，Bluestorm’s Space，平凡生活当然最主要的还是参考了课本《c语言程序设计》，毕竟课本比较专业先从基本概念入手机器语言：电脑CPU指令体系（系统）。CPU指令系统中包含大量的由0和1组成的指令码（序列），这些指令码可以直接由计算机识别。如下面两条指令码： 100000

2017-12-02 19:05:23 10192

原创 MCMC和贝叶斯统计在宇宙微波背景辐射（CMB）中应用

此文为笔者笔记，描述的是MCMC（马尔科夫链蒙特卡洛）方法和贝叶斯统计在CMB中的应用，描述的只是笔者的观点，有待修改（更关键的是，笔者主要描述的是CosmoMC程序相关的东西）。进入正题假设待测宇宙学参数只有六个，那么描述的则是六维空间。在参数ini文件中（params_CMB_defults.ini），每个参数都有一个事先估计的范围和最可几的值(路径为：~/programs/CosmoMC/ba

2017-12-02 16:00:20 1406

Mr_Cat123的博客