q703.m

MATLAB及在电子信息课程中的应用例题索引与程序

%《MATLAB在电子信息课程中的应用》第七章例7.3程序q703
% 序列x(n)的移位与周期延拓
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

%N: 观察窗口长度
%M: 序列x(n)长度
%m: 移位样点数
clear;close all
N=24;M=8;
m=input('输入移位值：m=');
if   (m<1|m>=N-M+1)			% 检验输入参数m是否合理
     fprintf('输入数据不在规定范围内！');
     break
  end
      n=0:N-1;
      x1=(0.8).^n; x2=[(n>=0)&(n<M)]; %产生x(n)
      xn=x1.*x2;
      xm=zeros(1,N);	% 设定xm的长度
      for k=m+1:m+M
         xm(k)=xn(k-m);
      end
      xc=xn(mod(n,8)+1); 		% 产生 x(n)的周期延拓
      xcm=xn(mod(n-m,8)+1); 	% 产生 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)的循环移位');axis([0,length(n),0,1]);
      ylabel('x(mod(n-m,8))');xlabel('n');