hc361.m

某本matlab书籍中的代码

%《数字信号处理教程——MATLAB释义与实现》第三章例3.6.1程序hc361
% 有限长正弦序列的频谱
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
x=[7,6,5,4,3,2];Nx=length(x);nx=0:Nx-1;		% 给出x序列		
nx1=-Nx:2*Nx-1;x1=x(mod(nx1,Nx)+1);	% 延拓为周期向量x1,注意mod用法
[y1,ny1]=sigshift(x1,nx1,2);			% 将x1右移两位，得到y1
RN=(nx1>=0)&(nx1<Nx);				% 在x1的位置向量nx1上设置主值窗
RN1=(ny1>=0)&(ny1<Nx); 				% 在y1的位置向量ny1上设置主值窗
subplot(4,1,1),stem(nx1,RN.*x1)		% 在子图上画出x1的主值
title('主值序列x')
axis([-6,15,0,10])
subplot(4,1,2),stem(nx1,x1) 		% 画出x1
title('周期序列x1')
axis([-6,15,0,10])
subplot(4,1,3),stem(ny1,y1)			% 画出y1
title('移位周期序列y1')
axis([-6,15,0,10])
subplot(4,1,4),stem(ny1,RN1.*y1) 	% 画出y1的主值
title('y1的主值序列y')
set(gcf,'color','w')              % 置图形背景色为白
axis([-6,15,0,10])