x=rand([2000,1]); % 离散数字信号,当然你也带入你的数字信号
fs=1024; % 对应于离散数字信号的采样速率
nfft=1024; % 做fft的点数
X=fftshift(fft(x,nfft));% 做fft变换
f=fs*(0:nfft-1)/nfft-fs/2;
figure;
subplot(211);
title('幅频特性');
plot(f,abs(X));
subplot(212);
plot(f,angle(X));
title('相频特性');