hc223.m

digital signal processing-MATLAB

%《数字信号处理教程——MATLAB释义与实现》第二章例2.2.3程序hc223
% 离散信号序列的移位、相加和相乘
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
x1=[0,1,2,3,4,3,2,1,0];ns1=-2;		% 给定x1及ns1
x2=[2,2,0,0,0,-2,-2]; ns2=2;		% 给定x2及ns2
nf1=ns1+length(x1)-1; nf2=ns2+length(x2)-1;
ny= min(ns1,ns2):max(nf1,nf2); 		% y(n)的位置向量
y1 = zeros(1,length(ny)); y2 = y1;       % 延拓序列初始化
y1(find((ny>=ns1)&(ny<=nf1)==1))=x1; % 给y1赋值x1 
y2(find((ny>=ns2)&(ny<=nf2)==1))=x2; %  给y2 赋值x2
ya = y1 + y2;                             	% 序列相加
yp = y1.* y2;                             	% 序列相乘
subplot(4,2,1), stem(ns1:nf1,x1,'.')		% 绘图
xlabel('nx1'),ylabel('x1'),axis([-5,10,0,4])
subplot(4,2,3), stem(ns2:nf2,x2,'.'),axis([-5,10,-2,2])
xlabel('nx2'),ylabel('x2')
subplot(4,2,2), stem(ny,y1,'.')		% 绘图
xlabel('ny'),ylabel('y1')
subplot(4,2,4), stem(ny,y2,'.')
xlabel('ny'),ylabel('y2')
line([ny(1),ny(end)],[0,0])		% 画x轴
subplot(4,2,6), stem(ny,ya,'.')
xlabel('ny'),ylabel('ya')
line([ny(1),ny(end)],[0,0])		% 画x轴
subplot(4,2,8), stem(ny,yp,'.')
xlabel('ny'),ylabel('yp')
line([ny(1),ny(end)],[0,0])		% 画x轴
set(gcf,'color','w')            % 置图形背景色为白