离散傅里叶变换和逆变换-matlab编程实现

最近写了一个压缩感知的程序，稀疏信号是一个频域稀疏信号，因此需要进行离散傅里叶变换。下面编程实现离散傅里叶变化:

一个周期序列的离散傅里叶变换如下：

其中：

一个N点的逆离散傅里叶变换给出为：

用matlab编程，我用矩阵乘法实现的

%% matlab离散傅里叶变换，不使用matlab函数，通过矩阵相乘的方法自己实现
%% Vsesion:1.0 Written by zhenhuaLiu@ 2021.11.15 HIT ATCI
function [Xk] = dft(xn,N)
% Computes Discrete Fourier Transform
%______________________________________________
% [Xk] = dft(xn,N)
% Xk = DFT coefficients array over 0 <= k <= N - 1
% xn = N-point finite - duration sequence 0 <= n <= N - 1
% N = Length of DFT
 
n = [0:1:N-1]; % row vector for n
k = [0:1:N-1]; % row vector for k
WN = exp(-1i*2*pi/N);%matlab提示要增强稳定性，使用li代替i
nk = n'*k;
WNnk = WN .^ nk;   %DFT matrix
Xk = xn * WNnk;

逆变换如下：

%% matlab离散傅里叶变换，不使用matlab函数，通过矩阵相乘的方法自己实现
%% Vsesion:1.0 Written by zhenhuaLiu@ 2021.11.15 HIT ATCI
function [xn] = idft(XK,N)
% Computes Discrete Fourier Transform
%______________________________________________
% [Xk] = dft(xn,N)
% Xk = DFT coefficients array over 0 <= k <= N - 1
% xn = N-point finite - duration sequence 0 <= n <= N - 1
% N = Length of DFT
 
n = [0:1:N-1]; % row vector for n
k = [0:1:N-1]; % row vector for k
WN = exp(1i*2*pi/N);
nk = k'*n;
WNnk = WN .^ nk;   %DFT matrix
xn = XK * WNnk;