假设原始的数据符合高斯分布,自由度为df,求卡方统计量\(\chi^2\)距离均值 \(\mu\) 多少个标准 […]
作者:打开君
使用EMCEE Python模块运行MCMC
In [1]: # 导入写好的似然函数 import sys sys.path.append('.. […]
✿ 阅读数:1,167  分类:文章使用EMCEE Python模块运行MCMC有1条评论
用python求不定积分和导数的解析方法和数值方法
解析求导 from sympy import diff x = sym.symbols('x […]
✿ 阅读数:1,051  分类:文章 给用python求不定积分和导数的解析方法和数值方法留言
scipy.integrate解ode的例子
下面是一个阻尼振动的例子 import numpy as np from scipy.integrate im […]
✿ 阅读数:812  分类:文章 给scipy.integrate解ode的例子留言
sympy中给算式赋值的例子
evalf是取浮点数,subs是替代。下面是两个例子。 import sympy omega, t = sym […]
✿ 阅读数:654  分类:文章 给sympy中给算式赋值的例子留言
用级数展开避免0处发散
当函数的分母随自变量变换,绝对值越来越小时,函数可能会不稳定。在这种情况下,用级数展开来近似这个函数,可以避免 […]
✿ 阅读数:750  分类:文章 给用级数展开避免0处发散留言
用GetDist分析马尔科夫链
一、样本的读取 1、从文件夹中读取 使用CosmoMC,cobaya,montepython等工具生成链时,输 […]
✿ 阅读数:2,124  分类:文章用GetDist分析马尔科夫链有4条评论
关于哈勃矛盾 (H_0 tension) 的梳理
测量H_0的方法有很多。比较精确的测量方法有两种,宇宙微波背景辐射 (CMB) + LCDM 和距离阶梯。二者 […]
✿ 阅读数:796  分类:文章 给关于哈勃矛盾 (H_0 tension) 的梳理留言
英文论文同义替换
英文论文忌讳一直用同一个词。所以在看论文时要留心别人在表达相同意思的时候采用了哪些表达方式,然后再自己借鉴过来 […]
✿ 阅读数:940  分类:文章 给英文论文同义替换留言