1 00:00:02,765 --> 00:00:04,488 大家好,我是 Shelly Kelly 2 00:00:04,488 --> 00:00:12,337 我来展示一下如何加载多元检测器采集的数据 3 00:00:12,337 --> 00:00:15,426 在 SSRL 采集,需要加载到 Athena 4 00:00:16,296 --> 00:00:20,068 以及如何简单地预处理 5 00:00:21,908 --> 00:00:25,665 第一步,打开 Athena 6 00:00:26,365 --> 00:00:30,358 到文件 (File),选择加载数据 7 00:00:32,212 --> 00:00:35,972 找到数据所在路径 8 00:00:41,761 --> 00:00:43,761 数据在这儿 9 00:00:45,153 --> 00:00:49,283 先打开第一个 10 00:00:52,193 --> 00:00:54,624 或者四个一起加载了 11 00:00:59,501 --> 00:01:02,508 这是多元检测器得到的数据 12 00:01:02,508 --> 00:01:07,064 这里是选择数据加载列的窗口 13 00:01:07,440 --> 00:01:11,348 这里是谱图预览窗口 14 00:01:11,825 --> 00:01:14,300 (调整下位置) 15 00:01:15,040 --> 00:01:17,713 这看着不像数据 16 00:01:17,713 --> 00:01:23,668 说明我们这么选不对 17 00:01:24,172 --> 00:01:30,733 你可以在右边看到纯文本数据文件长什么样 18 00:01:33,873 --> 00:01:37,576 每一列的名字位于这里(数字上方) 19 00:01:37,576 --> 00:01:40,494 请求能量的数值在第一列 20 00:01:41,501 --> 00:01:46,433 实时计数总和在下一列 21 00:01:46,638 --> 00:01:48,638 实时时钟数值在这里 22 00:01:48,907 --> 00:01:52,081 这之后的都是多元检测器得到的数据 23 00:01:52,474 --> 00:02:03,750 看起来想要的信号数据范围在... 24 00:02:04,217 --> 00:02:06,217 这是输入计数速率 25 00:02:08,025 --> 00:02:11,375 看下这一排排都是啥 26 00:02:11,861 --> 00:02:24,225 看起来我们需要的是 SCA1_2, 3, 4, 5, 6 27 00:02:24,225 --> 00:02:27,650 可以算下有多少列数据我们需要 28 00:02:27,650 --> 00:02:35,650 数下... (1, 2, 3, 4, 5, 6...) 29 00:02:35,650 --> 00:02:37,650 第五列是 I0 30 00:02:37,650 --> 00:02:42,250 从第六到第十一列(是信号) 31 00:02:43,150 --> 00:02:46,150 5 对应 I0 32 00:02:46,150 --> 00:02:48,650 也许是 11? 33 00:02:48,650 --> 00:02:53,425 这看起来像是反着的吸收边 34 00:02:54,022 --> 00:02:57,100 反一下,正常些了 35 00:02:57,100 --> 00:03:02,665 由于是荧光谱,不要勾选自然对数 36 00:03:03,125 --> 00:03:07,200 现在看起来是一个噪音较大的谱图 37 00:03:07,884 --> 00:03:12,500 数据类型是 μ(E),能量单位是 eV 38 00:03:12,500 --> 00:03:13,850 都是正确的 39 00:03:13,850 --> 00:03:18,725 如果想预览重新勾选的谱图,点击 replot 40 00:03:19,100 --> 00:03:26,500 来看下有多少检测器 41 00:03:26,500 --> 00:03:40,550 数数数到……30 42 00:03:40,853 --> 00:03:44,650 可以换到下一列数据预览 43 00:03:44,650 --> 00:03:45,825 最好检查下 44 00:03:45,825 --> 00:03:48,700 因为有时候检测器是坏的 45 00:03:48,700 --> 00:03:51,850 最好都每个看下 46 00:03:51,850 --> 00:03:57,475 确保每个输入的信号看起来都比较正常 47 00:03:57,475 --> 00:04:01,625 比如不能用第 15 列数据(看起来不像数据) 48 00:04:02,143 --> 00:04:04,143 那就跳过这列 49 00:04:04,143 --> 00:04:09,600 备注下: 跳过第 15 列 50 00:04:10,175 --> 00:04:13,600 第 16 列也跳过 51 00:04:16,675 --> 00:04:20,125 第 17 个还行 52 00:04:22,850 --> 00:04:24,850 第 18 个不行 53 00:04:27,195 --> 00:04:29,195 第 19 个还行 54 00:04:31,475 --> 00:04:36,975 第 20 个还行,第 21 个不行 55 00:04:39,116 --> 00:04:42,775 继续到 22,还行 56 00:04:43,250 --> 00:04:45,675 第 23 列还行 57 00:04:46,375 --> 00:05:00,450 (继续数下去……) 58 00:05:06,375 --> 00:05:11,125 第 30 个不行 59 00:05:12,850 --> 00:05:14,850 31 还行 60 00:05:22,800 --> 00:05:41,045 (继续数……) 61 00:05:41,045 --> 00:05:45,775 第 41 列看着不对,去掉 62 00:05:45,775 --> 00:05:51,875 我们从 5 开始,再到 10 63 00:05:52,525 --> 00:05:57,000 那就得到第 40 列 64 00:05:57,021 --> 00:06:00,175 之后的是不一样的数据 65 00:06:00,450 --> 00:06:04,698 回到最一开始 66 00:06:12,175 --> 00:06:15,325 那现在勾选上我们所需要的数列 67 00:06:15,325 --> 00:06:33,200 从 11 开始……(去掉之前看到不合格的数据) 68 00:06:33,307 --> 00:06:35,850 一直到 30 69 00:06:38,235 --> 00:06:44,875 可以看到信噪比有所提升 70 00:06:44,875 --> 00:06:49,650 下面的数学表达式也是同步更新的 71 00:06:49,650 --> 00:06:51,650 点击 OK 72 00:06:54,713 --> 00:06:58,800 这是谱图预览 73 00:06:58,800 --> 00:07:04,100 能量 (E),k,傅里叶变换(R),q 空间等等 74 00:07:04,100 --> 00:07:11,575 再把剩下三个加载进来 75 00:07:15,333 --> 00:07:20,450 之前我们选择的数据数列直接默认选择了 76 00:07:21,029 --> 00:07:26,100 还是有必要再确认下 77 00:07:26,100 --> 00:07:33,675 但现在暂时假定下这样选择是可以的 78 00:07:33,675 --> 00:07:37,250 要选上这些 79 00:07:45,208 --> 00:07:49,700 可以看到数据们是如何命名的 80 00:07:49,700 --> 00:07:52,225 这是锰 K 边数据 81 00:07:52,225 --> 00:07:55,525 看起来是个水溶液样品 82 00:07:55,525 --> 00:07:58,450 选了两个(光斑)位置来采集数据 83 00:07:58,675 --> 00:08:02,050 所以这里扫了两圈 84 00:08:02,350 --> 00:08:05,200 001 和 002 85 00:08:05,200 --> 00:08:09,700 还有两圈是在第二个位置上获得的 86 00:08:09,700 --> 00:08:11,700 001 和 002 87 00:08:11,700 --> 00:08:21,775 勾选上所有的数据,点击紫色 E 键展示所有数据 88 00:08:22,550 --> 00:08:25,575 第一眼看感觉它们相似 89 00:08:25,575 --> 00:08:31,079 放大到吸收边来看它们是否完全一致 90 00:08:36,395 --> 00:08:40,175 它们看起来确实很相近 91 00:08:41,014 --> 00:08:46,825 看不出明显差别 92 00:08:46,825 --> 00:08:54,175 有可能两个不同位置上的数据会有点不同 93 00:08:54,975 --> 00:08:58,100 看起来确实很像 94 00:08:58,100 --> 00:09:01,575 去 k 空间检查 95 00:09:01,575 --> 00:09:08,375 鉴于噪音较大,很难判断是否有区别 96 00:09:09,025 --> 00:09:11,025 再检查这两个 97 00:09:11,650 --> 00:09:14,675 很相似 98 00:09:14,675 --> 00:09:18,250 这里(信号)有点点升高 99 00:09:18,250 --> 00:09:22,913 也许有改变,但是不多 100 00:09:23,489 --> 00:09:27,275 看起来不错 101 00:09:29,850 --> 00:09:37,366 下一步是检查数据是否对齐 102 00:09:37,366 --> 00:09:43,350 吸收边看起来重叠得较好 103 00:09:43,350 --> 00:09:47,400 放大点观察 104 00:09:47,800 --> 00:09:54,525 看起来吸收边也没有偏移 105 00:09:54,525 --> 00:09:58,768 如果担心谱图有偏移 106 00:09:58,768 --> 00:10:02,225 你可以选择对齐数据 107 00:10:02,612 --> 00:10:08,150 到主窗口选项,点击,选择下方对齐数据 (Align data) 108 00:10:09,375 --> 00:10:22,300 我们现在是以 (35) 002 文件为参考,移动 (35) 001 文件使之对齐 109 00:10:22,300 --> 00:10:27,800 这是光滑后的一阶微分谱图 110 00:10:27,800 --> 00:10:30,375 两张谱图重叠程度高 111 00:10:30,375 --> 00:10:34,175 你可以选择自动对准 112 00:10:34,175 --> 00:10:38,275 试着让两个数据更好地对齐 113 00:10:38,975 --> 00:10:41,900 确实有点偏移,0.3 eV 114 00:10:41,900 --> 00:10:44,350 但数据噪音也大 115 00:10:44,350 --> 00:10:48,475 所以我更倾向于不移动 (0 eV) 116 00:10:48,475 --> 00:10:51,575 不使用对准这个功能 117 00:10:51,575 --> 00:10:55,825 你可以点击下一个数据 118 00:10:55,825 --> 00:11:01,850 需要对齐的数据名字就在这里了,文件名37点2 119 00:11:01,850 --> 00:11:05,550 你也可以调整第四个数据 120 00:11:05,800 --> 00:11:10,000 下面的按键是用来手动调整的 121 00:11:10,000 --> 00:11:13,550 如果想让谱图向左移动 0.5 eV 122 00:11:13,550 --> 00:11:18,100 那红色谱图就会向左移动 0.5 eV 123 00:11:18,700 --> 00:11:20,406 把它调整回来 124 00:11:20,406 --> 00:11:22,406 你都可以使用这些按键 125 00:11:22,406 --> 00:11:26,950 你也可以直接手动输入你已经知道的偏移数值 126 00:11:26,950 --> 00:11:28,950 重新加载 (replot) 127 00:11:28,950 --> 00:11:30,950 移动后的谱图就会出现 128 00:11:30,950 --> 00:11:34,425 调整回初始状态 129 00:11:34,425 --> 00:11:37,025 回到主菜单 130 00:11:37,025 --> 00:11:39,436 数据彼此都对齐 131 00:11:39,436 --> 00:11:41,825 把数据汇聚了 (merge) 132 00:11:41,825 --> 00:11:43,600 按键在上方 133 00:11:43,600 --> 00:11:45,600 汇聚它们的 μ(E) 数据 134 00:11:47,000 --> 00:11:50,650 蓝线是汇聚之后的谱图 135 00:11:50,650 --> 00:11:58,475 红线是基于被汇聚原始数据自动算出的标准方差 136 00:11:58,475 --> 00:12:04,225 勾选上新的谱图,点击紫色 E 键显示所有谱图 137 00:12:04,225 --> 00:12:06,225 或是是 k 空间 138 00:12:07,050 --> 00:12:09,050 看起来不错 139 00:12:09,050 --> 00:12:12,475 把 k 空间展示上限调到 12 140 00:12:13,375 --> 00:12:16,175 数据遍布整个窗口 141 00:12:16,175 --> 00:12:22,550 点击黄色 k 键,观察新的 k 空间谱图 142 00:12:22,550 --> 00:12:26,025 之前原始数据的噪音有所改善 143 00:12:26,300 --> 00:12:28,150 是好事 144 00:12:28,150 --> 00:12:34,275 现在来判断用于傅里叶变换需要的范围 145 00:12:34,275 --> 00:12:40,825 现在傅里叶变换的范围在这里可以看到 146 00:12:40,825 --> 00:12:43,725 从 3 到 10 147 00:12:44,250 --> 00:12:49,325 勾选窗函数 (window) 再展示 148 00:12:49,325 --> 00:12:57,075 就可以看到窗函数从 3 开始到 10 结束 149 00:12:57,075 --> 00:12:59,600 进行傅里叶变换 150 00:12:59,600 --> 00:13:02,157 如图所示 151 00:13:04,625 --> 00:13:06,975 来看下 152 00:13:06,975 --> 00:13:13,450 这是傅里叶变换的后强度 153 00:13:13,450 --> 00:13:18,625 来做个反傅里叶变换到 q 空间 154 00:13:18,625 --> 00:13:27,050 那从比 R = 1 略低的地方开始 155 00:13:27,500 --> 00:13:29,500 改成 0.8 156 00:13:30,100 --> 00:13:32,650 上限取到 3 合适 157 00:13:32,650 --> 00:13:35,900 (我们)试图把所有峰的信号都框进去 158 00:13:35,900 --> 00:13:37,900 点击 kq 159 00:13:38,375 --> 00:13:46,800 红线是反傅里叶变换的信号 160 00:13:46,800 --> 00:13:53,100 也就是从 R 空间中 0.8 到 3 之间的峰得到的 161 00:13:53,100 --> 00:13:56,954 看起来(反傅里叶变换)得到信号跟原始信号比较接近 162 00:13:57,400 --> 00:14:04,229 我会认为在 k 空间 8 到 10 之间的信号是真实的(而非噪音) 163 00:14:05,059 --> 00:14:10,450 复制数据 164 00:14:11,025 --> 00:14:15,125 右击数据,点击复制 165 00:14:15,125 --> 00:14:17,125 再复制一份 166 00:14:18,275 --> 00:14:22,500 如果傅里叶变换只到 8 167 00:14:23,120 --> 00:14:26,234 如果 8 - 10 之间的信号是真实的 168 00:14:26,234 --> 00:14:31,050 会看到不一样的傅里叶变换谱图(R空间) 169 00:14:31,050 --> 00:14:33,750 那就这个定为 8 170 00:14:34,500 --> 00:14:37,900 第二个设为 9 171 00:14:37,900 --> 00:14:40,425 第三个设为 10 172 00:14:41,125 --> 00:14:43,125 勾选这三个数据 173 00:14:43,700 --> 00:14:48,775 点击紫色 R 键 174 00:14:49,675 --> 00:14:57,300 你可以看到这三个数据都没有出现新的峰 175 00:14:57,300 --> 00:15:00,125 观察实部 176 00:15:00,650 --> 00:15:05,950 设 Rmax 为 10,可以看到后面的信号 177 00:15:05,950 --> 00:15:09,275 现在能看到噪音水平了 178 00:15:09,275 --> 00:15:11,800 绿线有一些噪音 179 00:15:11,800 --> 00:15:15,325 有一点信号一直到 R = 6 到 8 区间 180 00:15:15,325 --> 00:15:22,000 但是最左边的峰(R值小的地方)出现在三组数据 181 00:15:22,000 --> 00:15:25,950 如果做反傅里叶变换 182 00:15:25,950 --> 00:15:34,425 可以看到之前选择 8 的谱图在 8 - 10 之间几乎没有波动 183 00:15:34,425 --> 00:15:39,643 但是选 10 的谱图就会出现这里的信号 184 00:15:39,643 --> 00:15:41,867 选择更大的 k 值范围 185 00:15:41,867 --> 00:15:44,299 可以在 R 空间给出更多的信息 186 00:15:44,299 --> 00:15:47,213 更容易区分开的峰(代表傅里叶变换后强度) 187 00:15:47,213 --> 00:15:55,898 绿线相比蓝线的峰分开的更清楚,更窄一些 188 00:15:55,898 --> 00:15:58,500 但是并没有新的信息出现 189 00:15:58,500 --> 00:16:05,007 所以我认为 k 到 10 的数据就可以了 190 00:16:05,007 --> 00:16:07,415 来看下范围是否还可以再长一些 191 00:16:07,415 --> 00:16:11,994 展示 k 空间的数据 192 00:16:12,786 --> 00:16:16,435 看下能不能到 k = 12 193 00:16:16,566 --> 00:16:20,084 把这个从 8 改到 10 194 00:16:20,848 --> 00:16:24,413 11 还有 12 195 00:16:25,601 --> 00:16:28,492 去到 R 空间 196 00:16:28,492 --> 00:16:32,814 看起来不错 197 00:16:32,814 --> 00:16:34,814 没有新的信号 198 00:16:34,814 --> 00:16:38,698 k 选到 12 时可能会看到一些噪音 199 00:16:38,698 --> 00:16:43,083 (R 值高的区域)开始出现一些信号 200 00:16:43,083 --> 00:16:48,090 傅里叶变换的实部也能看到一些震荡 201 00:16:48,090 --> 00:16:57,680 可以使用 k 的三次方来强调这一部分的信号 202 00:16:57,680 --> 00:17:01,640 这里就能看到更明显的噪音 203 00:17:01,640 --> 00:17:06,707 使用 k 的三次方权重后,到 k 12 的噪音会更明显 204 00:17:06,707 --> 00:17:10,664 到反傅里叶变换 205 00:17:14,879 --> 00:17:19,547 (确认下参数是否正确……) 206 00:17:23,451 --> 00:17:25,451 把这个数据也加进来 207 00:17:25,451 --> 00:17:31,343 看起来这里也有一点震荡,不是很好 208 00:17:34,455 --> 00:17:39,264 把它的 q 空间和实验数据(k 空间)进行比较 209 00:17:39,264 --> 00:17:42,602 在 k = 12 附近确实信噪比差了 210 00:17:42,602 --> 00:17:45,488 但也不是特别差 211 00:17:46,648 --> 00:17:53,267 这就是如何决定用于傅里叶变换 k 的上限 (kmax) 212 00:17:53,267 --> 00:17:57,793 我们还没说 kmin (下限) 213 00:17:57,793 --> 00:17:59,793 仅仅使用默认值 3 214 00:17:59,793 --> 00:18:01,793 在这里 215 00:18:02,291 --> 00:18:07,072 我们可以利用指示 216 00:18:07,072 --> 00:18:10,155 在图中 k = 3 的地方 217 00:18:12,051 --> 00:18:14,051 大约是 3 218 00:18:18,642 --> 00:18:21,923 这里 标记出来 219 00:18:21,923 --> 00:18:25,374 但是现在是2.6,改成 3 220 00:18:25,892 --> 00:18:27,892 再次点击 k 键 221 00:18:27,892 --> 00:18:29,892 现在好了 222 00:18:30,240 --> 00:18:33,295 你可以让同样的数据在能量区间内展示 223 00:18:35,671 --> 00:18:39,321 这个黄圈是 E0,左上角这里 224 00:18:39,321 --> 00:18:41,321 也是定义了 k = 0 的地方 225 00:18:41,838 --> 00:18:48,628 k = 3 的位置离吸收边不远 226 00:18:48,628 --> 00:18:52,192 只是往 E0 后面移了一些 227 00:18:52,192 --> 00:18:55,728 甚至没有翻过吸收最强的区域 228 00:18:56,872 --> 00:18:59,208 这是个很低的能量 229 00:19:00,965 --> 00:19:06,676 试试 k = 0 230 00:19:06,676 --> 00:19:10,523 出于检查的目的,因为确实看起来太近了 231 00:19:11,015 --> 00:19:13,015 k = 0 232 00:19:16,690 --> 00:19:20,961 k = 3 233 00:19:22,178 --> 00:19:27,100 这好像看起来是能量 E0 + 3 而不是 k = 3 的地方 234 00:19:27,581 --> 00:19:30,551 看起来是个错误 235 00:19:30,551 --> 00:19:32,551 再来一次 236 00:19:36,633 --> 00:19:41,640 改回到能量区间 237 00:19:41,640 --> 00:19:43,640 上限 300 eV 238 00:19:43,640 --> 00:19:47,355 找下 k = 3 239 00:19:47,657 --> 00:19:49,657 这儿呢 240 00:19:50,363 --> 00:19:53,408 这是 k 空间开始的地方 241 00:19:53,408 --> 00:19:57,482 在吸收边之上,但是在这个剧烈震荡的区间 242 00:19:57,482 --> 00:19:59,482 所以在 k 空间 243 00:20:00,233 --> 00:20:02,233 就是这个向下的震荡 244 00:20:02,574 --> 00:20:08,118 在能量区间就在这里 245 00:20:08,882 --> 00:20:17,397 背景扣除麻烦的地方在于背景改变得很快 246 00:20:17,397 --> 00:20:19,397 在经过吸收边附近 247 00:20:19,397 --> 00:20:21,216 通常来说 248 00:20:21,216 --> 00:20:25,799 你希望傅里叶变换要从比较平的曲线开始 249 00:20:25,799 --> 00:20:28,712 这里的震荡不太容易受背景选择的影响 250 00:20:28,712 --> 00:20:36,803 尤其是在一开始想要搞清楚哪里是好的或是坏的信号 251 00:20:37,199 --> 00:20:42,206 想看到容易受背景方程影响的程度一个简单方法是 252 00:20:42,206 --> 00:20:46,449 就是改变 E0 值,也是 k = 0 253 00:20:46,449 --> 00:20:48,882 也是背景开始的地方 254 00:20:49,589 --> 00:20:53,437 我会复制数据成三份 255 00:20:53,437 --> 00:20:56,124 把所有参数设成一样 256 00:20:57,708 --> 00:21:00,254 傅里叶变换参数 257 00:21:00,254 --> 00:21:03,338 这些参数也是 258 00:21:05,403 --> 00:21:08,967 现在来改变 E0 的数值 259 00:21:08,967 --> 00:21:13,069 大概改变 0.5 eV 260 00:21:13,918 --> 00:21:17,114 点击 E 键 261 00:21:17,114 --> 00:21:20,311 看到 E0 想高能方向移动了一点 262 00:21:20,311 --> 00:21:22,687 背景方程也有点改变 263 00:21:23,091 --> 00:21:25,091 再来一次 264 00:21:25,091 --> 00:21:27,091 向高能移动 2 eV 265 00:21:27,091 --> 00:21:29,448 6552 eV 266 00:21:29,448 --> 00:21:31,448 点击 E 键 267 00:21:34,964 --> 00:21:37,482 我没看到有明显的改变 268 00:21:37,482 --> 00:21:42,432 现在同时展示三组 k 空间数据 269 00:21:43,423 --> 00:21:47,439 可以看到在 k 值低的区域有点不同 270 00:21:47,439 --> 00:21:50,976 背景移除方程和 χ(k) 数据 271 00:21:50,976 --> 00:21:52,927 之所以看起来(差别)没那么明显 272 00:21:52,927 --> 00:21:54,927 因为我们用的是 k 的二次方权重 273 00:21:54,927 --> 00:21:56,927 如果用 k 的一次方权重 274 00:21:56,927 --> 00:22:00,000 可以看到这个差异更加明显 275 00:22:00,792 --> 00:22:02,792 在我看来 276 00:22:02,792 --> 00:22:12,277 肯定不想在一开始处理数据的时候使用这部分 277 00:22:12,277 --> 00:22:13,663 因为它改变得很明显 278 00:22:13,663 --> 00:22:16,492 那还是可能把 kmin 改到 2.5 279 00:22:17,114 --> 00:22:23,819 我们来把 E0 调整的更多一些来看到更明显的差异 280 00:22:23,819 --> 00:22:25,819 现在这个差异不明显 281 00:22:25,819 --> 00:22:28,430 把这个改回来 282 00:22:29,307 --> 00:22:32,814 把这个改成 6552 eV 283 00:22:32,814 --> 00:22:35,700 把这个改成 6554 eV 284 00:22:36,718 --> 00:22:38,718 点击 k 键 285 00:22:38,718 --> 00:22:45,063 这下就可以直观地看到在 k 值低的区间差异了 286 00:22:45,063 --> 00:22:47,892 这也是背景方程差异比较大的地方 287 00:22:48,059 --> 00:22:53,182 点击 E 键观察背景方程是什么样子 288 00:22:54,200 --> 00:22:57,680 可以看到这里更窄了 289 00:22:59,858 --> 00:23:02,857 那个就像这样 290 00:23:04,724 --> 00:23:11,711 那么 kmin = 2.5 算是个不错的选择 291 00:23:11,711 --> 00:23:13,833 展示所有的数据 292 00:23:13,833 --> 00:23:16,662 设置三组数据的 k 都从 2.5 到 10 293 00:23:16,662 --> 00:23:18,662 点击 R 键 294 00:23:18,662 --> 00:23:20,961 傅里叶变换后的强度 295 00:23:21,386 --> 00:23:26,251 可以看到这里出于 E0 不同,有点偏移 296 00:23:26,561 --> 00:23:33,833 但总的来说峰形接近 297 00:23:33,833 --> 00:23:37,991 是个不错的开端 298 00:23:38,444 --> 00:23:40,452 查看实部谱图 299 00:23:40,452 --> 00:23:42,452 主要还是由不同 E0 带来的偏移 300 00:23:42,885 --> 00:23:44,752 这里有点强度的差异 301 00:23:44,752 --> 00:23:47,194 我好奇的一点是 302 00:23:47,581 --> 00:23:52,249 这组数据中有什么样的邻近原子 303 00:23:53,117 --> 00:24:00,735 我猜测有氧原子还有质量更大的邻近原子 304 00:24:02,489 --> 00:24:06,647 我们可以再对比一下 305 00:24:06,647 --> 00:24:12,899 设置所有参数对这三组数据一样 306 00:24:13,408 --> 00:24:16,690 在 R 空间中它们都一样了 307 00:24:17,058 --> 00:24:26,308 我们现在要以 k 不同的权重来对比谱图 308 00:24:26,308 --> 00:24:31,230 留意这些峰的高度 309 00:24:31,230 --> 00:24:34,427 选择 k 的二次方权重 310 00:24:35,898 --> 00:24:38,529 k 的三次方权重 311 00:24:38,896 --> 00:24:42,693 峰的相对高度几乎没有改变 312 00:24:42,693 --> 00:24:49,420 我没有看到在 k 的三次方权重下有任何一个峰的高度更突出 313 00:24:49,844 --> 00:24:53,210 选回 k 的二次方权重 314 00:24:53,748 --> 00:24:56,690 为了更好的对比 315 00:24:56,690 --> 00:25:01,075 可以把三个不同权重的图放在一起 316 00:25:01,075 --> 00:25:05,205 再缩放到第一个峰高度一样 317 00:25:05,205 --> 00:25:06,591 开整 318 00:25:06,591 --> 00:25:12,588 利用 kw (独立 k 权重)来展示 R 空间的图像 319 00:25:12,696 --> 00:25:22,913 这里我们使用任意 k 权重右边的方框,填 1 320 00:25:23,366 --> 00:25:25,366 下一组数据这里改成 2 321 00:25:25,912 --> 00:25:28,079 同理,最后个改成 3 322 00:25:28,543 --> 00:25:31,174 点击紫色 R 键 323 00:25:31,570 --> 00:25:38,331 需要点来点去下…… 324 00:25:39,377 --> 00:25:41,046 好了 325 00:25:41,046 --> 00:25:44,320 不同 k 权重的图放在了一起 326 00:25:44,320 --> 00:25:47,807 由于三者强度差距太大,几乎没法比较 327 00:25:48,203 --> 00:25:53,295 对最矮的那个数据乘以 30 328 00:25:53,295 --> 00:25:58,161 诶点错了 329 00:25:58,444 --> 00:26:01,527 对第一组数据乘以 30 330 00:26:01,527 --> 00:26:04,724 第二组数据乘以 5 331 00:26:05,007 --> 00:26:07,007 再点击紫色 R 键 332 00:26:07,411 --> 00:26:09,411 看起来接近了一些 333 00:26:10,466 --> 00:26:14,851 看起来第一组数据乘多了 334 00:26:14,851 --> 00:26:16,464 改成 25 335 00:26:16,464 --> 00:26:18,464 点击 R 键 336 00:26:19,321 --> 00:26:21,321 改成 20 337 00:26:21,321 --> 00:26:24,413 现在差不多了 338 00:26:24,413 --> 00:26:34,229 可以看到这三个峰的强度在不同 k 权重情况下相当接近 339 00:26:34,512 --> 00:26:39,462 那说明这三个峰所包含的邻近原子很可能是与第一壳层的邻近原子是一个类型 340 00:26:39,936 --> 00:26:43,196 如果说两个原子不一样,比如氧原子 341 00:26:43,196 --> 00:26:45,587 第二层是锰原子 342 00:26:45,587 --> 00:26:47,638 会看到第二壳层的峰强度变强 343 00:26:47,638 --> 00:26:51,711 这是因为锰原子在 k 值高的区域震荡更强烈 344 00:26:52,277 --> 00:26:58,981 所以这部分的信号会被锰原子主宰 345 00:26:59,406 --> 00:27:01,406 散射更多 346 00:27:01,782 --> 00:27:06,761 这个数据说明是有氧原子在不同的壳层,并且分得开 347 00:27:06,931 --> 00:27:12,333 数据中可以看到拍频的现象 348 00:27:12,673 --> 00:27:19,802 而不是一层氧原子和另一层金属原子 349 00:27:22,163 --> 00:27:28,514 这就是部分的 Athena 预数据处理 350 00:27:31,541 --> 00:27:33,541 感谢收听!