代码：

fs=100;N=128;
n=0:N-1;t=n/fs;
x=sin(2*pi*40*t)+sin(2*pi*15*t);
subplot(221);plot(n,x,'b');
xlabel('时间/s');ylabel('x');title('原始信号');
grid on;y=fft(x,N);
mag=abs(y);
f=n*fs/N;
subplot(222);plot(f(1:N/2),mag(1:N/2)*2/N,'b');
xlabel('频率/Hz');ylabel('振幅');title('原始信号的FFT变换');grid on;
xifft=ifft(y);
realx=real(xifft); %  快速傅里叶变换的实部
ti=[0:length(xifft)-1]/fs;
subplot(223),plot(ti,realx,'b');
xlabel('时间/s');ylabel('x');  title('利用傅里叶逆变换得到的信号') yif=fft(xifft,N);
mag=abs(yif);
f=[0:length(y)-1]'*fs/length(y);  %   ’  好像也可以不要
subplot(224);  plot(f(1:N/2),mag(1:N/2)*2/N);
xlabel('频率/Hz');ylabel('振幅');title('运用IFFT得到的信号的快速FFT变换');
grid on;