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)移位后的周期延拓
figure(2)
subplot(3,2,1); stem(n,xn,'.');         % 画出x(n)序列
set(gcf,'color','w')                    % 设置图形背景色为白色
title('x(n)=0.8^{n}R_{8}(n)');axis([0,length(n),0,1]);
ylabel('x(n)');xlabel('n');
subplot(3,2,2);stem(n,xm,'.');          % 画出xm(n)序列
title('x(n-m) 移位');axis([0,length(n),0,1]);
ylabel('x(n-m)');xlabel('n');
subplot(3,2,5); stem(n,xc,'.');         % 画出xc(n)序列
title('x(n)的周期延拓');axis([0,length(n),0,1]);
ylabel('x(mod(n,8))');xlabel('n');
subplot(3,2,6);stem(n,xcm,'.');         % 画出xcm(n)序列
title('x(n)的循环移位');axis([0,length(n),0,1]);
ylabel('x(mod(n-m,8))');xlabel('n');