fft ifft的Matlab实现
clc;
clear;
t=0:0.1:10;
x1=sin(t);
x2=sin(10*t);
x=x1+x2;
X_f=fft(x);
figure (1)
subplot(2,1,1);
plot(x);
xlabel(‘time’)
title(‘x1+x2 (time domain)’)
subplot(2,1,2);
plot(abs(X_f));
xlabel(‘frequency’)
title(‘x1+x2(frequency domian)’)
X1_f=X_f;
X1_f(10:95) = 0;
X2_f=X_f;
X2_f(1:9)=0;
X2_f(96:end)=0;
x1_reconstruc = ifft(X1_f);
x2_reconstruc = ifft(X2_f);
figure (2)
subplot(2,1,1);
plot(x1);
title(‘Original x1’)
subplot(2,1,2);
plot(real(x1_reconstruc));
title(‘Reconsturctured x1’)
figure (3)
subplot(2,1,1);
plot(x2);
title(‘Original x2’)
subplot(2,1,2);
plot(real(x2_reconstruc));
title(‘Reconsturctured x2’)
声明:本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人,不代表电子发烧友网立场。文章及其配图仅供工程师学习之用,如有内容侵权或者其他违规问题,请联系本站处理。 举报投诉
全部0条评论
快来发表一下你的评论吧 !