%《MATLAB在電子信息課程中的應用》第七章例7.6程序q706
% 序列卷積運算
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

% Nx：輸入序列x(n)長度
% Nh：h(n)長度
% m：x(n)的移位樣點數
clear;close all
Nx=20;Nh=10;m=5;
n=0:Nx-1;
x1=(0.9).^n;   %x1(n)
%產生x2(n)=x1(n-m)
x2=zeros(1,Nx+m);
for k=m+1:m+Nx
    x2(k)=x1(k-m);
end
nh=0:Nh-1;
h1=ones(1,Nh);   %h1(n)
h2=h1;             %h2(n)
y1=conv(x1,h1);  %y1(n)
y2=conv(x2,h2);  %y2(n)
Ny1=length(y1);  %y1(n)長度
Ny2=length(y2);  %y2(n)長度
subplot(3,2,1);stem(n,x1,'.');  
xlabel('n');ylabel('x1(n)');axis([0,Ny1,0,1]);
subplot(3,2,3);stem(nh,h1,'.');  
xlabel('n');ylabel('h1(n)');axis([0,Ny1,0,1]);
n=0:Ny1-1;
subplot(3,2,5);stem(n,y1,'.');  
xlabel('n');ylabel('y1(n)');axis([0,Ny1,0,max(y1)]);
n=0:length(x2)-1;
subplot(3,2,2);stem(n,x2,'.');
xlabel('n');ylabel('x2(n)');axis([0,Ny2,0,1]);
subplot(3,2,4);stem(nh,h2,'.');  
xlabel('n');ylabel('h2(n)');axis([0,Ny2,0,1]);
n=0:Ny2-1;
subplot(3,2,6);stem(n,y2,'.');
xlabel('n');ylabel('y2(n)');axis([0,Ny2,0,max(y2)])