下面是使用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('逆变换信号');
下面是显示结果:
下一篇我自己编程实现离散傅里叶变换和逆变换。