大毕设-MATLAB-FFT实现

引用来自:http://blog.csdn.net/sinwel/article/details/8115673

%仿真参数中的含义

% Ts 表示间隔Ts时间采样，这个越小越接近连续信号，而实际上不可能真正的仿真连续信号。

% df 表示变化后频率分量的最小间隔，即频谱中相邻的点间隔频率是多大，你也可以叫它 频率分辨力

% 为什么df = fs/N-1呢？
% 因为采样频率 fs=1/Ts，表示的意义是连续信号用fs去采样，得到了
% 离散的信号，因此实际上的连续信号在这里并没有给出，因为MATLAB不可能做的真正的
% 连续信号，而且既然采样间隔是Ts，所以时域信号也就只给出间隔Ts的信号，这样可以
% 刚好满足仿真需求，当然你也可以把时域的信号变为更细，把时间间隔变得更小，

% 也就是说信号的个数是 采样频率的整数倍，确保能够正常才够采样。

4.fs的选取，根据采样定义，fs>2fmax

如果信号最高频率分量为10HZ,则采样频率一定要大于或等于20HZ，最好选大点。这样仿真效果明显

clear all;

Ts=0.01;
    t=-1:Ts:1;
    %采样信号频率是fs=100；
    fs=1/Ts;
    %原信号频率是10；
    x=5*sin(2*pi*10*t);
    %N是总共的采样点数，N=101；
    N=length(x);
    fx=fft(x);
    %N个点有N-1个间隔；
    df=fs/(N-1);
    %让负频率也显示出来；
    n=-(N-1)/2:(N-1)/2;
    f=n*df;
    subplot(211);
    plot(t,x);
    grid;
    title('原始信号');
    xlabel('时间(t)');
    ylabel('幅度');
    subplot(212);
    plot(f,abs(fftshift(fx))*2/N); 
    grid;
    title('频谱');
    xlabel('频率(hz)');
    ylabel('频率分量');

% 此外MATLAB中FFT计算结果是一个N点的复数，且第n点表示的含义是
% 频率为 Fn = (n-1)*Fs/N 的复数值（cos 和 sin的分量）
% 它的模大小abs就是傅里叶的频谱幅度大小。

% 关于fftshift这个函数，它的作用是：
% shift zero-frequency component to center of spectrum.
% 因为FFT出来的结果是频率0~Fs分布的，如果想看双边频率谱，就需要将频率变为
% -Fs~0~FS。