最近為了研究期貨量化、以及不同策略對期望值的影響，我在 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...