最近为了研究期货量化、以及不同策略对期望值的影响，我在 Win11 上做了一轮安装测试，下面记录一下我的环境搭建过程 ✨📈 🐍 我的 Python 环境用 Miniforge 搭建安装时建议选择 “Just me” 👤✅这样更省心，能减少权限问题，也更不容易和系统里其他 Python 发生“默认解释器抢占”的混乱情况。 🧰 安装后配置（研究/回测环境）下面是我整理后的命...

最近在个人电脑上安装21cmFAST宇宙学演化模拟程序，这里简要记录安装过程以及遇到的问题： 首先安装conda环境12# 我的系统是Arch linuxyay -S miniconda3 安装后发现找不到conda命令，检查才知道环境变量设置不对，可以直接编辑.bashrc文件，我使用的是zsh终端，因此我修改.zshrc，文件最后加环境变量： 1export PATH="/o...

傅里叶变换(FFT)在各种巡天中经常被应用，其中，计算星系数密度的三维功率谱需要使用三维傅里叶变换，这其中需要注意归一化问题。 假设星系数密度涨落是$F(r)$，其中$r = (x,y,z)$ 傅里叶变换： 连续傅里叶变换对：\begin{align} F(\vec{k}) = \int F(\vec{r}) e^{-i\vec{k} \cdot \vec{r}} {\r...

在宇宙学中，星系数密度功率谱的计算方法一般是直接对三维坐标 $(x,y,z)$ 密度涨落进行FFT，这样可以得到三维 $(k_x, k_y, k_z)$ 傅里叶变换，但是为了更好的观察数据，一般会把三维FFT结果平均投影到以 $k = \sqrt{k_x^2 + k_y^2 + k_z^2}$ 为半径的球壳上，从而得到一维功率谱 $P(k)$ 即monopole 功率谱 $P_0(...

最近在使用 Jupyter Notebook 时遇到了一个奇怪的问题：更新 jupyter 和 getdist 后，matplotlib 的 plt 绘图突然无法显示。经过排查，发现是 Matplotlib 后端配置被 getdist 覆盖 导致的兼容性问题。 环境配置 GetDist：1.6.1 Matplotlib：3.3.4 Jupyter Notebook：1.0.0 问题复现12...

今天再补一个关于camb计算线性/非线性功率谱的坑。 camb在初始化时会有一个 pars.NonLinear 设置，设为 model.NonLinear_both 表示同时计算线性/非线性功率谱， model.NonLinear_none 表示计算线性功率谱，起初我以为设为 model.NonLinear_none 就可以得到的 results.get_matter_p...

今天利用camb计算宇宙的dark matter物质功率谱遇到一个很费解的问题，函数get_matter_power_interpolator()两个参数官方文档是这样写的： hubble_units – if true, output power spectrum in (({\rm Mpc}/h)^{3}) units, otherwise ({\rm Mpc}^{3}) k...

最近计算宇宙学的21cm 强度映射3D功率谱，我主要计算的多级功率谱。我使用numpy.fft计算，其中对于个数为偶数的数组tb，假设tb.shape=(Nx,Ny,Nz)，numpy.fft返回的数组tb_fft对应情况(以第一维度举例)： tb_fft[0]：0频率，表示输入信号的平均值，一般直接扔掉 tb_fft[1:Nx//2]：正频率 tb_ff...

最近在用matplotlib画图，需要在每个图上重新分成子图，也就是嵌套画图，记录一下 123456789101112131415161718192021222324252627282930313233343536373839404142434445from matplotlib import pyplot as pltfrom matplotlib.gridspec import Grid...