hc325.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第三章例3.2.5程序hc325
% 序列信号DTFT的对称性和时移性
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
x=[1,3,5,3,1];nx=[-1:3];
w=linspace(-8,8,1000);		% 设定频率向量
X=x*exp(-j*nx'*w);			% 用（3.2.1a）式计算DTFT
subplot(5,3,1),stem(nx,x,'.'),axis([-2,6,-1,6])		% 画序列图
title('原始序列')
ylabel('x(n)')
subplot(5,3,4),plot(w,abs(X)),	% 画幅频相频曲线
ylabel('幅度')
subplot(5,3,7),plot(w,angle(X))
ylabel('相角')
subplot(5,3,10),plot(w,real(X)), % 画实频虚频曲线
ylabel('实部')
subplot(5,3,13),plot(w,imag(X))
ylabel('虚部')
nx1=nx+2; X1=x*exp(-j*nx1'*w);			% 使x右移两位后，计算其DTFT
subplot(5,3,2),stem(nx1,x,'.'),axis([-2,6,-1,6])
title('右移两位')
subplot(5,3,5),plot(w,abs(X1)), subplot(5,3,8),plot(w,angle(X1))
subplot(5,3,11),plot(w,real(X1)), subplot(5,3,14),plot(w,imag(X1))
nx2=nx-1; X2=x*exp(-j*nx2'*w);		    % 使x左移一位后，计算其DTFT
subplot(5,3,3),stem(nx2,x,'.'),axis([-2,6,-1,6])
title('左移一位')
subplot(5,3,6),plot(w,abs(X2)), subplot(5,3,9),plot(w,angle(X2))
subplot(5,3,12),plot(w,real(X2)), subplot(5,3,15),plot(w,imag(X2))
set(gcf,'color','w')                        % 置图形背景色为白