%《MATLAB在電子信息課程中的應(yīng)用》第七章例7.3程序q703
% 序列x(n)的移位與周期延拓
% 電子工業(yè)出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

%N: 觀察窗口長度
%M: 序列x(n)長度
%m: 移位樣點(diǎn)數(shù)
clear;close all
N=24;M=8;
m=input('輸入移位值：m=');
if   (m<1|m>=N-M+1)			% 檢驗(yàn)輸入?yún)?shù)m是否合理
     fprintf('輸入數(shù)據(jù)不在規(guī)定范圍內(nèi)！');
     break
  end
      n=0:N-1;
      x1=(0.8).^n; x2=[(n>=0)&(n<M)]; %產(chǎn)生x(n)
      xn=x1.*x2;
      xm=zeros(1,N);	% 設(shè)定xm的長度
      for k=m+1:m+M
         xm(k)=xn(k-m);
      end
      xc=xn(mod(n,8)+1); 		% 產(chǎn)生 x(n)的周期延拓
      xcm=xn(mod(n-m,8)+1); 	% 產(chǎn)生 x(n)移位后的周期延拓
      subplot(3,2,1); stem(n,xn,'.');
      title('x(n)');axis([0,length(n),0,1]);
      ylabel('x(n)');xlabel('n');
      subplot(3,2,2);stem(n,xm,'.');
      title('x(n-m)');axis([0,length(n),0,1]);
      ylabel('x(n-m)');xlabel('n');
      subplot(3,2,5); stem(n,xc,'.');
      title('x(n)的周期延拓');axis([0,length(n),0,1]);
      ylabel('x(mod(n,8))');xlabel('n');
      subplot(3,2,6);stem(n,xcm,'.');
      title('x(n)的循環(huán)移位');axis([0,length(n),0,1]);
      ylabel('x(mod(n-m,8))');xlabel('n');