ft_spect(版本2.0)计算的输入信号的振幅和相位谱所需的频率分辨率和也为抑制滤波器的相位谱浮动四舍五入误差。

注意1:ft_spect不能消除频谱泄漏。

注意# 2:离散傅里叶变换(DFT)的输入信号的一个周期周期信号和周期信号的频谱可根据输入信号的长度。信号的采样频率Fs,在T = NΔt,垃圾箱(a.k.频率。频率分辨率的意思区分频率f1和f2)间隔Δf = 1 / T = Fs / N;因此,DFT的频率分辨率只取决于输入信号的长度(T),但补零不会增加频率分辨率和光谱没有透露更多信息;它只篡改箱之间的振幅。提高光谱分辨率,长时间的测量是必要的,因为DFT的输入信号为周期信号的一个周期;因此,重复输入信号是可以接受的,不产生任何工件。但是,在本例中,输入信号的长度增加,因此,光谱分辨率也会增加。

注意# 3:相位谱因为浮动四舍五入误差很吵。小四舍五入的“反正切”计算误差产生显著噪音相位谱的结果。为了抑制这种噪声,ft_spect使用一个阈值过滤。这意味着如果特定频率的振幅小于预定义的阈值,它把零代替它。

如果你在研究中使用这个程序,请引用以下文章:

Afshin Aghayan, Priyank贾斯瓦尔,哈米德Reza Siahkoohi (2016)。
“地震去噪使用冗余提升计划”。GEOPHYSICS, 81(3), V249-V260.
https://doi.org/10.1190/geo2015 - 0601.1

* *请改善ft_spect分享您的建议和想法
通过aghayan@okstate.edu或afshin.aghayan@gmail.com

1.0版(2017年春季)ft_spect v1.0是MATLAB R2013a编写和测试。

2.0版(2020年春季)应用以下更改:

1)你可以定义你想要的频率分辨率(Δf)
2)v1.0要快得多
3)添加一个例子来比较通常的FFT和这个程序的输出(只类型ft_spect演示;看看ft_spect_demo函数结束时程序的更多细节)

Afshin Aghayan
afshin.aghayan@gmail.com
405-334-7184