add: trapped electron fit

README.md

@@ -80,8 +80,32 @@ Dose Estimation by Simulation of China Space Station
 <div align=center><img src="docs/spectra-ep.png" style="max-width: 80%;"><br /><span style="color: gray;font-size: 12px;">AE8/AP8模型结果</span></div>
 <div align=center><img src="docs/spectra-ep-9.png" style="max-width: 80%;"><br /><span style="color: gray;font-size: 12px;">AE9/AP9模型结果</span></div>
 
+### 辐射环境拟合
+#### 银河宇宙射线
+[^4]
+
+#### 太阳辐射
+[^5]
+
+#### 俘获辐射带
+1. 电子
+    根据文献，在轨道较高（$L>2.5$）的内带中，电子能谱中有$\sim64\%$符合指数分布，$\sim4\%$符合幂律分布。[^6]
+
+    但在本例中，轨道较低（$L\sim1.06$），因此使用指数分布$j = j_0e^{-\frac{E}{E_0}}$和幂律分布$j=j_0E^{-\alpha}$进行拟合，随后使用`Origin Lab`对两者进行拟合比较，最终选用指数分布模型$j_0=1.333680\times10^6,\ E_0=0.0824$。
+
+    <div align=center>
+
+    | Model |   $R^2$  |    BIC    | BIC diff |
+    |:-----:|:-------:|:---------:|:--------:|
+    |   Ex  | 0.98107 | 402.50208 |     0    |
+    |  Pow  | 0.94279 | 423.51296 | 21.01088 |
+
+    </div>
+
+2. 质子
+
 ## 空间站结构
-1. 尺寸与分区[^4]  
+1. 尺寸与分区[^7]  
 <div align=center><img src="docs/size.webp" style="max-width: 50%;"></div>
 
 * 全长：16.6 m
@@ -101,8 +125,8 @@ Dose Estimation by Simulation of China Space Station
 2. 材料
 * 外壳： 3层 2mm 铝 + 10mm 芳纶 + 5mm 铝
         <div align=center><img src="docs/Shell.png" style="max-width: 40%;"></div>
-    * 5系铝合金[^5]
-    * 泰普龙[^6]
+    * 5系铝合金[^8]
+    * 泰普龙[^9]
         <div align=center><img src="docs/Taparan.png" style="max-width: 30%;"></div>
 * 填充：金属为主，使得总重大致相当
 * 内部：空气
@@ -126,6 +150,9 @@ Dose Estimation by Simulation of China Space Station
 [^1]: [中国空间站轨道参数](http://www.cmse.gov.cn/gfgg/zgkjzgdcs/)
 [^2]: 程彭超,闵锐.近地空间辐射环境与防护方法概述[J].辐射防护通讯,2017,37(01):14-21.
 [^3]: Bourdarie S, Xapsos M. The near-earth space radiation environment[J]. IEEE transactions on nuclear science, 2008, 55(4): 1810-1832.
-[^4]: [【知识点·航天】“天和”核心舱、“天宫”空间站和新一代载人飞船的最新知识点（干货版）](https://zhuanlan.zhihu.com/p/103709953)
-[^5]: [中铝造为“天和”号核心舱披“铠甲”壮“筋骨”](https://m.thepaper.cn/baijiahao_12484370)
-[^6]: [泰普龙产品相关知识](https://wenku.baidu.com/view/a2ecf93501d8ce2f0066f5335a8102d276a261cd.html)
+[^4]:
+[^5]:
+[^6]: Zhao H, Johnston W R, Baker D N, et al. Characterization and evolution of radiation belt electron energy spectra based on the Van Allen Probes measurements[J]. Journal of Geophysical Research: Space Physics, 2019, 124(6): 4217-4232.
+[^7]: [【知识点·航天】“天和”核心舱、“天宫”空间站和新一代载人飞船的最新知识点（干货版）](https://zhuanlan.zhihu.com/p/103709953)
+[^8]: [中铝造为“天和”号核心舱披“铠甲”壮“筋骨”](https://m.thepaper.cn/baijiahao_12484370)
+[^9]: [泰普龙产品相关知识](https://wenku.baidu.com/view/a2ecf93501d8ce2f0066f5335a8102d276a261cd.html)

assets/trapped-9.csv

@@ -0,0 +1,4 @@
+0.1, 0.2, 0.4, 0.6, 0.8, 1, 2, 4, 6, 8, 10, 15, 20, 30, 50, 60, 80, 100, 150, 200, 300, 400, 700, 1200, 2000
+170.46630985001116, 148.7169777853078, 116.02958215502514, 91.68016573576256, 75.4004563167008, 52.44979368353229, 27.39576993169936, 11.534263356467031, 4.927074386668564, 2.503935512743802, 1.2207035405229754, 0.6406616768835681, 0.32153683341049205, 0.18332837912910016, 0.11496359892828469, 0.08025880550325193, 0.05970056915795244, 0.0384373301355559, 0.022821770094163085, 0.011082462960424989, 0.0047728913059455905, 0.0007325376524683292, 5.947420364867388e-05, 1.918902440238389e-07, 0.0
+0.04, 0.07, 0.1, 0.25, 0.5, 0.75, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, 5.5, 6, 6.5, 7, 8.5, 10
+781425.9929592956, 679139.2146049731, 322379.40275874606, 56156.789327450584, 6157.486909985416, 494.7242655747493, 77.04500465644205, 17.196451096282715, 5.771096902928527, 1.655431341146984, 0.5327835060623752, 0.17776818060047872, 0.05960754201838149, 0.02284581379206805, 0.008191147518478015, 0.0034509863771389026, 0.0010065728536974858, 0.00025474598119156806, 4.308092673195076e-06, 0.0, 0.0

utils/ae9ap9.py

@@ -24,8 +24,6 @@ _, ax = plt.subplots(1, 2, dpi=150, figsize=(16, 6))
 
 ax1, ax2 = ax[0], ax[1]
 
-print(proton)
-
 ax1.plot(E_proton, proton)
 ax1.set_ylim([1e-3, 5 * 1e2])
 ax1.set_yscale('log')