How to plot ifft for a frequency signal?

7 次查看（过去 30 天）

Priyom Goswami 2021-8-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/893042-how-to-plot-ifft-for-a-frequency-signal

评论： Yazan 2021-8-6

FrequencyDATA.mat

I've a .dat file (attaching the .mat file of that .dat file after converting, since this platform doesn't allow to upload .dat file while asking question) which represent the frequency spectrum of a vibration signal from a bearing having a sampling frequency of (fs=32768). How to convert the frequency signal to time domain signal. Please help.

3 个评论
显示 1更早的评论隐藏 1更早的评论

Priyom Goswami 2021-8-6

that is what i'm unable to do so.

Scott MacKenzie 2021-8-6

编辑：Scott MacKenzie 2021-8-6

Have you tried? Please post the code you've attempted so far.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Yazan 2021-8-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/893042-how-to-plot-ifft-for-a-frequency-signal#answer_762302

编辑：Yazan 2021-8-6

在 MATLAB Online 中打开

clc, close all
fs = 32768;
fidx = z.Date23122013;
specOneSide = z.Time100351;
[fidx, idx] = unique(fidx);
specOneSide = specOneSide(idx);
fOneSide = (fidx./fidx(end))*fs/2;
subplot(1,2,1)
plot(fOneSide, specOneSide), xlabel('Frequency - Hz'), ylabel('Magnitude spectrum');
axis tight, title('One-sided spectrum');
specTwoSide = [flip(specOneSide(1:end)); specOneSide(2:end-1)];
N = length(specTwoSide);
f = -fs/2:fs/N:fs/2-fs/N;
subplot(1,2,2)
plot(f, specTwoSide), xlabel('Frequency - Hz'), ylabel('Magnitude spectrum');
axis tight, title('Two-sided spectrum');
sig = ifft(fftshift(specTwoSide));
figure,
t = (0:N-1)./fs;
plot(t, sig), axis tight, title('Signal');
xlabel('Time - sec'), ylabel('Amplitude');

5 个评论
显示 3更早的评论隐藏 3更早的评论

Priyom Goswami 2021-8-6

Thank You soo much. But i think there is something wrong with the data generated. Anyways, thank you once again.

Yazan 2021-8-6

Notice that the values of the Date23122013 (or your frequency axis points) are not unique. Also, for one value of Date23122013, you have sometimes multiple values of Time100351. My code uses the first appearing value of Time100351 for a given Date23122013 automatically, which is sometimes zero, but you might prefer otherwise. It is always a good idea to clean your data first.

请先登录，再进行评论。