exa2_8.m

新版Matlab 7辅助信号处理技术与应用

% exa2-8_fft_ifft.m for example 2-8
% to test fft.m and ifft.m

clear all;
 % 产生两个正弦加白噪声；
 N=256;
 f1=.1;
 f2=.2;
 fs=1;
 a1=5;a2=3;
 w=2*pi/fs;
 x=a1*sin(w*f1*(0:N-1))+a2*sin(w*f2*(0:N-1))+randn(1,N);
 
 % 应用FFT 求频谱； 
 subplot(3,1,1);
 plot(x(1:N/4));
 ylabel('原始信号')
 f=-0.5:1/N:0.5-1/N;
 X=fft(x);
 subplot(3,1,2);
 plot(f,fftshift(abs(X)));
 ylabel('FFT频谱')
 
 % 应用IFFT 求原信号； 
 y=ifft(X);
 subplot(3,1,3);
 plot(real(x(1:N/4))); 
 ylabel('IFFT信号')