如何将频域转换为时域

时间:2012-12-26 03:52:13

标签: matlab time fft ifft time-frequency

我创造了一个3个不同频率的信号并滤除了我不想要的信号。但是当我在matlab中使用ifft时,它会显示错误的图形。如何将我的频域频谱转换回我的3秒时域图?我的代码如下:

clc
clear all
Fs = 8192;
T = 1/Fs;                                             
%create tones with different frequency
t=0:T:1;
t2=1:T:2;
t3=2:T:3;

y1 = sin(2*pi*220*t); 
y2 = sin(2*pi*300*t2);
y3 = sin(2*pi*440*t3);                          

at=y1+y2+y3;


figure;
plot(t,y1,t2,y2,t3,y3),title('Tones with noise');

[b,a]=butter(2,[2*290/Fs,2*350/Fs],'stop');
e=filter(b,a,at);
et=(ifft(abs(e)));
figure,
plot(et)

1 个答案:

答案 0 :(得分:1)

就像现在一样,由于fft,et在频域中。你不需要fft。只需plot(e),您就可以获得时域滤波波形。 Yo可以检查频率中的过滤器性能。域名依据fft,只是

plot(abs(fftshift(fft(fftshift(e)))));
xlim([4000 5000])

修改

您在问题上写的代码有以下错误:at正好有1秒的信息(或8192个元素)。如果您plot(at),您会看到频率总和,但它们在同一时间happen。这是如何解决它:

clear all
Fs = 8192; % or multiply by 3 if needed
T = 1/Fs;                                             
%create tones with different frequency
t=0:T:3;

y1 = sin(2*pi*220*t).*(t<1); 
y2 = sin(2*pi*300*t).*(t<2 & t>=1);
y3 = sin(2*pi*440*t).*(t>=2);                          

at=y1+y2+y3;

[b,a]=butter(2,[2*290/Fs,2*350/Fs],'stop');
e=filter(b,a,at);
figure,
plot(t,e)

enter image description here

dt=t(2)-t(1);
N=length(at);
df=1/(N*dt);        % the frequency resolution (df=1/max_T)
if mod(N,2)==0
    f_vector= df*((1:N)-1-N/2);     % frequency vector for EVEN length vectors: f =[-f_max,-f_max+df,...,0,...,f_max-df]
else
    f_vector= df*((1:N)-0.5-N/2);   % frequency vector for ODD length vectors f =[-f_max,-f_max+fw,...,0,...,f_max]
end
freq_vec=f_vector;
fft_vec=fftshift(fft(e));

plot(freq_vec,abs(fft_vec))
xlim([0 1000])

enter image description here