hc362.m

这是数字信号处理——matlab释义与实现一书的源代码

%《数字信号处理教程——MATLAB释义与实现》第三章例3.6.2程序hc362
% 序列的循环折叠
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
clear
x = [3,4,5,6,8,9,10]; N = 10;nx=0:N-1;                              % 给出原始数据,假如nx从零开始
x1=[x,zeros(1,N-length(x))];                                        % 将x的长度通过补零扩展到N
y1 = x1(mod(-nx,N)+1);                                              % 把x循环折叠，求得y
subplot(2,2,1); stem([0:N-1],x1); title('初始序列')                 % 以下绘图
ylabel('x1(n)'); axis([-0.5,N,-1,10.5]);text(10.8,-1,'n')
subplot(2,2,3); stem([0:N-1],y1); title('循环折叠序列')
ylabel('x(mod(-n,N)+1)'); axis([-0.5,N,-1,10.5]);text(10.8,-1,'n')    
N1=7;x2(1:N1)=0;nx2=0:N1-1;
for i=0:N1-1 x2(mod(i,N1)+1)=x2(mod(i,N1)+1)+x(i+1);end
y2 = x2(mod(-nx2,N1)+1);                                            % 把x循环折叠，求得y
subplot(2,2,2); stem([0:N1-1],x2); title('初始序列')                % 以下绘图
ylabel('x2(n)');axis([-0.5,N1,-1,10.5]); text(10.8,-1,'n')
subplot(2,2,4); stem([0:N1-1],y2); title('循环折叠序列')
ylabel('x(mod(-n,N1)+1)');axis([-0.5,N1,-1,10.5]);text(10.8,-1,'n')     
set(gcf,'color','w')                                                % 置图形背景色为白