%《數字信號處理教程——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')                                                % 置圖形背景色為白