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
感谢收听！