0%

可选的练习

可选的练习

如果您有兴趣详细了解基带数据,则可以尝试以下内容。 这些主要是为了说明与相关处理无关的检查数据的一些工具。

制作基带通道的高分辨率频谱

这个实用程序m5spec可以读取基带数据并生成频谱。 该程序和包含频谱的输出文本文件具体信息可以参考DiFX参考手册。

这个程序可以在许多种VLBA基带数据(包括第一个数据集中的Mark5B数据)上工作,但不能在具有多个线程的VDIF数据上工作(例如,来自第二个数据集的数据)。 适用于数据集1数据的数据格式参数是Mark5B-2048-16-2,生成每个基带通道500个频谱点的完整命令是:

1
m5spec /home/difx_data/n6043/NRAO+301_0540_N6043_PT_No0008 Mark5B-2048-16-2 500 2000 spec.out

这里的数字2000是要处理的FFT帧数。 增加该数字将导致更多的平均值和信噪比更高的频谱。
可以更改这两个数字(频谱点数和要处理的FFT数),以了解这些参数如何更改输出。

gnuplot 程序可用于绘制数据。 首先在终端中启动gnuplot:

1
$ gnuplot

然后在提示符下,您可以发出绘图命令,例如:

1
p 'spec.out' u 1:2 w l

第1列(在1:2参数中标识)是频谱点频率偏移(从0开始,延伸到通道的带宽),第2列是第一个基带通道的幅度数据。 应该注意三个特征:

  • 带通朝着带宽的边缘下降,但不会降至零。 这意味着什么?
  • 数据中每1 MHz有尖峰。 这些是pcal信号。 随着频谱点数的增加,脉冲校准音的相对强度如何变化? 为什么?
  • RFI: 信号在带通顶部,与脉冲校准音无关,可能是RFI。 特别是S波段通道有RFI。

或许逐个查看16个基带通道会更有启发性,以了解真实数据的样子。

探测VDIF数据头

你可以使用名为printVDIFheader的程序解码并打印VDIF数据帧的帧头信息。
此程序接受VDIF数据文件(例如,来自ts036h观测)并打印每个帧的头的内容。
为了防止数据永远进行,您可能希望使用unix工具“ less”来控制输出。

1
2
cd /home/difx_data/ts036h
printVDIFheader NRAO+338_0039_TS036H_NL_No0001 | less

The header can be printed in binary with:
头可以以二进制形式打印:

1
printVDIFheader NRAO+338_0039_TS036H_NL_No0001 5032 | less

或者可以打印更多输出:

1
printVDIFheader NRAO+338_0039_TS036H_NL_No0001 5032 long | less

探测80 Hz噪声

一个名为m5fold的程序“折叠”给定频率处的功率,并产生该周期内相位的功率曲线。 在VLBA(以及其他几个VLBI天线)上,噪声功率是通过方波注入并由下游电子设备解码以确定系统温度的。

可以尝试查看n6043的数据:

1
2
cd /home/difx_data/n6043
m5fold NRAO+301_0540_N6043_PT_No0008 Mark5B-2048-16-2 50 20000 80 fold80

然后使用gnuplot绘制输出文件“ fold80”。 第一列是相位(以秒为单位表示自第一个相位bin以来的时间),第2列到第17列对应于16个基带通道。 上述命令中的数字“ 80”是要折叠的频率。

gnuplot

p 'fold80' u 1:5

试一试其他频率,例如与80 Hz不成比例的频率(例如79 Hz)。 如果您尝试40 Hz,您会期望什么?

你可以通过增加命令中的数字“ 50”来增加折叠轮廓中的样本数。 您可以通过更改“ 20000”来增加或减少要处理的数据量。 越大,输出信噪比越高。

研究其他人的相关处理

这个相关状态消息被发送到每个人都可以访问的本地网络上的多播地址。
因此,您可以在不干扰的情况下查看其他机器正在做什么。
要做到这一点,请启动2个以上的终端,并在每个终端中运行setup_difx setup脚本。 然后运行

1
2
$ errormon
$ statemon

这两个终端将显示与专用监视器终端相同的输出,但将包括房间中运行的所有作业的信息。

处无为之事,行不言之教;作而弗始,生而弗有,为而弗恃,功成不居!

欢迎关注我的其它发布渠道