ddfftt1.m

来自「这是一个经典的把matlab应用于DFT的例子,适合于初学者」· M 代码 · 共 30 行

%{ 
用FFT对序进行谱分析
x1(n)=R4(n),x2(n)=cos(n*pi/4)+cos(n*pi/8)
X1:存放x1(n)的向量,Y1:存放X1(k)
X2:存放x2(n)的向量,Y2:存放X2(k)
%}
echo on     %echo部分命令回显至command windows
x1=[1,1,1,1,0,0,0,0];
n=0:15;
x2=cos(pi*n/4)+cos(pi*n/8);
i=0:7;
subplot(3,2,1);stem(i,x1,'.');
axis([0 7 0 1]);       %规定x轴和y轴标值范围
xlabel('n');ylabel('x1(n)');
y1=fft(x1,8);          %x1(n)的8点FFT
subplot(3,2,3);stem(i,abs(y1),'.');
xlabel('(N=8 wk=2*pi*k/N)k');ylabel('|X1(k)|');
y1=fft(x1,64);          %x1(n)的64点fft
i=0:63;
subplot(3,2,5);stem(i,abs(y1),'.');
axis([0 63 0 4]);
xlabel('(N=64 wk=2*pi*k/N)k');ylabel('|X1(k)|');
y2=fft(x2);             %x2(n)的16点fft
                 %另绘一幅图
subplot(2,2,2);stem(n,x2);
title('x2(n)时域序列');xlabel('n');ylabel('x2(n)');
subplot(2,2,4);stem(n,abs(y2));
title('x2(n)的幅频特性');xlabel('(N=16 wk=2*pi*k/N)k');ylabel('|X2(k)|');
echo off