dft.m

离散傅里叶变换的程序

function [c,b]=DFT(x,N,str)
n=0:N-1;
k=-400:400;
w=(pi/100)*k;
b=w;
c=x*(exp(-j*pi/100)).^(n'*k);
magX=abs(c);
angX=angle(c);
subplot(2,1,1);
t=max(x);
n=0:N-1;
stem(n,x,'.');
xlabel('n');
ylabel('X(n)');
title('信号的原形');
text((0.3*N),(0.8*t),str);
hold on
n=0:N-1;
m=zeros(N);
plot(n,m);
subplot(2,1,2);
plot(w/pi,magX);
axis([-4 4 0 max(magX)*5/4]);
xlabel('w/pi');
ylabel('|Y(jw)|');
title('上图信号的傅氏变换|jw|');