hc323.asv

某本matlab书籍中的代码

%《数字信号处理教程——MATLAB释义与实现》第二章例3.2.3程序hc323
% 信号序列的DTFT
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
M=4;N=2*M+1;T=0.5;n=-16:16; 
x=[zeros(1,12),ones(1,N),zeros(1,12)];  % 给出输入序列
w=[-15:0.1:15]+1e-10;
X=sin(0.5*N*w*T)./sin(0.5*w*T);         % 给出频谱序列
subplot(1,3,1),stem(n,x,'.')            % 画出输入序列
axis([-20,20,-0.1,1.1]),grid on
xlabel('n'),title('(a) 序列幅度')
subplot(1,3,2),plot(w,X),grid on        % 画出频谱序列
xlabel('\Omega'),title('(b) 幅频特性')
subplot(1,3,3),plot(w,X),grid on        % 改变横轴比例，画出频谱序列
v=axis;axis([-pi/T,pi/T,v(3),v(4)])
xlabel('\Omega'),title('(c) 横轴放大后幅特性')
set(gcf,'color','w')                    % 置图形背景色为白