healpy函数合成中lmax增加的偏谱

Question

我最近使用python使用healpy的函数synfast。更具体地说，我使用了它的参数lmax，它设置了要计算的最大模式（默认值：lmax=3*nside-1）。正如您将看到的，我还使用了ud_grade函数。你知道吗

我研究了几个例子，在这些例子中，我用synfast从等于1的平坦光谱（温度TT和极化EE和BB）生成了非常简单的映射。然后使用anafast从这些映射中计算出光谱。我发现并不是所有情况下都能恢复初始输入频谱，如下图所示。你知道吗

在所有这些情况下，我将synfast的pixwin参数设置为True，然后使用函数pixwin(nside=128)校正anafast的输出信号。你知道吗

• 灰色虚线：输入光谱模型
• Blue：（正常默认情况下）synfast(nside=128, lmax=3*128-1)=>；anafast(lmax=3*128-1)
• 绿色：synfast(nside=256, lmax=3*256-1)=>；ud_grade(256->128)=>；anafast(lmax=3*128-1)
• 洋红色：synfast(nside=256, lmax=3*128-1 )=>；ud_grade(256->128)=>；anafast(lmax=3*128-1)
• 红色：synfast(nside=128, lmax=4*128-1 )=>；anafast(lmax=3*128-1)
• 青色：synfast(nside=128, lmax=10*128-1 )=>；anafast(lmax=3*128-1)

Plot : spectra estimations from different kind of map using syfast, ud_grade and anafast

• 正如你所看到的，在大多数情况下，我得到了有偏的光谱。你知道吗
• 即使是正常情况（蓝色）是稍微向下偏移高模式。这显然是正常的，因为healpy不能很好地处理模式>；2*nside-1的频谱估计。你知道吗
• 温度和偏振光谱（TT，EE和BB）的结果是相同的。你知道吗
• 如您所见，情况lmax=4*128-1和lmax=10*128-1给出了相同的结果。这来自synfast，它不允许lmax高于4*128-1，并且设置lmax=4*128-1，即使我们要求更高的lmax。你知道吗

有人对这种行为有什么解释吗？你知道吗

如何恢复初始输入频谱？你知道吗

先谢谢你

Answer 1

这是预期的影响，这是一个破坏性的行动，你消除了一些规模的信号，事实上，你恢复了偏见，这是一致的白噪声。你知道吗