%《MATLAB在電子信息課程中的應用》第七章例7.1程序q701
% 離散信號序列的移位、相加和相乘
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

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)的時間變量
xa1 = zeros(1,length(ny)); xa2 = xa1;       % 延拓序列初始化
xa1(find((ny>=ns1)&(ny<=nf1)==1))=x1; % 給xa1賦值x1 
xa2(find((ny>=ns2)&(ny<=nf2)==1))=x2; %  給xa2 賦值x2
ya = xa1 + xa2                              	% 序列相加
yp = xa1.* xa2                             	% 序列相乘
subplot(4,1,1), stem(ny,xa1,'.')		% 繪圖
subplot(4,1,2), stem(ny,xa2,'.')
line([ny(1),ny(end)],[0,0])		% 畫x軸
subplot(4,1,3), stem(ny,ya,'.')
line([ny(1),ny(end)],[0,0])		% 畫x軸
subplot(4,1,4), stem(ny,yp,'.')
line([ny(1),ny(end)],[0,0])		% 畫x軸