快速傅里叶变换实验教案分析.docxVIP

实验二：快速傅里叶变换 一．基本信号（函数）的FFT变换 1. 所利用的matlab代码基本如下： clc;clear; fs=input(选择采样频率：); N=input(选择采样点数：); n=0:N-1;t=n/fs; x=sin(2*pi*t+pi/6)+sin(4*pi*t)+cos(6*pi*t); subplot(3,1,1) stem(t,x) title([采样频率为 num2str(fs) 截断长度为 num2str(N) 的采样信号]); xlabel(t/T); ylabel(x); grid on; y=fft(x,N); y=fftshift(y)/N; mag=abs(y); phase=angle(y)*180/pi; f=n*fs/N-fs/2; subplot(3,1,2) stem(f,mag) title(快速傅里叶变换后的幅频图) xlabel(f/f0); ylabel(S); grid on; subplot(3,1,3) stem(f,phase) title(快速傅里叶变换后的相频图) xlabel(f/f0); ylabel(Φ); grid on; figure(1); （1） 采样频率fs=8f0，截断长度N=16 分析：信号中最高频率fc=3f0，采样频率为fs=8f0?，满足采样定理。采样点数N=16，根据?f=fsN=