fftanalysis.m

来自「数字信号处理程序——用FFT实现信号频谱分析的matlab程序」· M 代码 · 共 21 行

%已知序列 x(n) = 2*sin(0.48*pi*n) + cos(0.52*pi*n) ,0 =< n <=99，试绘制及它
%的离散傅里叶变换图。
%MATLAB实现程序：
clear all
N =100;
n =0:N -1;
xn = 2*sin(0.48*pi*n) + cos(0.52*pi*n);
XK = fft(xn,n);
magXK = abs(XK);
phaXK = angle(XK);
subplot(1,2,1)
plot(n,xn)
xlabel('n');
ylabel('x(n)');
title('x(n) N =100')
subplot(1,2,2)
k = 0:length(magXK) - 1;
stem(k,magXK,'.');
xlabel('k');
ylabel('X(K)');
title('X(k) N =100');