我的 FFT 方法有问题吗?

信息处理 fft 频率 自由度 信号分析 功率谱密度
2022-02-04 08:22:37

我正在尝试通过傅立叶变换分析系统。我的系统有 3 种主要模式。(1,0.05) (0.51,0.01),(0.46,0.01) 这些对被定义为 (freq,damping_ratio)
我已经通过特征值分析得到了这些结果。系统包含 3 个信号。所以这些频率 [1,0.51, 0.46]也必须在这三个信号中找到。所以我运行了 FFT,但我没有在图中看到这些频率。你心中的问题是什么?为什么这些频率在 FFT 图中不存在?
在这里你可以得到这个 txt 文件中包含 3 个信号的信号。我把最后一个信号 FFT(PSD) 放在这里。Tnx 在此处输入图像描述 在此处输入图像描述
如您所见,第三信号中不存在这些频率。其他 2 个信号也一样。我也会把我的代码放在这里。

t=linspace(0,30,3012);
fS=ceil(inv(t(2)-t(1))); 
x=data(:,3);
 Nfft = 10 * 2^nextpow2(length(t));
    psd = 20.*log10(fftshift(abs(fft(x,Nfft))));
    freqs=[0:Nfft - 1].*(fS/Nfft);
    freqs(freqs >= fS/2) = freqs(freqs >= fS/2) - fS;
    freqs=fftshift(freqs);
    figure(1);
    plot(freqs, psd); 
    xlim([-1.5  1.5]);
    xlabel('Frequency / Hz');
    title (sprintf('PSD'));
0个回答
没有发现任何回复~