离散傅里叶变化-matlab的fft
下面是使用matlab的自带的例子进行离散傅里叶变换的程序，非常的简单。
写下来主要是给自己看的。
%% matlab离散傅里叶变换 %% Vsesion:1.0 Written by zhenhuaLiu@ 2021.11.15 HIT ATCI clear clc close all fs = 1000; t = 0: (1/fs): (1-1/fs);%时间向量 x = sin(2*pi*100*t)+sin(2*pi*200*t);%多点频信号,分别是100Hz和200Hz N = length(x);%向量的长度 figure(1) plot(t, x); title(原始信号); X = fft(x);%DFT将输入信号当作是周期信号，所以需要除以N，又因为单边谱的幅度是双边谱的两倍，因此需要乘以2 X_shift = fftshift(X);%可以将将零频率的分量移到频谱的中心 index_X = (-N/2 : N/2-1)*fs;%生成以Hz为单位的频率下标 figure plot(index_X, abs(X_shift)*2/N);%对傅里叶变换幅度谱进行绘制 xlabel(HZ) ylabel(幅度) title(离散傅里叶变换的幅度谱); x_idft = ifft(X);%逆离散傅里叶变换， figure plot(t, x_idft);%对离散傅里叶变换的结果进行显示 title(逆变换信号);
下面是显示结果：
