我创建了一个简单的正弦振铃变迹函数,我在执行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维度中应用相同的过滤器,但我有一点心理障碍,并努力解决如何实现这一点。非常感谢您提供的任何帮助。
感谢。