dftmatlab.txt

来自「DFT：离散傅立叶变换 符运行后得频谱图」· 文本 代码 · 共 26 行

DFT：离散傅立叶变换

MATLAB代码：

clear; 
clc; 

N=12; % 数据长度
n=0:N-1; 
xn=sin(pi*n/6); % 以此序列为例，1/6是信号频率对1/2采样频率归一化的值
k=0:N-1; 
WN=exp(-j*2*pi/N); 
nk=n'*k; 
WNnk=WN.^nk; 
Xk=xn*WNnk; % DFT结果 

% 屏幕打印 
xn 
lXkl=abs(Xk) 


% 画图 
figure(1),plot(n,xn,'-O'); % 时域
figure(2),bar(k,abs(Xk)); % 频域