matlab对序列实现FFT与IFFT

第一部分

代码

%% 完成对给出实数序列的FFT与IFFT
N = 64;
n = 0:2*N-1;

% 定义一个函数句柄，表示自变量为 x，因变量为 y 的函数
f = @(n) cos(2*pi/N*7*n) + 1/2*cos(2*pi/N*19*n);

%计算FFT
y = fft(f(n),2*N);

%计算复数模
Y = abs(y);

%绘制离散序列图
stem(n, Y);
xlabel('k');
ylabel('|X(k)|');
title('离散序列的频谱图');

运行结果

第二部分

代码

%% 实现对离散频域序列的IFFT

N = 64;
k = 0 : N-1;

% 定义一个函数句柄，表示自变量为 n，因变量为 y 的函数 
f = @(k) 1/(1-0.6*exp(-1j*2*pi/N));

%进行IFFT
x = ifft(f(k), N);

%画出离散序列图
figure(1);
stem(k, real(x));
xlabel('n');
ylabel('x(n)');
title('离散频域序列反变换结果实部');

figure(2);
stem(k, imag(x));
xlabel('n');
ylabel('x(n)');
title('离散频域序列反变换结果虚部');

运行结果