我创建了一个简单的正弦钟变迹函数,在执行 2D FFT 并制作结果的等高线图之前将其应用于我的数据,如下所示。
P = 90
apd_fn = np.sin(((np.pi - (np.pi/180.0) * P) / 1024) * np.arange(0, 1024, 1) + (np.pi/180.0) * P)**2
fid_mat = ws_list * apd_fn
jres_spec = abs(fftshift(fft2(fid_mat)))
max_val = (jres_spec[:,400:]).max()
figure()
CS=contour(jres_spec,arange(0,max_val,max_val/30))
xlabel('frequency (ppm)')
show()
在 x 维度上开窗后,这会生成一个 2D 光谱图像,如下所示:
我现在想在 y 维度上应用相同的过滤器,但我有一点心理障碍,很难弄清楚如何实现它。您能提供的任何帮助将不胜感激。
谢谢。