hc431.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第四章例4.3.1程序hc431
% 用FFT计算离散时间傅立叶变换DTFT
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
% 程序hc431
x=[1,3,5,3,1];nx=0:4;T=0.5;      % 给定原始数据
N1=length(x); D1=2*pi/(N1*T);     % 求出序列长度及频率分辨率
k1=floor((-(N1-1)/2):((N1-1)/2));  % 求对称于零频率的FFT位置向量
X1=fftshift(fft(x,N1));           % 求对称于零频率的FFT序列值
N2=16; D2=2*pi/(N2*T);     % 求出序列长度及频率分辨率
k2=floor((-(N2-1)/2):((N2-1)/2));  % 求对称于零频率的FFT位置向量
X2=fftshift(fft(x,N2));           % 求对称于零频率的FFT序列值
w=-2*pi:0.01:2*pi;		% 生成连续频谱的细分频率 
X=x*exp(-j*0.5*nx'*w); %   计算连续频谱
subplot(1,2,1),plot(w,abs(X),k1*D1,abs(X1),'o:',k2*D2,abs(X2),'x-.')    % 画幅频特性图
xlabel('频率'),ylabel('幅度')
subplot(1,2,2),plot(w,angle(X),k1*D1,angle(X1),'o:',k2*D2,angle(X2),'x-.')   % 画相频特性图
xlabel('频率'),ylabel('相位')
set(gcf,'color','w')              % 置图形背景色为白