fft ifft的Matlab實現(xiàn)

　　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’）