一、按概率密度函数生成样本 参照python – How to simulate from an […]
post
已知\(\chi^2\),求与均值相差多少个标准差
假设原始的数据符合高斯分布,自由度为df,求卡方统计量\(\chi^2\)距离均值 \(\mu\) 多少个标准 […]
✿ 阅读数:2,394  分类:文章已知\(\chi^2\),求与均值相差多少个标准差有 1 条评论
使用EMCEE Python模块运行MCMC
In [1]: # 导入写好的似然函数 import sys sys.path.append('.. […]
✿ 阅读数:2,501  分类:文章使用EMCEE Python模块运行MCMC有 2 条评论
用python求不定积分和导数的解析方法和数值方法
解析求导 from sympy import diff x = sym.symbols('x […]
✿ 阅读数:1,971  分类:文章 给用python求不定积分和导数的解析方法和数值方法留言
scipy.integrate解ode的例子
下面是一个阻尼振动的例子 import numpy as np from scipy.integrate im […]
✿ 阅读数:1,409  分类:文章 给scipy.integrate解ode的例子留言
sympy中给算式赋值的例子
evalf是取浮点数,subs是替代。下面是两个例子。 import sympy omega, t = sym […]
✿ 阅读数:1,294  分类:文章 给sympy中给算式赋值的例子留言
用级数展开避免0处发散
当函数的分母随自变量变换,绝对值越来越小时,函数可能会不稳定。在这种情况下,用级数展开来近似这个函数,可以避免 […]
✿ 阅读数:1,333  分类:文章 给用级数展开避免0处发散留言
用GetDist分析马尔科夫链
一、样本的读取 1、从文件夹中读取 使用CosmoMC,cobaya,montepython等工具生成链时,输 […]
✿ 阅读数:4,251  分类:文章用GetDist分析马尔科夫链有 6 条评论
关于哈勃矛盾 (H_0 tension) 的梳理
测量H_0的方法有很多。比较精确的测量方法有两种,宇宙微波背景辐射 (CMB) + LCDM 和距离阶梯。二者 […]
✿ 阅读数:1,637  分类:文章 给关于哈勃矛盾 (H_0 tension) 的梳理留言